线性代数老师布置的求规范正交基,明天就要交了.实在是没办法,高手帮忙呀!我已经编出来了,感觉逻辑上是没有错误的,但不知道为什么无法实现.
以下是代码:
/*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编辑过]

谢谢,我感觉我的算法是正确的.
而且我无意中实现了，我把xishuc函数合并到了schmidt函数里，居然就可以正常运算了！我也不知道为什么。我一直以为Turbo C是莫名其妙的。但不管它了。
我还有一个问题，和线性代数无关的。我现在已经可以算出规范正交基了,但还想搞个小的容错.就是我想在用scanf函数输入了m后,加一个对m不是数字的容错，用printf输出Please input a number!
我原本打算用isdigit，但发现它是测试字符型数据的。对我的int型的m变量的测试结果根本是乱七八糟。我在网上查到有用字符串数组测试的，但如果那样，那就无法使用变量m了。所以现在可以用束手无策来形容我了。
所以我想请教高手，如何加这样的一个容错。给我指个方向好吗？提示一下，谢谢了！
以下是我改好的代码：
/*guifanzhengjiaoji*/
#include<stdio.h>
#include<math.h>
#include<ctype.h>
#include<stdlib.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("Row number is wrong!");
}

void schmidt(int g)
{
float w,sh,xia;
int i,j;
for(j=0;j<g;j++){
sh=vemulti(j,g);
xia=vemulti(j,j);
if(xia!=0){
w=sh/xia;
w=w*(-1);
for(i=0;i<m;i++)
b[i][j]=a[i][j]*w;}
else{
printf("Vectors like these is prohibited!\n");
exit(0);}
}
for(j=0;j<g;j++){
for(i=0;i<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);
if(y!=0){
y=1/y;
for(i=0;i<m;i++)
a[i][g]*=y;
}
else{
printf("Vectors like these is prohibited!\n");
exit(0);}
}

main()
{
int i,j;
printf("put in the number of every vector:\n");
scanf("%d",&m);

printf("how many vectors do you want:\n");
scanf("%d",&n);

if(m<=10&&m>=1&&n<=10&&n>1)
{
for(j=0;j<n;++j)
{
for (i=0;i<m;++i)
{
printf("input %d number of vector %d:",i+1,j+1);
scanf("%f",&a[i][j]);/*2*/
}
}
}
else if (n=1) {
printf("Vectors like this is prohibited!");
exit(0);}
else{
printf("Please input a number which is less than 11 and more than 0!");
exit(0);}
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("%-6.3f ",a[i][j]);}
printf("\n");}
}

[此贴子已经被作者于2007-1-3 21:09:08编辑过]

容错就是如果错误地输入了不该输的东西，可以正确的返回提示信息，而不是错误的继续运行。
就像我上边说的我想输入数字给m变量，但如果打入了一个字母，比如f，应该是程序运行中断，屏幕上显示Please input a number!而不是继续运行，得到乱七八糟的错误结果。
谢谢帮忙！

这两天太忙了,现在终于有空了.我把作业已经交上了,没提这个容错
谢谢7楼海蓝啸了!
逻辑上是个不错的解决方案.我的思路没发展到能把字符串再转化成数值.
在过两天,我先编一下试试.等OK了,我再把我的全部正常的代码贴上来.