DFT多相滤波器组信道化结构

描述了16通道DFT多相滤波器组信道化结构，输出为16个子信道的信号，程序包括滤波器设计

clear;
clc;
close all;
fid=fopen（‘filter.txt‘‘r‘）;
[fir_coffr_data]=fscanf（fid‘%f‘[1256]）;


f0=30e6;   %中心频率
T=1e-6;    %线性调频信号脉宽
B=16e6;     %线性调频信号带宽
K=B/T;      %调频斜率
Fs=1024e6;   %采样率
Ts=1/Fs; 
D=16;
N=1024;
t=linspace（0TN）;
st=cos（2*pi*f0*t+pi*K*t.^2）; %线性调频信号（实信号）
%  st=exp（j*（2*pi*f0*t+pi*K*t.^2））;%线性调频信号（复信号）

figure（1）
freq=linspace（0FsN）;
signal_am_freq=abs（fft（st））;
plot（tsignal_am_freq）;
xlabel（‘频率 Hz‘）;
ylabel（‘幅度‘）;

 %多相滤波器中抽取倍数为16倍
S=reshape（st1664）   %将总的信号分解成一个矩阵
S3 = flipud（S）; 
for m=1:64
  S（:m）= S3（:m）.*exp（-1*j*（pi/2）*（m））;
% S（:m）= S3（:m）.*（-1）^m;

end

 
% %==========================================================================
% %多相滤波器通道
% %==========================================================================
 n=256;      %滤波器系数个数
for p = 1:D
    for m = 1:n/D
        h（pm） = fir_coff（（m-1）*D+p）; 
    end
end
for K=1:D
    S4（K:）=filter（h（K:）1S（K:））;
end
for K=1:D
    S1（K:）= S4（K:）.*exp（1*j*（pi/（2*D）*（K-1）））;
end



% 以下是对每一列进行FFT转换，可得

[gh] = size（S1）;
for i=1:h
        output（:i）=D*ifft（S1（:i）16）;
end
 
 [ymyn] = size（output）;
 fsy = Fs/D;  %抽取后 采样点数的变化
 fy = linspace（-fsy/2fsy/2yn）;
 
 figure（2）
 %Y1=real（output（16:））
 ys（1:） = abs（fftshift（fft（output（1:））））;
 
    subplot（121）;
    plot（fyys（1:））;
%    ylim（[0 6e5]）;
    xlabel（‘频率 Hz‘）;
    ylabel（‘幅度‘）;
  axis（[-6.4e7 6.4e7 0 12]）;
    grid on;
    
   
    ys（2:） = abs（fftshift（fft（output（2:））））;
    
    subplot（122）;
    plot（fyys（2:））;
%    ylim（[0 6e5]）;
    xlabel（‘频率 Hz‘）;
    ylabel（‘幅度‘）;
 axis（[-6.4e7 6.4e7 0 12]）;
    grid on;
    
    figure（3）
for c=1:16
    ys（c:） = abs（fftshift（fft（output（c:））））;
    subplot（116c）;
    plot（fyys（c:））;
%    ylim（[0 6e5]）;
    xlabel（‘频率 Hz‘）;
    ylabel（‘幅度‘）;
    axis（[-6.4e7 6.4e7 0 20]）;
    grid on;
end

% 
% 
% M=output（1:）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5698  2015-04-07 18:38  filter.txt

     文件       2083  2018-06-07 17:57  channel.m

----------- ---------  ---------- -----  ----

                 7781                    2