原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
Matlab 数值积分与数值微分
Matlab 数值积分
一重数值积分的实现方法
变步长辛普森法、高斯 - 克朗罗德法、梯形积分法
变步长辛普森法
|精.
|品.
|可.
|编.
|辑.
|学.
|习.
|资.
|料.
*
|
*
|
*
|
*
|
|欢.
|迎.
|下.
|载.
Matlab 供应了 quad 函数和 quadl 函数用于实现变步长辛普森法求数值积分 . 调用格式为 :
[I,n]=Quad(fname,a,b,tol,trace) [I,n]=Quadl(fname,a,b,tol,trace)
Fname 是函数文件名， a,b 分别为积分下限、积分上限； tol 为精度掌握，默认为
1.0 × 10-6 ，trace 掌握是否绽开积分过程，假如为 0 如此不绽开，非 0 如此绽开， 默认不绽开 .
返回值 I 为积分数值； n 为调用函数的次数 .
例如：求 的值 .
先建立函数文件
fesin.m
function f=fesin(x) f=exp(-0.5*x).*sin(x+(pi/6));
再调用 quad 函数
[I,n]=quad(fesin,0,3*pi,1e-10) I =
0.9008
n =
365
例如：分别用 quad 函数和 quadl 函数求积分
|精.
|品.
|可.
|编.
|辑.
|学.
|习.
|资.
|料.
*
|
*
|
*
|
*
|
|欢.
|迎.
|下.
|载.
的近似值，比拟函数调用的次数 .
先建立函数文件
fesin.m
function f=fesin(x) f=exp(-0.5*x).*sin(x+(pi/6)); format long
[I,n]=quadl(fesin,0,3*pi,1e-10)
I =
0.900840787818886
n =
198
[I,n]=quad(fesin,0,3*pi,1e-10) I =
0.900840787826926
n =
365
可以发觉 quadl 函数调用原函数的次数比 quad 少，并且比 quad 函数求得的数值解更精确 .
|精.
|品.
|可.
|编.
|辑.
|学.
|习.
|资.
|料.
*
|
*
|
*
|
*
|
|欢.
|迎.
|下.
|载.
高斯- 克朗罗德法
Matlab 供应了自适应高斯 - 克朗罗德法的 quadgk 函数来求震荡函数的定积分， 函数的调用格式为：
[I,err]=quadgk(fname,a,b)
Err 返回近似误差 X 围，其他参数的意义与 quad 函数一样，积分上下限可以是 -Inf
或 Inf ，也可以是复数，假如为复数如此在复平面上求积分 .
例如：求积分
的数值 .
先编写被积函数的 m文件fsx.m
function f=fsx(x) f=x.*sin(x)./(1+cos(x).^2);
再调用 quadgk 函数
I=quadgk(fsx,0,pi) I =
2.4674
2.4674
例如：求积分
的值 .
|精.
|品.
|可.
|编.
|辑.
|学.
|习.
|资.
|料.
*
|
*
|
*
|
*
|
|欢.
|迎.
|下.
|载.
先编写被积函数的 m 文件fsx.m
function f=fsx(x) f=x.*sin(x)./(1+cos(x).^2);
再调用 quadgk 函数
I=quadgk(fsx,-Inf,Inf) I =
-9.0671e+017
梯形积分法
对于一些不知道函数关系的函数问题，只有试验测得的一组组样本点和样本值，由
表格定义的函数关系求定积分问题用梯形积分法， 其函数是 trapz 函数， 调用格式为：
I=Traps(X,Y)
X,Y 为等长的两组向量，对应着函数关系 Y=f(X)
X=(x 1,x 2, ?,x n)(x 1x 2 ?x n),Y=(y 1 ,y 2, ?,y n ) ，积分区间是 [x 1 ,x n]
例如：某次物理试验测得如下表所示的两组样本点 .
x
1.38
1.56
2.21
3.97
5.51
7.79
9.19
11.12
13.39
y
3.35
3.96
5.12
8.98
11.46
17.63
24.41
29,83
32.21
现变量 x 和变量 y 满意肯定的函数关系，但此关系未知，设 y=f(x) ，求积分的数值 .
X=[1.38,1.56,2.21,3.97,5.51,7.79,9.19,11.12,13.39];
Y=[3.35,3.96,5.12,8.98,11.46,17.63,24.41,29.83,32.21];
|精.
|品.
|可.
|编.
|辑.
|学.
|习.
|资.
|料.
*
|
文档评论（0）