view m-toolbox/classes/@ssm/ssmFromPzmodel.m @ 0:f0afece42f48
Import.
author |
Daniele Nicolodi <nicolodi@science.unitn.it> |
date |
Wed, 23 Nov 2011 19:22:13 +0100 (2011-11-23) |
parents |
|
children |
|
line source
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 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