《舍入误差及数值稳定实验(精)》由会员分享，可在线阅读，更多相关《舍入误差及数值稳定实验(精)（2页珍藏版）》请在金锄头文库上搜索。
1、实验报告一、实验目的： 舍入误差及数值稳定实验 ,误差在近似计算中的传播规律以及舍入误差所引起的数值不稳定性。对 n 0,1,2,3,20，计算定积分y n1xn0 xdx.5二、实验原理：根据递推方法， y(n-1)= 编写实验程序，找出计算数值稳定的计算方法，在MATLAB文件中，建立相应的 M文件，确立所需的函数表达式，通过函数调用，求出相应的实验结果。三、实验内容：建立相应的M文件，其函数名为f2 ，如下 :function y=f(n);format longy=zeros(21,1);y(1)=log(6.0)-log(5.0);for n=2:1:21y(n)=1.0/n-5.0
2、*y(n-1);end实验结果：ans =1.0e+012 *0.00000000000018-0.000000000000410.00000000000239-0.000000000011710.00000000005873-0.000000000293500.00000000146767-0.000000007338210.00000003669117-0.000000183455730.00000091727873-0.000004586393580.000022931967970.00057329919901-0.002866495994980.358311999373368.9577
3、9998433379建立相应的M文件，其函数名为f ，如下 :function y=f(n);y=zeros(3,7);y(21)=0.00873;for n=21:-1:2y(n-1)=1.0/(5.0*(n-1)-y(n)/5.0;end输入函数条件：f(21)实验结果 :ans =0.18230.04310.02430.01690.01300.01050.00890.08840.03430.02120.01540.01200.00990.00830.05800.02850.01880.01410.01120.00930.0087四、实验结果分析：对此数值方法，由于计算过程中误差不会扩大，且计算结果收敛，因此此数值方法是稳定的。五、实验体会通过做本实验，了解到了误差的重要性，并且对 MATLAB知识有了一定的了解，明白了如何定义 M文件以及定义向量形式得到结果。
