研究生數(shù)值分析課后題上機(jī)編程部分答案-文庫(kù)吧資料
2025-06-29 04:46本頁(yè)面
【正文】 0,k=1。第五題：用Romberg算法求積分近似值:依據(jù)Romberg法列出相應(yīng)的方程，先求出第一列的相應(yīng)元素，然后利用Romberg中三角法求出后面列的元素，將得到的第一行元素取差的絕對(duì)值，當(dāng)小于誤差要求的時(shí)候可以停止運(yùn)算，公式為： m= 1,2,…,l. k=1,2, …,lm+1 具體步驟為：（1） 先求出按梯形公式所得的積分值：（2） 把區(qū)間2等分，求出兩個(gè)小梯形面積之和，記為，即，由外推法得，（3） 把區(qū)間再等分，由Richardson外推法，構(gòu)造新序列，m= 1,2,…,l. k=1,2, …,lm+1（4） 或就停止計(jì)算，否則回到(3),計(jì)算（使用軟件：VC）：方法一：includeincludedefine f(x) pow(3,(x))*pow((x),)*(5*(x)+7)*sin((x)*(x))void main(){double t[9][9]。}（誤差分析、上機(jī)出現(xiàn)情況等）（1）在做除法運(yùn)算的時(shí)候，應(yīng)對(duì)除數(shù)是否等于零進(jìn)行判斷。}if(f1==0)printf(出錯(cuò)\n)。 if(f1==0)break。 f0=pow(x0,7)28*pow(x0,4)+14。fabs(x0x1)=。ⅳ. ，其中c是a,b中使達(dá)到的一個(gè).則對(duì)任意初始近似值， Newton迭代過(guò)程為，k=0，1，2…所生成的迭代序列{}平方收斂于方程在區(qū)間[a,b]上的唯一解а.對(duì)本題：所以：（使用軟件：VC）：includevoid main(){double x0,x1=,f0,f1。（2）編程的過(guò)程出現(xiàn)問(wèn)題如下：①第二個(gè)for循環(huán)i誤寫成i++，造成無(wú)限循環(huán)錯(cuò)誤②在乘法的運(yùn)算中漏掉了兩個(gè)小括號(hào)間的“”號(hào)第四題：牛頓法求方程近似根:依據(jù)牛頓迭代法進(jìn)行計(jì)算，設(shè)函數(shù)在有限區(qū)間[a,b]上二階導(dǎo)數(shù)存在，且滿足條件ⅰ. 。()近似值：%.16lf\n,s,z)。z=m[4]/2*pow(,2)+m[5]/2*pow(,2)(f[4]m[4]/6)+(f[5]m[5]/6)。i) m[i]=(g[i]x[i]*m[i+1])/p[i]。 for(i=9。 g[i+1]=d[i+1]l[i+1]*g[i]。i++) {l[i+1]=y[i]/p[i]。 for(i=1。 p[1]=2。 /*變量η*/ double x[10]={0,1,}。doublef[11]={0,0,}。 /*定義函數(shù)近似值，導(dǎo)數(shù)近似值*/ double b[11]={0,2,2,2,2,2,2,2,2,2,2},p[11],g[11],l[11]。 ，由，得出：（使用軟件：VC）：includeincludemain(){ int i。（2）對(duì)于列主元素消元法，個(gè)人感覺(jué)將a,b兩個(gè)矩陣，放在一個(gè)矩陣中進(jìn)行編程，較方便。 ，也就是j1應(yīng)該從j+1開始取值，如果從j開始取，會(huì)導(dǎo)致a[i][j1]=(a[i][j]/a[j][j])*a[j][j1]。j1=10。j1++)a[i][j1]=(a[i][j]/a[j][j])*a[j][j1]。} /*從x9到x1逐個(gè)解出結(jié)果*/}（誤差分析、上機(jī)出現(xiàn)情況等）（1）這道題目出現(xiàn)的最大問(wèn)題在 for(j1=j+1。n++) sum+=x[n]*a[m][n]。 for(n=m+1。m=1。} /*逐行消元，將矩陣轉(zhuǎn)化為上三角矩陣*/ printf(x9為：%lf\n,x[9]=a[9][10]/a[9][9])。j1=10。i=9。h++) a[g][h]=y[h]。 for(h=j。h=10。h++) y[h]=a[j][h]。} /*確定本列絕對(duì)值最大元素所在行：g行*/
化學(xué)相關(guān)推薦
鄂ICP備17016276號(hào)-1