北航數(shù)值分析大作業(yè)一(參考版)
2025-06-29 08:46本頁面
【正文】 。如果按摸第二大的特征值λ2對(duì)應(yīng)的特征向量x2的系數(shù)α2不為0，則求出該特征值。i取不同值時(shí)，即取不同的初始向量時(shí)，可能得到不同的特征值。以冪法為例（反冪法原理相同），選取初始迭代向量ui=ei(i=0,1，…500)；即u[j]=0,j≠i；u[j]=1,j=i。 getchar()。 if(((i+1)%3)==0) printf(\n)。 d1=mifa(u,MatrixC,0)。i=500。 assignment(MatrixC)。//測試不同迭代初始向量對(duì)λ1計(jì)算結(jié)果的影響。i++) { a*=MatrixC[2][i]。 for(i=0。 LU(MatrixC)。 printf(A的（譜范數(shù)）條件數(shù)condA2=%.12e\n,d)。 printf(與μ%02d=%+.12e最接近的特征值λi%02d=%+.12e\n,i,a,i,d)。 chuzhi(u)。i++) { a=d1+(i*(d501d1))/40。//用帶原點(diǎn)平移的反冪法求λik for(i=1。 ds=fmifa(u,MatrixC,0)。 printf(λ501=%.12e\n,d501)。 d1=a。 d501=a。 a=mifa(u,MatrixC,d)。 d=mifa(u,MatrixC,0)。//將矩陣A轉(zhuǎn)存為MatrixCassignment(MatrixC)。 printf( 《數(shù)值分析》計(jì)算實(shí)習(xí)題目第一題\n)。 double u[501]。 /*直接返回距離原點(diǎn)P最接近的A的特征值*/}//主函數(shù)main(){ int i。 c=b。i++) { b+=y[i]*u[i]。 }//求βk for(i=0。j++) { u[i]=array[ij+2][j]*u[j]。i) { for(j=i+1。 for(i=499。j++) { u[i]=array[ij+2][j]*u[j]。i++) { for(j=max2(0,(i2))。 } for(i=1。i=500。i++) { y[i]=u[i]/a。 //求yk1 for(i=0。i++) { a=a+u[i]*u[i]。 //求ηk1 for(i=0。) { a=0。 } //先將矩陣Doolittle分解 LU(array)。i=500。 double y[501]。
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1