matlab求解一元函数积分,matlab实现一元函数积分

一元函数的数值积分

● 函数1 trapz

功能 复合梯形求积公式

格式 T=trapz(Y) %用等距复合梯形求积公式近似计算Y 的积分。若Y 是一向

量，则trapz(Y)为Y 的积分；若Y 是一矩阵，则trapz(Y)

为Y 的每一列的积分。

T = trapz(X,Y) %用复合梯形求积公式计算Y 在X 点上的积分，X 和Y 可以

是负数。若X 为一列向量，Y 为矩阵，且size(Y,1) =

length(X)，则trapz(X,Y)通过Y 的第一个非单元集方向

进行计算。

T = trapz(…,dim) %沿着dim 指定的方向对Y 进行积分。若参量中包含X ，

则应有length(X)=size(Y,dim)。

例1 计算?π0sin xdx 的数值积分。计算可得2sin 0=?π

xdx 。 准备工作：

>>X = 0:pi/100:pi;

>>Y = sin(X);

函数调用：

>>T = trapz(X,Y)

可得

>>T=

1.9998

或者

>>T = pi/100*trapz(Y)

可得

>>T=

1.9998

例2 计算复积分

>>z = exp(i*pi*(0:100)/100);

>>trapz(z, 1./z)

可得

>>ans =

0.0000 + 3.1411i

produce

Z =

1.9998

● 函数2 quad 、quadl 、quad8

功能 数值定积分，自适应复合Simpson 求积公式。

格式 q = quad(fun,a,b) %近似地从a 到b 计算函数fun 的数值积分，误差为10-6。

若给fun 输入向量x ，应返回向量y ，即fun 是一单值函