hc741a.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第七章例7.4.1a程序hc741a
% 理想频率响应的频率样本设定
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=input('N= (N必须为奇数）');wc=input('wc=');   % 给定原始数据
N=N+mod(N+1,2);                                 % 使N成为奇数 
N1=fix(wc/(2*pi/N));N2=N-2*N1-1;                % 划分通带和阻带的样本点数
A=[ones(1,N1+1),zeros(1,N2),ones(1,N1)];	% 符幅特性样本序列
theta = - pi*[0:N-1]*(N-1)/N;			% 相位特性样本序列
wdl=[0,wc,wc,2*pi-wc,2*pi-wc,2*pi]/pi;Adl=[1,1,0,0,1,1];   % 绘制理想幅特性的频率和幅度数据
l = 0:N-1; wl = (2*pi/N)*l;             % 把样本下标化为实际频率值
subplot(2,1,1);plot(wl(1:N)/pi,A(1:N),'.',wdl,Adl); % 画图
axis([0,2,-0.1,1.1]),grid on
set(gcf,'color','w')
subplot(2,1,2);plot(wl(1:N)/pi,theta(1:N),'.'); 
set(gcf,'color','w');                % 置图形背景色为白
grid on