ex3040.m

the file contain many matlab signal source code and many example .

ws1=0.2*pi;wp1=0.35*pi;
wp2=0.65*pi;ws2=0.8*pi;
As=60;
tr_width=min([wp1-ws1,ws2-wp2])
M=ceil(10*pi/tr_width)+1
n=[0:1:M-1];
wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;
hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);
w_bla=(blackman(M))';
h=hd.*w_bla;
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
Rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))
As=-round(max(db(ws2/delta_w+1:1:501)))
% plots
figure(1);
subplot(221)
stem(n,hd);axis([0 M-1 -0.4 0.5]);
subplot(222)
stem(n,w_bla);axis([0 M-1 0 1.1])
subplot(223)
stem(n,h);axis([0 M-1 -0.4 0.5]);
subplot(224)
plot(w/pi,db);axis([0 1 -150 10])
grid on