Hsrh1=2614.86; %汽水分离再热器第一级抽气焓
Hzs1=1018.9756; %汽水分离再热器第一级疏水焓
Hsrh2=2783.348; %汽水分离再热器第二级抽气焓
Hzs2=1219.8773; %汽水分离再热器第二级疏水焓
Hfw=707.227; %分离器分离水比焓
h1=99.46; %高低压加热器中的平均焓升
h2=117.6; %再热器中的平均焓升
Hfwi=[146.6448 246.10 345.56 445.02 544.48 643.94 743.40]; %各级加热器进口焓（kj/kg）
Hfwo=[246.10 345.56 445.02 544.48 643.94 743.40 842.86]; %各级加热器出口焓(kj/kg) Hc=[2377.40 2525.49 2666.18 2803.096 2480.12 2448.787 2531.89]; %每一级的给水加热器的加热蒸汽比焓
Hw=[254.06 353.61 453.235 552.88 0 753.57 854.18]; %每一级给水加热器疏水比焓 H=100;
for Gslp=100:3000
if abs(H-Hfwo(5))/Hfwo(5)<0.005
H=(X*Hfwo(4)+(Ghes(1)+Ghes(2)+Gzc1+Gzc2)*Hw(6)+Gsdea*Hho+Hfw*Gslp*(xz1-xhz)/xhz)/Gfw; break
else
end
for Gfw=Gslp:3000;
Nfwpp=1000*Gfw*Hfwp/pw; %给水泵有效功率
Nfwpt=Nfwpp/(nfwpp*nfwptm*nfwptg); %给水泵实际功率
Gsfwp=Nfwpt/(Hfh-Hho); %给水泵耗气量
Ds=Gfw/(1+Ed); Gd=Gfw-Ds; %Ds为新蒸汽气量，Gd为排污量
X=Gslp+Gd+Gsfwp; %假设给水量
Gles(4)=h1*X/(Hc(4)-Hw(4))/nh;
Gles(3)=(h1*X-nh*Gles(4)*(Hw(4)-Hw(3)))/(Hc(3)-Hw(3))/nh;
Gles(2)=(h1*X-nh*(Gles(3)+Gles(4))*(Hw(3)-Hw(2)))/(Hc(2)-Hw(2))/nh;
Gles(1)=(h1*X-nh*(Gles(2)+Gles(3)+Gles(4))*(Hw(2)-Hw(1)))/(Hc(1)-Hw(1 ))/nh;
w1=(x-Gles(1)-Gles(2)-Gles(3)-Gles(4))*(Hli-Hlo)+Gles(1)*(Hli-Hc(1))+Gles(2)*(Hli-Hc(2))+Gles(3)*(Hli-Hc(3))+Gles(4)*(Hli-Hc(4));
W1=w1*nm*nge;
W2=10^6-W1; %上述程序Gles指的是低压给水再热器抽气，括号内数字为级数。W1和W2为高低压缸发电量
Gsrh1=Gslp*h2/(nh*(Hsrh1-Hzs1)); %第一级再热抽汽量
Gsrh2=Gslp*h2/(nh*(Hsrh2-Hzs2)); %第二级再热抽汽量
Ghes(2)=(Gfw*h1-nh*Gsrh2*(Hzs2-Hw(7)))/(nh*(Hc(7)-Hw(7)));
Ghes(1)=(Gfw*h1-nh*((Ghes(2)+Gsrh2)*(Hw(7)-Hw(6))+Gsrh1*(Hzs1-Hw(6))))/(nh*(Hc(6)-Hw(6)));
Gsdea=Gfw-X-Ghes(1)-Ghes(2)-Gsrh1-Gsrh2-Gslp*(xz1-xhz)/xhz;
Gt=Gsdea+Gslp*xz1/xhz; %上述Ghes为高压给水再热器抽气量，括号内数字为级数。Gsdea为除氧器耗气量
w2=Gt*(Hfh-Hho)+Ghes(1)*(Hfh-Hc(6))+Ghes(2)*(Hfh-Hc(7))+Gsrh1*(Hfh-Hzc1);