值，N(f0)为同频背景噪声平均功率，计 算过程如下，
[0050] (C1)对输入信号和输出信号采样得到的离散序列x(n)和y(n)做FFT，记为X( ω )和 Υ(ω)；
[0051] (C2)计算输出信噪比，设信号频率fo处谱峰值的序号为ω〇、采样频率为fs，则fo = (ω o/N-1 ))fs，S(fo) = |Υ( ω〇) 12 ,Ν为义样点数;N(fo)是 ω 0附近 ω 0+M,, ω 0+1, ω 〇-1,… ω ο-Μ处的平均功率，其计算公式为：
[0052]
[005引其中，Μ是根据采样点数Ν的大小适当选择的整数，Μ为^ ，
[0054] ，采样点数Ν越大Μ就越大。
[0055] 步骤(D)，对初始种群Q(to)再实施一次测量，得到一组状态R(t)，对其进行适应度 评估，依据调整策略对种群个体执行量子交叉操作，采用量子旋转口 U(t)更新种群为Q(t)， 利用量子非口执行量子变异操作，种群Q(t)经过运算后得到下一代群体Q(t+1);其中，量子 交叉操作:利用量子的想干特性构造量子间联合交叉，交叉算子的作用是实现个体的信息 交换，即充分利用染色体信息，一旦算法出现早熟现象，它就可W带来新个体产生，具体操 作为:①所有种群个体随机排序;②取当前个体的第一个基因作为新生个体的第一个基因， 依此取相邻个体的第二个基因为新个体第二个基因，循环往复直到新个体具有相同基因 数;③往返如此，直到新种群具有相同规模；
[0056] 量子旋转口更新种群操作:设f(xj)为当前测量值的适应度值，f(bj)为当前目标适 应度值，根据测量结果和调整策略，调整个体 <对应的量子比特，使概率幅对(〇1，01)朝着 有利于XI或bi的方向演化，此方法可W使算法收敛；
[0057] 量子变异操作:变异的作用既能使算法有局部捜索能力还能维持种群多样性，目 的是防止程序未成熟收敛，该操作实际上更改了该量子比特态叠加的状态，使得倾向于巧 缩到状态"Γ的变为倾向于巧缩到状态"0"，或者相反，具体步骤为:①根据变异概率化从种 群中选取若干个体;②对所选个体按确定的概率确定一个或多个变异位;③对选中若干的 量子位执行量子非口操作，使得该量子位的两个概率幅互相变换，运样的目的是使两条基 因链能够同时进行变异。显然，该操作保证了算法能捜索到空间每一点，从而使寻优结果具 有全局最优；
[0058] 步骤化），记录当前的最优解，其对应适应度值为下一次迭代的目标值；
[0059] 步骤(F)，将进化代数t = t+l，若t<T则跳到步骤(4);若*>1'，则转而判断SNRI是 否大于1，若SNRK1，则转步骤(2)，否则输出最大适应度值Ybest及其个体向量Xbest作为最优 解输出，个体向量Xbest为输入输出信噪比差值最大时a，b，k的值；
[0060] 其中，SNRI为随机共振信噪比增益，为输出信噪比和输入信噪比的比值，如公式 (3)所示，
[0061] C3)
[0062] 其中，門估633。。*为输出信噪比、門估633山为输入信噪比、5(时)。。*为输出信号功率 谱、N(f〇)Dut为输出噪声平均功率、S(f〇)in为输入信号功率谱、N(f〇)in为输入噪声平均功率；
[0063] 步骤(G)，将优化得到的最大适应度值Ybest及其个体向量Xbest，输入Duffing随机共 振模型，WIPIX雷达海杂波实测数据为干扰背景，实现此干扰背景中自适应检测微弱小目 标信号，包括W下步骤，
[0064] (Gl)Duffing系统随机共振模型，如公式(4)所示，
[00 财
（4)
[0066] 将公式(4)写成系统的形式，如公式(5)所示，
[0067]
巧）
[0068] 其中，激励信号e(t)为IPIX雷达海杂波实测数据（由McMaster大学提供的IPIX雷 达实测海杂波数据），其发射频率为9.39GHz，脉冲重复频率为1000化，采用VV极化方式，每 组数据包括131072个采样点；-ax+bx 3为非线性恢复力；:?为X的二次导数、i：为X的一次导 数、k是阻尼比、U为公式(4)转化成公式(5)的转换系数，&为U的一次导数；
[0069] (G2)根据公式(5)，得到系统势函数U(x)，如公式(6)所示，
[0070]
(6)
[0071] 其中，公式(6)描述了一个含有两个极小值，= 和一个极大值χ = 〇的双稳 系统，其的势垒高度为A U = a2/4b;
[0072] (G3)在激励信号e(t)输入时，得到系统势函数一临界值
当 信号幅值4<4。时，信号能量无法克服势垒的阻挡，系统势函数输出状态只能在某势阱附近 做局域周期运动；当信号幅值A>A。时，信号、噪声达到协同作用，一部分噪声能量将转移到 信号身上，系统势函数输出会在两稳态间大范围跃迁，进入SR状态，此时，部分海杂波能量 化为小目标信号能量，从而实现微弱信号检测。
[0073] W上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该 了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原 理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，运些变化和改进 都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界 定。
【主权项】
1.基于二维Duffing振子的海杂波微弱信号检测方法，其特征在于:包括以下步骤， 步骤(A)，QGA初始化，随机生成个体为η的种群种群中染色 体所有基因都被初始化为即染色体的全部可能状态等概率叠加；设置 Duff ing振子的参数范围；t为进化代数，初始化进化代数计数器t = 1，并设置最大进化代数 T； 其中，< 为第t代第j个染色体，如公式(1)所示，其中，m为染色体基因个数，g为编码每个基因量子比特数；炉和β分别为I 0>和11>的概率 幅，且满足归一化条件IfP+i PI2=I 表示分别为第m个染色体基因编码基因量子 比特数为g的I 〇>、11>状态的概率幅； Duff ing 振子的参数a，b，k范围设置分别为：ae [0.001,3· 5]，be [0.001，3.5]，ke [0.001，2.5]，a、b是Duff ing振子的双稳态结构参数，k是阻尼比，寻优精度为0.001; QGA初 始化参数设置为:种群数量G = 50，最大迭代步数Tmax = 50，维度d = 3; 步骤（B)，对初始种群中个体进行测量，由种群Q( to)量子坍塌获得一组状态其中，为第t代种群中第j个个体的测量值，测量过程为随机产 生一个0到1之间的数，若小于概率幅的平方，则测量结果为0,反之为1;<的表现形式为长 度是m的二进制串； 步骤(C)，计算个体适应度值，用适应度函数fitness对R(to)中每个个体进行评价，保留 最优个体状态及其适应度值作为初始目标值，所述适应度函数fitness是指Duffing振子随 机共振系统输出信噪比，如公式(2)所示，其中，S(fo)为信号功率谱在频率fo处的幅值，N(fo)为同频背景噪声平均功率； 步骤(D)，对种群Q(to)再实施一次测量，得到一组状态R(t)，对其进行适应度评估，依据 调整策略对种群个体执行量子交叉操作，采用量子旋转门U(t)更新种群为Q(t)，利用量子 非门执行量子变异操作，种群Q(t)经过运算后得到下一代群体Q(t+1); 步骤(E)，记录当前的最优解，其对应适应度值为下一次迭代的目标值； 步骤(F)，将进化代数t = t+l，若t < T则跳到步骤(4);若t>T，则转而判断SNRI是否大 于1，若SNRKl，则转步骤(2)，否则输出最大适应度值Ybest及其个体向量Xbest作为最优解输 出，个体向量Xb est为输入输出信噪比差值最大时a，b，k的值； 其中，SNRI为随机共振信噪比增益，为输出信噪比和输入信噪比的比值，如公式(3)所 示，其中，fitnesses为输出信噪比、fitnessinS输入信噪比、S(fo)cmt为输出信号功率谱、N (fo)cmt为输出噪声平均功率、S(f〇)in为输入信号功率谱、N(f〇)ir*输入噪声平均功率； 步骤(G)，将优化得到的最大适应度值Ybest及其个体向量Xbest,输入Duffing随机共振模 型，以IPIX雷达海杂波实测数据为干扰背景，实现此干扰背景中自适应检测微弱小目标信 号，包括以下步骤， (G1 )Duff ing系统随机共振模型，如公式(4)所示，I 将公式(4)写成系统的形式，如公式(5)所示，其中，激励信号e(t)为IPIX雷达海杂波实测数据，其发射频率为9.39GHz，脉冲重复频 率为1000Hz，采用VV极化方式，每组数据包括131072个采样点；-ax+bx3为非线性恢复力；X 为X的二次导数、i为X的一次导数、u为公式(4)转化成公式式(5)的转换系数，为u的一次 导数； (G2)根据公式(5)，得到系统势函数U(X)，如公式(6)所示，I 其中，公式(6)描述了一个含有两个极小值，和一个极大值X = O的双稳系统， 其的势皇高度为AU = a2/4b; (G3)在激励信号e(t)输入时，得到系统势函数一临界值当信号幅 值八<1时，信号能量无法克服势皇的阻挡，系统势函数输出状态只能在某势阱附近做局域 周期运动；当信号幅值A>A。时，信号、噪声达到协同作用，一部分噪声能量将转移到信号身 上，系统势函数输出会在两稳态间大范围跃迀，进入SR状态，此时，部分海杂波能量化为小 目标信号能量，从而实现微弱信号检测。2.根据权利要求1所述的基于二维Duffing振子的海杂波微弱信号检测方法，其特征在 于:步骤(C)，信号功率谱在频率fo处的幅值S(fo)和同频背景噪声平均功率N(fo)的计算过 程如下， (Cl)对输入信号和输出信号采样得到的离散序列x(n)和y(n)做FFT，记为Χ( ω )和Y (ω ); (C2)计算输出信噪比，设信号频率fo处谱峰值的序号为ω 〇、采样频率为fs，则&= ( ω 〇/ N-l))fs，S(fo)= |Υ( ω〇) |2，Ν 为采样点数;N(fo)是 ω〇 附近 ω〇+Μ，···，ω〇+1，ω〇-1，…ω〇-Μ 处的平均功率，其计算公式为：其中，M是根据采样点数N的大小适当选择的整数，M为# 2 Q
【专利摘要】本发明公开了一种基于二维Duffing振子的海杂波微弱信号检测方法，结合量子遗传算法和随机共振算法，增强海杂波信号有用目标信号，得到最优检测结果的方法，采用改进的遗传算法优化Duffing混沌振子，以输出信噪比为评价函数，以信噪比增益为判断指标，对系数联合编码智能获取最佳系统参数，根据所得参数对IPIX雷达实测海杂波数据进行自适应随机共振处理，在不损坏有用信号的同时，高效检测目标信号，本发明抗噪能力强，大幅度提高了系统输出信噪比，不但克服传统随机共振手动设置参数或单参数寻优检测效果不佳的缺陷，而且扩大了Duffing混沌振子随机共振在混沌领域检测小目标信号的应用范围，具有良好的应用前景。
【IPC分类】G01S7/41
【公开号】CN105487062
【申请号】CN201510970319
【发明人】行鸿彦, 卢春霞, 张强
【申请人】南京信息工程大学
【公开日】2016年4月13日
【申请日】2015年12月22日