数值计算方法实验五

1.用雅可比迭代法、高斯—塞德尔迭代法和松弛法求下列线性方程组的解，要求准确到1/2×10-3。

要求：

(1) 线性方程组的维数n、系数矩阵A的元素、列向量b的元素从键盘输入；

(2) 误差ε和松弛因子ω从键盘输入；

(3) 在同一个程序里将雅可比迭代法、高斯—塞德尔迭代法和松弛法编程实现；

(4) 绘制这三种算法求解线性方程组的算法框图；

(5) 改变ω的值，根据迭代计算结果比较ω取不同值的情况下迭代过程的收敛速度；

(6) 将这三种算法的每一步迭代计算结果打印输出；

(7) 根据计算结果，比较这三种算法的迭代收敛速度。

 

算法框图：

程序源代码：

程序源代码：
#include<iostream>
#include<cmath>
using namespace std;
int main(){
	int a[10][10],a1[10][10],a2[10][10];
	double b[10],b1[10],b2[10];
	double x[10]={0.0};
	double y[10]={0.0},ch[10] = {0.0};
	int n,i,j,k,N=50; 
	double e,w,sum;
	cout<<"请输入线性方程组的维数"<<endl;
	cin>>n;
	cout<<"请输入系数矩阵A的元素"<<endl;
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++){
		    cin>>a[i][j]; 
		    a1[i][j] = a[i][j];
		    a2[i][j] = a[i][j];
		} 
	}
	cout<<"请输入列向量b"<<endl;
	for(i=1;i<=n;i++){
		cin>>b[i];
		b1[i] = b[i];
		b2[i] = b[i];
	}
	cout<<"请输入误差e和松弛因子w"<<endl;
	cin>>e>>w;
	cout<<"雅可比迭代法："<<endl<<endl;
	k=1;
	do{
		for(i=1;i<=n;i++){
			sum=0;
			for(j=1;j<=n;j++){
				if(j!=i){
				    sum += a[i][j]*x[j]; 	
				}
			}
			y[i] = (b[i] - sum)/a[i][i];
			ch[i] = abs(x[i] - y[i]);
		}
		
		double max = ch[1];
	    for(i=2;i<=n;i++){
	    	if(ch[i] > max){
	    		max = ch[i];
			}
		}
		cout<<"第"<<k<<"次结果为： ";
	    for(i=1;i<=n;i++){
            printf("X%d:  %.6lf  ",i,y[i]);
		}		
		cout<<endl;
		if(max<e){
			break;
		}
		else{
			k++;
			if(k>N){
				cout<<"超过最大迭代次数，迭代失败"<<endl;
			} 
			for(i=1;i<=n;i++){
				x[i] = y[i];
			}
		}	
	}while(1);
	
	for(i=0;i<=10;i++){
	   x[i] = 0;
	   y[i] = 0;
	   ch[i] = 0;		
    }
	cout<<endl<<endl<<"高斯-赛德尔迭代法："<<endl;
	k=1;
	do{
		for(i=1;i<=n;i++){
			sum=0;
			for(j=1;j<=n;j++){
				if(j!=i){
				    sum += a1[i][j]*y[j]; 	
				}
			}
			y[i] = (b1[i] - sum)/a1[i][i];
			ch[i] = abs(x[i] - y[i]);
		}
		
		double max = ch[1];
	    for(i=2;i<=n;i++){
	    	if(ch[i] > max){
	    		max = ch[i];
			}
		}
		cout<<"第"<<k<<"次结果为： ";
	    for(i=1;i<=n;i++){
            printf("X%d:  %.6lf  ",i,y[i]);
		}	
		cout<<endl;	
		if(max<e){
			break;
		}
		else{
			k++;
			if(k>N){
				cout<<"超过最大迭代次数，迭代失败"<<endl;
			} 
			for(i=1;i<=n;i++){
				x[i] = y[i];
			}
		}	
	}while(1);

    for(i=0;i<=10;i++){
	   x[i] = 0;
	   y[i] = 0;
	   ch[i] = 0;		
	}
	cout<<endl<<endl<<"松弛法："<<endl;
	k=1;
	do{
		for(i=1;i<=n;i++){
			sum=0;
			for(j=1;j<=n;j++){
				if(j!=i){
				    sum += a2[i][j]*y[j]; 	
				}
			}
			y[i] = (b2[i] - sum)/a2[i][i];	
			y[i] = (1-w)*x[i] + w*y[i];
			ch[i] = abs(x[i] - y[i]);
		}
	
		double max = ch[1];
	    for(i=2;i<=n;i++){
	    	if(ch[i] > max){
	    		max = ch[i];
			}
		}
		cout<<"第"<<k<<"次结果为： ";
	    for(i=1;i<=n;i++){
            printf("X%d:  %.6lf  ",i,y[i]);
		}
		cout<<endl;		
		if(max<e){
			break;
		}
		else{
			k++;
			if(k>N){
				cout<<"超过最大迭代次数，迭代失败"<<endl;
			} 
			for(i=1;i<=n;i++){
				x[i] = y[i];
			}
		}	
	}while(1);

	return 0;
}

2.用高斯—约当消元法或全选主元高斯消元法求解下列线性方程组的解

要求：

(1) 线性方程组的维数n、系数矩阵A的元素、列向量b的元素从键盘输入；

(2) 用算法框图描述相应算法的实现过程；

(3) 将每一次消元过程得到的增广矩阵打印输出；

(4) 最后输出方程组的解。

算法框图：

程序源代码：

源代码：
#include<iostream>
using namespace std;
int main(){
	double a[10][10],b[10],x[10]={0};
	int n,i,j,count=1;
	cout<<"请输入线性方程组的维数"<<endl;
	cin>>n;
	cout<<"请输入系数矩阵A的元素"<<endl;
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++){
			cin>>a[i][j];
		}
	}
	cout<<"请输入列向量b"<<endl;
	for(i=1;i<=n;i++){
		cin>>b[i];
		a[i][n+1] = b[i];
	}
	int k,c,d;
	double ch;
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++){
			if(i!=j){
				ch = -a[j][i]/a[i][i];
				for(k=0;(i+k)<=n+1;k++){
				    a[j][i+k] = a[i][i+k]*ch+a[j][i+k];
			    }
			}	
		}
		cout<<endl;
		cout<<"第"<<count<<"次消元得到的增广矩阵"<<endl; 
		for(c=1;c<=n;c++){
			for(d=1;d<=n+1;d++){
                printf("%.3lf   ",a[c][d]);
			}
			cout<<endl;
		}	
	}
	cout<<endl<<"最后求得各x为"<<endl;
	for(i=1;i<=n;i++){
		cout<<"X"<<i<<": "<<a[i][n+1]/a[i][i]<<endl;
	}
	return 0;
}