hc761.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.6.1程序hc761
% 加窗后数字微分器的幅频特性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=6; tau=(N-1)/2; n=[0:N-1]+eps;          % 设定微分器长度
hd =-sin((n-tau).*pi)./(pi.*(n-tau).^2);    % 计算其矩形窗截断脉冲响应
hh=hd.*hamming(N)';                     % 加哈明窗后的系数向量
[Hd,wd]=freqz(hd,1);                    % 计算矩形窗截断微分器频率响应
[Hh,wh]=freqz(hh,1);                  % 计算哈明窗截断微分器频率响应
plot(wd,abs(Hd),wh,abs(Hh),':',wd,abs(wd),'--'),hold on  %绘制幅特性
legend('矩形窗','哈明窗','理想'),grid on
set(gcf,'color','w'),pause              % 置图形背景色为白
% 用信号处理工具箱中的最优化设计函数设计希尔波特变换器
F=[0,1],A=[0,pi];
b=firls(N-1,F,A,'differentiator')   % 进行最小二乘设计
b1=remez(N-1,F,A,'differentiator')  % 进行最优等波动设计