查看: 106|回复: 1|关注: 0

[已答复] if语句结果为真结果不执行 这是为什么呢

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 2020-7-1 16:25:34 | 显示全部楼层 |阅读模式
本帖最后由 yjy15232 于 2020-7-1 16:28 编辑
  1. clear allsigma=0.4;Vliquid=10^-3;n=1000;%%催化剂数量
  2. M=100;tmax=50000;N=100;rmax=0.0025;tao=tmax/M;l=rmax/N;k2=0.207;k=0.21*(1-sigma)+k2*sigma;
  3. cp2=2607.9;cp=1427.5*(1-sigma)+cp2*sigma;
  4. rou2=791.8;rou=2200*(1-sigma)+rou2*sigma;alpha=k/(rou*cp);
  5. d=2*rmax;u=0.5;miu=0.5525*10^-3;
  6. Re=rou2*d*u/miu;
  7. Pr=cp2*miu/k2;
  8. h=(k2/d)*(2+((0.4*(Re^0.5)+0.06*(Re^0.67))*(Pr^0.4)));
  9. t=0:tao:tmax;
  10. r=0:l:rmax;
  11. T=zeros(N+1,M+1);
  12. T(:,:)=40;a=1;p=zeros(1,101);
  13. %%初始条件,本行后循环开始
  14. while a<=101
  15.     if a>=2
  16. T0=p(a-1);
  17.     else
  18.         T0=40;
  19.     end
  20. Q=25.816*(1.634e-5*exp(-(a-1)*tao/1743)+3.314e-5*exp(-(a-1)*tao/3468.3)+4.743e-5*exp(-(a-1)*tao/7502.8)...
  21.     +5.016e-9);
  22. T(N+1,a)=T0+5*Q*(1-exp(-2*alpha*h.*(a-1)*tao/(rmax*k)))/(8*rmax^2*h*pi);
  23. b=N+1;
  24. q=Q/(4*pi*rmax^3/3);
  25. while b>=2
  26.         T(b-1,a)=T(b,a)+l*(rmax^2/((b-1)*l)^2)*(h*(T(N+1,a)-T0))/k-l*q*(4*pi*(rmax^3-((b-1)*l)^3)/3)/(4*k*pi*((b-1)*l)^2);
  27.         b=b-1;
  28. end
  29. for b=1:N+1
  30.         if T(b,a)>40+q*(a-1)*tao/(rou*cp)
  31.             T(b,a)=40+q*(a-1)*tao/(rou*cp);
  32.         end
  33. end
  34. QQ1=0;
  35. for b=1:N+1
  36.     QQ1=QQ1+rou*cp*4*pi*((b-1)*l)^2*l*(T(b,a)-40);
  37. end
  38. QQ1=n*QQ1;
  39. QQ=n*25.816*(0.5-0.02375*exp(-(a-1)*tao/1639.4)-0.3487*exp(-(a-1)*tao/7620.8)-0.1276*exp(-(a-1)*tao/3472.9));
  40. QQ2=QQ-QQ1;
  41. T1=40+QQ2/(rou2*cp2*Vliquid);
  42. if (abs(T0-T1)<0.01)
  43.     p(a)=T0
  44.     a=a+1;
  45. else T0=T1;
  46. end
  47. end
复制代码


if语句我单独验算过了 结果为真 最后T1和T0的结果也是收敛的 结果没有执行陷入了死循环,p数组也没有数据输出,求教这是为什么?万分感谢

回复主题 已获打赏: 0 积分

举报

新手

5 麦片

财富积分


050


0

主题

1

帖子

0

最佳答案
发表于 2020-7-2 12:45:43 | 显示全部楼层
你的if判断语句if (abs(T0-T1)<0.01),我这边运行的结果T0-T1=-0.0102, 执行else后面的语句,然后T0=T1,a值没有变,后面循环之后一直是同样a值,同样的T0-T1结果,当然就是死循环了。
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表