mc

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

%% Convert to partial fractions
pf = parfrac(pzm)

%% Convert back again

%% Check resp

pzmr = resp(pzm);
pfr  = resp(pf);

iplot(pzmr, pfr)

%% Make a filter

pfilt = miir(pzm)
filts = miir(pf)

fr = resp(filts)


