基于多项滤波的数字正交变换MATLAB仿真程序

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);