NBI消除窄带干扰消除

clc
close all
clear all
%% -初始化参数
fs=2;                                                   %采样频率
ns=1;
a=2.3;                      
tr=0.4*ns;                     
tf=2*ns; 
m=160;                                                    %序列延迟长度
nline=4096;                                             %序列点长（需要为偶数）
noise_n=14;                                             %单频噪声数量
fir_n=100;                                              %fir滤波器点数
dt=1/fs;
t=-nline/2*dt-m*dt:dt:nline/2*dt-dt;                    %离散时间点
t_temp=-nline/2*dt-dt-m*dt:dt:nline/2*dt-dt;            %离散时间点
nt=length（t）;
%% -sn vn yn信号序列产生
%-------------sn---------------
s=1./（exp（-a*t_temp/tr）+exp（a*t_temp/tf））;
sn=5*（s（2:end）-s（1:end-1））;                     %有用信号
%-------------yn---------------
freq=[0.611.82.134.85.25.76.16.46.777.89.3]‘*fs/20;
ampli=repmat（[0.5110.50.10.30.51110.50.31.50.5]‘1nt）;
% rand（‘seed‘1954）;
phase=repmat（2*pi*rand（noise_n1）1nt）;
yn=sum（ampli.*cos（2*pi*freq*t+phase））;          %单频噪声信号
%-------------vn---------------
vn=wgn（1n