首页 > 学术百科

MATLAB数值计算函数汇总

%% 向量最⼤值

% X=[1,2,3,4,5,7,6];

% y=max(X);%最⼤值

% [y,k]=max(X);%最⼤值及位置

%% 矩阵最⼤值

% X=[1,2,3;4,5,7;6,1,1];

% max(X)%求每列最⼤元素

% max(X,[],2)%求每⾏最⼤元素

% max(max(X))%求整个矩阵最⼤元素

%% 向量或矩阵对应元素⽐较取最⼤

% x=[4,5,6;1,4,8]

% y=[1,7,5;4,5,7]

% p=max(x,y)

湖北财税厅

%% 向量或矩阵与常数⽐较取最⼤

% x=[4,5,6;1,4,8]

% y=2

% p=max(x,y)

%% 求平均值和中值和标准⽅程

% x=[1,2,3,6,7,6];

% mean(x)

% median(x)

%% 矩阵元素求和和求积

% x=[4,5,6;1,4,8;2,5,6]

% sum(x)%列相加

% sum(x,2)%⾏相加

% prod(x)%列相乘

% prod(x,2)%⾏相乘

%% 矩阵元素累加和累乘

% x=[4,5,6;1,4,8;2,5,6]

% cumsum(x)%列累加

% cumsum(x,2)%⾏累加

% cumprod(x)%列累乘

% cumprod(x,2)%⾏累乘

%% 向量或矩阵⽅差

% x=[4,5,6;1,4,8;2,5,6]

% std(x,0,1)%列⽅差

% std(x,0,2)%⾏⽅差

%% 相关系数

% x=randn(10000,5);

% corrcoef(x)%把每列作为变量求他们的相关系数

% X=[1,2,3];Y=[3,2,1];

% corrcoef(X,Y)%求向量X,Y的相关系数

%% 元素排序不符合项报告

% x=[4,5,6;1,4,8;2,1,6]

% sort(x)%按列升序

% sort(x,2)%按⾏升序

% -sort(-x)%按列降序

% -sort(-x,2)%按⾏降序

% [a,b]=sort(x)%b为所在列号

%% ⼀维插值 t ⾃变量 f因变量 x为插值插值适⽤点相对准确⽆异常点

% t=0:5:30;

% f=[3,2,879,1835,2968,4136,5237];

% format bank

% interp1(t,f,x,'linear') %线性插值⽅法计算

% interp1(t,f,x,'nearest') %最近点插值⽅法计算

% interp1(t,f,x,'spline') %3次样条插值⽅法计算

% interp1(t,f,x,'PCHIP') %3次多项式插值⽅法计算

%% ⼆维插值 x,y⾃变量 z因变量

% x=0:0.1:1;

% y=0:0.1:2;

% [X,Y]=meshgrid(x,y);%划分成采样⽹格

% Z=X.^2+Y.^2;

% interp2(x,y,Z,0.5,0.5) %在(0.5,0.5) 插值

% interp2(x,y,Z,[0.5 0.6]',[0.4 0.5])

%% ⼆维插值 x,h⾃变量 T因变量

% x=0:2.5:10;

% h=[0:30:60]';

% T=[95,14,0,0,0;

% 88,48,32,12,6;

% 67,64,54,48,41];

% xi=[0:0.5:10];

% hi=[0:10:60]';

% temps=interp2(x,h,T,xi,hi,'spline');

% mesh(xi,hi,temps);

%% 曲线拟合和最⼩⼆乘法拟合适⽤点的误差可以⼤

% X=linspace(0,2*pi,50)%0-2pi分成50份

% Y=sin(X);

% P=polyfit(X,Y,3);%得到3次多项式的系数

% xi=linspace(0,2*pi,100);

% yi=sin(xi);

% y1=polyval(P,xi);%按照多项式系数计算拟合值

% plot(X,Y,xi,y1,':o')

%% 多项式系数加减乘除

% f=[3,-5,2,-7,5,6];

% g=[3,5,-3];

% g1=[0,0,0,g];

% f+g1

% f-g1

% conv(f,g)%乘法

% [Q,r]=deconv(f,g)%除法商送Q，余送r

%% 多项式的导函数

% P=[3,5,0,-8,1,-5];

% Q=[10,5,0,0,6,0,0,7,-1,0,-100];

% [p,q]=polyder(P,Q)%p为分⼦,q为分母

%% 多项式求值

% A=[1,8,0,0,-10];

% x=1.2;

% y1=polyval(A,x)%x=1.2时，多项式的值

% x=[-1,1.2,-1.4;2,-1.8,1.6;1,2,3];

% y2=polyval(A,x)%x=矩阵，各元素多项式的值

% y3=polyvalm(A,x)%x=矩阵(必须为⽅阵)，矩阵进⾏多项式计算的值

%% 多项式求根已知根求多项式

% A=[1,8,0,0,-10];

% x=roots(A)王治坪

% B=poly(x)

% f=inline('(x+1).^2')

% g=inline('2.*x+2')

% x=-3:0.01:3;

% p=polyfit(x,f(x),2);%求拟合多项式系数拟合求导

% dp=polyder(p);%求导数多项式的系数

% dpx=polyval(dp,x);%求导函数对应函数值

% plot(x,dpx)

%-----------------------------------------------

% dx=diff(f([x,3.01]))/0.01;% 直接对f(x)求数值导数微分求导% plot(x,dx);

%-----------------------------------------------

% gx=g(x) % 导函数求导数

% plot(x,gx);

%% 数值积分被积函数是⼀个表达式

% format long

% I=quad('exp(-x.^2)',0,1)

%% 数值积分被积函数是⼀个表格 (离散的点)

东京奥运会必将如期举行% x=0:0.01:1;

% y=exp(-x.^2);

% trapz(x,y)

%% ⼆重数值积分被积函数是⼀个表达式

% format long

% f=inline('exp(-x.^2/2).*sin(x.^2+y)','x','y');

% I=dblquad(f,-2,2,-1,1)

%% 给定函数进⾏傅⾥叶变换

% N = 128; %采样128个点

% T = 1; %⾃变量范围0-1

% t=linspace(0,T,N); %分割

% x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t);

凯恩斯革命% dt=t(2)-t(1) %采样周期

% f=1/dt; %采样频率山炮绿

% X=fft(x); %傅⾥叶变换

% f=f*(0:N-1)/N %频率轴

% plot(f,abs(X),'-*') %振幅，频率曲线

% xlabel('Frequency');

% ylabel('|F(k)|');

% ix=real(ifft(X)); %逆变换

% plot(t,x,'-*',t,ix,':')

% norm(x-ix)

作者：ChenBD

本文发布于:2024-09-22 19:30:12，感谢您对本站的认可！

本文链接：https://www.17tex.com/xueshu/4599.html

上一篇：eviews_基本命令

下一篇：基于有理插值样条的非线性回归算法

标签：函数矩阵元素计算变量

留言与评论（共有 0 条评论）