Mercurial > hg > ltpda

diff m-toolbox/test/test_pzmodel_fromParfrac.m @ 0:f0afece42f48

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_pzmodel_fromParfrac.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,27 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% TEST pzmodel/fromParfrac
+% 
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% HISTORY:     23-04-2009 L Ferraioli
+%                 Creation
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%% VERSION
+
+'$Id: test_pzmodel_fromParfrac.m,v 1.1 2009/10/21 12:50:43 luigi Exp $';
+
+%% test
+
+% build parfrac model
+pfm = parfrac([1 0.3+0.7i 0.3-0.7i],[-2 -3.1+0.35i -3.1-0.35i],0);
+
+% convert to pzmodel
+pzm = pzmodel(pfm);
+
+% get response of models
+f = [logspace(-5,1,300)]';
+rpfm = resp(pfm,plist('f',f));
+rpzm = resp(pzm,plist('f',f));
+
+% plot
+iplot(rpfm,rpzm)
+iplot(abs(rpfm)./abs(rpzm))
\ No newline at end of file