综述函数调用格式的应用实例1 :用梯形法求积分的实例2 :不同步骤对积分结果的影响
总的说来
数值积分问题是传统数值分析课程的重要内容。 如果被积函数公式未知,则需要根据实测数据用梯形算法求出积分的近似值。 介绍被积函数公式未知时数值积分问题的求解方法,即已知数据点求积分。
函数调用格式s=trapz(x,y ); 应用实例1 :梯形积分
x=[0:pi/30:pi] '; y=[sin(x ) cos (x ) sin ) x/2 ) ] s=trapz(x,y )的结果是s=[1. 99820.0000 ] 1.9995 ]
选择的步长为h=/30=0.1h=pi/30=0.1h=/30=0.1,因此结果存在较大误差。 其实积分问题与样条插值技术相结合,可以给出一个可以精确计算积分的MATLAB函数。 ((补充) ) ) ) ) ) )。
例2 :不同步长对积分结果的影响题目:用定步长法求解积分032 cos 15 xdxint _0^ {frac {3 pi } {2} { cos 15 xdx }
首先,创建被积函数的图像。 x=[ 0:0.01:3 * pi/2,3 * pi/2 ]; % //这样代入的话,可以确保包含3*pi/2点的y=cos(15*x )。 打印(x,y ) )。
从图像中观察到解区域内的被积函数有很强的振动。
不同步幅h=0.1、0.01、0.001、0.0001、0.00001、0.00001h=0.1、0.0001、0.00001、0.00001h=0.1、0.0001 3*pi/2 ) /求出理论值为1/15h0=[ 0.1,0.01,0.001,0.0001,0.0001,0.00001 ] v=[ ]; forh=h0x=[ 0: h :3 * pi/2,3 * pi/2 ]; y=cos(15*x ); I=trapz(x,y ); v=[v; h,I,a-I; 结束的结果如下。
可见,随着步幅h h h的减小,计算精度逐渐增加。