function filter_emu( Num )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
%num_size = sizeof(Num);

%利用多項濾波器的分支特性獲得兩個有半個取樣點時延的低通濾波器係數</span>

filt1 = zeros(1,8);
filt2 = zeros(1,8);
for m = 1:64
    if(mod(m,8) == 4)
        filt1(i) = Num(m);
    end
    if(mod(m,8) == 0)
        filt2(i) = Num(m);
        i = i + 1;
    end
end
%------生成窄帶訊號，中頻150MHz，頻寬不大於20MHz
%------模擬訊號x(t)=a(t)*cos[2*pi*f0*t+phi(t)]
f0 = 1.5e8;     %中心頻率
fs = 2e8;       %取樣頻率
N = 1600;       %取的樣本點數
n = 0:N-1;      %取的樣本序列
t = n/fs;       %獲得以1/fs為時間間隔的取樣序列
%a = 1+cos(2*pi*1000*t);      
                %獲取a(t)的取樣點
phi = 2*pi*2e6;
                %獲取phi
xt = cos(2*pi*f0*t+phi*t);
                %生成窄帶訊號並獲取其取樣點
%------進行2倍抽取並混頻
nt = 1:N/20;
xi = xt(2*nt).*((-1).^nt);  %獲取同相分量
xq = xt(2*nt-1).*((-1).^nt);%獲取正交分量

%------列印同相分量與正交分量
figure(1);
plot(nt,xi,'r');
hold on;
plot(nt,xq,'g');
grid on;
%------濾波
x1 = filter(filt1,1,xi);
x2 = filter(filt2,1,xq);
%------列印恢復後的同相正交分量
figure(2);
plot(nt,x1,'r');
hold on;
plot(nt,x2,'g');
grid on;

end