%半帶濾波器的設計
%半帶濾波器的系數滿足如下條件
%h(n)=0, when n=+2/-2,+4/-4,...,+2k/-2k, and k!＝0;
%設計內插因子為2的半帶濾波器
%截取符號數為10

clear;
clc;    
M=8;
w=blackman(17);                     %Blackman窗的數據
for n=0:2*M
    i=n+1;              
    H(i)=sinc((n-M)/2)*w(i);        %sinc=sin(pi*x)/(pi*x);
    if abs(H(i))<1e-10
        H(i)=0;
    end;
end;

% 量化比特數為8,最高位為符號位;
L=8;
Hb_8bit=round(H*2^(L-1));
save hb_8bit.mat Hb_8bit;
%給出頻率響應
Hb_8bit=Hb_8bit/sum(Hb_8bit);
freqz(Hb_8bit,1,'whole');

% 量化比特數為10,最高位為符號位;
L=10;
Hb_10bit=round(H*2^(L-1));
save hb_10bit.mat Hb_10bit;
%給出頻率響應
Hb_10bit=Hb_10bit/sum(Hb_10bit);
freqz(Hb_10bit,1,'whole');

% 量化比特數為12,最高位為符號位;
L=12;
Hb_12bit=round(H*2^(L-1));
save hb_12bit.mat Hb_12bit;
%給出頻率響應
Hb_12bit=Hb_12bit/sum(Hb_12bit);
%freqz(Hb_12bit,1,'whole');
%  0 1 0 -19 0 89 0 -300 0 1252 2048 1252 0 -300 0 89 0 -19 0 1 0