基于双稳态随机共振的微弱信号检测系统

资 源 简 介

该项目构建了一个基于非线性随机共振理论的微弱信号检测仿真平台，旨在强噪声背景下恢复和识别微弱的周期性特征信号。系统核心利用双稳态势阱模型，通过协同利用噪声能量实现信号能量的跃迁增强，从而显著提高输出信噪比。具体实现包含两个关键部分：1. 数值求解器（sr）：使用四阶龙格-库塔（Runge-Kutta）算法对含噪非线性朗之万方程进行离散化迭代求解，模拟粒子在双稳态势阱中的运动轨迹，这是实现随机共振效应的数学核心；2. 主控流程（rs）：作为系统的入口程序，负责生成测试数据（包括微弱正弦信号与高斯白噪声的叠加）、配置系统结构参数（如势垒高度参数a和b）、调用sr求解器进行非线性变换，以及执行后续的信号分析。通过对输入混合信号和经过随机共振系统处理后的输出信号进行时域对比及FFT频谱分析，项目能够清晰演示在特定噪声强度匹配下，微弱信号在频谱上的峰值显著增强现象，从而验证在低信噪比条件下检测弱信号的可行性。

详 情 说 明

基于随机共振的弱信号检测系统

项目简介

本项目构建了一个基于非线性随机共振（Stochastic Resonance, SR）理论的微弱信号检测仿真平台。系统利用双稳态势阱模型，在强噪声背景下通过非线性系统的协同作用，将噪声能量转移至信号能量，从而实现微弱周期性信号的增强与检测。

该仿真平台使用 MATLAB 开发，集成了信号生成、四阶龙格-库塔（Runge-Kutta）数值迭代求解、频谱分析及可视化功能，能够直观地演示在特定噪声强度下，微弱信号在频谱上的峰值显著增强现象，验证了低信噪比条件下检测弱信号的可行性。

功能特性

• 双稳态系统建模：核心采用朗之万方程描述的过阻尼双稳态系统，通过调节系统参数（线性项 $a$ 和非线性项 $b$）控制势阱结构。
• 高精度数值求解：内置四阶龙格-库塔（RK4）算法，对含噪非线性微分方程进行离散化迭代求解，确保仿真精度。
• 参数化仿真配置：支持自定义采样频率、仿真时长、信号频率/幅值、噪声强度 $D$ 以及系统结构能垒参数。
• 信号与噪声生成：能够生成微弱正弦信号与高斯白噪声的混合输入，模拟真实的低信噪比环境。
• 频域分析与评估：自动计算输入与输出信号的功率谱密度（PSD），并针对特征频率计算谱增益（dB），量化随机共振效果。
• 多维度可视化：提供时域波形对比、SR系统响应波形以及输入输出频谱对比图，清晰展示信号增强过程。

核心算法与实现细节

本项目的主要逻辑封装在单一入口函数及其辅助函数中，具体的代码实现逻辑如下：

1. 系统参数配置与初始化

程序首先进行环境清理，随后定义了仿真所需的关键参数：

• 仿真参数：设定采样频率为 50Hz（适用于SR处理低频信号的特性），总时长 1000秒，步长 $h$ 为采样周期的倒数。
• 信号参数：设定微弱信号幅值为 0.3，频率为 0.05Hz。该频率满足绝热近似条件。
• 双稳态参数：$a=1, b=1$，构建标准的双稳态势场 $V(x) = -a/2 x^2 + b/4 x^4$。
• 噪声参数：定义噪声强度 $D=0.4$，用于生成具有特定功率的随机干扰。

2. 混合信号生成

• 生成标准的正弦波作为纯净的微弱信号（Amplitude=0.3）。
• 生成高斯白噪声序列，噪声振幅按照 $sqrt{2D}$ 进行缩放，以匹配随机共振方程中的噪声项定义。
• 将纯净信号与噪声线性叠加，形成信噪比较低的输入信号 s_input。

3. 随机共振数值求解流程

这是系统的核心处理环节，通过内部函数 sr_runge_kutta 实现。程序并在控制台输出求解进度提示。

• 输入：含噪信号向量、系统参数 $a, b$、步长 $h$、数据长度 $L$。
• 动力学方程：求解 $frac{dx}{dt} = ax - bx^3 + Input(t)$。
• RK4 迭代：采用经典的四阶龙格-库塔法，在每一步迭代中计算四个斜率分量 ($k_1, k_2, k_3, k_4$)，加权平均后更新系统状态 $x$。
• 离散化处理：将输入的含噪信号视为每个时间步的驱动力，模拟粒子在双稳态势阱中的运动轨迹，最终输出系统响应信号 x_output。

4. 信号分析与性能评估

系统对输入和输出信号分别进行处理：

• 调用 calculate_psd 函数计算单边幅值谱。
• 特征提取：自动搜索距离目标频率 $f_0$ 最近的频点。
• 增益计算：提取特征频率处的幅值，计算输出相对于输入的频谱增益（$Gain = 10 log_{10}(Amp_{out}/Amp_{in})$）。
• 结果输出：在命令行打印特征频率、输入/输出幅值及具体的增益分贝数。

5. 结果可视化

程序会生成一个包含三个子图的综合图表：

1. 输入时域图：展示含噪输入信号与原始纯净信号的对比（为避免由于数据量大导致的显示混叠，仅绘制前2000个数据点）。
2. SR输出时域图：展示经过双稳态系统处理后的输出波形，观察是否存在周期性跃迁。
3. 功率谱密度对比图：在同一坐标系下绘制输入与输出的低频段频谱（$0 sim 5f_0$），并用虚线标记特征频率位置，直观展示随机共振带来的频谱峰值提升。

关键辅助函数说明

sr_runge_kutta

该函数实现了随机共振的数学核心。

• 算法：标准四阶龙格-库塔（RK4）。
• 逻辑：通过循环遍历整个输入信号序列。在每一步中，基于当前状态 $x_n$ 和当前时刻的外部驱动力（信号+噪声），计算下一时刻的状态 $x_{n+1}$。
• 方程形式：明确使用了 $dx/dt = f(x) + signal$ 的形式，其中 $f(x) = ax - bx^3$ 为系统自身的恢复力。

calculate_psd

该函数用于频域分析。

• 算法：快速傅里叶变换（FFT）。
• 处理：计算FFT模值，除以数据长度进行归一化，并转换为单边谱（即乘以2，且仅取前半部分频率），确保物理意义上的幅值正确性。

系统要求

• MATLAB R2016b 或更高版本。
• 不需要额外的工具箱（Toolbox），代码仅使用 MATLAB 基础数学和绘图函数。

使用方法

1. 将代码保存为 .m 文件。
2. 在 MATLAB 环境中打开该文件。
3. 直接运行文件。
4. 观察命令行输出的增益数据，以及弹出的波形和频谱对比图窗。
5. 若需研究不同噪声强度下的共振效果，可修改代码头部的 D 参数并重新运行。