| 查看: 738 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求大神帮忙看一下程序，为什么一直无法执行选择之后的循环已有2人参与
|
|||
|
编了一个程序如下，希望p<pjc时使用一个循环，p>=pjc时使用另一个循环。在编程过程中发现不管是使用if选择语句还是while循环语句都是只能执行p<pjc时的循环，当p>pjc时程序就停止了，请各位大神帮忙看一下是什么原因，谢谢大家。 clear all p0=101.3;kp=7.9;T0=300;md0=1.3;u0=0.5;kz=4; Tbv=2221;Tbp=1867; V=5000;kl=1.4;kr=1.25; pb=p0*kp; pjc=p0+5; t(1)=0;tj=0.001; Vr(1)=0; Vrj(1)=0; p(1)=p0; u(1)=u0; krj(1)=0; rr(1)=0; mdw(1)=md0; Tw(1)=T0; i=1; %开始循环计算 while p(i)<=pb & p(i)>=p0 while p<=pjc & i<319 i=i+1; t(i)=t(i-1)+tj; uz(i)=u0*(p(i-1)/p0)^0.25 Vrj(i)=(4*pi/3)*(((3*Vr(i-1)/(4*pi))^(1/3)+(kz*uz(i)*tj))^3)-Vr(i-1); krj(i)=(Vrj(i)/V)*((p(i-1)/p0)^(1/kl)); p(i)=p(i-1)+krj(i)*(pb-p0); Vr(i)=V-(V-Vr(i-1)-Vrj(i))*((p(i-1)/p(i))^(1/kl)); Vw(i)=V-Vr(i); rr(i)=((3*Vr(i))/(4*pi))^(1/3); u(i)=((rr(i)-rr(i-1))/tj); mdw(i)=md0*((p(i)/p0)^(1/kl)); Tw(i)=T0*((p(i)/p0)^((kl-1)/kl)); Tr(i)=Tbp+(Tbv-Tbp)*((p(i)-p0)/(pb-p0)); mdr(i)=(p(i)*Tbv*md0)/(pb*Tr(i)); end while p>pjc i=i+1; t(i)=t(i-1)+tj; uz(i)=u0*(p(i-1)/p0)^0.25 Vrj(i)=(4*pi/3)*(((3*Vr(i-1)/(4*pi))^(1/3)+(kz*uz(i)*tj))^3)-Vr(i-1); krj(i)=(Vrj(i)/V)*((p(i-1)/p0)^(1/kl)); p(i)=p(i-1)+krj(i)*(pb-p0); pb(i)=p(i) %pb(i)=p(i-1)+krj(i)*(pb-p0); Vr(i)=V-(V-Vr(i-1)-Vrj(i))*((p(i-1)/p(i))^(1/kl)); Vw(i)=V-Vr(i); rr(i)=((3*Vr(i))/(4*pi))^(1/3); u(i)=((rr(i)-rr(i-1))/tj); mdw(i)=md0*((p(i)/p0)^(1/kl)); %Tw(i)=T0*((p(i)/p0)^((kl-1)/kl)); %Tr(i)=Tbp+(Tbv-Tbp)*((p(i)-p0)/(pb-p0)); %mdr(i)=(p(i)*Tbv*md0)/(pb*Tr(i)); Vbw(i)=(V-Vr(i-1)-Vrj(i))*((p(i-1)/p(i))^(1/kl)); Vbr(i)=V-Vbw(i); mdbw(i)=mdw(i-1)*(Vw(i-1)-Vrj(i))/Vbw(i); mdbr(i)=(mdr(i-1)*Vr(i-1)+mdw(i-1)*Vrj(i))/Vbr(i); Sl=150;Sr=150; Vxr(i)=((Sr*tj)/mdbr(i))*((2*p(i-1)*mdr(i-1)*(((p0/p(i-1))^(2/kr))-((p0/p(i-1))^((kr+1)/kr))))^(1/2)); Vxl(i)=((Sl*tj)/mdbw(i))*((2*p(i-1)*mdw(i-1)*(((p0/p(i-1))^(2/kl))-((p0/p(i-1))^((kl+1)/kl))))^(1/2)); pr(i)=pb(i)*(((Vbw(i)-(Vxl(i)))/Vbw(i))^(kl)); pl(i)=pb(i)*(((Vbr(i)-(Vxr(i)))/Vbr(i))^(kr)); p(i)=((Vbw(i)*((pr(i))^(1/kl))+((pl(i)/((pr(i)+pl(i))/2))^((kl-kr)/(kl*kr)))*Vbr(i)*((pl(i))^(1/kl)))/V)^kl; end end plot(t,p) |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有23人回复
-
导师想让我从独立一作变成了共一第一
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19571.1
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4938
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
7楼2016-12-17 13:01:04