ex3043.m

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

M=60;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l;
Hrs=[ones(1,7),0.5925,0.1099,zeros(1,43),0.1099,0.5925,ones(1,6)];
Hdr=[1,1,0,0];wdl=[0,0.2,0.3,1];
k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;
angH=[-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2)];
H=Hrs.*exp(j*angH);
h=real(ifft(H,M));
[db,mag,pha,grd,w]=freqz_m(h,1);
[Hr,ww,a,L]=Hr_Type2(h);
% plots
subplot(221)
plot(wl(1:31)/pi,Hrs(1:31),'o',wdl,Hdr);
axis([0,1,-0.1,1.1]);grid
subplot(222)
stem(l,h);axis([-1,M,-0.1,0.3])
subplot(223)
plot(ww/pi,Hr,wl(1:31)/pi,Hrs(1:31),'o');grid
axis([0 1 -0.1 1.1])
subplot(224)
plot(w/pi,db);axis([0,1,-100,10]);grid