




已阅读5页,还剩3页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
实验报告实验项目名称 数值积分与数值微分 实验室 数学实验室 所属课程名称 数值逼近 实 验 类 型 算法设计 实 验 日 期 班 级 学 号 姓 名 成 绩 实验概述:【实验目的及要求】本次实验的目的是熟练数值分析第四章“数值积分与数值微分”的相关内容,掌握复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式。本次试验要求编写复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式的程序编码,并在MATLAB软件中去实现。【实验原理】数值分析第四章“数值积分与数值微分”的相关内容,包括:复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式的相应算法和相关性质。【实验环境】(使用的软硬件)软件:MATLAB 2012a硬件:电脑型号:联想 Lenovo 昭阳E46A笔记本电脑操作系统:Windows 8 专业版 处理器:Intel(R)Core(TM)i3 CPU M 350 2.27GHz 2.27GHz实验内容:【实验方案设计】第一步,将书上关于复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式的内容转化成程序语言,用MATLAB实现;第二步,分别用以上求积公式的程序编码求解不同的问题。【实验过程】(实验步骤、记录、数据、分析)实验的主要步骤是:首先分析问题,根据分析设计MATLAB程序,利用程序算出问题答案,分析所得答案结果,再得出最后结论。实验:用不同数值方法计算积分 (1)取不同的步长h.分别用复合梯形及复合辛普森求积计算积分,给出误差中关于h的函数,并与积分精确值比较两个公式的精度,是否存在一个最小的h,使得精度不能再被改善?(2)用龙贝格求积计算完成问题(1)。(3)用勒让德多项式确定零点,再代入计算高斯公式,使其精度达到10-4(1)在MATLAB的Editor中建立一个M-文件,输入程序代码,实现复合梯形求积公式的程序代码如下:function s=T(n)a=0.0000001;b=1;h=(b-a)/n;s=h*(f(a)+f(b)/2;if n1for k=1:n-1x=a+k*h;s=s+h*f(x);endE=s+4/9%复合梯形误差end在command Windows中输入命令:T(10),T(100)以及T(1000),得出的结果为:T(10)E = 0.0271ans = -0.4173 T(100)E = 0.0013ans = -0.4431 T(1000)E = 5.4375e-05ans = -0.4444建立一个新的M-文件,输入程序代码,实现切比雪夫多项式的程序代码如下: function t=S(n)a=0.0000001;b=1;h=(b-a)/n;t=h*(f(a)+f(b)/6;if n1for k=0:n-1 x0=a+(k+0.5)*h; x1=a+k*h; if k=0 t=t+4*f(x0)*h/6; else t=t+(4*f(x0)+2*f(x1)*h/6; endendE=t+4/9%复合辛普森误差end在command Windows中输入命令:S(10),S(100)以及S(1000),得出的结果为:S(10)E = 0.005ans = -0.4387 S(100)E = 2.4147e-04ans = -0.4442 S(1000)E = 9.1563e-06ans = -0.4444总结由结果(1)、(2)可知复合辛普森法求积分精度明显比复合梯形法求积的精度要高,且当步长取不同值时即n越大、h越小时,积分精度越高。实验结果说明不存在一个最小的h,使得精度不能再被改善。又两个相应的关于h的误差(余项) 其中属于a到b。可知h愈小,余项愈小,从而积分精度越高。(2)在MATLAB的Editor中建立一个M-文件,输入程序代码,实现龙贝格算法的程序代码如下:function q,n=Roberg(f,a,b)M=1;abs0=10;k=0;T=zeros(1,1);h=b-a;T(1,1)=(h/2)*(subs(f,a)+subs(f,b);while abs00.0001 k=k+1; h=h/2; p=0; for i=1:Mx=a+h*(2*i-1);p=p+subs(f,x); endT(k+1,1)=T(k,1)/2+h*p;M=2*M;for j=1:kT(k+1,j+1)=(4j)*T(k+1,j)-T(k,j)/(4j-1);endabs0=abs(T(k+1,j+1)-T(k,j);endq=T(k+1,k+1);n=k;在command Windows中输入命令:Fx,n=Roberg(sqrt(x)*log(x),10(-8),1),得出的结果为:Fx = -0.444387313932947n = 9(3)在MATLAB的Editor中建立一个M-文件,输入程序代码,实现龙贝格算法的程序代码如下:function ql,Ak,xk=guasslegendre(fun,a,b)% fun:被积函数% a,b:积分上下限% ql:积分结果% Ak:系数% xk:零点n=0;fun=(x)sqrt(x).*log(x);eps=1;while eps0.0001 syms x p=sym2poly(diff(x2-1)(n+1),n+1)/(2(n+1)*factorial(n+1); tk=roots(p); Ak=zeros(n+1,1); for i=1:n+1 xkt=tk; xkt(i)=; pn=poly(xkt); fp=(x)polyval(pn,x)/polyval(pn,tk(i); Ak(i)=quadl(fp,-1,1,0.0001);end xk=(b-a)/2*tk+(b+a)/2;fx=(b-a)/2).*fun(xk);ql=sum(Ak.*fx);eps=abs(ql+4/9);n=n+1;endE=eps在command Windows中输入命令:guasslegendre(sqrt(x)*log(x),10(-8),1),得出的结果为: guasslegendre(sqrt(x)*log(x),10(-8),1)E = 9.4667e-05ans = -0.4445【结论】(结果)复合求积法相比普通的求积公式而言精度要高,其中复合辛普森法求积分精度比复合梯形法求积的精度要高,龙贝格求积法使等距节点求积精度进一步提高。高斯求积公式具有最高代数精度。【小结】通过
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 针对于小样本数据的换热器污垢预测方法进行研究
- 2024年碳排放监测设备安装调试合规考核试卷
- 购买热镀锌合同(标准版)
- 简易道路施工合同(标准版)
- 2024年贵州电网招聘真题
- 厦门市湖里区产业投资集团有限公司人员招聘考试真题2024
- 城南旧事感想200个字
- 2025年电大房地产开发与管理专业《房地产估价》真题汇编及答案解析
- 考点解析-人教版八年级物理上册第5章透镜及其应用-透镜专项攻克练习题(含答案详解)
- 考点攻克人教版八年级物理上册第5章透镜及其应用-透镜专项训练试卷(含答案详解版)
- 未来工程师的展望-展示专业能力共创未来趋势
- 事业机关食堂承包合同6篇
- 整式(第一课时)课件冀教版七年级数学上册
- 《四川省电动汽车充电设施建设技术标准》
- 厂房办公室装饰施工方案
- 浙江省天域全国名校协作体2026届高三上学期10月联考政治试题(含答案)
- 医养结合政策课件
- 环卫机械安全培训内容课件
- 财务报表数据核对与审计清单
- 人工智能+智能制造人才培养策略研究报告产业人才需求分析
- 北师大版三年级上册第四单元《我们生活的空间(一)》单元检测卷(含答案)
评论
0/150
提交评论