clear all; clc
% Form the transfer function matrix
g11 = 10*tf([1 1],[1 0.2 100]); 
g12 = tf(1,[1 1]); 
g21 = tf([1 2],[1 0.1 10]); 
g22 = 5*tf([1 1],conv([1 2],[1 3])); 
G = [g11 g12;g21 g22]

% Form the state space representation
Gss = ss(G)

%Form the minimal realization
Gssm = minreal(ss(G))

%Showing data of the stat-space
[A,B,C,D] = ssdata(Gss)
[Am,Bm,Cm,Dm] = ssdata(Gssm)

%Plot the min-max singular values
sigma(A,B,C,D)