| 
 
	积分1860贡献 精华在线时间 小时注册时间2012-11-19最后登录1970-1-1 
 | 
 
| 
本帖最后由 爱的侍者 于 2013-1-29 21:06 编辑
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  
 用matlab计算结果正确
 用fortran计算结果不对
 贴出一段程序看一下为啥呢
 -51.76    6.7    6.6    0.6   -0.3
 这是fortran算出来的
 第一个为结果其他四个 分别为sumu(i,1)sumu(i,4)stmt(i,1)stmt(i,2)
 
 if(stmt(i,1).ne.99999.and.stmt(i,2).ne.99999.and.sumu(i,1).ne.99999.and.sumu(i,4).ne.99999)then
 ri(i)=(9.8/(stmt(i,1)+273.15))*((stmt(i,2)-stmt(i,1))/60+0.0065)/((sumu(i,4)-sumu(i,1))/60)**2
 else
 ri(i)=99999
 endif
 write(25,'(f20.1,1x,f7.1,1x,f7.1,1x,f7.1,1x,f7.1)')ri(i),sumu(i,2),sumu(i,4),stmt(i,1),stmt(i,2)
 
 这样出了问题
 
 正确结果
 >> (9.8/(0.6+273.15))*((-0.3-0.6)/60+0.0065)/(((6.6-6.7)/60).^2)
 ans =
 -109.5452
 应该是这个
 为啥呢
 
 | 
 |