走过路过的大神你们好，麻烦帮我看一下这个si指数的计算，弄了好长一段时间了 还...

秋刀鱼。

我算了一下前面都是正确的，就是假相当位温循环的地方有问题。下面是gs文件：
********假相当位温迭代循环
n=10000
while(n>=0)
'define t3=tlcl-0.0001* 'n''       tlcl是抬升凝结高度的温度。
'define tk=t3+273.15'
'define est=6.112*exp((17.67*t3)/(t3+243.5))'
'define etd=6.112*exp((17.67*td2)/(td2+243.5))'
'define e=etd'
'define u=(etd/est)*100'
'define tl=1/(1/(tk-55)-log(u/100)/2840)+55'
'define r=622*(e/(p2-e))'
*'define tl=55.0+2840.0/(3.5*log(tk)-log(e)-4.805)'
'define a2=tk*pow((1000/p2),(0.2854*(1.0-0.00028*r)))*exp((3.376/tl-0.00254)*r*(1+0.00081*r))'
'define a=a2-a1'
'define k=abs(a)'
  if(k<0.001);break;endif;
    n=n+1
  endwhile
'define si=t500-t3'
然后就弹出下面的问题：


求大神指导

男紫汗

对数计算的时候有负值

秋刀鱼。

男紫汗 发表于 2017-7-9 16:44
对数计算的时候有负值

可是etd和est都是正值的嘛 ，咋会算出来有负值？我之前也以为是不是log有问题 可是我后来特意去百度warning from log 好像是警告日志。如果n设置比较小 就没有这个问题 n比较大就有

男紫汗

秋刀鱼。 发表于 2017-7-10 09:21
可是etd和est都是正值的嘛 ，咋会算出来有负值？我之前也以为是不是log有问题 可是我后来特意去百度warni ...

我之前算位温也遇到过。。。怎么回事也忘了。

秋刀鱼。

男紫汗 发表于 2017-7-10 11:14
我之前算位温也遇到过。。。怎么回事也忘了。

那是什么问题啊  麻烦你想起了一定告诉我，谢谢。。哎、、、、

南南偶卖糕Oo

n=10000会不会不太对，你试试n=0开始