24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lichuanhang

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 729.1
帖子: 402
在线: 6.1小时
虫号: 360470
注册: 2007-05-03
性别: GG
专业: 生物化工

[交流] 【求助】哪位大哥能帮忙看看这个程序的错误啊【已完结】

#include
#include
#include
void main()
{
int i;
double a=1,b=3,c,d;
double K[3],C[3],E1=4.6447E+03,f1[3],f3[3];
double P=1.013,X;
X=a/b;
double F1=0,F2=0,F3=0,F4=0;
float t;
for(i=0;i<3;i++)
{C=((i+1)*4E+03)/P;
/*printf("%f\n",C);*/
}
printf("请输入温度初值t:";
scanf("%f",&t);
/*printf("%.1f\n",t);*/
do{

for(i=0;i<3;i++){
d=-E1/(1.8*t+492);
c=pow(2.718,d);
K=C*c;
/*printf("%f\n",K);*/
}
for(i=0;i<3;i++){
f3=K*X;
/*printf("%f\n",K);
printf("%f\n",X);*/
F3+=f3;
/*printf("%f\n",f3);*/
}
/*printf("%f\n",F3);*/
F2=F3-1;
F3=F4;
printf("%f\n",F2);

if (F2<0) break;

for(i=0;i<3;i++){
f1=K*X*1.8*E1/((1.8*t+492)*(1.8*t+492));
F1+=f1;

}
/*printf("%f\n",F1); */
t=t-F2/F1;
/*printf("%f\n",t); */
F1=F4;
}
while(1);
printf("%fPa时的泡点温度为：%0.1f",P,t);
}
谢谢各位，这个程序使用牛顿迭代求泡点温度，是我的VC有问题，结果没问题

[ Last edited by 余泽成 on 2009-10-15 at 11:07 ]

回复此楼

» 猜你喜欢

相逢一笑皆知己

1楼 2009-10-14 17:31:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

magic7004

金虫 (职业作家)

程序强帖: 2
应助: 0 (幼儿园)
金币: 444
红花: 4
帖子: 3278
在线: 156.2小时
虫号: 482895
注册: 2007-12-26
性别: GG
专业: 化学工程

★
余泽成(金币+1,VIP+0):辛苦，欢迎常来程序软件版！ 10-15 09:58

vc6.0试了下，语法错误：
1.  X=a/b这行报错，把它移动到float t;后面就OK了
2.  t=t-F2/F1;这行F2和F1是double，t是float，转换可能损失精度。

逻辑错误肯定有，搂住自己找吧。我试了下，结果是：
初值          泡点温度
0             0
1             1
3             3
5             5
200          死循环
300          死循环
500          死循环