%《數(shù)字信號處理教程——MATLAB釋義與實(shí)現(xiàn)》第三章演示程序fgp314
% 數(shù)字頻率與模擬頻率的關(guān)系
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
T=1;
Omega=[-2.5*pi:0.01:2.5*pi]/T;    % 設(shè)定比乃奎斯特頻率區(qū)間大2.5倍的模擬頻率區(qū)間
Omegad=mod(Omega-pi,2*pi)-pi;           % 映射到主值區(qū)間，得到數(shù)字頻率
plot(Omega,Omegad)                      % 畫出兩者的關(guān)系曲線
line([-2.5*pi,2.5*pi],[0,0]),hold on    % 繪制輔助線
line([0,0],1.2*[-pi,pi]),               % 繪制輔助線
plot([-2.5*pi,2.5*pi],[0.7,0.7],'-.')
set(gcf,'color','w')                        % 置圖形背景色為白
xlabel('模擬頻率\Omega'),title('\Omega0')   % 標(biāo)注
text(1,2,'數(shù)字頻率\omega=\OmegaT')
plot([0.7,0.7],[0,0.7],':')             % 補(bǔ)上三條豎的短虛線
plot([0.7,0.7]+2*pi,[0,0.7],':')
plot([0.7,0.7]-2*pi,[0,0.7],':')