mc

%% start from pzmodel
pzm = pzmodel(5, {[1 2]}, {[10 2]});
pzm.setName;


%% Convert to rational

rat = rational(pzm);

%% Convert back

pzm2 = pzmodel(rat)


%% Convert to parfrac

pf = parfrac(rat);
pf2 = parfrac(pzm);

%% resp

pzmr = resp(pzm);
ratr = resp(rat);
pfr  = resp(pf, plist('f', ratr));
pfr2  = resp(pf2, plist('f', ratr));

iplot(pzmr, ratr)
iplot(pzmr, pfr)
iplot(pzmr, pfr2)


