本发明专利技术公开了一种基于混沌量子遗传算法的随机共振微弱信号检测方法，将混沌量子遗传算法(Chaotic Quantum Genetic Algorithm,CQGA)运算得到的优化参数代入到由α稳定噪声驱动的Duffing双稳随机共振系统中，进行微弱信号检测，其中，CQGA是一种利用混沌搜索改进的量子遗传算法，以量子运算为基础，将量子比特的概率幅应用于遗传编码并通过量子的态矢量表达出来，执行混沌搜索操作更新个体染色体，本发明专利技术利用噪声能量增强信号能量，实现信噪比的大幅提高，不依靠更新表查找计算，满足高精度匹配要求，基于α稳定噪声背景拓宽了随机共振的应用范围，给非高斯噪声背景下的小目标检测提供了一种新方法。
【技术实现步骤摘要】
本专利技术属于微弱信号检测领域，具体涉及一种基于混沌量子遗传算法的随机共振微弱信号检测方法。
技术介绍
信号处理中，总是想方设法滤除或抑制噪声以保留有用信号，因此，强噪声背景下的微弱信号检测，是一种专门与噪声作斗争的技术，现代电子学领域，如通信、控制、广播、遥控遥测或其他电子系统，都存在处理微弱信号和噪声的问题。因此，有效提高信噪比是改善微弱信号检测技术的首要任务，进而满足现代科学研究和技术开发的需要。基本量子遗传算法(BasicQuantumGeneticAlgorithmBQGA)是一种改进的遗传算法，将量子计算引入到遗传算法中，结合两者优点使遗传算法具有量子的独特计算特性，在较小的规模内可以快速收敛到全局最优解，具有寻优能力强、计算时间短、收敛精度高等优点。虽然采用量子门更新种群操作在一定程度上改变了旋转角，但靠表查找计算复杂度高，不适合随机共振系统参数高精度匹配要求。以往的自适应随机共振方法中，其研究基本建立在高斯白噪声模型背景下。而基于高斯背景下的随机共振研究，一般假设高斯噪声理想化，但是在工程中经常遇到具有显著尖峰脉冲特性和拖尾特性的非高斯噪声(如：海杂波)。因此，拓宽随机共振的应用范围，对非高斯噪声背景下的小目标检测是急需解决的问题。
技术实现思路
针对上述所提出的问题，本专利技术为了提高随机共振在微弱信号检测领域的实用性以及增强信号检测精度和系统参数高精度匹配要求，以及拓宽随机共振的应用范围，提出了一种基于混沌量子遗传算法的随机共振微弱信号检测方法。为实现上述目的，本专利技术采用的技术方案为：一种基于混沌量子遗传算法的随机共振微弱信号检测方法，其特征在于：将经过混沌量子遗传算法(ChaoticQuantumGeneticAlgorithm,CQGA)运算得到的优化参数代入到由α稳定噪声驱动的Duffing双稳随机共振系统中，进行微弱信号检测，其中，CQGA是一种利用混沌搜索改进的量子遗传算法，以量子运算为基础，将量子比特的概率幅应用于遗传编码并通过量子的态矢量表达出来，执行混沌搜索操作更新个体染色体；其具体检测方法包括以下几个步骤：步骤(A)，量子混沌初始化及参数设置，在待寻优参数可行解空间中混沌产生个体为n的初始种群根据背景噪声确定Duffing随机共振系统参数范围，设置种群规模W，最大进化代数T，初始化进化代数计数器t＝1；其中，一个基因代表量子比特进行量子遗传算法计算，如公式(1)所示：式(1)中，和β分别为‘0’状态|0>和‘1’状态|1>的概率幅，且满足归一化条件染色体概率幅如公式(2)所示：式(2)中，为第t代第j个染色体，m为染色体基因个数，g为每个基因量子比特数，表示‘0’态概率幅具有m个染色体基因g个量子比特数，βmg表示‘1’态概率幅具有m个染色体基因g个量子比特数；步骤(B)，对种群实施一次测量，由初始种群量子坍塌获得一组状态其中，为第t代种群中第j个个体的测量值，其表现形式为长度是m的二进制串，测量过程为随机产生一个0到1之间的数，若小于概率幅的平方，则测量结果为0，反之为1；步骤(C)，将当前状态不同测量值分别载入由α稳定噪声驱动Duffing双稳随机共振系统计算输出信噪比SNR，比较不同测量值对应的信噪比，将信噪比最大的个体保留，作为初始目标向量；其中，所述由α稳定噪声驱动Duffing双稳随机共振系统计算输出信噪比SNR，如公式(3)所示：式(3)中，S(f0)为信号功率谱在频率f0处的幅值，N(f0)为同频背景噪声平均功率；步骤(D)，对Q(t0)实施一次测量并对其个体适应度评价，执行构造量子间联合交叉的量子交叉操作以及更改量子比特态叠加状态的量子变异操作，进而产生下一代个体Xinest；步骤(E)，对个体Xinest中各分量执行混沌搜索操作得到Xibest，计算其适应度值，如果Xibest位置优于Xinest，则替换Xinest，记录当前信噪比Yibest及其个体状态对应个体Xibest的参数值a,b,k；其中，混沌搜索是将混沌状态引入到优化变量中，混沌搜索中常用到的是Logistic映射，其迭代公式如下：zi+1＝μzi(1-zi),i＝0,1,2,…(4)式(4)中，μ∈[0,4]为控制参数，初值z0∈[0,1]；步骤(F)，若遗传代数t≤T则跳到步骤(D)，否则；将最大适应度值Ybest及其个体向量Xbest即信噪比最大时a,b,k的值作为最优解输出；步骤(G)，将优化后的最佳匹配参数，即个体向量Xbest包含的3个元素值a,b,k载入由α稳定噪声驱动的Duffing双稳随机共振系统，优化后的最佳匹配参数使非线性系统与激励信号实现了协同效应，产生共振输出，大幅提高输出信噪比，从而实现α稳定噪声背景下的自适应随机共振微弱信号检测；其中，Duffing随机共振系统模型为：式(5)中k为阻尼比，-ax+bx3为非线性恢复力，e(t)＝s(t)+n(t)表示信号和噪声的混合激励信号，由α稳定噪声驱动的Duffing双稳随机共振系统模型为：式(6)中，k为阻尼比，-ax+bx3为非线性恢复力，为x的二次导数、为x的一次导数，s(t)＝Asin(2πf0t)是驱动信号幅值为A频率为f0的待测信号，xDξ(t)为α稳定噪声，D代表α稳定噪声的强度放大系数；其双稳系统势函数为：式(7)中描述了一个含有两个极小值和一个极大值x3＝0的双稳系统，其势垒高度为ΔU＝a2/4b，c,d为稳态系统参数，以上x表示粒子的位置，即考虑一颗具有动力学行为的粒子，在具有双稳态势能的系统中运动过程；当Duffing系统有信号输入时，系统存在临界值当信号幅值A<Ac时，系统输出将在某势阱附近做局域周期运动；当A>Ac时，非线性系统和激励信号达到协同作用，一部分噪声能量向信号能量转移，系统输出在式(7)双稳系统的两稳态间大范围跃迁；其中，α稳定噪声数据通常由α稳定分布序列产生，α稳定分布通常用特征函数表达：φ(t)＝exp{jut-γ|t|α[1+jBsign(t)ω(t,α)]本文档来自技高网...
【技术保护点】
基于混沌量子遗传算法的随机共振微弱信号检测方法，其特征在于：将经过混沌量子遗传算法运算得到的优化参数a,优化参数b,优化参数k代入到由α稳定噪声驱动的Duffing双稳随机共振系统中，进行信号检测；其具体检测方法包括以下几个步骤：步骤(A)，量子混沌初始化及参数设置，在待寻优参数可行解空间中混沌产生个体为n的初始种群根据背景噪声确定Duffing随机共振系统参数范围，设置种群规模W，最大进化代数T，初始化进化代数计数器t＝1；其中，一个基因代表量子比特进行量子遗传算法计算，如公式(1)所示：式(1)中，和β分别为‘0’状态|0>和‘1’状态|1>的概率幅，且满足归一化条件染色体概率幅如公式(2)所示：式(2)中，为第t代第j个染色体，m为染色体基因个数，g为每个基因量子比特数，表示‘0’态概率幅具有m个染色体基因g个量子比特数，βmg表示‘1’态概率幅具有m个染色体基因g个量子比特数；步骤(B)，对种群实施一次测量，由初始种群量子坍塌获得一组状态其中，为第t代种群中第j个个体的测量值，其表现形式为长度是m的二进制串，测量过程为随机产生一个0到1之间的数，若小于概率幅的平方，则测量结果为0，反之为1；步骤(C)，将当前状态不同测量值分别载入由α稳定噪声驱动Duffing双稳随机共振系统计算输出信噪比SNR，比较不同测量值对应的信噪比，将信噪比最大的个体保留，作为初始目标向量；其中，所述由α稳定噪声驱动Duffing双稳随机共振系统计算输出信噪比SNR，如公式(3)所示：Y=SNR=10lgS(f0)N(f0)---(3)]]>式(3)中，S(f0)为信号功率谱在频率f0处的幅值，N(f0)为同频背景噪声平均功率；步骤(D)，对Q(t0)实施一次测量并对其个体适应度评价，执行构造量子间联合交叉的量子交叉操作以及更改量子比特态叠加状态的量子变异操作，进而产生下一代个体Xinest；步骤(E)，对个体Xinest中各分量执行混沌搜索操作得到个体Xibest，计算其适应度值，如果Xibest位置优于Xinest，则替换Xinest，记录当前信噪比Yibest及其个体状态对应个体Xibest的a,b,k；其中，混沌搜索是将混沌状态引入到优化变量中，混沌搜索中使用Logistic映射，其迭代公式如下：zi+1＝μzi(1‑zi),i＝0,1,2,… (4)式(4)中，μ∈[0,4]为控制参数，初值z0∈[0,1]；步骤(F)，若遗传代数t≤T则跳到步骤(D)，否则；将最大适应度值Ybest及其个体向量Xbest作为最优解输出，个体向量Xbest即信噪比最大时a,b,k的值；步骤(G)，将个体向量Xbest包含的3个优化参数a,优化参数b,优化参数k载入由α稳定噪声驱动的Duffing双稳随机共振系统，产生共振输出；其中，x表示粒子位置，原有Duffing随机共振系统模型为：x··+kx·-ax+bx3=e(t)---(5)]]>式(5)中k为阻尼比，‑ax+bx3为非线性恢复力，e(t)＝s(t)+n(t)表示信号和噪声的混合激励信号，由α稳定噪声驱动的Duffing双稳随机共振系统模型为：x··+kx·-ax+bx3=s(t)+xDξ(t)---(6)]]>式(6)中，k为阻尼比，‑ax+bx3为非线性恢复力，为x的二次导数、为x的一次导数，s(t)＝Asin(2πf0t)是驱动信号幅值为A频率为f0的待测信号，xDξ(t)为α稳定噪声，D代表α稳定噪声的强度放大系数；其双稳系统势函数为：U(x)=-c2x2+d4x4---(7)]]>式(7)中是一个含有两个极小值和一个极大值x3＝0的双稳系统，其势垒高度为△U＝a2/4b，c,d为稳态系统参数；当Duffing系统有信号输入时，系统存在临界值当信号幅值A<Ac时，系统输出将在某势阱附近做局域周期运动；当A>Ac时，非线性系统和激励信号达到协同作用，一部分噪声能量向信号能量转移，系统输出在式(7)双稳系统的两稳态间大范围跃迁；其中，α稳定噪声数据通常由α稳定分布序列产生，α稳定分布的特征函数表达为：φ(t)=exp{jut-γ|t|α[1+jBsign(t)ω(t,α)]}]]>ω(t,α)=2πlog|t|,α=1tan(πα2),α≠1---(8)]]>式(8)中，α∈(0,2]为特征指数，决定分布的脉冲特性和拖尾特性；γ为尺度参数；u为位置参数，表明分布的中心；B∈[...
【技术特征摘要】
1.基于混沌量子遗传算法的随机共振微弱信号检测方法，其特征在于：将经过混沌量子遗传算法运算得到的优化参数a,优化参数b,优化参数k代入到由α稳定噪声驱动的Duffing双稳随机共振系统中，进行信号检测；其具体检测方法包括以下几个步骤：步骤(A)，量子混沌初始化及参数设置，在待寻优参数可行解空间中混沌产生个体为n的初始种群根据背景噪声确定Duffing随机共振系统参数范围，设置种群规模W，最大进化代数T，初始化进化代数计数器t＝1；其中，一个基因代表量子比特进行量子遗传算法计算，如公式(1)所示：式(1)中，和β分别为‘0’状态|0>和‘1’状态|1>的概率幅，且满足归一化条件染色体概率幅如公式(2)所示：式(2)中，为第t代第j个染色体，m为染色体基因个数，g为每个基因量子比特数，表示‘0’态概率幅具有m个染色体基因g个量子比特数，βmg表示‘1’态概率幅具有m个染色体基因g个量子比特数；步骤(B)，对种群实施一次测量，由初始种群量子坍塌获得一组状态其中，为第t代种群中第j个个体的测量值，其表现形式为长度是m的二进制串，测量过程为随机产生一个0到1之间的数，若小于概率幅的平方，则测量结果为0，反之为1；步骤(C)，将当前状态不同测量值分别载入由α稳定噪声驱动Duffing双稳随机共振系统计算输出信噪比SNR，比较不同测量值对应的信噪比，将信噪比最大的个体保留，作为初始目标向量；其中，所述由α稳定噪声驱动Duffing双稳随机共振系统计算输出信噪比SNR，如公式(3)所示：Y=SNR=10lgS(f0)N(f0)---(3)]]>式(3)中，S(f0)为信号功率谱在频率f0处的幅值，N(f0)为同频背景噪声平均功率；步骤(D)，对Q(t0)实施一次测量并对其个体适应度评价，执行构造量子间联合交叉的量子交叉操作以及更改量子比特态叠加状态的量子变异操作，进而产生下一代个体Xinest；步骤(E)，对个体Xinest中各分量执行混沌搜索操作得到个体Xibest，计算其适应度值，如果Xibest位置优于Xinest，则替换Xinest，记录当前信噪比Yibest及其个体状态对应个体Xibest的a,b,k；其中，混沌搜...
【专利技术属性】
技术研发人员：行鸿彦，卢春霞，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。