北航數(shù)值分析大作業(yè)一-全文預覽
【正文】 ray[jk+2][k]=array[jt+2][t]*array[tk+2][k]。 } } if(k500) for(j=k+1。j=min2((k+2),500)。}void LU(double array[5][501]) /*對矩陣A進行Doolittle分解*/{ /*矩陣A轉存在C[5][501]中*/ int j,k,t。 for(i=0。 for(i=0。 c=b。 /*引入平移量*/ } //求βk for(i=0。 for(j=max2(i2,0)。i++) { y[i]=u[i]/a。i++) { a=a+u[i]*u[i]。) { a=0。 }}double mifa(double u[501],double array[5][501],double p) /*帶原點平移的冪法*/{ int i,j。j=500。 array[1][j]=。 j++。j=500。) { array[i][j]=0。//所有元素歸零 for(i=0。 else t=b。 }int min2(int a,int b) /*求兩個整型數(shù)最小值的子程序*/{ if(ab) return b。①cond（A）2=|λ1/λn|，其中λ1和λn分別是矩陣A的模最大和最小特征值。②使用帶原點平移的冪法（mifa（）函數(shù)），令平移量p=λmax，求出對應的按摸最大的特征值λ，max，如果λmax0,則λ1=λ，max+p；如果λmax0,則λ501=λ，max+p。λmin即為λs；如果λmax0,則λ501=λmax；如果λmax0,則λ1=λmax。（譜范數(shù)）條件數(shù)cond（A）2和行列式detA。 else return b。 if(ab) t=a。}void assignment(double array[5][501]) /*將矩陣A轉存為數(shù)組C[5][501]*/{ int i,j,k。j=500。 }//第0,4行賦值 for(
環(huán)評公示相關推薦
鄂ICP備17016276號-1