- 积分
- 5
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-7-20
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
%已知干球温度和相对湿度,求湿球温度
%该程序仅适用于常压、0-90℃范围内,通过迭代求得不是非常精确的湿球温度解
P=1000;%对应环境压力,hPa
Td=298;%干球温度,K
RH=0.5;%相对湿度,0-1
Psd=10^(10.11564-1687.54/(Td+230.23-273.15));%干球温度对应的饱和蒸汽压力,Pa
d=0.622*RH*Psd/(P*100-RH*Psd);%干球温度条件下含湿量,kg/kg(dry air)
h=1.01*Td+(2500+1.84*Td)*d;%干球温度对应焓值,KJ/kg
Tw=250;%假定一个湿球温度,K
hw=h;%湿球温度焓值=干球温度焓值
i=1;
Twk=Tw;
j=1;
while i>10^-7
Tw=Twk;
dw=(hw-1.01*Tw)/(2500+1.84*Tw);%湿球温度条件下含湿量,kg/kg(dry air)
Psw=dw*P*100/(dw+0.622);%湿球温度对应的饱和蒸汽压力,Pa
Twk=42.92+1687.54/(5.11564-log10(Psw*10^(-5)));%迭代出的新湿球温度,K
i=abs(Twk-Tw)/Tw;
j=j+1;
end
fprintf('湿球温度是 %g K\n',Tw);
|
|