| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 382 人关注过本帖
标题:高手们看看哪里出错了
只看楼主 加入收藏
lin5165
Rank: 1
等 级:新手上路
帖 子:14
专家分:0
注 册:2010-11-5
结帖率:0
收藏
已结贴  问题点数:10 回复次数:3 
高手们看看哪里出错了
#include"math.h"
 #include"stdio.h"
 main()
 {
 int i,j,k,n;
 static double a[5]={0.92,0.6,0.58,0.41,0.10};
 static double b[5]={0.7,0.58,0.43,0.35,0.228};
 static double d[5]={0.345,0.005,0.065,0.115,0.230};
  double t[7],xt[7],ta,aw,r[5],aaw,la[5],cc,eps,qloss;
  xt[0]=1720.0;  t[0]=1720; eps=0.0001; ta=293;
  aaw=2.56;
 for(i=1;i<7;i++)  t[i]=720;
 j=0;
 do{
 j=j+1;
 for(i=0;i<5;i++)
  la[i]=a[i]+b[i]*(t[i]+t[i+1])/2/1000;
   for(i=0;i<5;i++)
 r[i]=la[i]/d[i];
 for(i=1;i<5;i++)
xt[i]=r[i-1]/(r[i-1]+r[i])*xt[i-1]+r[i]/(r[i-1]+r[i])*t[i+1];
 aw=aaw*pow((t[5]-ta),0.25)+4.54*(pow((t[5]/100),4)-pow((ta/100),4))/(t[5]-ta
 xt[5]=r[4]/(aw+r[4])*xt[4]+aw/(aw+r[4])*ta;
 k=0; for(i=1;i<6;i++)
 for(i=1;i<6;i++)
{ cc=abs(xt[i]-t[i]);
 if(cc>eps) k=k+1;
  }
 for(i=0;i<6;i++)
 t[i]=xt[i];
 }while(k>0);
  qloss=la[0]/d[0]*(t[0]-t[1]);
  for(i=0;i<6;i++)
  printf("tt(%d)=%e\n",i,xt[i]);
  printf("qloss=%e,%d\n",qloss,j);
  }
高手们指点一下,谢谢啦
用这个程序应该能计算出
T[0]=1720
T[1]=1478
T[2]=1474
T[3]=1413
T[4]=1267
T[5]=371
搜索更多相关主题的帖子: double eps include 
2010-11-06 23:37
lin5165
Rank: 1
等 级:新手上路
帖 子:14
专家分:0
注 册:2010-11-5
收藏
得分:0 
没有人么
2010-11-07 12:16
vandychan
Rank: 15Rank: 15Rank: 15Rank: 15Rank: 15
等 级:贵宾
威 望:18
帖 子:2296
专家分:6418
注 册:2010-8-20
收藏
得分:5 
you a 有啊

到底是“出来混迟早要还”还是“杀人放火金腰带”?
2010-11-08 00:57
caokewen
Rank: 1
等 级:新手上路
帖 子:2
专家分:5
注 册:2010-11-8
收藏
得分:5 
回复 楼主 lin5165
aw=aaw*pow((t[5]-ta),0.25)+4.54*(pow((t[5]/100),4)-pow((ta/100),4))/(t[5]-ta);
2010-11-08 13:31
快速回复:高手们看看哪里出错了
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.011797 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved