热力系统计算模型

function a=atfun（）
%给水压力和凝结水压力以及水的比容
pfw=31.7;pcp=3.1;vv=0.0011;
%凝器参数
pc=0.0049;hc=2300.5;
[pctchcscvcxcrc]=WASPCNMEX（‘ph‘pchc）;
%主汽参数
p0=24.2;t0=566;
[p0t0h0s0v0x0r0]=WASPCNMEX（‘pt‘p0t0）;
%抽汽口参数
p=[6.687 4.434 2.133 1.069 0.363 0.197 0.099 0.045];
t=[370.9 315.9 468.8 367.5 234.6 168.8 104.1 0];
for i=1:8
    if i==8
        h（i）=2562.2;[p（i）t（i）h（i）s（i）v（i）x（i）r（i）]=WASPCNMEX（‘ph‘p（i）h（i））;
    else
        [p（i）t（i）h（i）s（i）v（i）x（i）r（i）]=WASPCNMEX（‘pt‘p（i）t（i））;
    end
end

%再热入口
prhi=p（2）;trhi=t（2）;
[prhitrhihrhisrhivrhixrhirrhi]=WASPCNMEX（‘pt‘prhitrhi）;
%再热出口
prho=3.99;trho=566;
[prhotrhohrhosrhovrhoxrhorrho]=WASPCNMEX（‘pt‘prhotrho）;
%中低压连通管参数
pil=p（4）;til=t（4）;
[piltilhilsilvilxilril]=WASPCNMEX（‘pt‘piltil）;
%上端差
dtu=[-1.7 0 0 0 2.8 2.8 2.8 2.8];
%下端差
dtd=[5.6 5.6 5.6 0 5.6 5.6 5.6 5.6];
%凝汽器出口，凝泵入口参数
[pcotcohcoscovcoxcorco]=WASPCNMEX（‘px‘pc0）;
%1H出口参数：给水参数
tfw=WASPCNMEX（‘p2t‘p（1））-dtu（1）;
[pfwtfwhfwsfwvfwxfwrfw]=WASPCNMEX（‘pt‘pfwtfw）;
%各给水入口压力
for i=1:8
    if i<4
        pw（i）=pfw;
    else
        pw（i）=pcp;
    end
end
%给水泵泵功，给水泵处给水焓升
tpu=vv*（pfw-p（4））*1000/0.9;
%凝结水泵泵功，凝结水泵处凝结水焓升
tcp=vv*（pcp-pco）*1000/0.9;
%各加热器入口给水参数
for i=1:8
    if i==8
        hw（i）=tcp+hco;
        [pw（i）tw（i）hw（i）sw（i）vw（i）xw（i）rw（i）]=WASPCNMEX（‘ph‘pw（i）hw（i））;
    elseif i==3
        hw（i）=WASPCNMEX（‘px2h‘p（i+1）0）+tpu;
        [pw（i）tw（i）hw（i）sw（i）vw（i）xw（i）rw（i）]=WASPCNMEX（‘ph‘pw（i）hw（i））;
    else
        tw（i）=WASPCNMEX（‘p2t‘p（i+1））-dtu（i+1）;
        [pw（i）tw（i）hw（i）sw（i）vw（i）xw（i）rw（i）]=WASPCNMEX（‘pt‘pw（i）tw（i））;
    end
end
%疏水参数
for i=1:8
    td（i）=tw（i）+dtd（i）;
end
for i=1:8
    if i==4
        hd（i）=hw（i）;
    else
        [pd（i）td（i）hd（i）sd（i）vd（i）xd（i）rd（i）]=WASPCNMEX（‘pt‘p（i）td（i））;
    end
end
%给水泵前给水参数
[pw（3）tw（3）hw（3）sw（3）vw（3）xw（3）rw（3）]=WASPCNMEX（‘px‘p（4）0）;
%确定tao值
for i=1:8
    if i==1
        tao（i）=hfw-hw（i）;
    else
        tao（i