资源描述:
《数值积分用matlab实现》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、东北大学秦皇岛分校数值计算课程设计报告数值积分及Matlab实现学院数学与统计学院专业信息与计算科学学号5133117姓名楚文玉指导教师张建波姜玉山成绩教师评语:指导教师签字:2015年07月14日数学与统计学院课程设计(实习)报告第14页1绪论在科研计算中,经常会碰到一些很难用公式定理直接求出精确解的积分问题,对于这类问题,我们一般转化为数值积分问题,用计算机来实现求解问题.1.1课题的背景对于定积分在求某函数的定积分时,在一定条件下,虽然有牛顿-莱布里茨公式可以计算定积分的值,但在很多情况下的原函数不易求出或非常
2、复杂.被积函数的原函数很难用初等函数表达出来,例如等;有的函数的原函数存在,但其表达式太复杂,计算量太大,有的甚至无法有解析表达式.因此能够借助牛顿-莱布尼兹公式计算定积分的情形是不多的.另外,许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数,对这类函数的定积分,也不能用不定积分方法求解,只能设法求其近似值.因此,探讨近似计算的数值积分方法是有明显的实际意义的,即有必要研究定积分的数值计算方法,以解决定积分的近似计算.而数值积分就是解决此类问题的一种有效的方法,它的特点是利用被积函数在一些节点上的信息求出
3、定积分的近似值.微积分的发明是人类科学史上一项伟大的成就,在科学技术中,积分是经常遇到的一个重要计算环节数值积分是数学上重要的课题之一,是数值分析中重要的内容之一.随着计算机的出现,近几十年来,对于数值积分问题的研究已经成为一个很活跃的研究领域.现在,数值积分在计算机图形学,积分方程,工程计算,金融数学等应用科学领域都有着相当重要的应用,所以研究数值积分问题有着很重要的意义.国内外众多学者在数值积分应用领域也提出了许多新方法.在很多实际应用中,只能知道积分函数在某些特定点的取值,比如天气测量中的气温、湿度、气压等,医
4、学测量中的血压、浓度等等.通过这个课题的研究,我们将会更好地掌握运用数值积分算法求出特殊积分函数的定积分的一些基本方法、理论基础;并且通过Matlab软件编程的实现,应用于实际生活中.1.2课题的主要内容框架1.2.1数值积分各求积公式简介数学与统计学院课程设计(实习)报告第14页简介牛顿-柯特斯求积公式及其辛普森求积公式,龙贝格求积公式,高斯求积公式的基本理论基础和方法.1.2.2求积公式的代码实现通过理解各种数值积分求积公式的原理方法,通过Matlab软件编程,实现以上求积公式.1.2.3应用举例通过简单举例,自
5、建一个相对简单和复杂的函数,用上面编写的Matlab源程序来解决实际问题,体会数值积分和Matlab的优势.2牛顿-柯特斯公式及Matlab实现2.1牛顿-柯特斯公式的基本原理方法设将积分区间[a,b]划分为n等分,步长为,选取等距节点构造出的差值型求积公式,(2.1)称为牛顿-柯特斯公式,式中称为柯特斯系数.根据,(2.2)引进变量代换,则有(2.3)当n=2时,此时柯特斯系数为,,,相应的求积公式就是辛普森求积公式:(2.4)2.2牛顿-柯特斯公式的Matlab实现function[C,g]=NCotes(a,b
6、,n,m)数学与统计学院课程设计(实习)报告第14页%a,b分别为积分的上下限;%n是子区间的个数;%m是被调用第几个被积函数;%当n=1时计算梯形公式;当n=2时计算辛普森公式,以此类推;I=n;h=(b-a)/i;z=0;forj=0:ix(j+1)=a+j*h;s=1;ifj==0s=s;elsefork=1:js=s*k;endendr=1;ifi-j==0r=r;elsefork=1:(I-j)r=r*k;endendifmod((I-j),2)==1q=-(I*s*r);elseq=i*s*r;endy=
7、1;fork=0:i数学与统计学院课程设计(实习)报告第14页ifk~=jy=y*(sym('t')-k);endendl=int(y,0,i);C(j+1)=l/q;z=z+C(j+1)*f1(m,x(j+1));endg=(b-a)*z3复合求积公式及Matlab实现3.1复合梯形公式的基本原理将区间[a,b]划分成n等分,分点,,在每个子区间[]()上采用梯形公式得:(3.1)记(3.2)称式(3.2)为复合梯形公式.3.2复合梯形公式的Matlab实现functions=trapr1(f,a,b,n)%f表示
8、被积函数;%a,b表示积分上下限;%n是子区间的个数;h=(b-a)/n;s=0;fork=1:(n-1)x=a+h*k;s=s+feval('f',x);数学与统计学院课程设计(实习)报告第14页endformatlongs=h*(feval('f',a)+feval('f',b))/2+h*s;3.3复合辛普森求积公式的基本原理将区间[a,