function?a=atfun（）
%給水壓力和凝結(jié)水壓力以及水的比容
pfw=31.7;pcp=3.1;vv=0.0011;
%凝器參數(shù)
pc=0.0049;hc=2300.5;
[pctchcscvcxcrc]=WASPCNMEX（‘ph‘pchc）;
%主汽參數(shù)
p0=24.2;t0=566;
[p0t0h0s0v0x0r0]=WASPCNMEX（‘pt‘p0t0）;
%抽汽口參數(shù)
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）;
%中低壓連通管參數(shù)
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];
%凝汽器出口，凝泵入口參數(shù)
[pcotcohcoscovcoxcorco]=WASPCNMEX（‘px‘pc0）;
%1H出口參數(shù)：給水參數(shù)
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;
%凝結(jié)水泵泵功，凝結(jié)水泵處凝結(jié)水焓升
tcp=vv*（pcp-pco）*1000/0.9;
%各加熱器入口給水參數(shù)
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
%疏水參數(shù)
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
%給水泵前給水參數(shù)
[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