LAB01_算法的数值稳定性实验.docVIP

Lab01．atlab语言的编程训练，初步体验算法的软件实现； 2．通过对稳定算法和不稳定算法的结果分析、比较，深入理解算法的数值稳定性及其重要性。 【实验内容】 1．用Matlab语言编写按递推公式 计算(n=0，1，2，……)的程序，并取I0=0.6321，计算n=0，1，2，……，9时In的值。 2. 用Matlab语言编写按递推公式 计算(n=0，1，2，……，9)的值程序。 3．分析比较两种算法的数值稳定性。 【实验仪器与软件】 1．CPU主频在1GHz以上，内存在128Mb以上的PC； 2．Matlab 6.0及以上版本。 实验讲评： 实验成绩： 评阅教师： 年 月 日 Lab01．，用分部积分可得计算的递推公式 其中 0计算时，可先计算出，然后再代入（1）式，逐次求出的值。 算法二： 对于积分，由（1）式可得： （2） 计算式，可先估计一个，再反推要求的 由于 可取 (3) 在试验中，取，由（3）式计算得,然后用递推公式（2）计算的值。 二、算法程序及计算结果 算法一程序及计算结果 clear all clc I(1)=0.6321; for n=2:12 I(n)=1-(n-1)*I(n-1); end I I = Columns 1 through 8 0.6321 0.3679 0.2642 0.2074 0.1704 0.1480 0.1120 0.2160 Columns 9 through 12 -0.7280 7.5520 -74.5200 820.7200 算法二程序及计算结果 clear all clc A(10)=0.0684; for n=10:-1:2 A(n-1)=1/(n-1)*(1-A(n)); end A A = Columns 1 through 8 0.6321 0.3679 0.2642 0.2073 0.1709 0.1455 0.1268 0.1121 Columns 9 through 10 0.1035 0.0684 三、算法分析 由： 知一切，且。 按算法一，从计算的结果中看到出现负值，和的误差更大。因此，当较大时，用近似显然是不正确的。这里计算公式与每步计算都是正确的，主要是初值有误差,由此引起以后各步计算的误差满足关系，容易推得 这说明有误差，则的误差就是的倍误差，这表明这种算法是数值不稳定的。 按算法二，虽然有误差，但由算出的误差就极小。实际上，当时。 反推计算中，误差是逐渐减小的。这表明这种算法是数值稳定的。 五、总结 一个算法如果输入数据有误差，而在计算过程中舍入误差不增长，则称此算法是数值稳定的，否则称此算法为不稳定的。在上述算法中，算法二是数值稳定的，算法一是数值不稳定的。同一个问题，当某一算法数值不稳定时，可以换一种思路对其修改，找到数值稳定的算法。 n=1，2， n=9，8，