9、MATLAB程序设计与应用刘卫国（第三版）课后实验九：数值微分与积分

目录

 一、

 二、

 三、

 四、

 五、

---

一、

产生具有10个元素的向量x，其元素是两位随机整数，求x的1~3阶差分。

---------------------------------------示例代码---------------------------------------------

X=unifrnd(10,99,1,10)
deltaX1=diff(X,1)
deltaX2=diff(X,2)
deltaX3=diff(X,3)

---------------------------------------运行结果---------------------------------------------

 二、

求矩阵在指定点的差分。

 ---------------------------------------示例代码---------------------------------------------

for x=1:3
    y=[x,x^2,x^3;1,2*x,3*x;0,2,6*x];
    diff(y)
end

---------------------------------------运行结果---------------------------------------------

 三、

求定积分。

  ---------------------------------------示例代码---------------------------------------------

%1
f=@(t)sqrt(cos(t.^2)+4*sin(2*t).^2+1);
I1=integral(f,0,2*pi)
%2
g=@(x)log(1+x)./(1+x.^2);
I2=integral(g,0,1)

---------------------------------------运行结果---------------------------------------------

 四、

求二重定积分。

   ---------------------------------------示例代码---------------------------------------------

%1
f1=@(x,y) 1./[sqrt(x.^2+y.^2)];
i1=integral2(f1,0,1,0,1);
%2
f2=@(x,y) y.*sin(x)+2.*cos(y)-1
I2=integral2(f2,0,2*pi,0,pi)

---------------------------------------运行结果---------------------------------------------

 五、

分别用矩形、梯形（trapz）公式计算由实验表3中数据给出的定积分

 

  ---------------------------------------示例代码---------------------------------------------

%矩形
x=[0.3:0.2:1.5]
L=ones(length(x),1)
y=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]
Q1=y*L*0.2 
%梯形
x=[0.3:0.2:1.5]
y=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]
I=trapz(x,y) 

---------------------------------------运行结果---------------------------------------------