MATLAB实验《数值微分与积分》

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

X=unifrnd(10,99,1,10)

deltaX1=diff(X,1)

deltaX2=diff(X,2)

deltaX3=diff(X,3)

 

2.求矩阵在指定点的差分。

for x=1:3

    y=[x,x^2,x^3;1,2*x,3*x;0,2,6*x];

    diff(y)

End

 

3.求定积分

%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)

 

 4.求二重定积分。

%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)

 

5，分别用矩形、梯形 (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)