线性方程组的求解急 - C语言论坛

tronlove

等　级：新手上路
帖　子：7
专家分：0
注　册：2005-12-20

楼主

问题点数：0 回复次数：3

线性方程组的求解急

1。方程组为齐次方程Ax=0,判断是否有零解，如果无零解，求出它的基础解系
2。方程组为非齐次方程Ax=b,判断是否有解，如果方程组无解，要给出相应的信息，如果有唯一解，求出它的解，如果有无穷多解，给出解的合理形式。
要求
1。矩阵A和矩阵B由键盘输入

小生先在此谢过

搜索更多相关主题的帖子: 线性方程　求解　

2005-12-20 00:47

303770957

等　级：小飞侠
威　望：6
帖　子：838
专家分：2125
注　册：2005-9-10

第 3 楼

得分:0

看完了写一份感想给我好么?
程序清单：
#include "stdlib.h"
#include "math.h"
#include "stdio.h"
#define MAX 255
int Guass(double a[],double b[],int n)
{int *js,l,k,i,j,is,p,q;
double d,t;
js=malloc(n*sizeof(int));
l=1;
for(k=0;k<=n-2;k++)
{d=0.0;
for(i=k;i<=n-1;i++)
for(j=k;j<=n-1;j++)
{t=fabs(a[i*n+j]);
if(t>d)
{d=t;js[k]=j;is=i;}
}
if(d+1.0==1.0)l=0;
else {if(js[k]!=k)
for(i=0;i<=n-1;i++)
{p=i*n+k;
q=i*n+js[k];
t=a[p];a[p]=a[q];a[q]=t;}
if(is!=k)
{for(j=k;j<=n-1;j++)
{p=k*n+j;
q=is*n+j;
t=a[p];a[p]=a[q];a[q]=t;}
t=b[k];b[k]=b[is];b[is]=t;}
}
if(l==0){free(js);
printf("该矩阵为奇异矩阵！\n");
return (0);}
d=a[k*n+k];
for(j=k+1;j<=n-1;j++)
{p=k*n+j;
a[p]=a[p]/d;}
b[k]=b[k]/d;
for(i=k+1;i<=n-1;i++)
{for(j=k+1;j<=n-1;j++)
{p=i*n+j;
a[p]=a[p]-a[i*n+k]*a[k*n+j];}
b[i]=b[i]-a[i*n+k]*b[k];
}
}
d=a[(n-1)*n+n-1];
if(fabs(d)+1.0==1.0)
{free(js);
printf("该矩阵为奇异矩阵！\n");
return (0);}
b[n-1]=b[n-1]/d;
for(i=n-2;i>=0;i--)
{t=0.0;
for(j=i+1;j<=n-1;j++)t=t+a[i*n+j]*b[j];
b[i]=b[i]-t;}
js[n-1]=n-1;
for(k=n-1;k>=0;k--)
if(js[k]!=k)
{t=b[k];b[k]=b[js[k]];b[js[k]]=t;}
free(js);
return (1);}
void main()
{int i,n;
double A[MAX],B[MAX];
printf("请输入要求解的齐次线性方程组的介数:n=");
scanf("%d",&n);
printf("请逐个输入%d介齐次线性方程组未知数的系数矩阵共%dX%d个:\n",n,n,n);
for(i=0;i<n*n;i++)
scanf("%lf",&A[i]);
printf("请逐个输入%d介齐次线性方程组右边的常量矩阵共%d个:\n",n,n);
for(i=0;i<n;i++)
scanf("%lf",&B[i]);
if(Guass(A,B,n)!=0)
{printf("要求解的%d介齐次线性方程组的解为:\n",n);
for(i=0;i<n;i++)
printf("X(%d)=%f\n",i+1,B[i]);
}
else printf("要求解的%d介齐次线性方程组无解.\n",n);
}

♂ 死后定当长眠，生前何须久睡。♀

2005-12-24 20:16