數(shù)值分析上機(jī)實(shí)驗(yàn)報(bào)告-文庫吧資料
2025-08-09 23:50本頁面
【正文】 ) Y[k][i]=Y[k][i]+(K[1][i]+2*K[2][i]+2*K[3][i]+K[4][i])/6。 for(i=1。i++) z[i]=Y[k][i]+e*K[3][i]。 for(i=1。i++) z[i]=Y[k][i]+e*K[2][i]/2。 for(i=1。i++) z[i]=Y[k][i]+e*K[2][i]/2。 for(i=1。i++) z[i]=Y[k][i]。m=m+e) {for(i=1。k++) {for(m=y[k1]。 for(k=1。j=3。i=4。i++) Y[1][i]=0。 for(i=1。 double y[5]={0,}。 y[3]=(10001000*x[2]100*x[2]100*x[3])*h。本題采用四階古典的RungeKutta公式： C語言程序代碼：includevoid fun(double x[4],double y[4],double h){y[1]=1*h。本程序采用的后置方式。，涉及函數(shù)定義，有循環(huán)，而且循環(huán)中又有判斷，編寫時(shí)需要注意該判斷條件是處于循環(huán)中，當(dāng)達(dá)到要求時(shí)跳出循環(huán)，終止運(yùn)算。[T,n]=mromb(f,1,3,) 問題討論：:把積分化為代數(shù)運(yùn)算,而實(shí)際上只需求T1(i),一般4或5次即能達(dá)到要求。T=R(J+1,J+1)。 n=2*n。 for k=1:J R(J+1,k+1)=(4^k*R(J+1,k)R(J,k))/(4^k1)。 S=S+feval(f,x)。S=0。while (erreps) J=J+1。J=0。R(1,1)=(h/2)*(feval(f,a)+feval(f,b))。} 運(yùn)行結(jié)果： MATLAB上機(jī)程序function [T,n]=mromb(f,a,b,eps)if nargin4,eps=1e6。} if(t[0][1]e) printf(t=%\n,t[0][1])。i=0。//判斷前后兩次所得的T(0)的差是否符合要求，如果符合精度要求則停止循環(huán) t[l][1]=(t[l1][1]+(ba)*s/pow(2,l1))/2。k=pow(2,l1)。fabs(t[0][j1]t[0][j2])=e。j=3。 //下為romberg算法 t[1][1]=(ba)*(f(b)+2*f((b+a)/2)+f(a))/4。 int i,j,l,k。 return(z)。 C語言程序代碼：includeincludedouble f(double x) //計(jì)算f(x)的值 {double z。（3）把區(qū)間再等分（即等分），得復(fù)化梯形公式，由與外推可得，如此，若已算出等分的復(fù)化梯形公式，則由Richardson外推法，構(gòu)造新序列， m=1,2,…,l, k=1,2,…,lm+1,最后求得。3. 編程過程中由于定義的數(shù)組比較多，需要仔細(xì)弄清楚各數(shù)組所代表的參數(shù)，要注意各下標(biāo)代表的含義，特別是在用追趕法計(jì)算的過程中。} 運(yùn)行結(jié)果： 問題討論1. 三次樣條插值效果比Lagrange插值好，沒有Runge現(xiàn)象，光滑性較好。g4=x[3]*pow(x0,2)/2+x[4]*pow(x1,2)/2(f[3]x[3]/6)+(f[4]x[4]/6)。 x1=。 printf(將x=：\n)。} printf(可得s(x)在區(qū)間[4,5]上的表達(dá)式。i) {x[i1]=(b[i1]s[i1][i]*x[i])/a[i1]。 for(i=9。}} x[9]=p10/q10。 if(i==8) {p10=b[i]。 //求d矩陣 a[i]=s[i][i]s[i1][i]*c[i]。i9。 b[0]=d[0]。}} printf(根據(jù)追趕法解三對角矩陣得：\n)。p++) //求r矩陣 { printf(%,s[i][p])。i++)
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1