MATLAB 水和水蒸气计算函数

matlab语言；水和水蒸气参数计算函数，多个m文件，可计算过热

function [X T V S] = PH（P H）
%***************************************
% 已知压力比焓求性质函数
%  Parameters: P      压力   MPa 
%              t      温度   ℃ 
%              V      比容   m^3/kg 
%              H      比焓   kJ/kg
%              S      比熵   kJ/（kg.℃ ）
%              x      干度  　<1（过热度℃）
%****作者：王雷 zrqwl2003@126.com*****************  
%***************************************
TS = TSK（P）;
HA=H;
[VG HG SG]=PTG（P TS）;
[VF HF SF]=PTF（P TS）;
ZHA=0;
 if （H >= HF） 
     if （H > HG）
        x = 1;
        n = 1;
        ZH = H / 4186.8;
        while（1）
            ZT = （sqrt（0.035081 + 0.085 * （ZH - 0.4949）） - 0.1873） / 0.085;
            T = 1000 * ZT - 273.15;
            if （T < TS）
                T = TS;
                HB = HG;
                V = VG;
                S = SG;
                ZT = （273.15 + TS） / 1000; 
                ZHB=0.08218+（0.085*ZT+0.1873）^2/0.085;
            else
                [V HB S]= PTG（P T）;
                ZHB = ZH;           
            end
          %  ZHA = ZHB;
            X = T - TS;
            if （X >= 1）
            else
                X = X + 1;
            end

            if （abs（H - HB） < 0.001）
                break;
            end
            if（n > 1）
                ZH = ZHB + （H - HB） * （ZHB - ZHA） / （HB - HA）;
            else
                ZH = 0.000238 * （H - HB） + ZHB;
            end
            n = n + 1;
            HA = HB;
            ZHA = ZHB;
            ZHB = ZH;
        end
     else
        T = TS;
        X = （H - HF） / （HG - HF）;
        V = VF + X * （VG - VF）;
        S = SF + X * （SG - SF）;
     end
 else
    X = 0;
    T = 0.2195 * HA;
    while（1）
        [V HB S]= PTF（P T）;
        X = T - TS;
        if （abs（H - HB） < 0.001）
            break;
        end
        T = T + 0.238846 * （H - HB）;
    end
 end
        if  X>1
            X=1;
        end
        if X<0
            X=0;
        end
end

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

     文件       2034  2017-05-15 15:47  新建文件夹 （2）\PH.m

     文件       2150  2017-03-31 21:07  新建文件夹 （2）\PS.m

     文件        569  2017-03-30 16:38  新建文件夹 （2）\PSK.m

     文件        270  2018-03-17 10:39  新建文件夹 （2）\PSU.m

     文件        763  2017-04-01 17:51  新建文件夹 （2）\PT.m

     文件       2186  2017-04-01 19:54  新建文件夹 （2）\PTF.m

     文件       4358  2017-03-30 16:37  新建文件夹 （2）\PTG.m

     文件       1197  2017-03-31 21:07  新建文件夹 （2）\PV.m

     文件        560  2017-03-30 16:37  新建文件夹 （2）\PX.m

     文件        456  2012-07-13 21:32  新建文件夹 （2）\TSK.m

     目录          0  2018-06-19 09:53  新建文件夹 （2）

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

                14543                    11