发明名称 一种低旁瓣随机跳频脉冲信号的调制方法
摘要 一种低旁瓣随机跳频脉冲信号的调制方法,涉及一种跳频脉冲信号的调制方法。解决了传统低旁瓣的随机跳频脉冲信号调制方法存在局部精确寻优能力差,收敛速度慢而导致随机跳频脉冲信号调制速度慢的问题,本发明先求出随机跳频脉冲信号的自相关函数,利用自相关函数计算随机跳频脉冲信号的偏差函数计算初始化频率组合,将初始化频率组合带入偏差函数获得最优初始频率组合再加入随机扰动,更新偏差函数,对更新的偏差函数进行检测,利用偏差函数值最大值点对应的时刻按梯度原理对频率组合进行a阶段频率调整,直到相邻两个阶段获得的偏差函数的旁瓣水平的差值小于设定门限时停止频率组合更新,完成随机跳频脉冲信号调制。本发明适用于信号调制领域。
申请公布号 CN103138799B 申请公布日期 2014.12.17
申请号 CN201310092791.0 申请日期 2013.03.21
申请人 哈尔滨工业大学 发明人 位寅生;赵德华;邓维波
分类号 H04B1/7176(2011.01)I;G01S7/282(2006.01)I 主分类号 H04B1/7176(2011.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 胡树发
主权项 一种低旁瓣随机跳频脉冲信号的调制方法,其特征在于,它的调制方法由以下步骤实现: 步骤一、待调制的随机跳频脉冲信号s(t)根据公式: <img file="FDA0000547960780000011.GIF" wi="1537" he="172" />求出待调制的随机跳频脉冲信号s(t)的自相关函数χ(τ),其中,T为脉冲宽度,N为随机跳频脉冲信号s(t)在一个信号周期内脉冲的个数,f<sub>n</sub>为在该信号周期内的第n个脉冲的频率,f<sub>m</sub>为在该信号周期内的第m个脉冲的频率; 步骤二、将步骤一获得的自相关函数χ(τ),根据公式: E(τ)=[χ(τ)‑d(τ)]<sup>2</sup>    (2) 获得待调制的随机跳频脉冲信号的偏差函数E(τ),其中,d(τ)为具有期望旁瓣结构特点的目标函数; 步骤三、利用公式: <img file="FDA0000547960780000012.GIF" wi="1251" he="155" />在[0,2]内依次以步进值为0.1选择21次k值,获得21组初始化频率组合;式中,N为一个信号周期内的脉冲个数,B为随机跳频脉冲信号的带宽,i为迭代次数; 步骤四、将获得的21组初始化频率组合分别代入步骤二中获得的待调制的随机跳频脉冲信号的偏差函数E(τ),共获得21个偏差函数E(τ),将其中旁瓣水平最低的偏差函数E(τ)对应的初始频率组合作为最优初始频率组合f<sub>0,1</sub>、f<sub>0,2</sub>、…、f<sub>0,N</sub>; 步骤五、向步骤四求出的最优初始化频率组合f<sub>0,1</sub>、f<sub>0,2</sub>、…、f<sub>0,N</sub>依次对应加入微小随机扰动c<sub>1</sub>、c<sub>2</sub>、…、c<sub>N</sub>,获得非对称的最优初始化频率组合;所述微小的随机扰动c<sub>1</sub>、c<sub>2</sub>、…、c<sub>N</sub>中的c<sub>2</sub>、…、c<sub>N‑1</sub>服从均值为零且均方差为B/(20N)的正态分布,c<sub>1</sub>和c<sub>N</sub>的值为零; 步骤六、将步骤五获得的非对称的最优初始化频率组合代入步骤二获得的待调制的随 机跳频脉冲信号的偏差函数E(τ),获得更新后待调制的随机跳频脉冲信号的偏差函数E(τ); 步骤七、检测步骤六中更新后待调制的随机跳频脉冲信号的偏差函数E(τ)输出的峰值点,得到除去主峰点后待调制的随机跳频脉冲信号的偏差函数E(τ)的峰值,同时检测偏差函数E(τ)在主瓣宽度预设值D处的输出值,比较该输出值与除去主峰点后待调制的随机跳频脉冲信号的偏差函数E(τ)的峰值获得最大值,从而获得最大值对应的时刻τ<sub>i,max</sub>; 步骤八、将偏差函数E(τ)的最大值对应的时刻τ<sub>i,max</sub>按梯度原理对频率组合进行a个阶段的频率调整,a为大于等于1的正整数,每个阶段获得一组调整后的频率组合,共获得a组调整后的频率组合; 获得a组调整后的频率组合的方法为: 步骤A、将偏差函数E(τ)的最大值对应的时刻τ<sub>i,max</sub>,根据公式: <img file="FDA0000547960780000021.GIF" wi="1611" he="263" />获得当前频率组合<img file="FDA0000547960780000022.GIF" wi="48" he="79" />的调整量<img file="FDA0000547960780000023.GIF" wi="690" he="98" />Δf<sub>i,j</sub>为第j个脉冲第i次迭代的频率调整量;式中η为步进因子,<img file="FDA0000547960780000024.GIF" wi="51" he="79" />为第i次迭代的频率组合;步骤B、利用步骤A获得的当前频率组合<img file="FDA0000547960780000025.GIF" wi="46" he="79" />的调整量<img file="FDA0000547960780000026.GIF" wi="99" he="79" />根据公式:<img file="FDA0000547960780000027.GIF" wi="298" he="79" />对频率组合进行L次更新,L=100,按最速下降法对步进因子η进行选择,通过公式: <img file="FDA0000547960780000028.GIF" wi="1053" he="163" />对频率组合<img file="FDA0000547960780000029.GIF" wi="46" he="79" />进行L次调整,获得L个频率组合的调整量<img file="FDA00005479607800000210.GIF" wi="101" he="80" />L=100,对每次调整获得的最大频率调整量<img file="FDA00005479607800000211.GIF" wi="326" he="93" />进行统计并求均值,获得L个最大频率调整量<img file="FDA00005479607800000212.GIF" wi="326" he="93" />的均值M;步骤C、以<img file="FDA00005479607800000213.GIF" wi="519" he="93" />作为最大调整量,将<img file="FDA00005479607800000214.GIF" wi="496" he="93" />带入公 式(4),对频率组合<img file="FDA0000547960780000031.GIF" wi="44" he="79" />进行Q次梯度调整,Q≥100,对应获得Q个频率组合的调整量的<img file="FDA0000547960780000032.GIF" wi="99" he="80" />并利用公式<img file="FDA0000547960780000033.GIF" wi="262" he="78" />对频率组合进行更新,获得Q个频率组合;步骤D、将步骤C获得的Q个频率组合分别代入待调制的随机跳频脉冲信号的偏差函数E(τ),获得Q个偏差函数E(τ),记录Q个偏差函数E(τ)中出现旁瓣水平的最优的结果,获得偏差函数E(τ)的精英保留曲线,直至连续W(W=100)次未出现水平度更低的结果时此阶段调整结束,选择偏差函数E(τ)中出现旁瓣水平的最优结果的频率组合作为第a阶段调整后的频率组合,步骤九、将步骤八获得的a组调整后的频率组合代入步骤二获得的待调制的随机跳频脉冲信号的偏差函数E(τ),获得a个偏差函数E(τ); 步骤十、判断步骤九获得的a个偏差函数E(τ)相邻两个偏差函数E(τ)的旁瓣水平的差值是否小于设定门限ΔPSL,如果判断结果为是,则停止频率组合更新,将更新后的频率组合作为调制结果,完成随机跳频脉冲信号调制;否则,令a=a+1,进行下一阶段调整,返回执行步骤八。 
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号