fig723.m

来自「高西全2009年最新版<数字信号处理>源码程序」· M 代码 · 共 32 行

%《数字信号处理--原理、实现及应用》第六章 图7.2.2绘图程序:fig722.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月


close all;clear all
w=[-pi,-pi/2,-pi/2,pi/2,pi/2,pi];
m=[0,0,1,1,0,0];
subplot(3,1,1);plot(w/pi,m);axis([-1,1,-0.3,2])
wk=[-512:512]*2*pi/1024;
Wrg=sin(wk*31/2)./sin(wk/2);
subplot(3,1,2);plot(wk/pi,1.2*Wrg/32);%axis([-1,1,-0.3,2])
wk=[-512:512]*2*pi/1024;
Wrg=sin((wk-pi/2)*31/2)./sin((wk-pi/2)/2);
subplot(3,1,3);plot(wk/pi,1.2*Wrg/32);%axis([-1,1,-0.3,2])
hold on;plot(w/pi,m,':');axis([-1,1,-0.3,2])
figure(2)
Wrg=sin((wk-pi/2+2*pi/31)*31/2)./sin((wk-pi/2+2*pi/31)/2);
subplot(3,1,1);plot(wk/pi,1.2*Wrg/32);%axis([-1,1,-0.3,2])
hold on;plot(w/pi,m,':');axis([-1,1,-0.3,2])
Wrg=sin((wk-pi/2-2*pi/31)*31/2)./sin((wk-pi/2-2*pi/31)/2);
subplot(3,1,2);plot(wk/pi,1.2*Wrg/32);%axis([-1,1,-0.3,2])
hold on;plot(w/pi,m,':');axis([-1,1,-0.3,2])
hn=fir1(30,1/2,boxcar(31));
wk1=[0:1023]*2*pi/1024;
Hk=fft(hn,1024);Hk1=Hk.*exp(j*wk1*15);
Hks=fftshift(Hk1);
wk=[-511:512]*2*pi/1024;
subplot(3,1,3);plot(wk/pi,Hks);%axis([-1,1,-0.3,2])
figure(3)
subplot(2,1,1);plot(wk/pi,Hks)
hold on;plot(w/pi,m,':');axis([-1,1,-0.3,2])