线性代数老师布置的求规范正交基,明天就要交了.实在是没办法,高手帮忙呀!我已经编出来了,感觉逻辑上是没有错误的,但不知道为什么无法实现.
以下是代码:
/*guifanzhengjiaoji*/
#include<stdio.h>
#include<math.h>
float a[10][10];
float b[10][10];
int m,n;
float vemulti(int g,int l) /*vector multiply*/
{
int i;
float r=0;
if(g>=0&&g<n&&l>=0&&l<n){
for(i=0;i<m;i++)
r+=a[i][g]*a[i][l];
return r;}
else printf("liehaocuowu!");
}
int xishuc(int h,int f)
{
float q;
float sh,xia;
int i;
sh=vemulti(h,f);
xia=vemulti(h,h);
q=sh/xia;
q=q*(-1);
for(i=0;i<m;i++)
b[i][h]=a[i][h]*q;
return 0;
}
void schmidt(int g) /*正交化函数*/
{
int i,q,j;
q=0;
while(q<g){
if(xishuc(q,g)==0);
q++;}
for(j=0;j<g;j++){
for(i=0;j<m;i++)
a[i][g]=a[i][g]+b[i][j];
}
}
void biaozhun(int g)
{
float w;
float y;
int i;
w=vemulti(g,g);
y=sqrt(w);
y=1/y;
for(i=0;i<m;i++)
a[i][g]*=y;
}
main()
{
int i,j;
printf("put in the number of every vector:\n");
scanf("%d",&m); /*1*/
printf("how many vectors do you want:\n");
scanf("%d",&n);
if(m<10&&m>=0&&n<10&&n>=0)
{
for(j=0;j<n;++j)
{
for (i=0;i<m;++i)
{
printf("input the %dth number of vector %d:",i+1,j+1);
scanf("%f",&a[i][j]);/*2*/
}
}
}
else
printf("too many you want!");
for(j=1;j<n;j++){
schmidt(j);
}
for(j=0;j<n;j++)
biaozhun(j);
printf("the result is:\n");
for(i=0;i<m;i++){
for(j=0;j<n;j++){
printf("%f ",a[i][j]);}
printf("\n");}
}
高手帮忙呀!
[此贴子已经被作者于2007-1-3 21:50:17编辑过]