【正文】 《數(shù)值方法》實驗報告 32 } 。 } deltax=abs(X[N1]x)。 k++。 m++。j=i+deltan。 } m++。j++) { if(j=0) b=bX[j]*[m]。 for(j=ideltan。i++) //使用高斯 — 賽德爾迭代計算方程組的解 X { b=[k]。 for(i=0。 k=0。 bool run=true。 double delta=1e6。i++) //初始化數(shù)組 X X[i]=0。 //生成用于存儲方程組解的數(shù)組 X for(i=0。 deltan=int()。 } //用于計算方程組的解 double *save(typeA A,typeA B,int N,int n1) { int i,j,k,m,deltan。 system(pause)。iN。 //求解方程組 cout方程組的解為： endl。i++) cin[i]。 for(i=0。 //存儲帶狀方程組的值 cin。i++) cin[i]。 for(i=0。 cin。 =new double [n1]。 cout請輸入矩陣的每行最大元素數(shù)： 。 double *X=new double [N1]。 cout請輸入矩陣的階數(shù) N： 。 //實驗的主函數(shù) int main() { int N,n1,i。 int b。 } 4 includeiostream 《數(shù)值方法》實驗報告 30 include using namespace std。 } deltax=abs(X[N1]x)。 k++。 m++。j=i+deltan。 } m++。j++) { if(j=0) b=bX[j]*[m]。 for(j=ideltan。i++) //使用高斯 — 賽德爾迭代計算方程組的解 X { b=[k]。 for(i=0。 k=0。 bool run=true。 double delta=1e6。i++) //初始化數(shù)組 X X[i]=0。 //生成用于存儲方程組解的數(shù)組 X for(i=0。 deltan=int()。 } //用于計算方程組的解 double *save(typeA A,typeA B,int N,int n1) { int i,j,k,m,deltan。 system(pause)。iN。 //求解方程組 cout方程組的解為： endl。i++) cin[i]。 for(i=0。 //存儲帶狀方程組的值 cin。i++) cin[i]。 for(i=0。 cin。 =new double [n1]。 cout請輸入矩陣的每行最大元素數(shù)： 。 double *X=new double [N1]。 cout請輸入矩陣的階數(shù) N： 。 //實驗的主函數(shù) int main() { int N,n1,i。 int b。 } 3 includeiostream include using namespace std。 x[i]=y[i]/LU[i][i]。ji。i=0。 y[i]=B[i]。ji。i=N。 else B[i]=0。i=N。 float *B=new float[N]。 float *y=new float[N]。 } } return LU。k++) LU[j][k]=LU[j][k]c*LU[i][k]。 for(k=i。j=N。 k++。 B[i]=B[k]。 LU[k]=b。k=N) { b=LU[i]。 while(LU[i][i]==0amp。i=N。 //生成二維數(shù)組 LU，用于存放 L矩陣及 U矩陣 LU=A。 float *b=new float [N]。 } //求解 LU矩陣 的函數(shù)， A為系數(shù)矩陣 ， N為矩陣階數(shù) float **luchange(float **A,int *B,int N) { int i,j,k,a。 } system(pause)。j++) coutstd::leftsetw(15)invA[j][i]。i++) { for(j=0。 //輸出矩陣 A的逆矩陣 for(i=0。i++) //若存在求解逆矩陣 invA[i]=LUsave(LU,B[i],i,N)。 } } for(i=0。 《數(shù)值方法》實驗報告 26 system(pause)。i=N。 LU=luchange(A,B,N)。i=N。j++) cinA[i][j]。i++) for(j=0。 //輸入矩陣 A的值 for(i=0。i++) A[i]=new float [N]。 //生成用于存放矩陣 A的二維數(shù)組 for(i=0。 int *B=new int [N]。 N=N1。 cout請輸入矩陣的階數(shù)： 。 float **luchange(float**,int*,int)。 //實驗的 main函數(shù) int main() { int N,i,j。 return B。i=N。 《數(shù)值方法》實驗報告 25 B[0]=N。 double *B=new double [N]。 return A。j=N。i=N。i++) A[i]=new double [N]。 //生成二維動態(tài)數(shù)組A[N1][N1] for(i=0。 } //用于產(chǎn)生階數(shù)為 N的矩陣 A的函數(shù) double **buildA(int N) { int i,j。 x[i]=y[i]/U[i][i]。ji。i=0。 y[i]=B[i]。ji。i=N。 U[j][i]=c。k=N。j++) { c=U[j][i]/U[i][i]。i++) //將 A轉(zhuǎn)換為 LU矩陣 { 《數(shù)值方法》實驗報告 24 for(j=i+1。 for(i=0。 //生 成 保存結(jié)果的列矩陣 X double *y=new double [N]。 //生成二維數(shù)組 U，用于儲存 Lamp。 double c。 return 0。i++) //輸出 AX=B解的列 矩陣 X coutx[i+1]=x[i]endl。 for(i=0。 x=lufact(A,B,N)。 A=buildA(N)。 //輸入矩陣的階數(shù) ， 用于生成動態(tài)矩陣 cinN。 double *lufact(double**,double*,int)。 double **buildA(int)。 double **A。 //實驗的 main函數(shù) int main() { int N,i。 } return x。j) y[i]=y[i]x[j]*U[i][j]。i) //計算解 X的值 { for(j=N。 } for(i=N。j++) B[i]=B[i]y[j]*U[i][j]。i++) //計算中間矩陣 Y的值 { for(j=0。 } } for(i=0。k++) U[j][k]=U[j][k]c*U[i][k]。 for(k=i。j=N。i=N。 //生產(chǎn)用于保存中間值的列矩陣Y U=A。 //生成二維數(shù)組 U float *x=new float [N]。 float c。 return 0。i++) //輸出 AX=B解的列向量 X coutx[i+1]=x[i]endl。 for(i=0。 //輸入矩陣 B的值 《數(shù)值方法》實驗報告 22 x=lufact(A,B,N)。i=N。 cout請輸入矩陣 B的值： endl。k=N。i=N。 cout請輸入矩陣 A的值： endl。i=N。 float **A=new float*[N]。 float *B=new float