fgp851.m

来自「matlab数字信号处理代码」· M 代码 · 共 19 行

%《数字信号处理教程——MATLAB释义与实现》第八章图8.5.1生成程序fgp851
% 数字低通滤波器变换为其它数字频带的映射
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
r=1;theta=linspace(0,2*pi,501);             % 将圆周角均分为501份
subplot(2,2,1)
set(gcf,'color','w')                        % 置图形背景色为白
polar(theta,r*ones(1,501),':'),hold on      % 画出整个单位圆
polar(theta([[452:501],[1:50]]),r*ones(1,100),'.'),hold on % 加深低通频率范围的圆周
subplot(2,2,2)
polar(theta,r*ones(1,501),':'),hold on      % 画出整个单位圆
polar(theta([180:321]),r*ones(1,[321-180+1]),'.') % 加深高通频率范围的圆周
subplot(2,2,3)
polar(theta,r*ones(1,501),':'),hold on      % 画出整个单位圆
polar(theta([[80:120],[381:421]]),r*ones(1,82),'.') % 加深带通频率范围的圆周
subplot(2,2,4)
polar(theta,r*ones(1,501),':'),hold on      % 画出整个单位圆
polar(theta([[1:80],[120:381],[422:501]]),r*ones(1,[80+262+80]),'.') % 加深带阻频率范围的圆周