Mercurial > hg > ltpda
view m-toolbox/classes/@rational/fromPzmodel.m @ 13:e05504b18072 database-connection-manager
Move more functions to utils.repository
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % FUNCTION: fromPzmodel % % DESCRIPTION: Construct a rational from a pzmodel % % CALL: r = fromPzmodel(a, pl) % % PARAMETER: pl - plist % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function r = fromPzmodel(r, pli) VERSION = '$Id: fromPzmodel.m,v 1.11 2011/08/15 12:41:45 hewitson Exp $'; % get pzmodel info ii = rational.getInfo('rational', 'From Pzmodel'); % Set the method version string in the minfo object ii.setMversion([VERSION '-->' ii.mversion]); % Combine input plist with default values pl = applyDefaults(ii.plists, pli); % Set fields pzm = find(pl, 'pzmodel'); %--- Convert to rational % get normalising gains pg = 1; for jj=1:numel(pzm.poles) if pzm.poles(jj).q > 0.5 pg = pg .* 4*pi*pi*pzm.poles(jj).f*pzm.poles(jj).f; else pg = pg .* 2*pi*pzm.poles(jj).f; end end zg = 1; for jj=1:numel(pzm.zeros) if pzm.zeros(jj).q > 0.5 zg = zg .* 4*pi*pi*pzm.zeros(jj).f*pzm.zeros(jj).f; else zg = zg .* 2*pi*pzm.zeros(jj).f; end end % construct rational terms if isempty(pzm.poles) % in case on no poles den = 1; else den = poly(vertcat(pzm.poles(:).ri))./pg; end if isempty(pzm.zeros) % in case on no zeros num = pzm.gain; else num = pzm.gain.*poly(vertcat(pzm.zeros(:).ri))./zg; end r.num = num; r.den = den; % Override some properties from the input pzmodel if isempty(pl.find('ounits')) pl.pset('ounits', pzm.ounits); end if isempty(pl.find('iunits')) pl.pset('iunits', pzm.iunits); end if isempty(pl.find('name')) pl.pset('name', sprintf('rational(%s)', pzm.name)); end if isempty(pl.find('description')) pl.pset('description', pzm.description); end % Add history r.addHistory(ii, pl, [], pzm.hist); % Set object properties r.setObjectProperties(pl); end