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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 【求助】哪位英雄帮我看一下这个小程序【已完结】
北京石油化工学院2026年研究生招生接收调剂公告
171/2返回列表12下一页
查看: 1081  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂求院校招收 +6 鹤鲸鸽 2026-03-28 6/300 2026-03-29 08:15 by fmesaito
[考研] 0703化学 +11 妮妮ninicgb 2026-03-27 11/550 2026-03-29 06:45 by 544594351
[考研] 070305高分子化学与物理 304分求调剂 +7 c297914 2026-03-28 7/350 2026-03-29 01:17 by 544594351
[硕博家园] 求调剂 330分 085600材料与化工 +3 gqhhh 2026-03-22 3/150 2026-03-29 00:52 by 544594351
[考研] 一志愿华理，数一英一285求A区调剂 +8 AZMK 2026-03-25 12/600 2026-03-28 18:15 by AZMK
[考研] 0703化学求调剂 +9 奶油草莓. 2026-03-22 10/500 2026-03-28 13:30 by 唐沐儿
[考研] 调剂 +3 好好读书。 2026-03-28 3/150 2026-03-28 12:04 by 王保杰33
[考研] 291求调剂 +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[考研] 311求调剂 +3 希望上岸阿小杨 2026-03-23 3/150 2026-03-28 07:57 by 热情沙漠
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 0703化学求调剂，各位老师看看我！！！ +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 085600材料与化工调剂 +10 A-哆啦Z梦 2026-03-23 16/800 2026-03-27 15:13 by caszguilin
[考研] 305求调剂 +5 哇卢卡库 2026-03-26 5/250 2026-03-27 14:01 by laoshidan
[考研] 276求调剂。有半年电池和半年高分子实习经历 +10 材料学257求调剂 2026-03-23 11/550 2026-03-27 10:13 by YCIT- LHL
[考研] 求调剂 +6 林之夕 2026-03-24 6/300 2026-03-27 08:38 by hypershenger
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
[考研] 351求调剂 +4 麦克阿磊 2026-03-24 4/200 2026-03-27 00:32 by wxiongid
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 296求调剂 +4 汪！？！ 2026-03-25 7/350 2026-03-25 16:41 by 汪！？！
信息提示
关闭
请填处理意见
关闭