水和水蒸气热力性质计算程序代码matlab中调用

该程序组共包含10个子程序，全部为计算水和水蒸汽性质的子程序，采用的是国际公式化委员会制定的水和水蒸气热力性质（IFC67）公式。各子程序分别是： TSK.m 求某下压力饱和温度。 PSK.m 某温度下饱和压力。 HS.m 已知比焓、比熵，求其它性质。 PX.m 已知压力、干度，求其它性质。 PV.m 已知压力、比熵，求其它性质。 PTG.m 已知压力、温度，求饱和汽、过热蒸汽的性质。 PTF.m 已知压力、温度，求饱和水、过冷水的性质。 PT.m 已知压力、温度，求其它性质。 PS.m 已知压力、比熵，求其它性质。 PH.m 已知压力、比焓，求其它性质 该程序是汽轮机设计和热力系统设计的好帮手，减少了查表的麻烦。

function [X PT V] = HS（H S）
%************************************************
% 已知比焓、比熵求过热蒸汽、饱和蒸汽、湿蒸汽性质函数
%  Parameters: P      压力 MPa 
%              T      温度  
%              V      比容  
%              H      焓  
%              S      熵  
%              X      干度（过热度）
% ****作者：王雷 zrqwl2003@126.com*****************
% ************************************************
G=[-32230.91316763.0198-2684.70102126.0153826.08691173-0.515969833];
w=[0.08275963 -1.2121626 -73.614096.6074894-0.4065599];
HG = 0;
for i = 1 :1: 6
    HG = HG + G（i） * S ^ （i - 1）;
end
if （H > HG） 
    X = 1;
    ZH = H / 4186.8;
    ZS = 0.238846 * S;
    n = 1;
    while（1）
        ZT = （sqrt（0.035081 + 0.085 * （ZH - 0.4949）） - 0.1873） / 0.085;
        ZP = exp（13.064 - 9.072 * ZS + 3.3982 * log（ZT））;
        T = 1000 * ZT - 273.15;
        P = 98.0665 * ZP;
        TS = TSK（P）;
        if （T >= TS） 
            ZHB = ZH;
            ZSB = ZS;
        else
            T = TS;
            ZT = （T + 273.15） / 1000;
            ZHB = （（0.085 * ZT + 0.1873） ^ 2 - 0.035081） / 0.085 + 0.4949;
            ZSB = （13.064 - log（ZP） + 3.3982 * log（ZT）） / 9.072;
        end
        X = T - TS;
        if （X >= 1）
        else
            X = X + 1;
        end
        [V HB SB]=PTG（P T ）;
        if （（abs（H - HB） < 0.0005） && （abs（S - SB） < 0.000001））
            break;
        end
        if （n > 1）
            ZH = ZHB + （H - HB） * （ZHB - ZHA） / （HB - HA）;
            ZS = ZSB + （S - SB） * （ZSB - ZSA） / （SB - SA） ;
        else
            ZH = ZHB + 0.0002388 * （H - HB）;
            ZS = ZSB + 0.2388 * （S - SB）;
        end
        n = n + 1;
        HA = HB;
        SA = SB;
        ZHA = ZHB;
        ZSA = ZSB;
        ZHB = ZH;
        ZSB = ZS;
    end
else
    TB = H / S;
    while（1）
        GF = 0;
        for i = 1 :1: 5
            GF = GF + w（i） * （TB * 0.01） ^ （i - 1）;
        end
        T = （H - GF） / S - 273.15;
        if （abs（T - TB） < 0.001）
            break;
        end
        TB = T;
    end
    P = PSK（T）;
    [VG HG SG]=PTG（PT）;
    [VF HF SF]=PTF（PT）;
    X = （H - HF） / （HG - HF）;
    V = VF + X * （VG - VF）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2276  2013-03-11 18:28  IFC67\HS.m

     文件       1901  2012-07-13 21:31  IFC67\PH.m

     文件       2050  2012-07-13 21:31  IFC67\PS.m

     文件        564  2012-07-13 21:31  IFC67\PSK.m

     文件        692  2012-07-13 21:32  IFC67\PT.m

     文件       2181  2012-07-13 21:32  IFC67\PTF.m

     文件       4353  2012-07-13 21:32  IFC67\PTG.m

     文件       1097  2012-07-13 21:32  IFC67\PV.m

     文件        555  2012-07-13 21:32  IFC67\PX.m

     文件        456  2012-07-13 21:32  IFC67\TSK.m

     目录          0  2014-04-29 15:35  IFC67

----------- ---------  ---------- -----  ----

                16125                    11