MATLAB数值积分与微分.docVIP

1. 数值微分 MATLAB没有直接求数值导数的函数，diff函数可求数值微分，其调用格式为： DX=diff(X)：DX(i)=X(i+1)-X(i)，i=1,2,…,n-1。 DX=diff(X,n)：n阶 DX=diff(A,n,dim)：计算矩阵A，dim=1时(缺省状态)，按列；dim=2，按行。 V=vander(1:4) DV=diff(V) gradient函数可按行或列求矩阵的数值微分 [dx,dy]=gradient(V) %除了第一行(列)与最末一行(列)按前向公式，其余的dx=(y(:,i+1)-y(:,i-1))/2; dy=(y(i+1,:)-y(i-1,:))/2; gradient(V) 编程实现以上求一阶或二阶近似的算法(diffc.m) h=0.05;x=0:h:pi; y=sin(x)./(x+cos(x)); y1=diffc(y,h,1,2); %n=1 y2=diffc(y,h,2,2); %n=2 y3=diffc(y,h,3,2); %n=3 y4=diffc(y,h,4,2); %n=4 subplot(2,2,1);plot(x,y,x(1:length(x)-4),y1); subplot(2,2,2);plot(x,y,x(1:length(x)-4),y2);subplot(