for(i=1:size(y3,2))
y4(i)=N3*(y3(i)*K(i)/(A+K(i)))/N4; y5(i)=N3*(y3(i)/(A+K(i)))/N4; end
%disp('help:'); %disp(y4(1)); % disp(y4(2)); % disp(y4(3)); % disp(y4(4)); % disp(y4(5)); %disp('over'); N4=N3/(1+A); N5=N3-N4;
y4(1)=N3*(y3(1)*K(1)/(A+K(1)))/N4; %disp('y41'); %disp(y4(1));
y4(2)=N3*(y3(2)*K(1)/(A+K(2)))/N4; y4(3)=N3*(y3(3)*K(1)/(A+K(3)))/N4; y4(4)=N3*(y3(4)*K(1)/(A+K(3)))/N4; y4(5)=N3*(y3(5)*K(1)/(A+K(3)))/N4; y5(1)=N3*(y3(1)/(A+K(1)))/N4; y5(2)=N3*(y3(2)/(A+K(2)))/N4; y5(3)=N3*(y3(3)/(A+K(3)))/N4; y5(4)=N3*(y3(4)/(A+K(3)))/N4; y5(5)=N3*(y3(5)/(A+K(3)))/N4;
function fun=forigin(A,y31,y32,y33,y34,y35,k1,k2,k3,k4,k5);
fun=y31*(1-k1)/(A+k1)+y32*(1-k2)/(A+k2)+y33*(1-k3)/(A+k3)+y34 *(1-k4)/(A+k4)+y35*(1-k5)/(A+k5);
function fun1=fnext(A,y31,y32,y33,y34,y35,k1,k2,k3,k4,k5)
fun1=-y31*(1-k1)/(A+k1)^2-y32*(1-k2)/(A+k2)^2-y33*(1-k3)/(A+k3)^2-y34 *(1-k4)/(A+k4)^2-y35*(1-k5)/(A+k5)^2;
%分离器衡算（假设使用理想分离器计算） function [y4 N4 y5 N5]=f3(y3,N2,N3); %syms y31 y32 y33 y34 y35; y5=[1,0,0,0,0]; N5=N3*y3(1); N4=N3-N5; y4(1)=0;
y4(2)=N3*y3(2)/N4; y4(3)=N3*y3(3)/N4;
21
y4(4)=N3*y3(4)/N4; y4(5)=N3*y3(5)/N4; disp('N4') disp(N4) disp('y4') disp(y4)
%%%%放空%%%%%
function [y9 N9 y10 N10 y6 N6 y7 N7 y8 N8]=f4(y4,N4) N10=N4*12.5/827.1; y10=y4;
%N7=N10*y10(2); %N6=N10-N7; N7=0.16;
y7=[0,1,0,0,0]; N6=N10-N7; y6(1)=0; y6(2)=0;
y6(3)=N10*y10(3)/N6; y6(4)=N10*y10(4)/N6; y6(5)=N10*y10(5)/N6; N8=N4-N10; y8=y4;
N9=N8+N7; y9(1)=0;
y9(2)=(N8*y8(2)+N7*y7(2))/N9; y9(3)=N8*y8(3)/N9; y9(4)=N8*y8(4)/N9; y9(5)=N8*y8(5)/N9; disp('N7'); disp(N7);
%%%%提氢放空%%%%%
function [y9 N9 y10 N10 y6 N6 y7 N7 y8 N8]=f4(y4,N4) N10=N4*12.5/827.1; y10=y4;
%N7=N10*y10(2); %N6=N10-N7; N7=0.16;
y7=[0,1,0,0,0]; N6=N10-N7; y6(1)=0; y6(2)=0;
y6(3)=N10*y10(3)/N6; y6(4)=N10*y10(4)/N6; y6(5)=N10*y10(5)/N6;
22
N8=N4-N10; y8=y4;
N9=N8+N7; y9(1)=0;
y9(2)=(N8*y8(2)+N7*y7(2))/N9; y9(3)=N8*y8(3)/N9; y9(4)=N8*y8(4)/N9; y9(5)=N8*y8(5)/N9; disp('N7'); disp(N7);
23