实验一-用递推公式计算定积分

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业实验一 用递推公式计算定积分09信息 符文飞 1、实验目的：由于一个算法是否稳定，十分重要。如果算法不稳定，则数值计算的结果就会严重背离数学模型的真实结果，因此，在选择数值计算公式来进行近似计算时，我们应特别注意选用那些在数值计算过程中不会导致误差迅速增长的公式。体会稳定性在选择算法中的地位误差扩张的算法是不稳定的，是我们所不期望的；误差衰竭的算法是稳定的是我们努力寻求的，这是贯穿本课程的目标通过上机计算，了解舍入误差所引起的数值不稳定性。2、实验题目： 对=0,1,2

2、,20,计算定积分3、实验原理由于y(n)=01xnx+5 = 在计算时有两种迭代方法，如下：方法一： y(n)= 1 n 5*y(n-1),n=1,2,3,取y(0)= 011方法二：利用递推公式：y(n-1)= 1 5n-1 5而且，由1 126 = 1 6* 01x20dx01可取：y(20)1 2*(1 105+4、实验内容：算法1的程序： y0=log(6.0)-log(5.0);y1=0;n=1;while n=30 y1=1/n-5*y0; fprintf(y%d=%-20f,n,y1); y0=y1; n=n+1; if mod(n,1)=0; fprintf(n) enden

3、d算法2的程序：y0=(1/105+1/126)/2;y1=0;n=1;while n=30 y1=1/(5*n)-y0/5; fprintf(y%d=%-20f,n,y1) y0=y1; n=n+1; if mod(n,1)=0 fprintf(n) endend5、实验结果对于算法1：y1=0. y2=0. y3=0. y4=0. y5=0. y6=0. y7=0. y8=0. y9=0. y10=0. y11=0. y12=0. y13=0. y14=0. y15=0. y16=0. y17=0. y18=0. y19=0. y20=0. y21=-0. y22=0. y23=-0. y

4、24=2. y25=-14. y26=73. y27=-368. y28=1841. y29=-9209. y30=46049. 对于算法2：y1=0. y2=0. y3=0. y4=0. y5=0. y6=0. y7=0. y8=0. y9=0. y10=0. y11=0. y12=0. y13=0. y14=0. y15=0. y16=0. y17=0. y18=0. y19=0. y20=0. y21=0. y22=0. y23=0. y24=0. y25=0. y26=0. y27=0. y28=0. y29=0. y30=0. 6、实验结果分析：由实验结果可以看到，算法1在计算过程中

5、误差会增长，所以算法1不稳定。算法2在计算过程中误差逐渐减小，所以算法2稳定。为了不影响数值计算结果的精确度与真实性，在实际应用中，我们应选用数值稳定的算法2，尽量避免使用数值不稳定的算法1. 兰亭序永和九年，岁在癸丑，暮春之初，会于会稽山阴之兰亭，修禊事也。群贤毕至，少长咸集。此地有崇山峻岭，茂林修竹；又有清流激湍，映带左右，引以为流觞曲水，列坐其次。虽无丝竹管弦之盛，一觞一咏，亦足以畅叙幽情。是日也，天朗气清，惠风和畅，仰观宇宙之大，俯察品类之盛，所以游目骋怀，足以极视听之娱，信可乐也。夫人之相与，俯仰一世，或取诸怀抱，晤言一室之内；或因寄所托，放浪形骸之外。虽取舍万殊，静躁不同，当其欣于所遇，暂得于己，快然自足，不知老之将至。及其所之既倦，情随事迁，感慨系之矣。向之所欣，俯仰之间，已为陈迹，犹不能不以之兴怀。况修短随化，终期于尽。古人云：“死生亦大矣。”岂不痛哉！每览昔人兴感之由，若合一契，未尝