ex3039.m

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

wp=0.2*pi;
ws=0.3*pi;
tr_width=ws-wp
M=ceil(6.6*pi/tr_width)+1
n=[0:M-1];
wc=(ws+wp)/2
hd=ideal_lp(wc,M);
w_ham=(hamming(M))';
h=hd.*w_ham;
[db,mag,pha,grd,w]=freqz_m(h,[1])
delta_w=2*pi/1000;
Rp=-(min(db(1:wp/delta_w+1)))
As=-round(max(db(ws/delta_w+1:501)))
% plots
subplot(221)
stem(n,hd)
axis([0 M-1 -0.1 0.3])
subplot(222)
stem(n,w_ham)
axis([0 M-1 0 1.1])
subplot(223)
stem(n,h)
axis([0 M-1 -0.1 0.3])
subplot(224)
plot(w/pi,db);grid on
axis([0 1 -100 10])