%《數(shù)字信號(hào)處理教程——MATLAB釋義與實(shí)現(xiàn)》第七章例7.6.1程序hc761
% 加窗后數(shù)字微分器的幅頻特性
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
N=6; tau=(N-1)/2; n=[0:N-1]+eps;          % 設(shè)定微分器長(zhǎng)度
hd =-sin((n-tau).*pi)./(pi.*(n-tau).^2);    % 計(jì)算其矩形窗截?cái)嗝}沖響應(yīng)
hh=hd.*hamming(N)';                     % 加哈明窗后的系數(shù)向量
[Hd,wd]=freqz(hd,1);                    % 計(jì)算矩形窗截?cái)辔⒎制黝l率響應(yīng)
[Hh,wh]=freqz(hh,1);                  % 計(jì)算哈明窗截?cái)辔⒎制黝l率響應(yīng)
plot(wd,abs(Hd),wh,abs(Hh),':',wd,abs(wd),'--'),hold on  %繪制幅特性
legend('矩形窗','哈明窗','理想'),grid on
set(gcf,'color','w'),pause              % 置圖形背景色為白
% 用信號(hào)處理工具箱中的最優(yōu)化設(shè)計(jì)函數(shù)設(shè)計(jì)希爾波特變換器
F=[0,1],A=[0,pi];
b=firls(N-1,F,A,'differentiator')   % 進(jìn)行最小二乘設(shè)計(jì)
b1=remez(N-1,F,A,'differentiator')  % 進(jìn)行最優(yōu)等波動(dòng)設(shè)計(jì)