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;
}