数值分析之MATLAB实验一误差分析

《数值分析之MATLAB实验一误差分析》由会员分享，可在线阅读，更多相关《数值分析之MATLAB实验一误差分析（2页珍藏版）》请在装配图网上搜索。

1、实验 第一次：误差分析对，按照下面两种算法计算定积分.算法1：利用递推公式，取.算法2：利用递推公式注意到,取.思考：从计算结果看，哪个算法是不稳定的，哪个算法是稳定的。法一程序： t=log(6.0)-log(5.0); n=0;y=zeros(1,21); y(1)=t;for k=2:21 y(k)=1/k-5*y(k-1); n=n+1; endy(1:6)y(7:21)法二程序： clear allclcy=zeros(21,1);n=1;y1=(1/105+1/126)/2;y(21)=y1;for k=21:-1:2 y(k-1)=1/(5*(k-1)-y(k)/5; n=n+1

2、;endy法一结果： ans = 0.1823 -0.4116 2.3914 -11.7069 58.7343 -293.5049ans = 1.0e+012 * Columns 1 through 6 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 Columns 7 through 12 0.0000 -0.0001 0.0006 -0.0029 0.0143 -0.0717 Columns 13 through 15 0.3583 -1.7916 8.9578法二结果： 10.0884110.014120.058120.01330.0431130.01240.0343140.011250.0285150.010560.0243160.009970.0212170.009380.0188180.008990.0169190.0083100.0154200.0087