數(shù)值分析第二次大作業(yè)-全文預(yù)覽
【正文】 k。 for(k = 0。j10。}//構(gòu)造函數(shù)translation 用于進(jìn)行原點平移void translation (double x) { int i=0,j=0。 } } } //sp 5 } }double min(double v1,double v2){ return((v2v1)?v1:v2)。 i++) { for(j = 0。 i = m。 i = m。 pr[i] = tmp_pr / hr。 j++) { tmp_vr = tmp_vr + mk[j][i] * ur[j]。 tmp_pr = 0。 ur[r] = mk[r][r]cr。 i++) ur[i] = 0。 cr = * sgn(mk[r][r]) * dr。 for(i = r。i=m。for (r=0。 i++) mk[i][i]+= mt。 k++) tmp_sum+= A[i][k] * A[k][j]。 j = m。 double tmp_sum。 qmk(mk,m,ms,mt)。 else continue。 m。 m。amp。 tzr[0]= s1_r。 s2_r = ( *fb sqrt(temp))/ 2。 s2_r = *fb。 temp = fb * fb * detdk。 } else if(m==1) return 0。 flag=1。amp。 //step5 double detdk,fb,temp,s1_r,s1_i,s2_r,s2_i。 else return 1。 j++) if (A[i][j] amp。 } } } } //以上的A[][]存在一定的誤差，消除誤差 for(i=0。i 10。 for(i= 0。 for(i= 0。 sum= sum+A[i][j]* u[j]。 sum= 0。 u[r+1]= A[r+1][r] c。i++) u[i]= 0。 c= 1*sgn(A[r+1][r])*d。 for (i=r+1。i10。 //開始循環(huán) for (r=0。}//定義函數(shù)對矩陣A進(jìn)行擬上三角化void nssj(){ int r,i,j。j++) A[i1][j1]=sin(*i+*j)。 for(i = 1。根據(jù)數(shù)值分析課本的相關(guān)知識，步驟3中帶雙步位移的QR分解法的流程圖如下：本次作業(yè)所使用的編譯環(huán)境為：Visual c++ 三、 程序源代碼includeincludedouble A[10][10],rr[10][10],qq[10][10],rq[10][10],uk[10][10]。已知： (i,j=1,2,……,10)二、 算法設(shè)計1. 按照題目給出的矩陣定義對矩陣A賦初值：對應(yīng)的函數(shù)為initA()；2. 為了減少求特征值和特征向量過程中的計算量，在對矩陣進(jìn)行QR分解前先進(jìn)行擬上三角化：對應(yīng)的函數(shù)為nssj()； 3. 對擬上三角化后的矩陣A使用帶雙步位移的QR分解法逐次迭代（最大迭代次數(shù)L=500），逐個求出其特征值，對應(yīng)的函數(shù)為tezhengzhi()；這個過程中包含著兩個子程序：QR()和qmk()，分別用來對矩陣Mk進(jìn)行QR分解并得到Ak+1和計算mk的值；4. 使用帶原點平移的反冪法求出其對應(yīng)的特征向量，對應(yīng)的函數(shù)為：fmifa()，這個過程中求解線性方程時用到列主元的高斯消元法，對應(yīng)的函數(shù)為：gauss(double )。 //定義矩陣的特征值數(shù)組,r實部、i虛部//對矩陣A進(jìn)行初始化，賦值void initA(){ int i,j。 j= 10。i++) A[i 1][i 1]=cos(i+*i)*。 double sum=0,tmp=0。 for(i= r+2。 else
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1