数学实验报告(样例)(数值积分)

《数学实验报告(样例)(数值积分)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数学与物理系实验报告（2011级）课程名称：数值分析实验题目：数值积分专业：2011数学与应用数学班级：姓名：学号：任课教师：成绩评定：填写日期：年月日合肥学院数学与物理系制数学实验报告实验序号：日期：20年月日班级姓名学号实验名称数值积分问题背景描述：利用牛顿莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适用于被积函数的原函数能用初等函数表达出来的情形。如果这点办不到或者不容易办到，这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是—条实验记录曲线或者是一组离

2、散的采样值，这时只能应用近似方法去计算相应的定积分。实验目的：本实验将主要研究定积分的三种近似计算算法：复化（梯形、抛物线、牛顿柯特斯）算法，在不同精度要求下，看等分点数与运行时间，从而提出一些想法。同时也研究一下外推法。实验原理与数学模型：这里还应把牛顿柯特斯公式打上。实验所用软件及版本：Matlab7.0主要内容（要点）：对三种数值积分，采用不同的精度要求，计算对分点数与耗时，采用拟合曲线的方法，去推想更高精度要求时所花时间。改进的方法有外推法。从60年代开始，中国计算科学家们为中国原子弹的研制成功、

3、氢弹原理的突破和发展作出了巨大的贡献。中国当年在核武器研制上达到与美、苏基本相抗衡的水平时，前后只进行过338次核试验。而两个超级大国则分别进行了936次和716次。这其中的奥密从某种程度上讲就是由于我国的研制人员在经费和物质条件的限制下，更多地依靠了数值模拟实验手段：试想一想，这为国家节省了多少试验物资和经费。邓稼先在学海茫茫的大千世界中，终于找到了天体物理中的托马斯——费米理论。这个理论过去是用到中子星球上的，拿原子弹和中子星相比，温度大致差不多，都有极高的温度，但密度却没有中子星那么高。所以必须修正

4、托马斯——费米理论来推导出原子弹所需的高温高压下的状态方程，他们创造性地利用外推法，求出了极高温高压下的核材料的状态方程，并巧妙地与低压状态方程连接，得出了相当大区域内完整的状态方程，满足了理论设计要求。（这里写上是进行爱国主义教育。正式报告可以不写）也用了误差估计的方法推导出所需要的对分点数。实验过程记录（含：基本步骤、主要程序清单及异常情况记录等）：%%%fhtxgs_symbol.mclearall;formatlongg;symsfxwcnabss1s2;f=input('f=');a=input

5、('a=');b=input('b=');jieguo=[];forwc=1:9tic;n=1;s1=(subs(f,'x',a)+subs(f,'x',b))*(b-a)/n/2;n=2;s2=(subs(f,'x',a)+subs(f,'x',b))*(b-a)/n/2+subs(f,'x',a+(b-a)/n)*(b-a)/n;whileabs(s2-s1)>0.5*10^(-wc)s1=s2;n=n*2;s=0;h=(b-a)/n;fori=1:n-1s=s+subs(f,'x',a+i*h);e

6、nd;s=h*s+(subs(f,'x',a)+subs(f,'x',b))*h/2;s2=s;end;jieguo=[jieguo;[vpa(toc,10)vpa(0.5*10^(-wc),wc)nvpa(s2,wc)vpa(pi,wc+1)vpa(int('4/(1+x*x)',0,1)-s,20)]];end;jieguo%%%fxspsgs_symbol.mclearall;formatlong;symsfxabns1s2swchi;f=input('f=');a=input('a=');b=in

7、put('b=');jieguo=[];forwc=1:9%%%wc=15?tic;n=1;s1=(subs(f,'x',a)+4*subs(f,'x',(a+b)/2)+subs(f,'x',b))*(b-a)/n/6;n=2;s2=s1+1;whileabs(s2-s1)>0.5*10^(-wc)s1=s2;n=n*2;s=0;h=(b-a)/n;fori=1:ns=s+4*subs(f,'x',a+(i-1/2)*(h));end;fori=1:n-1s=s+2*subs(f,'x',a+i*h);

8、end;s=(subs(f,'x',a)+subs(f,'x',b)+s)*h/6;s2=s;s2-s1;end;jieguo=[jieguo;[vpa(toc,10)vpa(0.5*10^(-wc),wc)nvpa(s2,wc)vpa(pi,wc+2)vpa(vpa(int('4/(1+x*x)',0,1),wc)-s,wc+2)]];end;jieguo%newton_cotes_symbol.mclearall;clc;