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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 【求助】哪位英雄帮我看一下这个小程序【已完结】
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1072  |  回复: 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的回帖

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的回帖
查看全部 17 个回答

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(金币+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的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 求调剂 330分 085600材料与化工 +3 gqhhh 2026-03-22 3/150 2026-03-29 00:52 by 544594351
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 317求调剂 +6 十闲wx 2026-03-24 6/300 2026-03-28 13:27 by Iveryant
[考研] 339求调剂，想调回江苏 +6 烤麦芽 2026-03-27 8/400 2026-03-28 10:40 by 烤麦芽
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +8 曼111 2026-03-24 9/450 2026-03-28 07:58 by YYYYX1234
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 348求调剂 +4 小懒虫不懒了 2026-03-27 5/250 2026-03-27 12:47 by 果果妈咪
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 329求调剂 +7 钮恩雪 2026-03-25 7/350 2026-03-27 04:28 by wxiongid
[考研] 351求调剂 +4 麦克阿磊 2026-03-24 4/200 2026-03-27 00:32 by wxiongid
[考研] 333求调剂 +6 wfh030413@ 2026-03-23 6/300 2026-03-26 22:45 by 学员8dgXkO
[考研] 321求调剂 +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
[考研] 材料与化工328分调剂 +6 。，。，。，。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
信息提示
关闭
请填处理意见
关闭