源程序编写在附件中,
有时程序会运行成功,但是有时会有如下错误,
"Floatint point error:Domain"
请问一下大侠,这是什么原因,应该怎么样修改程序,使运行不再出这种错误!
不晓得这个原因啊,期待达人解答.
贴一下我的程序,写的很粗糙,扩展性做的不好,也没有考虑溢出.见谅哈~
你修改下,加上clrscr()和getch(),还有注释就OK了.
#include<stdio.h>
#include<stdlib.h>
// 向量空间维数
#define VR 3
// 向量长度
#define VN 3
// 求内积
float neiji(float x1[], float x2[], int n)
{
float neiji = 0.0;
for(int i=0; i<n; i++)
{
neiji += (x1[i] * x2[i]);
}
return neiji;
}
// 向量相减
void xiangjian(float x1[], float x2[], int n)
{
for(int i=0; i<n; i++)
{
x1[i] -= x2[i];
}
}
// 数乘向量
void xiangcheng(float x1[], float x2[], float m, int n)
{
for(int i=0; i<n; i++)
{
x1[i] = m * x2[i];
}
}
// 施密特正交化
void zhengjiaohua(float (*b)[VN], float (*a)[VN])
{
int i,j;
float temp;
float tempNeiji[VR];
float tempXiangcheng[VN];
for(i=0; i<VR; i++)
{
for(j=0; j<VN; j++)
{
*(b[i]+j) = *(a[i]+j);
}
}
for(i=1; i<VR; i++)
{
tempNeiji[i-1] = neiji(b[i-1], b[i-1], VN);
for(j=0; j<i; j++)
{
temp = neiji(b[j], a[i], VN);
xiangcheng(tempXiangcheng, b[j], temp/tempNeiji[j], VN);
xiangjian(b[i], tempXiangcheng, VN);
}
}
}
int main()
{
int i,j;
float a[VR][VN];
float b[VR][VN];
for(i=0; i<VR; i++)
{
printf("please input a[%d] which has %d numbers: \n", i, VN);
for(j=0; j<VN; j++)
{
scanf("%f", &a[i][j]);
}
}
printf("\nthe numbers of inputed arr: \n");
for(i=0; i<VR; i++)
{
for(j=0; j<VN; j++)
{
printf("%.2f ", a[i][j]);
}
printf("\n");
}
zhengjiaohua(b, a);
printf("\nthe numbers of formated arr: \n");
for(i=0; i<VR; i++)
{
printf("\n");
for(j=0; j<VN; j++)
{
printf("%.2f ",b[i][j]);
}
}
}