微弱信号检测方法课件

微弱信号检测方法第10章肖海林 微 弱 信 号 检 测 方 法 第10章肖海林随机共振检测方法混沌检测方法粒子滤波检测方法压缩感知检测方法本章内容第第1010章章 微弱信号检测方法微弱信号检测方法 随机共振检测方法 本章内容 第10章 微弱信号检测方法10.2 随机共振检测方法10.2 随机共振检测方法10.2.1 随机共振系统随机共振系统随机共振系统稳态随机共振系统稳态随机共振系统单稳态模型单稳态模型双稳态模型双稳态模型耦合双隐态模型耦合双隐态模型双隐态阵列模型双隐态阵列模型多隐态模型多隐态模型阈值随机共振系统阈值随机共振系统超阈值模型超阈值模型亚阈值模型亚阈值模型阈值阵列模型阈值阵列模型混沌随机共振系统混沌随机共振系统Chuas电路模型电路模型Duffing阵子模型阵子模型神经随机共振系统神经随机共振系统HHHH神经模型神经模型FNFN神经模型神经模型人工神经模型人工神经模型整合放电模型整合放电模型广义随机共振系统（广义随机共振系统（Noise Enhanced System）10.2.1 随机共振系统随机共振系统稳态随机共振系统单10.2.1 随机共振系统n 双稳态系统模型双稳态系统模型周期变化周期变化n 阈值系统模型阈值系统模型n混沌系统模型混沌系统模型郎之万方程郎之万方程10.2.1 随机共振系统 双稳态系统模型周期变化 阈值10.2.1 随机共振系统n 神经系统模型神经系统模型n神经系统模型神经系统模型广义的随机共振系统原理广义的随机共振系统原理10.2.1 随机共振系统 神经系统模型 神经系统模型广10.2.1 随机共振度量方法n 信噪比和处理增益信噪比和处理增益n 互相关函数互相关函数10.2.1 随机共振度量方法 信噪比和处理增益 互相关10.2.1 随机共振度量方法n 互信息互信息n 误码率和信道容量误码率和信道容量10.2.1 随机共振度量方法 互信息 误码率和信道容量n 双稳态随机共振系统原理双稳态随机共振系统原理系统势函数系统势函数10.2.2 双稳态随机共振系统郎之万方程郎之万方程逃逸速率逃逸速率双稳态系统的临界值双稳态系统的临界值 设计最佳双稳态随机共振的关键就是寻找合适的系统参数，设计最佳双稳态随机共振的关键就是寻找合适的系统参数，使得输出信噪比达到最佳值也就是该极大值点使得输出信噪比达到最佳值也就是该极大值点 双稳态随机共振系统原理系统势函数10.2.2 双稳态随n 基于双稳态随机共振的能量检测流程基于双稳态随机共振的能量检测流程10.2.3 基于双稳态随机共振系统能量检测算法 基于双稳态随机共振的能量检测流程10.2.3 基于双稳10.2.4 广义随机共振系统广义的随机共振原理广义的随机共振原理基于广义随机共振系统的能量检测算法流程基于广义随机共振系统的能量检测算法流程 相比于传统的能量检测算法，基于广义随机共振的能量检测相比于传统的能量检测算法，基于广义随机共振的能量检测还需要添加随机共振噪声这一步骤还需要添加随机共振噪声这一步骤10.2.4 广义随机共振系统广义的随机共振原理基于广义10.3 混沌检测方法10.3 混沌检测方法 n 混沌的定义混沌的定义10.3.1 非线性动力学系统中的混沌闭区间闭区间I I上的上的上的连续自映射上的连续自映射 若满足下列条件若满足下列条件：P P的的周周期期点点的的周周期期无无上上界界；P P具具有有任任意意正正整整数数周周期期点点，即即对对于于任任意意，有有 ，使，使 （非不动点的（非不动点的n n周期点）周期点）闭区间闭区间I I上存在不可数子集上存在不可数子集S S，满足，满足对对 ，当，当 时，有时，有 ；对对 ，有，有 ；对对 和和 P P 的任一周期点的任一周期点y,y,有有 ；则称则称P P是混沌的。是混沌的。一个非线性动力学系统是否会产生混沌要取决于系统参数的选择。一个非线性动力学系统是否会产生混沌要取决于系统参数的选择。混沌的定义10.3.1 非线性动力学系统中的混沌闭（1 1）自时域直接观察法）自时域直接观察法（2 2）自）自Poincare截面法截面法（3 3）频谱分析法）频谱分析法（4 4）Lyapunov指数法指数法（5 5）Kolmogorov 法法 （6 6）Melnikov 法法（7 7）Shilnikov 法法（8 8）分型理论分析）分型理论分析（9 9）分维数计算法）分维数计算法10.3.2 混沌运动的分析方法（1）自时域直接观察法10.3.2 混沌运动的分析方法10.3.3 Duffing振子的运动特性研究无阻尼、无驱动力的无阻尼、无驱动力的DuffingDuffing方程相平面轨迹图方程相平面轨迹图一对内轨道的参数方程：一对内轨道的参数方程：k k趋向于趋向于1 1时，内轨道趋向于同宿轨道时，内轨道趋向于同宿轨道同轨道参数方程：同轨道参数方程：10.3.3 Duffing振子的运动特性研究无阻尼、无10.3.3 Duffing振子的运动特性研究外轨的时间参数方程：外轨的时间参数方程：k k趋向于趋向于1 1时，轨道趋向于同宿轨道时，轨道趋向于同宿轨道Hamilton Hamilton 函数：函数：10.3.3 Duffing振子的运动特性研究外轨的时间10.3.3 Duffing振子的运动特性研究Duffing产生混的条件DuffingDuffing振子周期轨道存在的条件振子周期轨道存在的条件 当当 ，不存在周期闭轨，而当，不存在周期闭轨，而当 时，存在一条稳定的周期时，存在一条稳定的周期闭轨和一条不稳定的周期闭轨。闭轨和一条不稳定的周期闭轨。10.3.3 Duffing振子的运动特性研究Duffi10.3.4 参数对混沌振子运动的影响混沌质子的参数敏感性混沌质子的参数敏感性参数的敏感性与对初值的敏感性是等价的参数的敏感性与对初值的敏感性是等价的稳定周期轨道的选择稳定周期轨道的选择 设设 选择合适的选择合适的 和和 的临界值，的临界值，系统才能够混沌向周期运系统才能够混沌向周期运动状态转变的分岔值动状态转变的分岔值10.3.4 参数对混沌振子运动的影响混沌质子的参数敏感10.3.4 参数对混沌振子运动的影响混沌振子的噪声统计特性混沌振子的噪声统计特性系统周期轨道上的噪声统计特性输输出出噪噪声声在在一一个个周周期期内内的的各各轨轨迹迹点点上上统统计计特特性性都都不不同同，为为非非平平稳稳过过程程。但但在在若若干干周周期期之之间间，呈呈现现为为循循环环平平稳稳过过程程。可可见见输出噪声均方值最小点在输出噪声均方值最小点在A A或或E E附近附近10.3.4 参数对混沌振子运动的影响混沌振子的噪声统计10.4 粒子滤波检测法10.4 粒子滤波检测法 n 基于粒子滤波基于粒子滤波算法算法的存在概率检测算法的存在概率检测算法。u 连续的目标状态和离散的目标存在状态连续的目标状态和离散的目标存在状态。一维一维连续状态变量根据状态空间模型预测和更新连续状态变量根据状态空间模型预测和更新 存在变量依据转移概率公式实现预测和更新存在变量依据转移概率公式实现预测和更新n 基于粒子滤波算法的似然比检测方法。基于粒子滤波算法的似然比检测方法。u 当近似似然比大于给定门限时，认为目标存在当近似似然比大于给定门限时，认为目标存在u 当近似似然比小于给定门限时，认为目标不存在当近似似然比小于给定门限时，认为目标不存在10.4.1 粒子滤波背景知识 基于粒子滤波算法的存在概率检测算法。连续的目标状态和 n 状态空间模型状态空间模型 系统方程系统方程 观测方程观测方程n 后验概率密度函数后验概率密度函数u根据根据观测观测方程方程及其概率形式及其概率形式10.4.2 状态空间模型和后验概率密度函数u根据根据观测观测方程方程及其概率形式及其概率形式 状态空间模型 系统方程 后验概率密度函数根据观测方程10.4.3 卡尔曼滤波和扩展卡尔曼滤波卡尔曼滤波卡尔曼滤波u后验概率的递推过程后验概率的递推过程卡尔曼滤波器框图卡尔曼滤波器框图10.4.3 卡尔曼滤波和扩展卡尔曼滤波卡尔曼滤波后验概10.4.3 卡尔曼滤波和扩展卡尔曼滤波其他最优滤波器其他最优滤波器uK K时刻更新的方差时刻更新的方差uK K时刻更新的均值时刻更新的均值uK K时刻预测后验概率密度函数时刻预测后验概率密度函数10.4.3 卡尔曼滤波和扩展卡尔曼滤波其他最优滤波器K10.4.3 卡尔曼滤波和扩展卡尔曼滤波uK K时刻后验概率密度函数时刻后验概率密度函数扩展卡尔曼滤波器扩展卡尔曼滤波器扩展卡尔曼滤波器框图扩展卡尔曼滤波器框图10.4.3 卡尔曼滤波和扩展卡尔曼滤波K时刻后验概率密10.4.3 卡尔曼滤波和扩展卡尔曼滤波u预测均值与协方差预测均值与协方差u更新均值与协方差更新均值与协方差其他次优滤波器其他次优滤波器u预测后验概率密度函数预测后验概率密度函数10.4.3 卡尔曼滤波和扩展卡尔曼滤波预测均值与协方差10.4.4 粒子滤波算法蒙特罗积分蒙特罗积分u预测后验概率密度函数预测后验概率密度函数 蒙特卡罗近似的收敛速度只样本维数蒙特卡罗近似的收敛速度只样本维数N N有关，与积分位数有关，与积分位数 无关无关重要性采样重要性采样10.4.4 粒子滤波算法蒙特罗积分预测后验概率密度函数10.4.4 粒子滤波算法序贯重要性采样序贯重要性采样u后验概率密度函数后验概率密度函数 10.4.4 粒子滤波算法序贯重要性采样后验概率密度函数10.4.5 各种粒子滤波算法序贯重要性采样粒子滤波算法序贯重要性采样粒子滤波算法uSIRSIR算法中的权值：算法中的权值：SIRSIR算法收敛较慢。另外，算法收敛较慢。另外，SIRSIR算法在每个观测时刻都执行重采样过算法在每个观测时刻都执行重采样过程，必然会导致样本大量重复，损失样本多样性。程，必然会导致样本大量重复，损失样本多样性。10.4.5 各种粒子滤波算法序贯重要性采样粒子滤波算法10.4.5 各种粒子滤波算法其他其他粒子滤波算法粒子滤波算法 规划粒子滤波算法规划粒子滤波算法基于基于MCMCMCMC步进步进与与SIRSIR算算法法相相比比，APFAPF算算法法利利用用k k时时刻刻的的观观测测信信息息从从k-1k-1时时刻刻的的更更新新样样本本中中直直接接采采样样，因因此此可可更更接接近近真真实实的的目目标标状状态态。当当过过程程噪噪声声较较小小时时，APFAPF算算法法的的性性能能优优于于SIRSIR。当当过过程程噪噪声声较较大大时时，此时此时APFAPF性能较性能较SIRSIR差差辅助粒子滤波算法辅助粒子滤波算法10.4.5 各种粒子滤波算法其他粒子滤波算法 规划粒10.4.6 代价参考粒子滤波算法代价参考代价参考粒子滤波算法粒子滤波算法 初始化初始化重采样重采样样本更新及其代价样本更新及其代价状态估计状态估计10.4.6 代价参考粒子滤波算法代价参考粒子滤波算法 10.4.6 代价参考粒子滤波算法 前向前向-后向代价参考粒子滤波算法后向代价参考粒子滤波算法步骤一：步骤一：代价参考粒子滤波算法估计目标在各个时刻的状态，得代价参考粒子滤波算法估计目标在各个时刻的状态，得到各个时刻的粒子到各个时刻的粒子-代价集合代价集合 步骤二：步骤二：然后以然后以K K的粒子的粒子-代价集合代价集合 ，作为初始时刻作为初始时刻的样本的样本-代价集合代价集合，估计，估计 ，并以此估计目标在各并以此估计目标在各个时刻的状态个时刻的状态。10.4.6 代价参考粒子滤波算法 前向-后向代价参考粒10.5 压缩感知检测方法10.5 压缩感知检测方法u定义定义10.5.2 压缩感知理论的基本框架压缩感知理论框架压缩感知理论框架信号的稀疏变换信号的稀疏变换 信号信号X X在正交基在正交基 下的变换系数下的变换系数 假如对于假如对于0P20P0,R0,这些系数满足这些系数满足:则则说明在某种意义下系数向量说明在某种意义下系数向量 是稀疏的是稀疏的。压缩感知问题描述压缩感知问题描述定义10.5.2 压缩感知理论的基本框架压缩感知理论框架10.5.3 压缩感知的核心问题压缩感知的稀疏字典设计压缩感知的稀疏字典设计K-SVDK-SVD算法实现的迭代步骤为算法实现的迭代步骤为：（1 1）求）求x x的稀疏编码的稀疏编码（MPMP或或OMPOMP算法）算法）：（2 2）更新字典更新字典D D。过完备字典能够由多个正交基组成过完备字典能够由多个正交基组成，且过完备字典包含的原子最大个数且过完备字典包含的原子最大个数为为 ，并且其相干系数仍为并且其相干系数仍为单一正交基单一正交基10.5.3 压缩感知的核心问题压缩感知的稀疏字典设计K10.5.3 压缩感知的核心问题压缩感知的测量矩阵设计压缩感知的测量矩阵设计限制等距性质限制等距性质对于任意给定的整数对于任意给定的整数 ，矩阵矩阵 的限制等距常数的限制等距常数 对于任意对于任意K K稀疏向量稀疏向量 均满足均满足则称矩阵则称矩阵 满足满足K K阶阶RIPRIP性质性质。10.5.3 压缩感知的核心问题压缩感知的测量矩阵设计限10.5.3 压缩感知的核心问题压缩感知的重构算法设计压缩感知的重构算法设计步骤步骤一：一：基于采样值基于采样值 ，对测量矩阵对测量矩阵 ，求得估计值求得估计值 ；步骤二：步骤二：利用利用 求重构信号求重构信号稀稀疏疏度度为为K K的的x x信信号号在在变变换换基基下下的的系系数数向向量量为为 ,对对于于由由感感知知 矩矩阵阵 获得的获得的M M个采样个采样数据，如果数据，如果 ，则，则 优优化化问问题题将将以以极极大大概概率率精确恢复精确恢复 ，式中式中 是一个与恢复精度相关的常数。是一个与恢复精度相关的常数。定理定理10-5-310-5-310.5.3 压缩感知的核心问题压缩感知的重构算法设计步10.5.3 压缩感知的核心问题n 重构算法的两大类：重构算法的两大类：以匹配追踪为代表的贪婪类算法以匹配追踪为代表的贪婪类算法 以基追踪（以基追踪（Basic Pursuit BP）算法为代表的凸优化类算法）算法为代表的凸优化类算法 范数优化图范数优化图10.5.3 压缩感知的核心问题 重构算法的两大类：以 谢谢！谢谢！