hc743.m

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

%《数字信号处理教程——MATLAB释义与实现》第七章例7.4.3程序hc743
% 用fir2函数生成限带波形
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=input('N=');M=N-1;  % 给定序列长度
f=[0,0.45,0.5,1];   % 给定预期频谱特性
m=[0,1,0,0];
h=fir2(M,f,m,boxcar(N));      % 生成脉冲响应序列
[H,w]=freqz(h,1);    % 计算脉冲序列的频谱
subplot(3,1,1),plot(f,m)    % 画出设定频谱形状
title('设定频谱'),
xlabel(''),ylabel('m')
subplot(3,1,2),stem([0:N-1],h,'.')    % 画出生成的脉冲序列形状
title('用fir2函数生成的脉冲序列'),
xlabel(''),ylabel('h(n)')
subplot(3,1,3),plot(w/pi,abs(H))    % 画出实际频谱形状
title('实际频谱'),
xlabel('\omega/\pi'),ylabel('abs(H)')
set(gcf,'color','w')                % 置图形背景色为白