北航數值分析1-jacobi法計算矩陣特征值-全文預覽
【正文】 le m[(N+1)*N/2+1] = {}。 } putchar(10)。j++) { if(fabs(muk m[(j+1)*j/2]) tempabsmin) { lumdak = m[(j+1)*j/2]。 double lumdak = 。 printf( Det(A)=%.11e\n\n,deta)。 } if(fabs(minlumda) fabs(maxlumda)) conda = fabs(maxlumda) / absminlumda。i++) { if(m[(i+1)*i/2] maxlumda) maxlumda = m[(i+1)*i/2]。 double maxlumda = pow((double),12)。 printf(結果輸出:\n\n)。 i = N 。 for(int i = 2 。 } } } MAX = max。qq)。 for(int i = 1 。 sinz = sin(atan(2 * apq / (app aqq))/2)。}//void calCosSin(double *m){ double app = m[(p+1)*p/2]。 } m[(ipr1)*ipr/2+ipc] = m[(ipr1)*ipr/2+ipc] * cosz + m[(iqr1)*iqr/2+iqc] * sinz。 } if(i q) { iqr = i。 i != q) { if(i p) { ipr = i。 for(int i = 1。 break。}elsevoid PTS(double *m){}endifvoid recounti(int i , int *pp, int *qq){ for(int j = 0 。 for(int i = 1 。 } putchar(10)。 } putchar(10)。 i = PFrols 。int kk。int q。但由于Jacobi法在迭代過程中會破壞矩陣的形態(tài)，所以原來為零的元素可能會變?yōu)榉橇悖@就導致原來的二維數組無法存儲迭代后的矩陣。所以該題可以用Jacobi法求解。數值分析 2015/11/10準備工作216。由以上分析可知，用冪法和反冪法無法完成所有問題的求解，而用Jacobi法求得矩陣所有特征值后可以求解題目中所給的各個問題。 數據結構設計由于矩陣是對稱陣，上下帶寬均為2，所以可以考慮用二維數組壓縮存儲矩陣上半帶或下半帶。//include includeincludeincludedefine N 501define V (N+1)*N/2+1define e define a(i) ( * (i)) * sin( * (i)) * pow(e , / (i
高考資料相關推薦
鄂ICP備17016276號-1