数值分析第一章数值试验.docx

数值分析实验报告(一)一、实验名称误差传播与算法稳定性二、实验目的1．理解数值计算稳定性的概念。2．了解数值计算方法的必要性。3．体会数值计算的收敛性与收敛速度。三、实验内容计算，四、算法描述由，知 则 可得递推关系 1． ， 2． ， 下面分别以1，2递推关系求解方案1 ， 当时，, 递推公式为（1）方案2 ， 当时， 则 即取递推初值 递推公式为（2）取递推公式（1）中的初值 ，得取递推公式（2）中的初值 ，得五、程序流程图由于实验方案明显、简单，实现步骤及流程图省略。六、实验结果n0123456789100.4023594781085250.1494101304728690.0876474673817830.0614214664878880.0471446333780280.0382138416554930.0321132062527930.0276859841510870.0243285039622280.0216956517872210.0195760870531952.113096392496391e+0035.280240981240977e+0021.321310245310244e+002-32.9494227994227708.299855699855693-2.0249639249639230.547907647907647-0.1012626262626260.0565656565656570.0136363636363640.021590909090909七、实验结果分析由递推公式（1）知当时，应当为精确解，递推公式的每一步都没有误差的取舍，且由 得，误差随递推公式逐步缩小。由递推公式（2）知 ，为估计值，并不精确，且计算结果< ，出现负值。由此看出，当减小时，用递推公式（2）中的近似是不正确的。主要原因是初值不是精确值。设有误差，由递推公式（1）知，则有。误差随的减小而迅速增加，增加到的倍。由此可见，递推公式计算的误差不仅取决于初值的误差，公式的精确性，还依赖于误差的传递即递推计算的稳定性。综上所述，在递推计算中，数值计算方法是非常重要的，误差估计、误差传播及递推计算的稳定性都会直接影响递推结果。八、附录(程序)方案1syms I1 I2I1=(1/4)*log(5)for n=1:10 I2=1/(4*n)-(1/4)*I1 I1=I2;end I1 = 0.402359478108525I2 = 0.149410130472869I2 = 0.087647467381783I2 = 0.061421466487888I2 = 0.047144633378028I2 = 0.038213841655493I2 = 0.032113206252793I2 = 0.027685984151087I2 = 0.024328503962228I2 = 0.021695651787221I2 = 0.019576087053195方案2syms I1 I2I2=19/880for n=10:-1:1I1=1/n-4*I2I2=I1;end I2 = 0.021590909090909I1 = 0.013636363636364I1 = 0.056565656565657I1 = -0.101262626262626I1 = 0.547907647907647I1 = -2.024963924963923I1 = 8.299855699855693I1 = -32.949422799422770I1 = 1.321310245310244e+002I1 = -5.280240981240977e+002I1 = 2.113096392496391e+003