Mercurial > hg > ltpda
comparison m-toolbox/test/test_convert_pzmodel_parfrac.m @ 0:f0afece42f48
Import.
| author | Daniele Nicolodi <nicolodi@science.unitn.it> |
|---|---|
| date | Wed, 23 Nov 2011 19:22:13 +0100 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:f0afece42f48 |
|---|---|
| 1 mc | |
| 2 | |
| 3 %% start from pzmodel | |
| 4 pzm = pzmodel(1, {[1 2], 10}, [3 2]) | |
| 5 pzm.setName; | |
| 6 | |
| 7 %% Convert to partial fractions | |
| 8 pf = parfrac(pzm) | |
| 9 | |
| 10 %% Convert back again | |
| 11 | |
| 12 %% Check resp | |
| 13 | |
| 14 pzmr = resp(pzm); | |
| 15 pfr = resp(pf); | |
| 16 | |
| 17 iplot(pzmr, pfr) | |
| 18 | |
| 19 %% Make a filter | |
| 20 | |
| 21 pfilt = miir(pzm) | |
| 22 filts = miir(pf) | |
| 23 | |
| 24 fr = resp(filts) | |
| 25 | |
| 26 |