[发明专利]基于二维Duffing振子的海杂波微弱信号检测方法

摘要

本发明公开了一种基于二维Duffing振子的海杂波微弱信号检测方法，结合量子遗传算法和随机共振算法，增强海杂波信号有用目标信号，得到最优检测结果的方法，采用改进的遗传算法优化Duffing混沌振子，以输出信噪比为评价函数，以信噪比增益为判断指标，对系数联合编码智能获取最佳系统参数，根据所得参数对IPIX雷达实测海杂波数据进行自适应随机共振处理，在不损坏有用信号的同时，高效检测目标信号，本发明抗噪能力强，大幅度提高了系统输出信噪比，不但克服传统随机共振手动设置参数或单参数寻优检测效果不佳的缺陷，而且扩大了Duffing混沌振子随机共振在混沌领域检测小目标信号的应用范围，具有良好的应用前景。

主权项

1.基于二维Duffing振子的海杂波微弱信号检测方法，其特征在于：包括以下步骤，步骤(A)，QGA初始化，随机生成个体为n的种群种群中染色体所有基因都被初始化为即染色体的全部可能状态等概率叠加；设置Duffing振子的参数范围；t为进化代数，初始化进化代数计数器t＝1，并设置最大进化代数T；其中，为第t代第j个染色体，如公式(1)所示，其中，m为染色体基因个数，g为编码每个基因量子比特数；和β分别为|0>和|1>的概率幅，且满足归一化条件β_mg表示分别为第m个染色体基因编码基因量子比特数为g的|0>、|1>状态的概率幅；Duffing振子的参数a，b，k范围设置分别为：a∈[0.001，3.5]，b∈[0.001，3.5]，k∈[0.001，2.5]，a、b是Duffing振子的双稳态结构参数，k是阻尼比，寻优精度为0.001；QGA初始化参数设置为：种群数量G＝50，最大迭代步数Tmax＝50，维度d＝3；步骤(B)，对初始种群中个体进行测量，由种群Q(t₀)量子坍塌获得一组状态其中，为第t代种群中第j个个体的测量值，测量过程为随机产生一个0到1之间的数，若小于概率幅的平方，则测量结果为0，反之为1；的表现形式为长度是m的二进制串；步骤(C)，计算个体适应度值，用适应度函数fitness对R(t0)中每个个体进行评价，保留最优个体状态及其适应度值作为初始目标值，所述适应度函数fitness是指Duffing振子随机共振系统输出信噪比，如公式(2)所示，其中，S(f0)为信号功率谱在频率f0处的幅值，N(f0)为同频背景噪声平均功率；步骤(D)，对种群Q(t0)再实施一次测量，得到一组状态R(t)，对其进行适应度评估，依据调整策略对种群个体执行量子交叉操作，采用量子旋转门U(t)更新种群为Q(t)，利用量子非门执行量子变异操作，种群Q(t)经过运算后得到下一代群体Q(t+1)；步骤(E)，记录当前的最优解，其对应适应度值为下一次迭代的目标值；步骤(F)，将进化代数t＝t+1，若t≤T则跳到步骤(4)；若t＞T，则转而判断SNRI是否大于1，若SNRI＜1，则转步骤(2)，否则输出最大适应度值Ybest及其个体向量Xbest作为最优解输出，个体向量Xbest为输入输出信噪比差值最大时a，b，k的值；其中，SNRI为随机共振信噪比增益，为输出信噪比和输入信噪比的比值，如公式(3)所示，其中，fitnessout为输出信噪比、fitnessin为输入信噪比、S(f0)out为输出信号功率谱、N(f0)out为输出噪声平均功率、S(f0)in为输入信号功率谱、N(f0)in为输入噪声平均功率；步骤(G)，将优化得到的最大适应度值Ybest及其个体向量Xbest，输入Duffing随机共振模型，以IPIX雷达海杂波实测数据为干扰背景，实现此干扰背景中自适应检测微弱小目标信号，包括以下步骤，(G1)Duffing系统随机共振模型，如公式(4)所示，将公式(4)写成系统的形式，如公式(5)所示，其中，激励信号e(t)为IPIX雷达海杂波实测数据，其发射频率为9.39GHz，脉冲重复频率为1000Hz，采用VV极化方式，每组数据包括131072个采样点；‑ax+bx³为非线性恢复力；为x的二次导数、为x的一次导数、u为公式(4)转化成公式式(5)的转换系数，为u的一次导数；(G2)根据公式(5)，得到系统势函数U(x)，如公式(6)所示，其中，公式(6)描述了一个含有两个极小值，和一个极大值x＝0的双稳系统，其的势垒高度为ΔU＝a²/4b；(G3)在激励信号e(t)输入时，得到系统势函数一临界值当信号幅值A＜A_c时，信号能量无法克服势垒的阻挡，系统势函数输出状态只能在某势阱附近做局域周期运动；当信号幅值A＞A_c时，信号、噪声达到协同作用，一部分噪声能量将转移到信号身上，系统势函数输出会在两稳态间大范围跃迁，进入SR状态，此时，部分海杂波能量化为小目标信号能量，从而实现微弱信号检测；其中，所述步骤(C)，信号功率谱在频率f0处的幅值S(f0)和同频背景噪声平均功率N(f0)的计算过程如下，(C1)对输入信号和输出信号采样得到的离散序列x(n)和y(n)做FFT，记为X(ω)和Y(ω)；(C2)计算输出信噪比，设信号频率f0处谱峰值的序号为ω0、采样频率为fs，则f0＝(ω0/N‑1)fs，S(f0)＝|Y(ω0)|2，N为采样点数；N(f0)是ω0附近ω0+M，…，ω0+1，ω0‑1，…ω0‑M处的平均功率，其计算公式为：其中，M是根据采样点数N的大小适当选择的整数，M为