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

AAAAAA实验一用递推公式计算定积分09信息符文飞09211210071、实验目的：由于一个算法是否稳定，十分重要。如果算法不稳定，则数值计算的结果就会严重背离数学模型的真实结果，因此，在选择数值计算公式来进行近似计算时，我们应特别注意选用那些在数值计算过程中不会导致误差迅速增长的公式。体会稳定性在选择算法中的地位．误差扩张的算法是不稳定的，是我们所不期望的；误差衰竭的算法是稳定的．是我们努力寻求的，这是贯穿本课程的目标．通过上机计算，了解舍入误差所引起的数值不稳定性。2、实验题目：对n=0,1,2,…,20,计算定积分dxxxynn1053、实验原理由于y(n)==–在计算时有两种迭代方法，如下：方法一：y(n)=–5*y(n-1),n=1,2,3,…,20;取y(0)==ln6-ln5≈0.182322方法二：AAAAAA利用递推公式：y(n-1)=-*y(n),n=20,19,…,1.而且，由=*≤≤*=可取：y(20)≈*()≈0.008730.4、实验内容：算法1的程序：y0=log(6.0)-log(5.0);y1=0;n=1;whilen=30y1=1/n-5*y0;fprintf('y[%d]=%-20f',n,y1);y0=y1;n=n+1;ifmod(n,1)==0;fprintf('\n')endend算法2的程序：y0=(1/105+1/126)/2;AAAAAAy1=0;n=1;whilen=30y1=1/(5*n)-y0/5;fprintf('y[%d]=%-20f',n,y1)y0=y1;n=n+1;ifmod(n,1)==0fprintf('\n')endend5、实验结果对于算法1：y[1]=0.088392y[2]=0.058039y[3]=0.043139y[4]=0.034306y[5]=0.028468y[6]=0.024325y[7]=0.021233y[8]=0.018837AAAAAAy[9]=0.016926y[10]=0.015368y[11]=0.014071y[12]=0.012977y[13]=0.012040y[14]=0.011229y[15]=0.010519y[16]=0.009904y[17]=0.009304y[18]=0.009035y[19]=0.007457y[20]=0.012713y[21]=-0.015946y[22]=0.125183y[23]=-0.582439y[24]=2.953862y[25]=-14.729311y[26]=73.685015y[27]=-368.388036y[28]=1841.975892y[29]=-9209.844979AAAAAAy[30]=46049.258229对于算法2：y[1]=0.198254y[2]=0.060349y[3]=0.054597y[4]=0.039081y[5]=0.032184y[6]=0.026897y[7]=0.023192y[8]=0.020362y[9]=0.018150y[10]=0.016370y[11]=0.014908y[12]=0.013685y[13]=0.012648y[14]=0.011756y[15]=0.010982y[16]=0.010304y[17]=0.009704y[18]=0.009170y[19]=0.008692AAAAAAy[20]=0.008262y[21]=0.007871y[22]=0.007517y[23]=0.007192y[24]=0.006895y[25]=0.006621y[26]=0.006368y[27]=0.006134y[28]=0.005916y[29]=0.005713y[30]=0.0055246、实验结果分析：由实验结果可以看到，算法1在计算过程中误差会增长，所以算法1不稳定。算法2在计算过程中误差逐渐减小，所以算法2稳定。为了不影响数值计算结果的精确度与真实性，在实际应用中，我们应选用数值稳定的算法2，尽量避免使用数值不稳定的算法1.兰亭序永和九年，岁在癸丑，暮春之初，会于会稽山阴之兰亭，修AAAAAA禊事也。群贤毕至，少长咸集。此地有崇山峻岭，茂林修竹；又有清流激湍，映带左右，引以为流觞曲水，列坐其次。虽无丝竹管弦之盛，一觞一咏，亦足以畅叙幽情。是日也，天朗气清，惠风和畅，仰观宇宙之大，俯察品类之盛，所以游目骋怀，足以极视听之娱，信可乐也。夫人之相与，俯仰一世，或取诸怀抱，晤言一室之内；或因寄所托，放浪形骸之外。虽取舍万殊，静躁不同，当其欣于所遇，暂得于己，快然自足，不知老之将至。及其所之既倦，情随事迁，感慨系之矣。向之所欣，俯仰之间，已为陈迹，犹不能不以之兴怀。况修短随化，终期于尽。古人云：“死生亦大矣。”岂不痛哉！每览昔人兴感之由，若合一契，未尝不临文嗟悼，不能喻之于怀。固知一死生为虚诞，齐彭殇为妄作。后之视今，亦犹今之视昔。悲夫！故列叙时人，录其所述，虽世殊事异，所以兴怀，其致一也。后之览者，亦将有感于斯文。