就随便举个例子好了
比如一个4*4的矩阵,如何表示成上三角矩阵。
不好意思,前几天我论坛登不上,无法及时给你们回复.
我两个都有做但是都有一些问题.下面是我的做法:
如果是简单地把主对角线下的数去掉:
#include<stdio.h>
void main()
{
int N[3][3];
int i,j;
printf("请输入3*3的距阵\n");
for(i=0;i<3;i++)
{
for(j=0;j<3;j++)
{
scanf("%d",&N[i][j]);
printf("%5d",N[i][j]);
}
printf("\n");
}
printf("该距阵的上三角距阵为:\n");
for(i=0;i<3;i++)
{
for(j=0;j>=i&&j<3;j++)/*这两句有问题,运行结果只显示第一行,不知道错在哪里?*/
printf("%5d",N[i][j]);
printf("\n");
}
}
如果是线形代数上定义的上三角距阵:
我的做法:(题目算是难的,很难把算法完全写清楚。我们线形代数上的距阵化成上三角距阵是不能用小数表示的,最多是用分数,不知道C能不能表示分数?我用int型的,除非每一次运算都能整除,否则,结果就是不正确的)不知道怎么做才能正确表示出来。麻烦大家啦^-^
#include<stdio.h>
void main()
{
int N[3][3];
int i,j,k,l,m,n,o,p;
printf("请输入3*3的距阵\n");
for(i=0;i<3;i++)
{
for(j=0;j<3;j++)
{
scanf("%d",&N[i][j]);
printf("%5d",N[i][j]);
}
printf("\n");
}/*原样输出数组*/
for(i=0;i<3;i++)
if(N[i][0]!=0)
break;
for(j=0;j<3;j++)
{
m=N[i][j]/N[i][0];
N[i][j]=m;
k=N[0][j];
N[0][j]=N[i][j];
N[i][j]=k;
}/*化第一列的N[0][0]=1*/
for(i=1;i<3;i++)
{
for(j=0;j<3;j++)
{ n=N[i][j]-N[i][0]*N[0][j];
N[i][j]=n;
}
}/*把第一列除N[0][0]外的元素化为0*/
if(N[1][1]==0)
for(j=0;j<3;j++)
{ l=N[2][j];
N[2][j]=N[1][j];
N[1][j]=l;
}
else
for(j=0;j<3;j++)
{
o=N[1][j]/N[1][1];
N[1][j]=o;
p= N[2][j]-N[1][j]*N[2][1];
N[2][j]=p;
}/*化N[2][1]为0*/
printf("该距阵的上三角距阵为:\n");
for(i=0;i<3;i++)
{
for(j=0;j<3;j++)
printf("%5d",N[i][j]);
printf("\n");
}/*输出上三角距阵*/
}