資源描述：
《計算方法源程序》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、計算方法實驗報告實驗一舍入誤差與數(shù)值穩(wěn)定性班級：計算機cO96姓名：郭睿學號:088115實驗一舍入誤差與數(shù)值穩(wěn)定性目的與要求：1、通過上機編程，復(fù)習鞏固以前所學程序設(shè)計語言；2、通過上機計算，了解舍入誤差所引起的數(shù)值不穩(wěn)定性。3、通過上機計算，了解運算次序?qū)τ嬎憬Y(jié)果的影響，從而盡量避免大數(shù)吃小數(shù)的現(xiàn)象。實驗內(nèi)容：?通過止反兩個實例的計算，了解利用計算機進行數(shù)值計算中舍入誤并所引起的數(shù)值不穩(wěn)定性，深入理解初始小的舍入誤差可能造成誤差積累從而對計算結(jié)果的巨大影響。?通過實際編程，了解運算次序?qū)τ嬎憬Y(jié)果的影響，了解實數(shù)運算符合的結(jié)合律和分配律在計算機里不一定
2、成立。1對n=0,1,2,...,20計算定積分算法1利用遞推公式nln6-ln5?0.182322算法2利用遞推公式人n=20,19,...,l注意到110511112011nr訂x%*占注和x%算法一：#include#include#include#includeHstdlib.hHvoidmain(void){floaty_0=log(6.0)-log(5.0),y_l;intn=l;system(”cis”);/*清屏*/printf(”y[0]=%?20f”,y_0);while(l){y_
3、l=1.0/n-5*y_0;printf(Ky[%d]=%-20f,y_l);/*輸出*/if(n>=20)break;y_0=y_l;n++;if(n%3==0)printf(Hn);}getch();/*保持用戶屏幕*/運行結(jié)杲:Y[0]=0.182322Y[l]=0.088392Y⑵=0.058039Y[3]=0.043128Y[41=0.034310Y[5]=0.028448Y[6]=0.024428Y[7]=0.020719Y[8]=0.021407Y[9]=0.004076Y[10]=0.079618Y[ll]=-0.307181Y[
4、12]=1.619237Y[13]=-8.019263Y[14]=40.167744Y[15]=-200.772069Y[16]=1003.922729Y[17]=-5019.554688Y[18)=25097.828125Y[19]=-125489.085938Y(20
5、=627445.500000算法二#include#include#include#include"stdlib.hHmain()floaty_0=(l/105.0+l/l26.0)/2,y_l;intn=20;system(ncls
6、");/*清屏*/printf(Hy[201=%-20fy_0);while(l){y_l=l/(5.0*n)-y_0/5.0;printf(Ky[%d]=%-20f-1,y_l);if(n<=l)break;y_O=y_l;n-;if(n%3==0)printf(nn);}getch();運行結(jié)杲:Y[20]=0.008730Y[19]=0.008254Y[18]=0.008876Y[17]=0.009336Y[16]=0.009898Y[151=0.010520Y[14]=0.011229Y[13]=0.012040Y[12]=0.0129
7、77Y[ll]=0.014071Y[10]=0.015368Y[9]=0.016926Y⑻=0.018837Y[7]=0.021233Y[6]=0.024325Y[5]=0.028468Y[4]=0.034306Y[3]=0.043139Y(2]=0.058039Y(l]=0.088392Y[0
8、=0.182322說明：從計算結(jié)果可以看出，算法1是不穩(wěn)定的，而算法2是穩(wěn)定的。計算方法實驗報告實驗二拉格朗日插值與牛頓插值班級：計算機c096姓名:學號：088115實驗二拉格朗日插值與牛頓插值一、目的與要求：熟悉拉榕朗日插值多項式和牛頓插值多項式，注意其不同
9、特點；二、實驗內(nèi)容：通過拉格朗口插值和牛頓插值多項式的兩個實例的計算，了解兩種求解方法,分析他們的優(yōu)缺點。三、程序與實例算法1.輸入X,.,y;.(i=0,l,2,...,n),令L(xM)=0;2.對=0,1,2,.,?n計算?1/W=nn/=0j右X-Xjxi-XjL,kL“+l,(x)y,程序與實例例1已知函數(shù)農(nóng)Xj0.561600.562800.564010.56521y,-0.827410.826590.825770.82495用三次拉格朗日多項式求x=0.5635的函數(shù)近似值。拉格朗日插值：#include#include<
10、conio.h>#include#include"stdli
