学数学也要学点偷懒的技术
2008-03-08 12:25:41 来源:WEB开发网核心提示: 我做这个微小软件其实是想在做线性代数作业偷懒,她很快帮你解出N阶行列式的值,学数学也要学点偷懒的技术,N阶线性方程的解,还能进行矩阵的相乘,不能进行运算!"); exit(0); }for(i=0;i<m;i++) { for(j=0;j<n;j++) c[i][j]=0; }print
我做这个微小软件其实是想在做线性代数作业偷懒,她很快帮你解出N阶行列式的值,N阶线性方程的解,还能进行矩阵的相乘。我拿它出来其实是让大家在学习其他课程是也能学C语言。以后我还会对它升级。大家多多留意喔!
#include<stdio.h>
#define N 10
#define M 10
#define S 10
void njie_valu();
void njief_value();
float Dvalue(float a[N][N],int n);
void juzhench(); void juzhench()
{
int i,j,k,m,n,s,s1;
float a[M][S],b[S][N],c[M][N];
PRintf("矩阵A是m*n矩阵。 ") ;
printf("请分别输入m,s的值:");
scanf("%d,%d",&m,&s);
printf("矩阵B是s*n矩阵 ");
printf("请分别输入s,n的值: ");
scanf("%d,%d",&s1,&n);
if(s!=s1) {
printf("两次输入的s的值不一样,不能进行运算! ");
exit(0);
}
for(i=0;i<m;i++) {
for(j=0;j<n;j++)
c[i][j]=0;
}
printf("请输入矩阵A ");
for(i=0;i<m;i++) {
for(j=0;j<s;j++)
scanf("%f",&a[i][j]);
}
printf("请输入矩阵B: ");
for(i=0;i<s;i++) {
for(j=0;j<n;j++)
scanf("%f",&b[i][j]);
}
for(i=0;i<m;i++) {
for(j=0;j<n;j++) {
for(k=0;k<s;k++)
c[i][j]=c[i][j]+a[i][k]*b[k][j];
}
}
printf("相乘后的结果是: ");
for(i=0;i<m;i++) {
for(j=0;j<n;j++)
printf("%f ",c[i][j]);
printf(" ");
}
getch();
} void njie_valu()
{
int i,j,n,k,h=0;
float x[N],a[N][N],D=1.0;
printf("说明:输入的n阶行列式,必须保证主对角线上的前n-1个元素都不为零! ");
printf("请输入是几阶行列式: ");
scanf("%d",&n);
printf("请以行列式的正确顺序输入个元素: ");
for(i=0;i<n;i++){
for(j=0;j<n;j++)
scanf("%f",&a[i][j]);
}
for(i=0,j=i;i<n-1,j<n-1;i++,j++)
if(a[i][j]==0)
h++;
if(h==n-1) {
printf("请确保主对角线上前n-1个元素都不为零,请您调整!☆ ");
printf("");
exit(0);
}
for(k=0;k<n-1;k++) {
for(i=k+1;i<n;i++)
x[i]=a[i][k]/a[k][k];
for(i=k+1;i<n;i++) {
for(j=0;j<n;j++)
a[i][j]=a[i][j]-a[k][j]*x[i];
}
}
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
printf("%f ",a[i][j]);
printf(" ");
}
printf(" ");
for(i=0,j=0;i<n,j<n;i++,j++)
D=D*a[i][j];
printf("D= %f",D);
getch();
}
void njief_value()
{
float D;
float Dvalue(float [N][N],int n);
int i,j,k,n,h,m;
float a[N][N],x[N],b[N],D1[N],a1[N][N],vable[N];
printf("说明::输入的n阶行列式,必须保证n阶线性方程的系数及系数和常数项构成的行列式主对角线上的前n-1个元素都不为零! ");
printf("请输入是几阶线性方程: ");
scanf("%d",&n);
printf("请以行列式的正确顺序输入个元素: ");
for(i=0;i<n;i++){
for(j=0;j<n;j++)
scanf("%f",&a[i][j]);
}
for(i=0;i<n;i++){
for(j=0;j<n;j++)
printf("%f ",a[i][j]);
printf(" ");
}
printf("按正确顺序输入常数项: ");
for(i=0;i<n;i++)
scanf("%f",&b[i]);
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
a1[i][j]=a[i][j];
}
D=Dvalue(a1,n);
printf("D= %f",D);
getch();
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
a1[i][j]=a[i][j];
}
for(i=0;i<n;i++) {
for(k=0;k<n;k++)
a1[k][i]=b[k];
for(m=0;m<n;m++) {
for(h=0;h<n;h++)
printf("%f ",a1[m][h]);
printf(" ");
}
printf(" ");
D1[i]=Dvalue(a1,n);
printf("D%d= %f ",i+1,D1[i]);
getch();
for(h=0;h<n;h++) {
for(m=0;m<n;m++)
a1[h][m]=a[h][m];
}
}
for(i=0;i<n;i++)
vable[i]=D1[i]/D;
for(i=0;i<n;i++){
printf("x%d=%f",i+1,vable[i]);
printf(" ");
}
getch();
} float Dvalue(float a[N][N],int n)
{
int i,j,k;
float x[N],D=1.0;
for(k=0;k<n-1;k++) {
for(i=k+1;i<n;i++)
x[i]=a[i][k]/a[k][k];
for(i=k+1;i<n;i++) {
for(j=0;j<n;j++)
a[i][j]=a[i][j]-a[k][j]*x[i];
}
}
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
printf("%f ",a[i][j]);
printf(" ");
}
printf(" "); for(i=0,j=0;i<n,j<n;i++,j++)
D=D*a[i][j];
return(D);
}
void main()
{
int i,j=1;
do {
clrscr();
printf("1.解n阶行列式 2.解n阶线性方程 3.求矩阵相乘 0.退出 ");
printf("请输入所需功能的序号:");
scanf("%d",&i);
switch(i){
case 1: njie_valu();break;
case 2: njief_value();break;
case 3: juzhench();break;
case 0: j=0;break;
}
}while(j);
}
#define N 10
#define M 10
#define S 10
void njie_valu();
void njief_value();
float Dvalue(float a[N][N],int n);
void juzhench(); void juzhench()
{
int i,j,k,m,n,s,s1;
float a[M][S],b[S][N],c[M][N];
PRintf("矩阵A是m*n矩阵。 ") ;
printf("请分别输入m,s的值:");
scanf("%d,%d",&m,&s);
printf("矩阵B是s*n矩阵 ");
printf("请分别输入s,n的值: ");
scanf("%d,%d",&s1,&n);
if(s!=s1) {
printf("两次输入的s的值不一样,不能进行运算! ");
exit(0);
}
for(i=0;i<m;i++) {
for(j=0;j<n;j++)
c[i][j]=0;
}
printf("请输入矩阵A ");
for(i=0;i<m;i++) {
for(j=0;j<s;j++)
scanf("%f",&a[i][j]);
}
printf("请输入矩阵B: ");
for(i=0;i<s;i++) {
for(j=0;j<n;j++)
scanf("%f",&b[i][j]);
}
for(i=0;i<m;i++) {
for(j=0;j<n;j++) {
for(k=0;k<s;k++)
c[i][j]=c[i][j]+a[i][k]*b[k][j];
}
}
printf("相乘后的结果是: ");
for(i=0;i<m;i++) {
for(j=0;j<n;j++)
printf("%f ",c[i][j]);
printf(" ");
}
getch();
} void njie_valu()
{
int i,j,n,k,h=0;
float x[N],a[N][N],D=1.0;
printf("说明:输入的n阶行列式,必须保证主对角线上的前n-1个元素都不为零! ");
printf("请输入是几阶行列式: ");
scanf("%d",&n);
printf("请以行列式的正确顺序输入个元素: ");
for(i=0;i<n;i++){
for(j=0;j<n;j++)
scanf("%f",&a[i][j]);
}
for(i=0,j=i;i<n-1,j<n-1;i++,j++)
if(a[i][j]==0)
h++;
if(h==n-1) {
printf("请确保主对角线上前n-1个元素都不为零,请您调整!☆ ");
printf("");
exit(0);
}
for(k=0;k<n-1;k++) {
for(i=k+1;i<n;i++)
x[i]=a[i][k]/a[k][k];
for(i=k+1;i<n;i++) {
for(j=0;j<n;j++)
a[i][j]=a[i][j]-a[k][j]*x[i];
}
}
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
printf("%f ",a[i][j]);
printf(" ");
}
printf(" ");
for(i=0,j=0;i<n,j<n;i++,j++)
D=D*a[i][j];
printf("D= %f",D);
getch();
}
void njief_value()
{
float D;
float Dvalue(float [N][N],int n);
int i,j,k,n,h,m;
float a[N][N],x[N],b[N],D1[N],a1[N][N],vable[N];
printf("说明::输入的n阶行列式,必须保证n阶线性方程的系数及系数和常数项构成的行列式主对角线上的前n-1个元素都不为零! ");
printf("请输入是几阶线性方程: ");
scanf("%d",&n);
printf("请以行列式的正确顺序输入个元素: ");
for(i=0;i<n;i++){
for(j=0;j<n;j++)
scanf("%f",&a[i][j]);
}
for(i=0;i<n;i++){
for(j=0;j<n;j++)
printf("%f ",a[i][j]);
printf(" ");
}
printf("按正确顺序输入常数项: ");
for(i=0;i<n;i++)
scanf("%f",&b[i]);
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
a1[i][j]=a[i][j];
}
D=Dvalue(a1,n);
printf("D= %f",D);
getch();
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
a1[i][j]=a[i][j];
}
for(i=0;i<n;i++) {
for(k=0;k<n;k++)
a1[k][i]=b[k];
for(m=0;m<n;m++) {
for(h=0;h<n;h++)
printf("%f ",a1[m][h]);
printf(" ");
}
printf(" ");
D1[i]=Dvalue(a1,n);
printf("D%d= %f ",i+1,D1[i]);
getch();
for(h=0;h<n;h++) {
for(m=0;m<n;m++)
a1[h][m]=a[h][m];
}
}
for(i=0;i<n;i++)
vable[i]=D1[i]/D;
for(i=0;i<n;i++){
printf("x%d=%f",i+1,vable[i]);
printf(" ");
}
getch();
} float Dvalue(float a[N][N],int n)
{
int i,j,k;
float x[N],D=1.0;
for(k=0;k<n-1;k++) {
for(i=k+1;i<n;i++)
x[i]=a[i][k]/a[k][k];
for(i=k+1;i<n;i++) {
for(j=0;j<n;j++)
a[i][j]=a[i][j]-a[k][j]*x[i];
}
}
for(i=0;i<n;i++) {
for(j=0;j<n;j++)
printf("%f ",a[i][j]);
printf(" ");
}
printf(" "); for(i=0,j=0;i<n,j<n;i++,j++)
D=D*a[i][j];
return(D);
}
void main()
{
int i,j=1;
do {
clrscr();
printf("1.解n阶行列式 2.解n阶线性方程 3.求矩阵相乘 0.退出 ");
printf("请输入所需功能的序号:");
scanf("%d",&i);
switch(i){
case 1: njie_valu();break;
case 2: njief_value();break;
case 3: juzhench();break;
case 0: j=0;break;
}
}while(j);
}
更多精彩
赞助商链接