单效溴化锂吸收式制冷机热力仿真matlab程序

function h = EqLiBrH2O(x,t)
% 
if nargin<2||isempty(t),error('输入参数：缺少温度');end
if nargin<1||isempty(x),error('输入参数：缺少浓度');end

if (x>40 && x<70) %---检查浓度范围
    if (t>15 && t<165) %---检查温度范围
        %---基本参数
        A0 = -2024.33; A1 = 163.309; A2 = -4.88161; A3 = 6.302948*10^(-2); A4 = -2.913705*10^(-4);
        B0 = 18.2829; B1 = -1.1691757; B2 = 3.248041*10^(-2); B3 = -4.034184*10^(-4); B4 = 1.8520569*10^(-6);
        C0 = -3.7008214*10^(-2); C1 = 2.8877666*10^(-3); C2 = -8.1313015*10^(-5); C3 = 9.9116628*10^(-7); C4 = -4.4441201*10^(-9);
        %---创建多项式系数
        P0 = A0 + t*B0 + t*t*C0;
        P1 = A1 + t*B1 + t*t*C1;
        P2 = A2 + t*B2 + t*t*C2;
        P3 = A3 + t*B3 + t*t*C3;
        P4 = A4 + t*B4 + t*t*C4;
        %---计算焓值
        P = [P4 P3 P2 P1 P0];
        h = polyval(P,x);
    else
        disp('错误：温度超过范围: t>15 && t<165');
    end
else
    disp('错误：溴化锂浓度超出范围: x>40 && x<70');
end