北航數(shù)值分析大作業(yè)一-文庫(kù)吧資料
2025-07-02 08:46本頁(yè)面
【正文】 } }}double fmifa(double u[501],double array[5][501],double p) { /*帶原點(diǎn)平移的反冪法*/ int i,j。t++) { array[jk+2][k]=array[jt+2][t]*array[tk+2][k]。j++) { for(t=max3(0,k2,j2)。 } } if(k500) for(j=k+1。t=(k1)。j=min2((k+2),500)。k=500。}void LU(double array[5][501]) /*對(duì)矩陣A進(jìn)行Doolittle分解*/{ /*矩陣A轉(zhuǎn)存在C[5][501]中*/ int j,k,t。i++) { a[i]=0。 for(i=0。i++) { a[i]=(*rand()/RAND_MAX)。 for(i=0。 /*直接返回A的特征值*/ }void chuzhi(double a[]) /*用隨機(jī)數(shù)為初始迭代向量賦值*/{ int i。 c=b。i++) { b+=y[i]*u[i]。 /*引入平移量*/ } //求βk for(i=0。j++) { u[i]+=array[ij+2][j]*y[j]。 for(j=max2(i2,0)。i=500。i++) { y[i]=u[i]/a。 //求yk1 for(i=0。i++) { a=a+u[i]*u[i]。 /*選用第一種迭代格式*///求ηk1for(i=0。) { a=0。 /* array[5][501]為矩陣A的轉(zhuǎn)存矩陣*/ double y[501]。 }}double mifa(double u[501],double array[5][501],double p) /*帶原點(diǎn)平移的冪法*/{ int i,j。 j++。j=500。 j+
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1