%《數字信號處理教程——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')  % 進行最優等波動設計