射频直采多相抽取信道化滤波器设计仿真

% author：张宗卫

% description：多相抽取滤波器设计仿真

% date：2020.2.2

%使用带通采样定律对射频信号直接采样，fs=6.2M,载波9.8M，经过采样后频谱搬移至2.6m，

%针对目前比较流行的204B接口，数据随路时钟clk=fs/4,一个时钟周期传输四个采样点，特别适合%多相抽取滤波器设计，该仿真使用此滤波器结构设计了带通滤波器实现载波9.8M和9.79M信号的

%分离，适合用于信道化滤波器设计。

clc

clear all

load('Bpf2600Coe.mat')%导入滤波器参数，该滤波器为fc1经过带通采样后

fs=620*10e3;%采样频率

f1=980*10e3;

f2=979*10e3;

% step1 产生脉冲

L=600000;

t=1/fs:1/fs:L/fs;

am=zeros(1,L);

TPulse=125;

t1=(4/1000)*fs;

t2=400000;

t3=424800;

for i=1:(4/1000)*fs

am(i)=sin(2*pi*TPulse*t(i));

end

for i=1:t1

ts(i)=sin(2*pi*TPulse*t(i));

end

j=1;

for i=t2+1:t3

am(i)=ts(j);

j=j+1;

end

figure(1)

plot(t,am,'r');

title('脉冲调制信号');

%step2 将脉冲信号加调制

fc1=sin(2*pi*f1*t)*2^14;% 调制1

fc2=sin(2*pi*f2*t)*2^14;% 调制2

fs1=am.*fc1;

fs2=am.*fc2;

%第一个脉冲为fc1频率的调制，第二个脉冲为fc2频率的调制

for i=1:L

if(i

Rec(i)=fs2(i);

else

Rec(i)=fs1(i);

end

end

figure(2)

plot(t,Rec)

title('射频脉冲信号');

%频谱分析

NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(Rec,NFFT)/L;

f = fs/2*linspace(0,1,NFFT/2+1);

figure(3)

% Plot single-sided amplitude spectrum.

plot(f,2*abs(Y(1:NFFT/2+1)))

title('Single-Sided Amplitude Spectrum of y(t)')

xlabel('Frequency (Hz)')

ylabel('|Y(f)|')

%step3 该信号通过设计好的带通滤波器

fRec=filter(Num,1,Rec);

figure(4)

plot(t,fRec)

title('信号过带通滤波器')

%step4 多相滤波器设计仿真：思想，将抽取放在滤波器前，降低运算速率%滤波器分解

b1=Num(1:4:1503);

b2=Num(2:4:1503);

b3=Num(3:4:1503);

b4=Num(4:4:1503);

%数据分解

Rec1=Rec(1:4:L);

Rec2=Rec(2:4:L);

Rec3=Rec(3:4:L);

Rec4=Rec(4:4:L);

%四路数据分别过四个滤波器

frec1=filter(b1,1,Rec1);

frec2=filter(b2,1,Rec2);

frec3=filter(b3,1,Rec3);

frec4=filter(b4,1,Rec4);

%合成滤波数据

ffrec=frec1+frec2+frec3+frec4;

tt=t(1:4:L);

figure(5)

plot(tt,ffrec)

title('信号过多相抽取带通滤波器')