fgp314.m

digital signal processing-MATLAB

%《数字信号处理教程——MATLAB释义与实现》第三章演示程序fgp314
% 数字频率与模拟频率的关系
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
T=1;
Omega=[-2.5*pi:0.01:2.5*pi]/T;    % 设定比乃奎斯特频率区间大2.5倍的模拟频率区间
Omegad=mod(Omega-pi,2*pi)-pi;           % 映射到主值区间，得到数字频率
plot(Omega,Omegad)                      % 画出两者的关系曲线
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')   % 标注
text(1,2,'数字频率\omega=\OmegaT')
plot([0.7,0.7],[0,0.7],':')             % 补上三条竖的短虚线
plot([0.7,0.7]+2*pi,[0,0.7],':')
plot([0.7,0.7]-2*pi,[0,0.7],':')