| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 948 人关注过本帖
标题:关于施密特正交
只看楼主 加入收藏
caoaihua3
Rank: 1
等 级:新手上路
帖 子:38
专家分:0
注 册:2006-4-20
收藏
 问题点数:0 回复次数:3 
关于施密特正交
最近在做一个关于施密特正交的问题:
源程序编写在附件中,

有时程序会运行成功,但是有时会有如下错误,
"Floatint point error:Domain"
请问一下大侠,这是什么原因,应该怎么样修改程序,使运行不再出这种错误!
搜索更多相关主题的帖子: 施密特正交 Domain point 
2006-10-27 17:08
caoaihua3
Rank: 1
等 级:新手上路
帖 子:38
专家分:0
注 册:2006-4-20
收藏
得分:0 
hgtfwseR.txt (1.15 KB) 关于施密特正交



希望大家帮忙啊
2006-10-27 17:13
C语言学习者
Rank: 4
等 级:贵宾
威 望:13
帖 子:1278
专家分:0
注 册:2006-9-26
收藏
得分:0 
虽然我不明白施密特正交,但我可以建议逐个调试一下,多插入printf()去看看。

谁有强殖装甲第二部,可以Q我460054868
2006-10-27 23:50
fans
Rank: 1
等 级:新手上路
帖 子:17
专家分:0
注 册:2006-9-4
收藏
得分:0 

不晓得这个原因啊,期待达人解答.
贴一下我的程序,写的很粗糙,扩展性做的不好,也没有考虑溢出.见谅哈~
你修改下,加上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]);
}
}
}

2006-10-28 19:13
快速回复:关于施密特正交
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.033396 second(s), 10 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved