clc    ;
echo on ;
% ****************************************************************************
%  Demo 941.m  This program builds a state space model  of the Shell Column
%  Copyright 1992 by Babu Joseph
%  Washington University in St.Louis
%******************************************************************************
% All time delays have been replaced by (1,1) Pade approximations
%
% We start by entering the transfer functions and getting step response coeffs
% input 1 is top draw
% input 2 is side draw
% input 3 is bottoms duty
% input 4 is the intermediate reflux duty
% input 5 is the upper reflux duty
% output 1 is top end point
% output 2 is side end point
% output 3 is bottoms reflux temp
% the program returns a state space model of the form [ac,bc,cc,dc]
pause % press any key
clc ; echo off
g11num=[ -54.68 4.05];
g11den=[ 675 63.5 1];
g12num=[-48.51 5.39];
g12den=[450 59 1];
g13num=[-43.8 4.38];
g13den=[330 43 1];
g21num=[-24.78 1.77];
g21den=[840 74 1];
g22num=[-40 5.72];
g22den=[420 67 1];
g23num=[-48 4.42];
g23den=[484 55 1];
g31num=[-79.38 5.88];
g31den=[675 63.5 1];;
g32num=[-51.75 6.9];
g32den=[300 47.5 1];
g33num=[0 7.2];
g33den=[19 1];
g41num=[-16.25 1.2];
g41den=[607.5 58.5 1];
g42num=[-11.4 1.52];
g42den=[187.5 32.5 1];
g43num=[0 1.14];
g43den=[27 1];
g51num=[-19.44 1.44];
g51den=[540 53.5 1];
g52num=[-13.73 1.83];
g52den=[150 27.5 1];
g53num=[0 1.26];
g53den=[32 1];
[a1,b1,c1,d1]=tf2ss(g11num,g11den);
[a2,b2,c2,d2]=tf2ss(g12num,g12den);
[a3,b3,c3,d3]=tf2ss(g13num,g13den);
[a4,b4 ,c4 ,d4 ]=tf2ss(g21num,g21den);
[a5 ,b5 ,c5 ,d5 ]=tf2ss(g22num,g22den);
[a6 ,b6 ,c6 ,d6 ]=tf2ss(g23num,g23den);
[a7 ,b7 ,c7 ,d7 ]=tf2ss(g31num,g31den);
[a8 ,b8 ,c8 ,d8 ]=tf2ss(g32num,g32den);
[a9 ,b9 ,c9 ,d9 ]=tf2ss(g33num,g33den);
[a10,b10,c10,d10]=tf2ss(g41num,g41den);
[a11,b11,c11,d11]=tf2ss(g42num,g42den);
[a12,b12,c12,d12]=tf2ss(g43num,g43den);
[a13,b13,c13,d13]=tf2ss(g51num,g51den);
[a14,b14,c14,d14]=tf2ss(g52num,g52den);
[a15,b15,c15,d15]=tf2ss(g53num,g53den);
ginum=[1];giden=[1];
[a16,b16,c16,d16]=tf2ss(ginum,giden);
[a17,b17,c17,d17]=tf2ss(ginum,giden);
[a18,b18,c18,d18]=tf2ss(ginum,giden);
[a19,b19,c19,d19]=tf2ss(ginum,giden);
[a20,b20,c20,d20]=tf2ss(ginum,giden);
[a21,b21,c21,d21]=tf2ss(ginum,giden);
[a22,b22,c22,d22]=tf2ss(ginum,giden);
[a23,b23,c23,d23]=tf2ss(ginum,giden);
nblocks=23;
blkbuild;
inputs=[16 17 18 19 20];
outputs=[21 22 23];
q=[1 16 0 0 0 0
   2 16 0 0 0 0
   3 16 0 0 0 0
   4 17 0 0 0 0
   5 17 0 0 0 0
   6 17 0 0 0 0
   7 18 0 0 0 0
   8 18 0 0 0 0
   9 18 0 0 0 0
  10 19 0 0 0 0
  11 19 0 0 0 0
  12 19 0 0 0 0
  13 20 0 0 0 0
  14 20 0 0 0 0
  15 20 0 0 0 0
  21 1 4 7 10 13
  22 2 5 8 11 14
  23 3 6 9 12 15];
[ac,bc,cc,dc]=connect(a,b,c,d,q,inputs,outputs) ;