导航:网站首页 >
C语言 解线性方程组 用c语言编程求线性方程组的解
C语言 解线性方程组 用c语言编程求线性方程组的解
相关问题:
匿名网友:
#include
#include
#include
#include
#define N 20
/*以下程序为不选主元的三角分解法(Doolittle)*/
main()
{
int i,j,k,s;
float a[N][N]={0},L[N][N]={0},U[N][N]={0},sigma1,sigma2,b[N],y[N],x[N];
/*为L主对角线元素赋1*/
for(i=0;i
{
L[i][i]=1;
}
/*输入矩阵a*/
for(i=0;i
{
printf("请输入矩阵第%d行元素:\n",i+1);
for(j=0;j
scanf("%f",&a[i][j]);
}
/*计算U第一行的元素和L第一列的元素*/
for(i=0;i
{
U[0][i]=a[0][i];
L[i][0]=a[i][0]/U[0][0];
}
for(k=1;k
{
/*计算矩阵U*/
for(j=k;j
{
sigma1=0;
for(s=0;s<=k-1;s++)
sigma1+=L[k][s]*U[s][j];
U[k][j]=a[k][j]-sigma1;
}
/*计算矩阵L*/
for(i=k;i
{
sigma2=0;
for(s=0;s<=k-1;s++)
sigma2+=L[i][s]*U[s][k];
L[i][k]=(a[i][k]-sigma2)/U[k][k];
}
}
/*输出矩阵L、U*/
printf("a矩阵为:\n");
for(i=0;i
{
for(j=0;j
printf("%5.1f ",a[i][j]);
printf("\n");
}
printf("L矩阵为:\n");
for(i=0;i
{
for(j=0;j
printf("%5.1f ",L[i][j]);
printf("\n");
}
printf("U矩阵为:\n");
for(i=0;i
{
for(j=0;j
printf("%5.1f ",U[i][j]);
printf("\n");
}
printf("请输入b矩阵\n",i+1);
for(i=0;i
scanf("%f",&b[i]);
/*回代法求解方程组Ly=b*/
for(i=0;i
{
sigma1=0;
for(k=0;k<=i-1;k++)
sigma1+=L[i][k]*y[k];
y[i]=b[i]-sigma1;
}
for(i=N-1;i>=0;i--)
{
sigma2=0;
for(k=i+1;k
sigma2+=U[i][k]*x[k];
x[i]=(y[i]-sigma2)/U[i][i];
}
printf("x为:\n");
for(i=0;i
printf("%5.1f ",x[i]);
printf("\n");
}
匿名网友:
1.C语言 解线性方程组
问:求C语言解线性方程组的代码。 另外,在网上找了不少编好的程序,运行总...
2.用c语言编程求线性方程组的解
问:例如 A={1 1}B={1} 求 AX=B的 x的解 1 2 2 在线等
3.C语言线性方程组求解
问:功能需求 1)输入线性方程组的系数矩阵和常数列 2)判断系数矩阵是否可逆 ...
4.c语言解线性方程组的编程题 【做的好会多给分】
问:这么一道矩阵题,(跟在字母后面的数为下标,因为打起来费事,所以偷懒这...
5.怎样用C语言实现线性方程组的求解
答:输入:系数矩阵 或者 字符串 后者要对字符串处理获得系数 求解:可以采用高斯消元法等方法求解,可以参考网页链接
6.求线性方程组C语言
问:已知多项式f(x)的系数数组f【】,求g(x)的系数数组g【】使得f(x)*...
7.有关C语言解线性方程组问题,求解。
问:这段代码 我用高斯消元解线性方程组 运行到回代输出一次 回代后又输出一...
8.用C语言做Gmres算法,求解线性方程组,各位有代码吗
答:#include#include#include#include#defineN20/*以下程序为不选主元的三角分解法(Doolittle)*/main(){inti,j,k,s;floata[N][N]={0},L[N][N]={0},U[N][N]={0},sigma1,sigma2,b[N],y[N],x[N];/*为L主对角线元素赋1*/for(i=0;i
9.C语言实现cholesky法解任意线性方程组
问:C语言实现cholesky法解任意线性方程组网上找了好多都不能用,要不就是不...
10.c语言 方程组的解数
问:给出一个方程组:a的平方+b=n a+b的平方=m和n、m。 计算有多少对整数(a,...
问题推荐