數(shù)值代數(shù)上機報告-文庫吧資料
2024-08-03 10:44本頁面
【正文】 =n。 l[i][k]=(R[i][k]q)/r[k][k]。pk。i=n。 r[k][i]=R[k][i]q。pk。i=n。k=n。i++) y0[i]=1。 /*第一步迭代，用Doolittle計算R*y1=y0中的y1*/ for(i=1。ji。i=n。 else L[i][j]=0。j=n。i=n。 L[i][k]=(B[i][k]q)/R[k][k]。pk。i=n。 R[k][i]=B[k][i]q。pk。i=n。k=n。 else B[i][j]=A[i][j]。j=n。i=n。lm)。 printf(\n請輸入近似特征值：)。j++) scanf(%lf,amp。i++) for(j=1。 for(i=1。n)。 printf(請輸入矩陣的階數(shù)：n=)。 double A[N][N],B[N][N],L[N][N],R[N][N],l[N][N],r[N][N]。三、 源程序：includeincludedefine e define N 100void main(){ int i,j,k,n,p,s,kk=1,flag=1。 } return 0。 //return 0。i++) printf(%.4f\n,nx[i])。 for(i=0。 if(fabs(xmaxoxmax)epsilon) { printf(按模最大特征值=%.4f\n,xmax)。jn。j++) nx[j]/=xmax。 for(j=0。jn。 } xmax=。kn。j++) { nx[j]=0。i++) { for(j=0。 for(i=0。i++) x[i]=1。 for(i=0。j++) scanf(%lf,amp。i++) for(j=0。 for(i=0。 return 1。 return 1。n)。 printf(輸入矩陣維數(shù)n: )。 static double a[N][N]。 int i,j,k。二、 算法：1) 輸入系數(shù)矩陣 A；2) 取出事向量為；3) 利用乘冪法的迭代格式 計算得出原矩陣按模最大特征值和相應(yīng)特征的向量。 getch()。in。 printf(\n原線性方程組的解為：\n)。 } x=(float *)malloc(3*sizeof(float))。j++) scanf(%f,amp。i++) { for(j=0。 for(i=0。n)。 printf(請輸入系數(shù)矩陣的階數(shù)：n=)。 float a[N][N]。i++) x[i]=y[i]。 } for(i=0。 } if(k=MAX) { printf(The Method is disconvergent)。 } if(epsilon EPS) { printf(\n迭代次數(shù)為：%d次\n,k)。 } y[i]=(a[i][n]s)/a[i][i]。j++) { if(j==i) continue。 for(j=0。in。 k++。x[i])。in。 printf(請輸入迭代的初始向量：\n)。 x=(float *)malloc(n*sizeof(float))。 double epsilon。 }}四、 計算結(jié)果與分析：分析：運行結(jié)果與預(yù)想的結(jié)果相近，誤差對結(jié)果的影響不是很大，比較理想五、 參考文獻：[1]刑志棟. 矩陣數(shù)值分析. 陜西： 陜西科學(xué)技術(shù)出版社， 2005[2]譚浩強. C語言程序設(shè)計. 北京：清華大學(xué)出版社，2005Jacobi迭代法報告3一、 目的意義：通過有限次的算術(shù)運算得到線性方程組的近似值二、 算法：1) 輸入系數(shù)矩陣 A；2) 輸入迭代的初始向量；3) 利用公式 ，計算得出原線性方程組的近似解。 j++) { sum = a[i][j]*x[j]。 for (int j=i+1。 i=0。 } }}void huidai(){ x[n1] = b[n1]/a[n1][n1]。 jn。 b[i] = tmp*b[row]。 i++) { int j=row。 }}void xiaoyuan(int row){ for (int i=row+1。 a[row1][j] = a[row2][j]。 jn。 b[row2] = tmp。 tmp = b[row1]。 } return rowmark。 rowmark = i。 in。 int rowmark = row。 system(pause)。 printf(\n)。 in。 printf(\n線性方程組的解為:\n )。 } xiaoyuan(i)。 return。 if (rowmark == 1) { printf(多解!)。 in1。 } getchar()。 i++) { scanf(%lf, amp。 for (i=0。 getchar()。 j++) scanf(%lf, amp。 i++) { for (int j=0。 for (int i=0。 getchar()。 scanf(%d, amp。 /*消元*/void huidai()。 /* 選主元*/void hangjiaohuan(int row1, int row2)。int
研究報告相關(guān)推薦
鄂ICP備17016276號-1