{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

BME_LAB_03_B

# BME_LAB_03_B - F_1=0.6 F_2=0.2 The mass flow rate of mixer...

This preview shows pages 1–2. Sign up to view the full content.

F_1=0.6; F_2=0.2; % The mass flow rate of mixer 1 is 1 % x is the mass flow rate of feed 1 % y is the mass flow rate of feed 2 % 0.6x+0.2y=0.5 % x+y=1 a=[F_1 F_2; 1 1]; b=[0.5;1]; c=inv(a)*b; fprintf('The mass flow rate of alloy feed 1 is: %4.3f\n',c(1,:)) fprintf('The mass flow rate of alloy feed 2 is: %4.3f\n',c(2,:)) d=[0.6 0.2 ;0.2 0.6;0.2 0; 0 0.2]; e=d*c; fprintf('The mass flow rate of F in mixer 1 is: %4.3f\n',e(1,:)) fprintf('The mass flow rate of G in mixer 1 is: %4.3f\n',e(2,:)) fprintf('The mass flow rate of H in mixer 1 is: %4.3f\n',e(3,:)) fprintf('The mass flow rate of K in mixer 1 is: %4.3f\n',e(4,:)) % G % 0.3+0 = P_m4; % H % 0.15+0.5*x_m5= P_m4; % K % 0.05+(m4-1)-1.5x_m5=m4*0.6-2p_m4 % % p_m4 = 0.3 % p_m4-0.5x_m5=0.15 % 2p_m4+0.4m4-1.5x_m5=0.95 coe=[1 0 0;1 -0.5 0;2 -1.5 0.4]; res=[0.3;0.15;0.95]; unk=inv(coe)*res; u fprintf('P_m4, the mass of compound G at the outflow of mixer 2 is: %4.3f\n',unk(1,:)) fprintf('X.m5, the mass inflow of compound F in alloy feed 5 is: %4.3f\n',unk(2,:)) fprintf('m4, the outflow mass rate of mixer 2 is: %4.3f\n',unk(3,:)) T=unk(3,1)-2*unk(1,1)-0.4*unk(3,1); fprintf('k, the mass of compound K at the outflow of mixer 2 is: %4.3f\n',T) m5=unk(3,1)-1;

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}