应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 【求助】哪位英雄帮我看一下这个小程序【已完结】
171/2返回列表12下一页
查看: 1205  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

li343

至尊木虫 (职业作家)

[交流] 【求助】哪位英雄帮我看一下这个小程序【已完结】

程序编译正常，有结果（不知道是否正确），但发现v一直为1，无论温度T怎么给，v都不变，这就是个问题了，请高手帮看一下是啥原因。
    先谢谢了！

#include
#include

void main()
{
   int i,j;
   double t,a1,a2,a3,P,Pc,Tc,T,Tr;
   double c0,c1,c2,c3,Mole,a_T,a1_T,v,R,X0,w;
   double m,ac,b,Afa_T;
   double h1,h2,H;

   a1=6.6400537;
   a2=8.3412027;
   a3=27.498677;

   if(-100    {
          c0=3.089329; c1=0.1020401; c2=-2.400100e-4; c3=4.126521e-7;
   }
   if(298<=T<1500)
   {
          c0=-2.207131; c1=0.1002983; c2=-0.5530711e-4; c3=0.1187821e-7;
   }
   
   X0=249.4758;

   Tc=408.16;  //临界温度/K
   Pc=3684e+3;    //临界压力/kPa
   Mole=58.12; //分子量
   w=0.1756;   //偏心因子
   R=8.314;    //气体常数

   T=338.15;  
   Tr=T/Tc;
   t=1-Tr;
   
   m=0.37464+1.54226*w-0.22992*w*w;
   ac=0.45724*pow((R*Tc),2)/Pc;
   b=0.07780*R*Tc/Pc;
   Afa_T=pow((1+m*(1-pow(Tr,0.5))),2);
   a_T=ac*Afa_T;
   a1_T=-m*a_T*pow(Tr,0.5)*pow(Afa_T,-0.5)/T;

   double V[1000];
   V[0]=0.7;
   for(i=0;i<=1000;i++)
   {
           V[i+1]=pow(((R*T-P*b)*pow(V,2)+(2*R*T*b+3*P*b-a_T)*V-(P*pow(b,3)+R*T*pow(b,2)-a_T*b))/P,(1/3));
v=V;
   }
   
   printf("%f\n",m);
   printf("%f\n",ac);
   printf("%f\n",b);
   printf("%f\n",Afa_T);
   printf("%f\n",a_T);
   printf("%f\n",a1_T);

   printf("\n";

   P=Pc*exp((a1+a2*pow(t,1.78)+a3*pow(t,6.06))*log(Tr));  //饱和蒸汽压,Pa

   h1=c0*T+c1*pow(T,2)/2+c2*pow(T,3)/3+c3*pow(T,4)/4;
   h2=(a_T-T*a1_T)*log((v-0.4142*b)/(v+2.4142*b))/(2*sqrt(2)*a1_T*b);
   H=4.1868*h1/Mole+h2+P*v-R*T+X0; //焓值,J/(kg*K)

   printf("%lf\n",P/1000);
   printf("%lf\n",v);
   printf("%lf\n",H/1000);

}

[ Last edited by 余泽成 on 2009-11-26 at 20:25 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-07-30 10:35:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sitonyl

铁杆木虫 (文坛精英)
★
li343(金币+1,VIP+0): 7-30 13:50
v=V;
应该是这个赋值语句有问题，我觉得应该是v=*(V+i)；或v=V；这样的。
赞一下(2人) 回复此楼
zz
2楼2009-07-30 10:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li343

至尊木虫 (职业作家)
引用回帖:
Originally posted by sitonyl at 2009-7-30 10:48:
v=V;
应该是这个赋值语句有问题，我觉得应该是v=*(V+i)；或v=V；这样的。

哪个语句原来就是你写的那样的，不知道咋不显示的！

[ Last edited by li343 on 2009-7-30 at 11:41 ]
赞一下 回复此楼
3楼2009-07-30 11:32:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

羽寒月

银虫 (小有名气)
★
li343(金币+1,VIP+0): 7-30 13:50
为什么我编译就通不过捏~~
test.c: In function `main':
test.c:48: error: incompatible type for argument 1 of `pow'
test.c:48: error: invalid operands to binary *
test.c:49: error: incompatible types in assignment
test.c:5: warning: return type of 'main' is not `int'
赞一下(1人) 回复此楼
4楼2009-07-30 12:35:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

羽寒月

银虫 (小有名气)
★ ★ ★ ★ ★ ★ ★
li343(金币+2,VIP+0): 7-30 13:30
li343(金币+5,VIP+0): 7-30 13:50
找到两个错误：
1、 变量P没有赋值
2、1/3要改成1.0/3.0，整型运算1/3=0，所以你的v一直是1.0
赞一下(3人) 回复此楼
5楼2009-07-30 12:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ms98

★ ★
li343(金币+1,VIP+0): 7-30 13:50
rjjy(金币+1,VIP+0):呵呵！ 7-30 17:15
if(-100 我看到这儿就不想往下看了。
话说来提问题是好事，但是提问之前不动手调试就不太好了。
另外，编译器不报错只表明程序中没有编译器能找到的错，主要还是靠人来看。

[ Last edited by ms98 on 2009-7-30 at 13:03 ]
赞一下(11人) 回复此楼
6楼2009-07-30 13:02:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

richsuper

银虫 (小有名气)
★
wangen994(金币+1,VIP+0):谢谢参与讨论 7-30 21:54
同志啊，多写代码，多多调试，看这是if(-100 ？
赞一下(10人) 回复此楼
7楼2009-07-30 21:14:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjpm

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
强贴留名！
引用回帖:
Originally posted by richsuper at 2009-7-30 21:14:
同志啊，多写代码，多多调试，看这是if(-100 ？

回复此楼
8楼2009-07-30 22:31:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shipsw

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wangen994(金币+2,VIP+0):感谢你参与讨论 8-10 08:26
循环体 for 内部出了个赋值语句 v=V意外。其他的没有调用变量v。所以这个语句应该移到循环体外部。
而且v 与V 是不同类型的。所以也不对。
赞一下(11人) 回复此楼
选择了这条路也就放弃了其他的路了。
9楼2009-08-07 09:59:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magic7004

金虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
搂住的代码很强大
赞一下(1人) 回复此楼
流氓不可怕，可怕的是流氓有文化，有文化又BH的流氓无敌~~！
10楼2009-08-09 21:10:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li343 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了，不知道这次命运如何 +8 雨打竹帘 2026-08-30 12/600 2026-09-03 06:24 by sincosx
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 6/300 2026-09-02 23:59 by rM1TE0WVDIIY
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +14 Kittylucky 2026-08-27 15/750 2026-09-02 19:35 by liumei0601
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗？ +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
信息提示
关闭
请填处理意见
关闭