适合参数估计的认知雷达最优波形设计方法
【专利摘要】本发明公开了一种适合参数估计的认知雷达最优波形设计方法，其利用更加符合实际情况的实际目标冲激响应构建信道模型，根据该信道模型，利用注水法得到所需波形的能量谱密度，可从回波中得到最大的目标信息量；以能量谱密度和时域复波形幅度恒定作为约束条件，利用迭代的加权最小二乘法得到恒幅的时域复波形，不会浪费发射机功率；雷达发射设计的时域复波形，得到目标回波，利用目标回波进行参数估计；一个参数估计积累时间达到后，利用当前的目标回波更新目标的冲激响应，并重新进行波形设计和发射，使当前发射的波形随目标和环境的变化而变化，始终保证当前得到的目标回波中包含最大的目标信息量，从而利于目标参数估计。
【专利说明】适合参数估计的认知雷达最优波形设计方法
【技术领域】
[0001]本发明涉及认知雷达波形设计研究领域，尤其涉及一种认知雷达最优估计波形设计方法。
【背景技术】
[0002]认知雷达赋予了雷达系统感知环境、理解环境、学习、推理并判断决策的能力，使雷达系统能够适应复杂多变的电磁环境，从而提升雷达系统的性能。认知雷达的主要特点是引入雷达闭环系统，包括两个部分:知识辅助和自适应波形设计。雷达通过目标回波信号处理来获取目标信息，而不同的发射波形会使目标回波中携带的目标信息量不同，对认知雷达而言体现在对不同目标采用与观测视角匹配的最优波形，从而提高雷达性能。因此波形设计对于认知雷达进行目标信息提取具有重要意义。
[0003]波形设计包括最优估计波形设计和最优检测波形设计，分别适用于参数估计和目标检测。基于信息论的波形设计方法体现了认知雷达感知环境的特点，使回波中包含的目标信息最大，以最大程度获取目标信息，更适合用于目标参数估计，进而达到目标识别的目的。本发明针对最优估计波形设计进行方案设计。
[0004]现有认知雷达最优估计波形设计方法的其中一种方案是:
[0005]1、假设已知目标冲激响应h(t)，构建如图1所示的信道模型；
[0006]2、根据信道模型，利用注水法得到所需设计波形的能量谱密度；
[0007]3、根据所述能量谱密度得到时域波形幅度。
[0008]上述估计波形设计方案的缺陷是:
[0009]1、通常情况下目标的冲激响应是未知的，而上述方案根据给定目标冲激响应或假定目标的目标冲激响应构建信道模型，不符合实际情况；
[0010]2、最终合成的时域复波形不仅包括幅度，还包括相位，但是目前公开的文件中或是未提及如何分配相位(一般未提及)；或是参照选择的参考信号的相位来分配相位，但波形性能与时域波形的恒幅性能需要折衷(仅一人)。
【发明内容】
[0011]有鉴于此，本发明提供了一种适合参数估计的认知雷达最优估计波形设计方法，该方法得到了最优估计雷达时域波形。
[0012]本发明的适合微动参数估计的认知雷达最优波形设计方法包括:
[0013]一种适合参数估计的认知雷达最优波形设计方法，包括如下步骤:
[0014]步骤1、雷达发射线性调频信号得到目标回波，进而得到目标的冲激响应；
[0015]步骤2、基于目标的冲激响应构建雷达发射信号经由目标反射到雷达接收机的信道模型；根据该信道模型，利用注水法得到所需设计波形的能量谱密度；
[0016]步骤3、在初始随机分配波形幅度谱序列中各元素相位的基础上，以所述能量谱密度和时域复波形幅度恒定作为约束条件，利用迭代的加权最小二乘法得到恒幅的时域复波形，且每次迭代后均由该时域复波形的傅里叶变换得到波形幅度谱序列各元素新的相位；所述波形幅度谱为能量谱密度开根号计算得到；
[0017]步骤4、雷达发射设计的时域复波形，得到目标回波，利用目标回波进行参数估计；
[0018]步骤5、一个参数估计积累时间达到后，利用当前的目标回波更新目标的冲激响应，返回步骤2。
[0019]优选地，步骤3中，所述利用迭代的加权最小二乘法得到恒幅的时域复波形的步骤具体为:
[0020]设，给定的波形能量谱为P(k)，其中，k = 1，2，3，…，K_l，k为频域采样点，K为频域总采样点数，则真实的幅度谱为
【权利要求】
1.一种适合参数估计的认知雷达最优波形设计方法，其特征在于，包括: 步骤1、雷达发射线性调频信号得到目标回波，进而得到目标的冲激响应； 步骤2、基于目标的冲激响应构建雷达发射信号经由目标反射到雷达接收机的信道模型；根据该信道模型，利用注水法得到所需设计波形的能量谱密度； 步骤3、在初始随机分配波形幅度谱序列中各元素相位的基础上，以所述能量谱密度和时域复波形幅度恒定作为约束条件，利用迭代的加权最小二乘法得到恒幅的时域复波形，且每次迭代后均由该时域复波形的傅里叶变换得到波形幅度谱序列各元素新的相位；所述波形幅度谱为能量谱密度开根号计算得到； 步骤4、雷达发射设计的时域复波形，得到目标回波，利用目标回波进行参数估计； 步骤5、一个参数估计积累时间达到后，利用当前的目标回波更新目标的冲激响应，返回步骤2。
2.如权利要求1所述的方法，其特征在于，步骤3中，所述利用迭代的加权最小二乘法得到恒幅的时域复波形的步骤具体为: 设，给定的波形能量谱为P(k)，其中，k = 1,2,3,…，K-l，k为频域采样点，K为频域总采样点数，则真实的幅度谱为
假设设计的时域复波形为￡/(/?) = η = 0,1，…，JV-1，N≤K, η为时域采样点，N为时域总采样点数，M表示该时域复波形的模，且根据帕塞瓦尔定理可得
采样点η处时域复波形的相位； A(k) (k = 0，1，…，K-ι)是a (η)的K点傅里叶变换，即
指傅里叶变换，令 A (k) = A(k) eJ$(k), s =
[S(0)eJ$(0),S(l)eJ$(1), -,Sa-De^^]1, a = [a(0)，a(l)，…，a(N_l) ]T，A = [A(O),
A(I)，…，A(K-1)]T，KXN 维的 DFT 矩阵为:
爾的上角标表示指数kn,则A = Da ; 幅度谱估计的加权最小二乘误差为J(9) = (s-Da)Hff(s-Da) #为权矩阵，并规定为正定矩阵；所以a的线性加权最小二乘估计为:a = (D1WD) 1DllWs ； 利用迭代法求解a，迭代的收敛条件设定为
上角标⑴和(1-1)均为迭代次数，常数C为设定的门限； 迭代过程具体步骤如下: 步骤S1、在O~2π范围内随机选取数值给ΦΟΟ赋初值为
,上标为迭代次数，当前为O ;迭代次数加I,即令i = I ;
步骤 S2、计算 s = [S (0)eJ4> (。)，S (l)eJ4> ⑴，,S (N) eJφ 妇)],利用 a = (D11WD) 1 D11Ws得到第i次时域复信号的估计Sui ； 步骤S3、把I?中的每一个样本值除以它们本身的模再乘以M，得到具有模为M的时域复信号的估计I5、 步骤S4、利用计算得到由于A(k) = |A(k)|ej<Hk)，由可得到Mk)(k = O, I,…，K-1)的第 i 次估计值 Φ (1) (k) (k = O, I,…，K-1)； 步骤S5、判断所述收敛条件是否满足；若不满足，则令i自加1，重复步骤S2?S5 ;若满足收敛条件，则停止迭代，得到最终的时域复波形估计-- 0
【文档编号】G01S7/02GK104198993SQ201410367813
【公开日】2014年12月10日 申请日期:2014年7月29日 优先权日:2014年7月29日
【发明者】傅雄军, 李婷, 高梅国, 王才, 赵会朋, 李文静 申请人:北京理工大学