e271.m

matlab算法集 matlab算法集

%-----------------------------------------------------------------------
% Example 2.7.1: Chemical Absorption Process 
%-----------------------------------------------------------------------

% Initialize 

   clc                 % clear command window
   clear               % clear variables
   n = 5;              % number of stages
   a = 0.72;           % slope 	
   b = 0.0;            % intercept 
   K = 1.63;           % gas to liquid ratio  
   xf = 0.01;          % liquid feed concentration  
   yf = 0.06;          % gas feed concentration 
   delta = a*K;        % stripping factor 
   A = zeros (n,n);
   c = zeros (n,1);
   
% Compute A and c 

   fprintf ('Example 2.7.1: Chemical Absorption Process\n');
   for i = 1 : n
      A(i,i) = -(1 + delta);
      if i < n
         A(i,i+1) = 1;
         A(i+1,i) = delta;
      end
   end
   c(1) = K*(b - yf);
   c(n) = -xf;
   show ('A',A)
   show ('c',c)
   
% Find liquid concentrations 

   show ('K(A)',condnum(A,0))
   x = gauss (A,c); 
   show ('x',x)
   y = a*x + b;
   show ('y',y)   
   show ('||r(x)||',residual(A,c,x)) 
%-----------------------------------------------------------------------