MATLAB程序设计与应用95311.docVIP

《Matlab程序设计与应用》课 程 论 文 定积分的近似计算与Matlab实现 姓名：陈林 学号：110314233 班级：统计学二班 摘要： 定积分计算的基本公式是Newton-Leibniz公式，但对于较复杂的函数要求其精确值是比较困难的，而且很多形式较简单的初等函数的原函数也不是初等函数，如∫sinxx dx，∫e-x2dx，∫dxlnx等等，其原函数都不是初等函数，无法利用Newton-Leibniz公式计算。因此，研究定积分的近似计算是非常必要的，而且，可以通过计算机编程来方便计算定积分的近似值。本文考虑了定积分近似计算的矩形法、梯形法、抛物线法和对这三种方法进行修改后的定积分近似计算公式,并给出了每种方法的Matlab程序。 关键词： Matlab；定积分；近似计算；矩形法；梯形法；抛物线法；修改公式 引言： 定积分计算的基本公式是牛顿－莱布尼兹公式。但被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线法。 实验内容： 1、矩形法 根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即 在几何意义上，这是一系列小形面积近似曲边梯形的结果，所以把这个近似计算方法称为矩形法。不过，只有当积分区间被分割的很细时，矩形法才有一定的精确度。 2、梯形法 等分区间 x0=ax1…xi=a+[(b-a)/n]*i…xn=b，x=(b-a)/n 相应函数值为 y0，y1，…，yn 曲线上相应的点为 P0，P1，…，Pn 将曲线的每一段弧用过点，pi的弦（线性函数）来代替，这使得每个上的曲边梯形成为真正的梯形，其面积为 于是各个小梯形面积之和就是曲边梯形面积的近似值， 即 称此式为梯形公式。 3、抛物线法 由梯形法求近似值，当为凹曲线时，它就偏小；当为凸曲线时，它就偏大。若每段改用与它凸性相接近的抛物线来近似时，就可减少上述缺点，这就是抛物线法。 4、改写矩形法和抛物线法 现给出利用求和函数，避免了for循环的梯形法程序如下： format long n=100;a=0;b=1; syms x fx fx=1/(1+x^2); i=1:n; xj=a+(i-1)*(b-a)/n; %所有左点的数组 xi=a+i*(b-a)/n; %所有右点的数组 fxj=subs(fx,x,xj); %所有左点值 fxi=subs(fx,x,xi); %所有右点值 f=(fxi+fxj)/2*(b-a)/n; %梯形面积 inum=sum(f); %加和梯形面积求解 integrate=int(fx,0,1); integrate=double(integrate) fprintf(The relative error between inum and real-value is about: %g\n\n,... abs((inum-integrate)/integrate)) 学习此程序设计方法，用函数sum改写矩形法和抛物线法的程序，避免for循环。改写后的矩形法程序如下： format long n=100;a=0;b=1; syms x fx fx=1/(1+x^2); i=1:n xj=a+(i-1)*(b-a)/n; xi=a+i*(b-a)/n; fxj=subs(fx,x,xj); fxi=subs(fx,x,xi); f=(fxi+fxj)/2*(b-a)/n; inum=sum(f); integrate=int(fx,0,1); integrate=double(integrate) fprintf(The relative error between inum1 and real-value is about: %g\n\n,... abs((inum1-integrate)/integrate)) fprintf(The relative error between inum2 and real-value is about: %g\n\n,... abs((inum2-integrate)/integrate)) f