算法数值稳定性MATLAB实验报告递推对比含源码.docxVIP

实验算法的数值稳定性实验

【实验目的】

1．掌握用MATLAB语言的编程训练，初步体验算法的软件实现；

2．通过对稳定算法和不稳定算法的结果分析、比较，深入理解算法的数值稳定性及其重要性。

【实验内容】

计算积分(n=0,1,2......，10)

其中a为参数，分别对a=0.05及a=15按下列两种方案计算，列出其结果，并对其可靠性，说明原因。

方案一用递推公式 (n=1,2,......,10)

递推初值可由积分直接得

方案二用递推公式 (n=N,N-1,......,1)

根据估计式 当

或

当

取递推初值为

当

或

当

计算中取N=13开始

【解】:手工分析怎样求解这题。

【计算机求解】:怎样设计程序？流程图？变量说明？能否将某算法设计成具有形式参数的函数形式？

【程序如下】:

%myexp1_1.m---算法的数值稳定性实验

%见P11实验课题(一)

%

functionyyjj

globalna

N=20;%计算N个值

a=0.05;%或者a=15

%%--------------------------------------------

%%[方案I]用递