數(shù)值分析上機(jī)題參考答案-wenkub
2023-06-22 19:16:31 本頁面
【正文】 39。i++)coutx1[i]39。w=w39。}l++。}x1[i]=(1w)*x0[i]+w*(b[i]sum)/a[i][i]。j++){ sum=sum+a[i][j]*x1[j]。i9。 do { for(i=0。t++) { l=0。 float b[9]={15,27,23,0,20,12,7,7,10}。 float w,t。} //輸出解向量}結(jié)果：方程的解為：x1= ，x2= ，x3= ，x4= ，x5= ，x6= ，x7= ，x8= ，x9= 。x39。x[i]=(a[i][n+1]sum3)/a[i][i]。float sum3=0。//float x[10]。39。39。39。39。39。39。39。39。39。39。39。39。39。\n39。j=(n+1)。a[i][k]=(a[i][k]sum2)/(a[k][k])。i++)//算出第k列列元素的值{float sum2=0。q++)sum1+=a[k][q]*a[q][j]。j=n+1。j++) //交換兩行元素{float s2=a[k][j]。i++)//比較第k步分解的第k列值的大小{if(fabs(s1[i])m){m=fabs(s1[i])。}int l=k。for(q=1。k=n。39。39。39。39。39。39。39。39。39。39。39。39。39。39。39。39。39。\t39。i=n。} //進(jìn)行第一列主元互換for(i=2。j=(n+1)。i=n。i=n。39。39。39。39。39。39。\n39。39。39。39。39。39。39。cout39。i++){for(j=1。j=(n+1)。cout輸入數(shù)組a:endl。（1）編制解n階線性方程組Ax=b的列主元三角分解法的通用程序；（2）用所編制的程序解線性方程組RI=V，并打印出解向量，保留五位有效數(shù)；（3）本編程之中，你提高了哪些編程能力？程序?yàn)椋篿ncludeincludevoid main(void){int i,j,n,k,q。分別為方程的精確解。 } cout輸出方程根收斂的區(qū)間值:\n。 while(flag==1) { cout方程的根為:39。 delay(30000)。type1=type139。 //調(diào)試值用coutcount=count39。 x0=z。 //f(x)導(dǎo)函數(shù)的表達(dá)式。 //f(x)的表達(dá)式。n0。 coutk39。 x1=x0+a。 int k=0。 df=x*x1。 f=x*x*x/3x。1．由Newton方法的局部收斂性可知存在，當(dāng)時(shí)，Newton迭代序列收斂于根。從大到小的順序計(jì)算得到的結(jié)果的有效位數(shù)少。j){s=1/(j*j1)。 } return sum。 for(j=2。（1）編制按從大到小的順序，計(jì)算的通用程序。（2）編制按從小到大的順序，計(jì)算的通用程序。j=N。} 按從小到大的順序計(jì)算的通用程序?yàn)椋篿ncludefloat sum(float N){float j,s,sum=0。sum+=s。計(jì)算機(jī)在進(jìn)行數(shù)值計(jì)算時(shí)會(huì)出現(xiàn)“大數(shù)吃小數(shù)”的現(xiàn)象，導(dǎo)致計(jì)算結(jié)果的精度有所降低，我們?cè)谟?jì)算機(jī)中進(jìn)行同號(hào)數(shù)的加法時(shí)，采用絕對(duì)值較小者先加的算法，其結(jié)果的相對(duì)誤差較小。試確定盡可能大的。 //f(x)的表達(dá)式。 //f(x)導(dǎo)函數(shù)的表達(dá)式。 cout請(qǐng)輸入初值x0:。 k++。\t39。n)。 return(f)。 return (df)。 while(count0) { x1=x0f(x0)/df(x0)。\t39。\n39。 //調(diào)試值用 } return 0。\n39。 coutdeltaeps。分析：對(duì)于不同的初值，迭代序列會(huì)收斂于不同的根，所以在某個(gè)區(qū)間內(nèi)求根對(duì)于初值的選取有很大的關(guān)系。float a[10][11],s[10],s1[10]。for(i=1。j++)cina[i][j]。j=(n+1)。\n39。39。39。39。39。39。39。//進(jìn)行第一行和第一列元素的求取39。39。39。39。39。39。39。i++){s[i]=a[i][1]。i++)if(fabs(s[i])max){max=fabs(s[i])。j++){float b=a[1][j]。i=n。i++){for(j=1。cout39。39。39。39。39。39。39。\n39。39。39。39。39。39。39。39。39。39。39。k++){for(i=k。qk。float m=fabs(s1[k])。l=i。a[k][j]=a[l][j]。j++)//算出第k行行元素的值{float sum1=0。a[k][j]=a[k][j]sum1。for(q=1。}} //第k步分解結(jié)束for(i=1。j++)couta[i][j]39。} //輸出改變后的數(shù)組//輸出解39。39。39。39。39。39。39。39。39。39。39。39。39。39。for(i=n1。for(j=i+1。} //回代過程for(i=1。i39。分析：我感覺是提高了查錯(cuò)誤點(diǎn)的能力和編循環(huán)語句的能力，即利用很規(guī)整的迭代公式進(jìn)行編程。 float m[9]。 float max(float m[9])。 float x0[9]={1,1,1,1,1,1,1,1,1}。i9。i++){sum=0。}for(j=i+1。} //解出九個(gè)解for(i=0。}while(max(m)=eps)。\n39。=39。\n39。i9。不是合理的迭代值。用所編制程序求車門的3次樣條插值函數(shù)S(x),并打印出S(i+)(i=0,1,…9)。//存放數(shù)組h[j] float u[11]。//存放數(shù)組M[j] float b[11]。 float xx。 for(i=0。 for(i=0。 cout輸入兩個(gè)端點(diǎn)值:\n。 //利用書本上的算法求出所需要的值 //求出h[j]的值 for(j=0。h39。]39。\t39。 for(j=1。 } //求出d[j]的值 for(j=1。 d[n]=6*(df[1](y[n]y[n1])/h[n1])/h[n1]。u39。]39。\t39。j++) { cout39。j39。v[j]39。j=n。[39。=39。 //利用書本上的追趕法求解方程組 for(i=0。 t[0]=b[0]。i++) { l[i]=u[i]/t[i1]。 for(i=n1。i=n。[39。=39。j++) { s[0]=y[j]。 cout當(dāng)x的值在區(qū)間39。j39。39。時(shí),輸出插值多項(xiàng)式的系數(shù):\n。s39。]39。\t39。應(yīng)用復(fù)化Simpson公式，采用逐次二分步長(zhǎng)的方法，編制計(jì)算I(f)的通用程序。程序如下：includeincludedefine PI define error double f(double x,double y){ double answ。 int i,j,m=1,n=1。// double k=(dc)/(2*m)。 answ+=f(a,c)。i++) temp+=
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1