北航數(shù)值分析1-jacobi法計(jì)算矩陣特征值(編輯修改稿)
2024-09-01 03:44 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 osz + m[(iqr1)*iqr/2+iqc] * sinz。 m[(iqr1)*iqr/2+iqc] = m[(ipr1)*ipr/2+ipc] * sinz + m[(iqr1)*iqr/2+iqc] * cosz。 } } m[(p1)*p/2+q] = 0。 PTS(m)。}//void calCosSin(double *m){ double app = m[(p+1)*p/2]。 double aqq = m[(q+1)*q/2]。 double apq = m[(p1)*p/2+q]。 cosz = cos(atan(2 * apq / (app aqq))/2)。 sinz = sin(atan(2 * apq / (app aqq))/2)。}//void find_pq(double *m){ double max = 。 int pp = 0。 int qq = 0。 for(int i = 1 。 i = V 。 i++) { if(fabs(m[i]) max) { recounti(i,amp。pp,amp。qq)。 if(pp != qq) { max = fabs(m[i])。 p = pp。 q = qq。 } } } MAX = max。}void init(double *m){ for(int i = 1 。 i = N 。i++) m[(i+1)*i/2] = a(i)。 for(int i = 2 。 i = N 。 i++) m[(i1)*i/2+i1] = b。 for(int i = 3 。 i = N 。 i++) m[(i1)*i/2+i2] = c。 PTS(m)。}void calFinal(double *m){ printf(\n)。 printf(結(jié)果輸出:\n\n)。 double conda。 double deta = 。 double minlumda = pow((double),12)。 double maxlumda = pow((double),12)。 double absminlumda = pow((double),12)。 for(int i = 1 。 i =N 。i++) { if(m[(i+1)*i/2] maxlumda) maxlumda = m[(i+1)*i/2]。 if(m[(i+1)*i/2] minlumda) minlumda = m[(i+1)*i/2]。 if(fabs(m[(i+1)*i/2]) absminlum
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1