數(shù)值逼近》課程設計報告-文庫吧
2025-01-03 23:44 本頁面
【正文】 td。const int ARRAY_LEN (10000)。class Comsim{protected:double a,b,h,n,x[ARRAY_LEN],y[ARRAY_LEN]。char f[ARRAY_LEN]。public:void getab(){ cout請輸入該積分的上下限(即區(qū)間):。 cinab。 coutendl。}void cal_nh(){ int c。 cout幾點的復化拋物線公式？。 cinc。 n=2*c1。 coutn的值為：n。 h=(ba)/(2*n)。 couth的值為：h。}void cal_x(){ int i=0。 double temp=0。 for(i。i2*n+1。i++) { temp=i*h。 x[i]=a+temp。 } /*coutx的值為：endl。 for(i=0。i2*n+1。i++) { coutx[i]=x[i] 。 if(i%5==0) coutendl。 } coutendlendl。*/}void get_f(){ char temp[ARRAY_LEN]。 cout請輸入f(x)的表達式。 cintemp。 strcpy(f,temp)。}void cal_y(){ int i=0。 double temp=0。 for(i=0。i2*n+1。i++) { temp=13*(x[i]x[i]*x[i])*exp(*x[i])。 y[i]=temp。 } /*couty的值為：endl。 for(i=0。i2*n+1。i++) { couty[i]=setprecision(8)y[i] 。 if(i%5==0) coutendl。 } coutendlendl。*/}double result(){ double temp1=y[0],temp2=0,temp3=0,sum1=0,sum=0。 int i=1。 for(i=1。i2*n+1。i++) { if(i==2*n) { temp1=temp1+y[i]。/////當有判斷條件時，要先進行判斷，不滿足時才進行原始計算 } else if(i%2==0) { temp2=temp2+2*y[i]。 } else temp3=temp3+4*y[i]。 } sum1=temp1+temp2+temp3。 sum=h*sum1/3。 return sum。}void display(){ double m。 m=result()。 cout積分區(qū)間為[a b],被積函數(shù)為f的積分endl。 cout用n+1點復化拋物線公式endl計算結果為:setprecision(8)mendl。 }}。int main(){ Comsim sim。 ()。 ()。 ()。 ()。 ()。 ()。 return 0。} 五、計算結果與分析分別運用復化梯形公式和復化拋物線公式計算該積分的結果如下圖：。由上面計算結果可以看出運用復化拋物線公式的計算量比較大，但是利用復化梯形公式要將積分區(qū)間等分的個數(shù)要大于復化拋物線公式。并且由該程序對于已知結果的積分進行運算，可知復化simpson積分公式計算精度要高一些。六、參考文獻[1] 譚浩強. C語言程序設計[M]. 北京：清華大學出版社，2005.一、目的意義 (1)機械設計目的是要設計出一種能達到預定功能要求、具有性能好、成本低、價值最優(yōu)，能滿足市場需求的機械產品。 (2) 掌握牛頓插值和分段線性插值多項式的算法實現(xiàn)。 (3) 學會分析誤差和精度，熟悉運用變成語句（4）提高用數(shù)值算法解決實際問題的能力二、內容要求 機械設計問題：萬能拉撥機中有一個圓柱形凸輪(見圖1)，其底圓半徑R=30cm,凸輪的上端面不在同一平面上，要根據從動桿位移變化的需要進行設計制造。將底圓周長36等分為xi (i=0,1, …， 36)，每一圓弧段長為h=,對應于每 一分點的柱高為yi (i=0,1, …， 36)。為方便，將圓柱展開為平面，柱面的的頂端成為圖2所示的平面曲線，并已知該曲線上的37個點的坐標（表1）。 y C A
高考資料相關推薦
鄂ICP備17016276號-1