- 积分
- 1834
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-11-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 爱的侍者 于 2013-1-29 21:06 编辑
用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
应该是这个
为啥呢
|
|