OSR = 32;
H = synthesizeNTF(5, OSR, 10)
N = 8192;
f = 85;
t = 0:N-1;
u = 0.5*sin(2*pi*f/N*t);
[a g b c] = realizeNTF(H, 'CRFB')
ABCD = stuffABCD(a, g, b, c, 'CRFB')
[v xn xmax y] = simulateDSM(u, ABCD, 10);
poles = H.p{:}
zeros = H.z{:}
