clear ring12000d;

transplot(spec12000d);

pause;

twotransplot(spec12000d/10, spec12000);

pause;

ring12000da = filtpz([pzpoint(12000, 1-1/sqrt(10)), pzpoint(12000, 1-1/sqrt(10))], [], stdimp);

spec12000da = fft(ring12000da);

clear ring12000da;

twotransplot(spec12000da, spec12000);

pause;

ring12000ta = filtpz([pzpoint(12000, 0.75), pzpoint(12000, 0.75), pzpoint(12000, 0.75)], [], stdimp);

spec12000ta = fft(ring12000ta);

clear ring12000ta;

twotransplot(spec12000ta/6.4, spec12000);

pause;

ring440c12000 = filtpz([pzpoint(440, 0.9), pzpoint(12000, 0.9)], [], stdimp);

spec440c12000 = fft(ring440c12000);

clear ring440c12000;

transplot(spec440c12000);

pause;

ring440p12000 = filtpz([pzpoint(440, 0.9), pzpoint(12000, 0.9)], ...
   [(pzpoint(440, 0.9) + pzpoint(12000, 0.9))/2], stdimp);

spec440p12000 = fft(ring440p12000);

clear ring440p12000;

transplot(spec440p12000);

pause;

impz440 = filtpz([], [pzpoint(440, 0.9)], stdimp);

specz440 = fft(impz440);

clear impz440;

transplot(specz440);

pause;

impz6220 = filtpz([], [pzpoint(6220, 0.9)], stdimp);

specz6220 = fft(impz6220);

clear impz6220;

transplot(specz6220);

pause;

impez6220 = filtpz([], [pzpoint(6220, 1.1)], stdimp);

specez6220 = fft(impez6220);

clear imipez6220;

transplot(specez6220);

pause;

impez6220a = filtpz([], [pzpoint(6220, 10)], stdimp);

specez6220a = fft(impez6220a);

clear impez6220a;

transplot(specez6220a);

pause;
