北航大學(xué)數(shù)值分析計算實(shí)習(xí)報告(編輯修改稿)
2024-08-18 01:55 本頁面
【文章內(nèi)容簡介】 ,)。 for(i=0。i501。i++) u[i]=i+1。 iterations = amp。 = Get_Fabs_Eigenvalue(C[2],u,iterations)。//將求得的絕對值最大特征值放到min變量中 += 。 if() { temp = 。 = 。 = temp。 temp1 = 。 = 。 = temp1。 } printf(最小特征值為：%.12e,迭代次數(shù)為：%d\n,)。 printf(最大特征值為：%.12e,迭代次數(shù)為：%d\n,)。 /**********************************/ /*以下利用反冪法求解模最小的特征值*/ /**********************************/ InitC(C)。//初始化C Initu(u)。 DoolittleC(C,501,2,2)。 = detA(C)。 iterations = amp。 = Get_min_Eigenvalue(C,u,501,2,2,iterations)。 = fabs()。 printf(絕對值最小特征值為：%.12e,迭代次數(shù)為：%d\n,)。 printf(A的行列式為：%.12e；,)。 printf(A的條件數(shù)cond(A)2為：%.12e\n,)。 /**********************************************/ /*以下利用反冪法求與數(shù)列Uk中元素最相近的特征值*/ /**********************************************/ double Uk[39]。//保存Uk的值并保存與Uk[i]最接近的λ double B[39]。 int diedai。 iterations = amp。diedai。 for( i=1。i=39。i++) { Uk[i1] = + i*( )/40。 InitC(C)。 Initu(u)。 for(int j =0。j501。j++) C[2][j] = Uk[i1]。 DoolittleC(C,501,2,2)。 B[i1] = Get_min_Eigenvalue(C,u,501,2,2,iterations)。 B[i1] +=Uk[i1]。 printf(μ%2d=%,與μ%2d最相近的特征值λ=%\n,i,Uk[i1],i,B[i1])。 } return 0。}double detA(double C[5][501]){ int i =0。 double e =1。 for(i =0。i501。i++) e*=C[2][i]。 return e。}void InitMatrix(double *p){ for(int i=1。i=501。i++) { *p=(*i)*sin(*i)*exp()。 p++。 }}void Initu(double *p){ for(int i=1。i=501。i++) { if(i=500) *p=0。 else *p=1。 p++。 } // *p=1。}void A_sub_minI(double *a,double min){ for(int i=1。i=501。i++) { *a=*amin。 a++。 }}double NeiJi(double *a,double *b){
試題試卷相關(guān)推薦
鄂ICP備17016276號-1