c--解線性方程組的幾種方法-資料下載頁
2024-08-04 10:39本頁面
【正文】 =a[n][k]*a[k][j]。 a[n][j]=temp。 } for(i=n+1。ilenth。i++) //求第n+1層的下三角矩陣部分即L { for(k=0,temp=0。kn。k++) temp+=a[i][k]*a[k][n]。 a[i][n]=(a[i][n]temp)/a[n][n]。 } }}int Doolittle_check(double temp_a[][Number],double temp_b[Number]) //若行列式不為零,將系數(shù)矩陣調(diào)整為順序主子式大于零{ int i,j,k,maxi。double lik,temp。 for(k=0。klenth1。k++) { j=k。 for(maxi=i=k。ilenth。i++) if(temp_a[i][j]temp_a[maxi][j]) maxi=i。 if(maxi!=k) { exchange_hang(k,maxi)。 for(j=0。jlenth。j++) { temp=temp_a[k][j]。 temp_a[k][j]=temp_a[maxi][j]。 temp_a[maxi][j]=temp。 } } for(i=k+1。ilenth。i++) { lik=temp_a[i][k]/temp_a[k][k]。 for(j=k。jlenth。j++) temp_a[i][j]=temp_a[i][j]temp_a[k][j]*lik。 temp_b[i]=temp_b[i]temp_b[k]*lik。 } } if(temp_a[lenth1][lenth1]==0) return 0。 return 1。}void exchange_hang(int m,int n) //交換a[][]中和b[]兩行{ int j。 double temp。 for(j=0。jlenth。j++) { temp=a[m][j]。 a[m][j]=a[n][j]。 a[n][j]=temp。 } temp=b[m]。 b[m]=b[n]。 b[n]=temp。}void exchange(int m,int i) //交換A_y[m],A_y[i]{ int temp。 temp=A_y[m]。 A_y[m]=A_y[i]。 A_y[i]=temp。}void exchange_lie(int j) //交換未知量b[]和第i列{ double temp。int i。 for(i=0。ilenth。i++) { temp=a[i][j]。 a[i][j]=b[i]。 b[i]=temp。 }}void exchange_a_lie(int m,int n) //交換a[]中的兩列{ double temp。int i。 for(i=0。ilenth。i++) { temp=a[i][m]。 a[i][m]=a[i][n]。 a[i][n]=temp。 }}void exchange_x(int m,int n) //交換未知量x[m]與x[n]{ char temp。 temp=x[m]。 x[m]=x[n]。 x[n]=temp。}void recovery() //用其中一種方法求解后恢復(fù)數(shù)據(jù)以便用其他方法求解{ for(int i=0。ilenth。i++) for(int j=0。jlenth。j++) a[i][j]=copy_a[i][j]。 for(i=0。ilenth。i++) b[i]=copy_b[i]。 for(i=0。ilenth。i++) x[i]=39。a39。+i。 a_sum=0。 lenth=copy_lenth。}
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號-1