0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044 function varargout = reduce_model(varargin)
0045
0046
0047 utils.helper.msg(utils.const.msg.MNAME, ['running ', mfilename]);
0048
0049
0050 if utils.helper.isinfocall(varargin{:})
0051 varargout{1} = getInfo(varargin{3});
0052 return
0053 end
0054
0055
0056 if ~nargin==2, error('wrong number of inputs!'), end
0057
0058 if ~isequal(class(varargin{1}),'ssm'), error(['argument is not a ssm but a ', class(varargin{1})]),end
0059 if ~isequal(class(varargin{2}),'plist'), error(['argument is not a plist but a ', class(varargin{2})]),end
0060
0061 sys = varargin{1};
0062 options = varargin{2};
0063 Nsys = length(sys);
0064 Noptions = length(options);
0065
0066
0067 inputs = find(options, 'inputs');
0068 outputs = find(options, 'outputs');
0069 states = find(options, 'states');
0070
0071
0072 if Noptions==1
0073 sys_out(:,1) = reshape(sys,Nsys,1);
0074 else
0075 for j=1:Noptions
0076 sys_out(:,j)=copy(reshape(sys,Nsys,1),true);
0077 end
0078 end
0079 for i_sys=1:Nsys
0080 for i_options=1:Noptions
0081
0082
0083
0084
0085 [As] = ssm.cell_select(sys.amats, sys.sssizes, sys.sssizes, states, states, sys.ssvarnames, sys.ssvarnames);
0086 [Ms] = ssm.cell_select(sys.mmats, sys.sssizes, sys.sssizes, states, states, sys.ssvarnames, sys.ssvarnames);
0087 [Bs, sssizes] = ssm.cell_select(sys.bmats, sys.sssizes, sys.inputsizes, states, inputs, sys.ssvarnames, sys.inputvarnames);
0088 [Cs] = ssm.cell_select(sys.cmats, sys.outputsizes, sys.sssizes, outputs, states, sys.outputvarnames, sys.ssvarnames);
0089 [Ds, outputsizes, inputsizes] = ssm.cell_select(sys.dmats, sys.outputsizes, sys.inputsizes, outputs, inputs, sys.outputvarnames, sys.inputvarnames);
0090
0091 [ssini] = ssm.cell_select(sys.ssini, sys.sssizes, 1, states, {true}, sys.ssvarnames, {' '});
0092 [inputvarnames2] = ssm.cell_select(sys.inputvarnames, 1, sys.inputsizes, {true}, inputs, {' '}, sys.inputvarnames);
0093 [ssvarnames2] = ssm.cell_select(sys.ssvarnames, 1, sys.sssizes, {true}, states, {' '}, sys.ssvarnames);
0094 [outputvarnames2] = ssm.cell_select(sys.outputvarnames, 1, sys.outputsizes, {true}, outputs, {' '}, sys.outputvarnames);
0095
0096 sys_out(i_sys, i_options).amats = As;
0097 sys_out(i_sys, i_options).mmats = Ms;
0098 sys_out(i_sys, i_options).bmats = Bs;
0099 sys_out(i_sys, i_options).cmats = Cs;
0100 sys_out(i_sys, i_options).dmats = Ds;
0101 sys_out(i_sys, i_options).inputvarnames = inputvarnames2;
0102 sys_out(i_sys, i_options).ssvarnames = ssvarnames2;
0103 sys_out(i_sys, i_options).outputvarnames = outputvarnames2;
0104 sys_out(i_sys, i_options).inputsizes = inputsizes;
0105 sys_out(i_sys, i_options).sssizes = sssizes;
0106 sys_out(i_sys, i_options).outputsizes = outputsizes;
0107 sys_out(i_sys, i_options).ssini = ssini;
0108
0109 sys_out(i_sys, i_options).name = [sys_out(i_sys, i_options).name, ' reduced'];
0110 sys_out(i_sys, i_options).description = [sys_out(i_sys, i_options).description, ' reduced'];
0111 sys_out(i_sys, i_options).addHistory(ssm.getInfo(mfilename), options(i_options) , {''}, sys_out(i_sys, i_options).hist );
0112 end
0113 end
0114 sys_out = validate(sys_out);
0115 varargout = {sys_out};
0116 end
0117
0118 function ii = getInfo(varargin)
0119 if nargin == 1 && strcmpi(varargin{1}, 'None')
0120 sets = {};
0121 pls = [];
0122 elseif nargin == 1 && ~isempty(varargin{1}) && ischar(varargin{1})
0123 sets{1} = varargin{1};
0124 pls = getDefaultPlist(sets{1});
0125 else
0126 sets = {'Default'};
0127 pls = [];
0128 for kk=1:numel(sets)
0129 pls = [pls getDefaultPlist(sets{kk})];
0130 end
0131 end
0132
0133 ii = minfo(mfilename, 'ssm', '', 'STATESPACE', '$Id: resp.m,v 1.17 2008/07/22 10:22:38 ingo Exp $', sets, pls);
0134 end
0135
0136 function plo = getDefaultPlist(set)
0137 switch set
0138 case 'Default'
0139 plo = plist('inputs','ALL','outputs','ALL', 'states', 'ALL');
0140 otherwise
0141 error('### Unknown parameter set [%s].', set);
0142 end
0143 end