+
−%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
−% 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))