數(shù)值逼近課程設(shè)計(jì)報(bào)告-wenkub
2023-04-07 03:20:24 本頁面
【正文】 m。in+1。 } coutendlendl。 } /*couty的值為：endl。 for(i=0。 strcpy(f,temp)。 } coutendlendl。 } /*coutx的值為：endl。 for(i。 h=(ba)/n。 cinc。 cinab。const int ARRAY_LEN (10000)。數(shù)值逼近課程設(shè)計(jì)數(shù)值逼近課程設(shè)計(jì)報(bào)告 21 一、目的意義 (1)進(jìn)一步熟悉掌握復(fù)化梯形和復(fù)化拋物線公式(2)學(xué)會比較復(fù)化梯形公式和復(fù)化拋物線公式如何達(dá)到所要求的精度(3) 提高編程能力(4)通過數(shù)值方法求出很難求得原函數(shù)的積分和解析表達(dá)是沒有明確的給出積分的近似值二、內(nèi)容要求 積分計(jì)算問題：分別用復(fù)化梯形和復(fù)化Simpsom求積公式計(jì)算積分，并比較計(jì)算量（精度為108）。class Comt{protected:double a,b,h,n,x[ARRAY_LEN],y[ARRAY_LEN]。 coutendl。 coutendl。 couth的值為：hendl。in+1。 for(i=0。*/}void get_f(){ char temp[ARRAY_LEN]。 coutendl。in+1。 for(i=0。*/}double result(){ double temp=y[0],sum。i++) { if(i==n) temp=temp+y[i]。}void display(){ double m。 }}。 ()。 return 0。char f[ARRAY_LEN]。}void cal_nh(){ int c。 coutn的值為：n。 double temp=0。 x[i]=a+temp。i++) { coutx[i]=x[i] 。 cout請輸入f(x)的表達(dá)式。 double temp=0。 y[i]=temp。i++) { couty[i]=setprecision(8)y[i] 。 int i=1。/////當(dāng)有判斷條件時(shí)，要先進(jìn)行判斷，不滿足時(shí)才進(jìn)行原始計(jì)算 } else if(i%2==0) { temp2=temp2+2*y[i]。 return sum。 cout用n+1點(diǎn)復(fù)化拋物線公式endl計(jì)算結(jié)果為:setprecision(8)mendl。 ()。 ()。并且由該程序?qū)τ谝阎Y(jié)果的積分進(jìn)行運(yùn)算，可知復(fù)化simpson積分公式計(jì)算精度要高一些。將底圓周長36等分為xi (i=0,1, …， 36)，每一圓弧段長為h=,對應(yīng)于每 一分點(diǎn)的柱高為yi (i=0,1, …， 36)。三、問題解決的方法與算法方法：Newton差值和分段線性插值 Newton插值設(shè)y=f(x)在節(jié)點(diǎn)處的函數(shù)值已知，則關(guān)于點(diǎn)的一階差商為：則關(guān)于點(diǎn)的二階差商為：則關(guān)于點(diǎn)的k階差商為：一階差商二階差商三階差商X0f(X0)X1f(X1)f[X0, X1]X2f(X2)f[X1, X2]f[X0,X1, X2]X3f(X3)f[X2, X3]f[X1, X2,X3]f[X0,X1, X2 X3] 插值公式 插值余項(xiàng) 插值函數(shù)設(shè)y=f(x)在節(jié)點(diǎn)處的函數(shù)值為為了提高近似程度，可以考慮用分段線性插值來逼近原函數(shù)，這時(shí)的插值函數(shù)為分段函數(shù)：在區(qū)間上的線性函數(shù)為： 插值公式 插值余項(xiàng)第一步：輸入要進(jìn)行計(jì)算的x的值第二步：輸入各組y[i]及x[i]的值，i=0,1,2,…,n第三步：判斷條件：if(x=0amp。x=)則用分段線性插值進(jìn)行計(jì)算，如果兩種情況均不是則輸出輸入錯(cuò)誤第四步：計(jì)算出柱高S第五步：輸出S四、計(jì)算程序// : 定義控制臺應(yīng)用程序的入口點(diǎn)。 double h=。 for(i=0。 } printf(請輸出所有x的值：\n)。 if(i%3==0) printf(\n)。jARRAY_LEN。 ///////////////////求差商 //printf(各階插商值為：\n)。ji。jARRAY_LEN。i++) { mul=mul*(ax[i])。 int j=0。 printf(請輸出所有x的值：\n)。 if(i%2==0) printf(\n)。j2。i2。 double N=0。 printf(請輸入所要求的x的值
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1