Matlab期末考 第五章数值积分和微分

（一）数值积分

积分函数

（1）sum(y*dx)，y为积分函数，dx为x的微分

（2）trapz(y*dx)

trapz使用梯形法对函数进行近似积分，sum是用长方形法对函数近似积分，trapz函数的误差比sum的误差更小

例子

求y=x^2的积分，积分上限是3，下限是0

dx=0.1
x=0:dx:3;
y=x.^2;
y1=sum(y*dx);
y2=trapz(y*dx);
y1
y2
>>
y1=9.4550
y2=9.0050

积分的答案是9，trapz更接答案。

补充

二重积分：trapz(trapz(z*dx*dy))

例子

对z=sin(x).^2+cos(y).^2，进行二重积分，-10*pi<=(x,y)<=10*pi

x=linspace(-10*pi,10*pi,200);
y=linspace(-10*pi,10*pi,200);
dx=pi/100
dy=pi/100
[X,Y]=meshgrid(x,y);%二维网格化
Z=sin(X).^2+cos(Y).^2;
i=trapz(trapz(Z*dx*dy));
i

（二）数值微分

微分命令

一阶导，