《matlab数值积分实例.docx》由会员分享，可在线阅读，更多相关《matlab数值积分实例.docx（4页珍藏版）》请在金锄头文库上搜索。
1、matlab数值积分实例matlab数值积分实例 编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（matlab数值积分实例）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快 业绩进步，以下为matlab数值积分实例的全部内容。 数值积分1。求积分，在积分区间中，点与点之间的间隔取为0。1.解:(一）用MATBLE编写复化梯形求积函数:func
2、tion I=T_quad（x，y）n=length（x)；m=length(y）；if n =m errorendh=(x(n）-x(1)）/(n-1）;a=1 2ones（1，n2) 1;I =h/2sum（a。*y);输入：x=1：0。1:1；y=exp（-x.2）;I =T_quad(x,y)运行得到：I = 1.4924(二）用MATBL编写复化Simpson求积函数：function I=S_quad(x，y)；n=length（x）;m=length（y)；if n =m errorendif rem（n1,2) =0 I=T_quad（x,y)； return；endN=（n1
3、）/2;h=(x（n)-x（1）)/N;a=zeros(1,n)；for k=1：N a（2k-1)=a（2k1)+1；a(2k）=a(2k)+4； a（2*k+1)=a(2*k+1）+1;end I=h/6sum（a。*y）；输入：x= 1:0。1:1;y=exp（-x。2）;I= S_quad(x，y)运行得到:I = 1。4936(三)用MATBL编写复化Cotes求积函数：function I=C_quad(x,y)；n=length(x)；m=length（y）;if n =m errorendif rem（n-1，4) =0 I=S_quad（x，y）； returnendN=(n
4、1）/4；h=(x(n）-x（1)）/N；a=zeros（1，n)；for k=1:N a（4k3）=a(4*k3）+7； a(4*k-2）=a(4k-2)+32; a（4k1）=a（4k-1）+12； a(4k）=a（4*k）+32； a（4*k+1)=a(4*k+1）+7；endI=h/90*sum(a。*y);输入:x= 1:0.1:1；y=exp(x。2）;I= C_quad（x，y)运行得:I = 1。4936(四）利用trapz(）函数,采用复化梯形公式求积分输入:x=-1：0。1:1;y=exp（x。2）;I=trapz（x,y)输出：I = 1.49242。 求积分,取精度要求
5、（一）用MATLAB编写自适应步长的梯形公式function I=T_quad_iter（fun,a，b，ep）if nargin4 ep=1e-5;endN=1;h=baT=h/2*（feval（fun，a)+feval(fun,b)）;while 1 h=h/2；I=T/2; for k=1：N I=I+h*feval(fun,a+（2k1）h）; end if abs(I-T)ep break end N=2N; T=IEnd输入:I=T_quad_iter(（x)exp(x。2）,-1,1）输出：I = 1。4936(二）用MATLAB编写自适应步长的Simpson公式function
6、 I=S_quad_iter（fun，a,b,ep)if nargin4 ep=1e5endN=1；h=baT1=h/2*（feval(fun,a)+feval(fun,b)）；S0=T1;while 1 h=h/2;T2=T1/2; for k=1:N T2=T2+h*feval(fun,a+（2*k1）*h） end I=(4*T2T1)/3; if abs(IS0）ep break end N=2N;T1=T2; S0=IEnd输入:I=S_quad_iter(（x） exp（x.2），1，-1)输出：I = 1。4936(三）用quad（）函数,采用自适应步长的Simposon求积分输入：I=quad（(x） exp（x。2) ，-1，1）输出：I = 1.4936
