1.一种基于压缩感知的抗间歇式主瓣干扰方法,其特征在于,所述方法包括:对间歇式干扰回波信号建模,获得雷达间歇式干扰回波信号模型,所述雷达间歇式干扰回波信号模型为:
R(n)=S(n)+Z(n)+E(n)其中,R(n)为实测回波信号,E(n)为雷达接收回路噪声,S(n)为原始发射信号由目标反射的理想回波,n为回波信号采样点数序号,Z(n)为间歇式干扰项;
利用自适应数字波束形成技术对所述雷达间歇式干扰回波信号模型中的副瓣干扰进行消除后,建立间歇式干扰稀疏回波信号模型,所述间歇式干扰稀疏回波信号模型为:其中,R为自适应数字波束形成技术处理后的雷达间歇式干扰回波信号, 为稀疏矢量,s为目标回波信号,z为干扰信号,A=[Φ,Id]为稀疏变换矩阵,Φ为傅里叶反变换基,Id为单位矩阵,E为对目标回波信号进行自相关运算;
将求解所述间歇式干扰稀疏回波信号模型转换为1‑范数规则化的最优化问题,获得转化后的间歇式干扰稀疏回波信号模型;
利用快速阈值迭代稀疏重构算法求解对所述转化后的间歇式干扰稀疏回波信号模型进行求解分析,消除间歇式主瓣干扰信号,输出间歇式主瓣干扰抑制后的目标回波信号。
2.根据权利要求1所述的方法,其特征在于,所述原始发射信号由目标反射的理想回波为:
S(n)=Asexp(2πnfc/fs)其中,As为回波信号幅度值,n为回波信号采样点数序号,fc为雷达工作频率,fs为采样频率。
3.根据权利要求1所述的方法,其特征在于,所述间歇式干扰项为:其中,P(n)为服从高斯分布的高斯白噪声,n为回波信号采样点数序号,Tp={[ti,ti+1]}i∈N为间歇式干扰时间区间,ti为第i个采样时间点,ti+1为第i+1个采样时间点,N为干扰信号采样点数, 为高斯白噪声, 为噪声方差,反映了雷达回波中除了干扰以外的其他噪声。
4.根据权利要求1所述的方法,其特征在于,所述自适应数字波束形成技术的处理公式为:
H
s.t.wa(θ0)=1
其中,y(t)为目标回波信号,E为对目标回波信号进行自相关运算,Pout为目标回波信号自相关运算结果,w为导向矢量加权系数,H表示转置计算,a(θ0)为目标导向矢量,θ0为目标方向角,s.t.为约束条件。
5.根据权利要求1所述的方法,其特征在于,所述1‑范数规则化的最优化问题求解公式为:
其中,L(s,z)为1‑范数规则化结果,λ1为回波信号范数加权系数,λ2为干扰信号范数加权系数,A为稀疏变换矩阵,R为自适应数字波束形成技术处理后的雷达间歇式干扰回波信号,x为稀疏矢量,s为目标回波信号,z为干扰信号。
6.根据权利要求1所述的方法,其特征在于,所述利用快速阈值迭代稀疏重构算法求解对所述转化后的间歇式干扰稀疏回波信号模型进行求解分析,消除间歇式主瓣干扰信号,获得间歇式主瓣干扰抑制后的目标回波信号的步骤,包括:T T
(1)初始参数设置,s0=[0,...,0] ,z0=[0,...,0];
(2)前向预测,Rapp=ifft(s0)+z0;
(3)残差计算,r=Rapp‑R;
(4)后向预测,
其中,s0为回波信号初始值,z0为干扰信号初始值,ifft(s0)为对回波信号初始值s0进行逆傅里叶变化,Rapp为前向预测结果,sapp为目标回波信号后向预测结果,zapp为干扰信号后向预测结果,r为残差,R为自适应数字波束形成技术处理后的雷达间歇式干扰回波信号,表示第j次迭代估计的回波信号, 表示第j+1次迭代估计的干扰信号,K为预测阈值系数,fft(r)表示对残差信号进行傅里叶变换;
(5)阈值处理,
t1=median(|sapp|)λ1t2=median(|zapp|)λ2其中,t1为回波信号阈值,t2为干扰信号阈值,median表示求中值运算,λ1为回波信号范数加权系数,λ2为干扰信号范数加权系数,sth为目标回波信号阈值处理结果,zth为干扰信号阈值处理结果,α表示 或 大于0时, 否则为0;
(6)更新系数,
其中,a(j)为第j次迭代时的阈值系数,a(j+1)为第j+1次迭代时的阈值系数。
(7)预测更新,
(8)回到第(2)步,直到r<σ为止,r<σ,σ为噪声标准方差;
(9)输出 为间歇式主瓣干扰抑制后的目标回波信号,为估计的间歇式主瓣干扰信号;
(10)去除估计的间歇式主瓣干扰信号 获得间歇式主瓣干扰抑制后的目标回波信号