误差传播与算法稳定性实验报告_第1页
误差传播与算法稳定性实验报告_第2页
误差传播与算法稳定性实验报告_第3页
误差传播与算法稳定性实验报告_第4页
误差传播与算法稳定性实验报告_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上数值分析实 验 分析报 告 姓名: 学号: S 日期: 2016.10.15 班级: 1602 一、实验名称误差传播与算法稳定性二、实验目的体会稳定性在选择算法中的地位。误差扩张的算法是不稳定的，是我们所不期望的；误差衰竭的算法是稳定的，是我们努力寻求的，这是贯穿本课程的目标。三、实验内容计算 四、算法描述利用分部积分可得：可得递推公式为：1、2、下面分别以1，2递推关系求解方案一：当时 （保留六位有效数字）方案二；当时，当时，这里取（保留六位有效数字） 五、程序流程图由于实验方案明显、简单，实现步骤及流程图省略。六、实验结果计算结果如表1-1：1-1 计算结果表10

2、.0.20.0.30.0.40.0.50.0.60.0.70.0.80.0.90.0.100.0.110.0.120.0.130.0.140.0.150.0.160.0.17-1.0.1819.0.19-376.0.207532.790.七、实验结果分析1、通过表1-1可以看出，算法一在前15项中迭代值基本保持一致，但是从16项开始就有了较大的差别。同时：又当时，恒成立。但是算法一中当n=17、19时明显不满足。算法二相比较于算法一在目前的迭代次数下比较精确。2、设算法一中的计算误差为，由递推计算到的误差为；算法二中的计算误差为，由向前递推计算到（）的误差为。设算法一中的精确值为，计算值为；算

3、法二中的精确值为，计算值为根据题设可得：算法一：递推可得：算法二：同理可得：3、根据算法一的公式可得：当时，由此可见最后的结果误差不仅仅取决于初始值的误差，还会随着递推公式的不断运行，误差在不断的增大。当时，由此可见最后的结果误差不仅仅取决于初始值的误差，还会随着递推公式的不断运行，误差在不断的减小。4、通过前三问，可以发现算法二比算法一更加稳定。八、附录（程序）算法一：function ek= jifeng(e0,n ) e(1)=vpa(e0,6); B=zeros(n,2); B(:,1)=1:n;for i=1:n e(i+1)=1-(i+1)*e(i) ; ek=e(i+1); B(i,2)=vpa(ek,6); end xlswrite(算法一,B,sheet1)算法二：function i,ek = jifeng1( e0,n )e(n+1)=e0;B=zeros(n,2);B(:,1)=1:n;for i=n:-1:1 e(i)=(1-e(i+1)/(i+1); ek=e(i); B(i,2)=vpa(ek,6

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论