Subversion Repositories sciir

close all
clear
clc
 
f = [0 0.25 0.35 1.0];      % Corner frequencies
m = [1 1 0 0];              % Filter magnitudes        
n = 3;                      % Filter order
 
[b,a] = yulewalk(n, f, m);  % Design IIR Filter
 
[h,w] = freqz(b,a,128);     % FFT
 
figure
plot(f,m,w/pi,abs(h))       % Plot frequency response
 
figure
step(filt(b,a))             % Plot step response
 
figure
impulse(filt(b,a))          % Plot impuse response
 
figure
zplane(b,a)                 % Plot poles/zeros
 
[sos, g] = tf2sos(b, a);    % Second Order Sections
 
afx = num2fixpt(a, sfix(8), 2^-6);
bfx = num2fixpt(b, sfix(8), 2^-6);
fvtool(bfx, afx)
fvtool(b, a)