首页 > 编程知识 正文

rand函数 matlab(matlab实现数值积分 【一】(trapz函数))

时间:2023-05-06 15:30:11 阅读:121347 作者:3667

综述函数调用格式的应用实例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的减小,计算精度逐渐增加。

版权声明:该文观点仅代表作者本人。处理文章:请发送邮件至 三1五14八八95#扣扣.com 举报,一经查实,本站将立刻删除。