微多普勒处理雷达数据

微多普勒处理雷达数据 已通过测试

% 变步长频域对消
function [e W mu] = blmsfft_vss（xdLNmu0mu_maxmu_minlampdabetaPianYiW）
% e 对消输出
% x input reference signal 
% d 回波
% L filter length
% N FFT length
% mu0 初始步长
% mu_max最大步长
% mu_min最小步长
%lambda= 0.99; 系数遗忘因子
%beta= 0.01; 步长调整系数
% PianYi  直达波提前 点数
% W filter coeff
ntr = length（x）;            %  temporary number of iterations 
ntrB = floor（（ntr-N）/（N-L））;  %  temporary number of block iterations
e = zeros（size（x））;                      %  initialize error signal vector
ee = zeros（1N）;
power=sum（abs（ x（1+L:N）.^2 ））;   % 能量归一化用，power可设为常数
mu（1）=mu0;
%

%  Main loop 
for k=0:ntrB 
    xx = x（1+k*（N-L）+PianYi:N+k*（N-L）+PianYi）;
    X=fft（xx）;       % 直达波信号组成N点向量，每次步进N-L点N-L即为块长度
    y=ifft（X.*W）;                           % 滤波预测值，取y的后N-L点
    ee（1+L:N）=d（1+L+k*（N-L）:N+k*（N-L））-y（1+L:N）;  % 回波- 预测值=对消中间向量ee，回波也每次步进N-L点
    e（1+L+k*（N-L）:N+k*（N-L））=ee（1+L:N）;    % ee的后N-L（块长度）点作为对消输出
    xe=ifft（conj（X）.*fft（ee））;               % xe的作用是更新频域滤波器系数
    xe（1+L:N）=0;                             % xe前L点有用，后N-L点置为0
    W=W+fft（xe）*mu（k+1）/power;               % 归一化LMS即NLMS更新公式，power做为能量归一化用，power可设为常数以减小运算量
    mu（k+2）=lampda*mu（k+1）+beta*sum（abs（ ee（1+L:N）））/sum（abs（ d（1+L+k*（N-L）:N+k*（N-L））））;
    if mu（k+2）>=mu_max
       mu（k+2）=mu_max;
    elseif mu（k+2）<=mu_min
       mu（k+2） = mu_min;
    end;
end
% coeff = W（1:L）.‘;       %  save final coefficient vector

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-03-21 09:23  micro\
     文件        1346  2011-01-13 19:14  micro\CAF_CIC_FIR.m
     文件        4533  2013-05-16 08:49  micro\CAF_mDgen.m
     文件        5140  2012-07-24 17:03  micro\DBF_cancellation_CAF.m
     文件        1605  2011-04-08 21:44  micro\blmsfft_vss.m
     文件        1230  2013-05-15 21:22  micro\mD_CAF_CIC_FIR.m
     文件        4373  2013-05-07 09:52  micro\tfrspwv.m
     文件        3069  2013-05-07 09:52  micro\tfrstft.m
     文件        5356  2005-10-13 23:15  micro\tftb_window.m