程序代码 用c语言程序实现改进的平方根法解方程组
余年寄山水
972次浏览
2020年07月29日 00:33
最佳经验
本文由作者推荐
游大理日记-呼幺喝六
#include "stdio.h"
#include "math.h"
#define N 20
int main()//主函数
{
int i,j,k,n;
float m; //定义变量
float a[N][N],l[N][N],lt[N][N],t[N][N],d[N][N];
float b[N],x[N],y[N];//定义矩阵
printf("ttt改进的平方根法（推理法）解方程组
");
printf("请输入矩阵的阶数：");
scanf("%d",&n);
printf("
");
printf("请输入系数矩阵a:
");
for(i=0;i
for(j=0;j
scanf("%f",&a[i][j]);//输入方程组系数矩阵a[][]
}
}
printf("
请输入常数项矩阵 b:
");
for(i=0;i
for(i=0;i
lt[i][i]=1;// 令lt[i][i]=1
}
for(i=0;i
t[i][j]=0;//令t[i][j]=0
}
for(j=0;j
lt[i][j]=0;// 令lt[i][j]=0
}
t[0][0]=a[0][0];
for(i=1;i
t[i][0]=a[i][0];//计算第一行的t[][]
lt[0][i]=a[0][i]/t[0][0];//计算第一列的lt[][]
}
for(i=1;i
for(j=1;j<=i;j++)//计算第2行到第n-1行的t[][]
{
t[i][j]=a[i][j];
for(k=0;k
t[i][j]=t[i][j]-t[i][k]*lt[k][j];
}
}
printf("
");
for(j=i+1;j
lt[i][j]=a[i][j];
for(k=0;k<=i-1;k++)
{
lt[i][j]=lt[i][j]-t[i][k]*lt[k][j];
}
lt[i][j]=lt[i][j]/t[i][i];
}
printf("
");
}
for(j=1;j
t[n-1][j]=a[n-1][j];
for(k=0;k<=j-1;k++)
{
t[n-1][j]=t[n-1][j]-t[n-1][k]*lt[k][j];
}
}
printf("
");
printf("输出矩阵t[i][j]
");
for(i=0;i
printf("%f",t[i][j]); printf(" ");//输出矩阵t[][]
}
printf("
");
}
printf("输出矩阵lt[i][j]
");
for(i=0;i
printf("%f",lt[i][j]); printf(" ");//输出矩阵lt[][]
}
printf("
");
}
printf("输出矩阵l[i][j]
");
for(i=0;i
l[i][j]=lt[j][i];
printf("%f",l[i][j]); printf(" ");//输出矩阵l[][]
}
printf("
");
}
printf("输出矩阵d[i][j]
");
for(i=0;i
if(i!=j)
d[i][j]=0;
else
d[i][j]=t[i][i];
printf("%f",d[i][j]); printf(" ");//输出矩阵d[][]
}
printf("
");
}
y[0]=b[0];
for(i=1;i
m=0.0;
for(k=0;k
m+=lt[k][i]*y[k];
y[i]=b[i]-m;
}
}// 计算yi
x[n-1]=y[n-1]/t[n-1][n-1];
for(i=n-1;i>=0;i--)
{t
m=0.0;
for(k=i+1;k
m+=lt[i][k]*x[k];
x[i]=y[i]/d[i][i]-m;
}
} //计算xi
for(i=0;i
",i,x[i]); //输
出xi
}