数值计算方法实验报告(一)算法的数值稳定性实验.docVIP

数值分析实验报告 姓 名 学 号 系 别 数学系 班级 主讲教师 指导教师 实验日期 专业 课程名称 数值计算方法 同组实验者 无 一、实验名称： 实验一、算法的数值稳定性实验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)*y(i-1); n=n+1; end n 算法二的程序如下： function [ y,n ] = Myfun_2() syms k y_9; y=zeros(10,1); n=1; y