Mercurial > hg > ltpda
view m-toolbox/classes/@parfrac/fromRational.m @ 4:e3c5468b1bfe database-connection-manager
Integrate with LTPDAPreferences
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: fromRational % % DESCRIPTION: Construct a parfrac from a rational TF % % CALL: pf = fromRational(a, pl) % % PARAMETER: pl - plist % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function pf = fromRational(pf, pli) VERSION = '$Id: fromRational.m,v 1.12 2011/08/15 12:31:59 hewitson Exp $'; % get pzmodel info ii = parfrac.getInfo('parfrac', 'From Rational'); % 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 rat = find(pl, 'rational'); % if the denominator is not a polynomial but a single number, then the % model is not really a rational model. It is practically a polynomial % model and cannot be converted in partial fractions. Output PF model % will be a zero residues and poles. Dterms will contain the numerator % polynomial if numel(rat.den)==1 pf.res = 0; pf.poles = 0; pf.dir = (rat.num)./(rat.den); pf.pmul = 0; else % convert a true rational model in partial fractions [res, poles, dterms, pmul] = utils.math.cpf('INOPT', 'RAT', ... 'NUM', rat.num, ... 'DEN', rat.den, 'MODE', 'SYM'); pf.res = res; pf.poles = poles; pf.dir = dterms; pf.pmul = pmul; end % Set properties from rational object if isempty(pl.find('ounits')) pl.pset('ounits', rat.ounits); end if isempty(pl.find('iunits')) pl.pset('iunits', rat.iunits); end if isempty(pl.find('name')) pl.pset('name', sprintf('parfrac(%s)', rat.name)); end if isempty(pl.find('description')) pl.pset('description', rat.description); end % Add history pf.addHistory(ii, pl, [], rat.hist); % Set object properties pf.setObjectProperties(pl); end