--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/@ssm/ssmFromPzmodel.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,67 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% FUNCTION:    ssmFromPzmodel
+%
+% DESCRIPTION: Construct a statespace model from a pzmodel
+%
+% CALL:        see ssm, this function is private
+%
+% TODO :  Modify using ss2zp/zp2ss
+%
+% VERSION :    '$Id: ssmFromPzmodel.m,v 1.26 2011/04/08 08:56:23 hewitson Exp $';
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function ssmout = ssmFromPzmodel(varargin)
+  
+  VERSION = '$Id: ssmFromPzmodel.m,v 1.26 2011/04/08 08:56:23 hewitson Exp $';
+  utils.helper.msg(utils.const.msg.MNAME, ['running ', mfilename]);
+  % get info
+  ii = ssm.getInfo('ssm', 'from pzmodel');
+  % Set the method version string in the minfo object
+  ii.setMversion([VERSION '-->' ii.mversion]);
+  
+ if nargin ~=2
+    error('ssmFromPzmodel need 2 inputs : (obj, plist) ')
+  elseif isa(varargin{1}, 'pzmodel') && isa(varargin{2}, 'plist')
+    pl = combine(varargin{2}, ii.plists);
+    pzm = varargin{1};
+  else
+    error('### Please input  (<object>,<plist>)');
+  end
+  
+  ssmout = ssm.initObjectWithSize(size(pzm,1),size(pzm,2));
+  
+  for i_pzms=1:numel(pzm)
+    [A,B,C,D] = utils.math.pzmodel2SSMats(pzm(i_pzms));
+
+    ssmout(i_pzms).dmats = {D};
+    ssmout(i_pzms).amats = {A};
+    ssmout(i_pzms).bmats = {B};
+    ssmout(i_pzms).cmats = {C};
+    
+    ssmout(i_pzms).name  = pzm(i_pzms).name;
+    ssmout(i_pzms).timestep = 0;
+    
+    inputstr = 'input';
+    outputstr = 'output';
+    ssstr =  'state';
+    Nss = size(A,1);
+    
+    ssmout(i_pzms).inputs = ssmblock.makeBlocksWithData({inputstr}, [], {{inputstr}}, {pzm(i_pzms).iunits}, []);
+    ssmout(i_pzms).outputs = ssmblock.makeBlocksWithData({outputstr}, [], {{outputstr}}, {pzm(i_pzms).ounits}, []);
+    ssmout(i_pzms).states = ssmblock.makeBlocksWithSize(Nss, ssstr);
+    
+    ssmout(i_pzms).addHistory(ii, pl, {''}, pzm(i_pzms).hist);
+    
+    if ~strcmp(pl.find('name'),'None')
+      ssmout(i_pzms).name = pl.find('name');
+    end
+    if ~strcmp(pl.find('description'),'')
+      ssmout(i_pzms).description = pl.find('description');
+    end
+
+  end
+  
+  
+end