数值计算方法实验报告(一)算法的数值稳定性实验
数值分析实验报告
姓 名学 号系 别数学系班级主讲教师指导教师实验日期专业课程名称数值计算方法同组实验者无一、实验名称:
实验一、算法的数值稳定性实验1.进行Matlab语言的编程训练,初步体验算法的软件实现;
2.通过对稳定算法和不稳定算法的结果分析、比较,深入理解算法的数值稳定性及其重要性。1.用Matlab语言编写按递推公式
计算(n=0,1,2,……)的程序,并取I0=0.6321,计算n=0,1,2,…,9时In的值。
2. 用Matlab语言编写按递推公式
计算(n=0,1,2,……,9)的值程序。
3.分析比较两种算法的数值稳定性。,要运用下面两种算法进行计算。 .
算法一:
由积分估值得:
取y(0)=0.6321时,运用递推公式为:
进行计算。
算法二:
取n=9时,得:
所以y(9)=0.0684时,运用递推公式为:
进行计算。
算法一的程序如下:
function [y,n] =myfun()
y_0=0.6321;
n=1;
y(1)=y_0;
for i=2:10
y(i)=1-(i-1)*