双稳态随机共振系统信噪比增益计算工具

资 源 简 介

本项目主要用于计算双稳态随机共振系统的输入与输出信噪比（SNR）增益，旨在量化分析非线性系统中随机共振现象的显著性。随机共振是一种在特定非线性系统、弱周期信号与背景噪声共同作用下，通过适度增加噪声强度反而使系统输出信噪比得到提升的奇特物理现象。

详 情 说 明

双稳态随机共振系统输入输出信噪比增益计算项目 README

项目介绍

本项目是一个基于 MATLAB 开发的数值仿真工具，专门用于研究和量化双稳态系统中发生的随机共振（Stochastic Resonance, SR）现象。随机共振是一种反直觉的物理过程，在非线性系统、微弱周期信号与背景噪声的共同作用下，通过调整或增加噪声强度，反而能够增强系统对信号的响应能力，提高输出端的信噪比。

本项目通过模拟双稳态势能模型，计算不同噪声强度下的系统响应，并从频域角度分析输入与输出信号的功率谱，最终给出信噪比增益曲线，为非线性信号处理和弱信号检测研究提供定量分析依据。

功能特性

1. 自动化的噪声强度扫描：程序能够在预设的噪声强度范围内（0.05 至 1.2）进行步进式仿真，自动记录每个强度下的系统表现。
2. 高精度数值求解：采用四阶龙格库塔算法（RK4）迭代求解随机二阶微分方程，确保非线性系统状态转移计算的准确性。
3. 全面的信噪比量化计算：内置输入/输出功率谱密度（PSD）分析模块，能够精准提取特定频率下的信号功率与邻域噪声功率。
4. 多维度可视化展示：提供时域波形对、功率谱对比图以及随机共振特性曲线（SNR 与 Gain 随噪声强度的变化趋势）。

实现逻辑与功能详细说明

本程序的运行流程严格遵循以下几个逻辑阶段：

1. 参数初始化

在程序开始阶段，定义了双向稳态系统的核心物理参数。包括：

• 势能参数 a 和 b：用于构建系统的双稳态势能函数。
• 信号参数：设置待检测周期信号的幅值 A 为 0.15，频率 f0 为 0.01Hz。
• 采样参数：设置采样频率为 5Hz，总仿真时长为 2000 秒，以保证足够的频率分辨率。
• 噪声范围：设定噪声强度 D 从 0.05 递增至 1.2，步长为 0.05。

2. 信号合成与系统仿真

针对每一个噪声强度 D，程序执行以下操作：

• 信号生成：生成纯净的余弦信号，并叠加服从标准正态分布的高斯白噪声。
• 状态演化：使用四阶龙格库塔法求解非线性方程 dx/dt = a*x - b*x^3 + s(t) + n(t)。该算法通过在每个时间步长内计算四个预测斜率（k1-k4），极大地提高了系统运动轨迹在双稳态势垒间跳变模拟的稳定性。

3. 频域分析与评价指标

• 功率谱计算：对原始输入信号和系统输出响应进行快速傅里叶变换（FFT）。在计算过程中，程序会自动去除直流分量以消除零频干扰，并计算单边功率谱密度。
• 信噪比（SNR）提取：在目标频率 f0 附近寻找峰值作为信号功率。同时，选取频率峰值邻域（排除中心点及其紧邻点，取附近的 10 个频率点）的平均功率作为噪声基底。
• 增益计算：通过输出信噪比减去输入信噪比（均为 dB 单位），得到系统的信噪比增益值。

4. 结果展示

程序最后生成三组图表：

• 时域分析：对比原始信号、含噪输入信号与经过双稳态系统处理后的输出波形，直观展示系统在特定噪声强度下的协同作用。
• 频域分析：展示输入与输出信号的功率谱，体现信号能量在 f0 处的集中度变化。
• 特性曲线：描绘输入 SNR、输出 SNR 随噪声强度变化的轨迹，以及 Gain 曲线。通常在特定噪声强度下，Gain 曲线会出现峰值，这正是随机共振现象发生的有力判据。

关键函数与算法细节

四阶龙格库塔算法 (RK4)

在每个时间步长 dt 内，程序通过计算四个分段斜率来更新位移 x 的值。这种方法比简单的欧拉法更适合处理具有陡峭势垒跳变的双稳态方程，能够捕捉到噪声诱导的系统跨势垒跃迁过程。

功率谱密度 (PSD) 计算逻辑

采用基于 FFT 的周期图法。核心细节包括：

• 使用 nextpow2 优化 FFT 计算长度。
• 对信号进行去均值处理，消除直流项对低频信号检测的干扰。
• 将双边谱转换为单边谱，并根据能量守恒原则对非直流分量乘以 2。

信噪比计算逻辑

信噪比的计算并非简单的全局信噪比，而是针对特定频率 f0 的窄带信噪比。程序通过：

• 定位：在频率轴寻找最接近 f0 的采样点。
• 隔离：将该点功率定义为 Ps。
• 采样：提取目标点周围一定范围内的功率值并求均值得到 Pn。
• 转换：使用 10 * log10(Ps / Pn) 公式转换为分贝。

使用方法

1. 启动 MATLAB 环境。
2. 将包含程序代码的项目文件夹设置为当前工作路径。
3. 在命令行窗口直接运行计算脚本或点击编辑器中的运行按钮。
4. 程序将自动进行循环迭代计算，并在完成后弹出三个结果分析窗口。
5. 观测“随机共振特性曲线”图，寻找增益曲线的最高点，该点对应的噪声强度即为该系统参数下的最佳共振噪声强度。

系统要求

• 软件环境：MATLAB R2016b 及以上版本（涉及绘图函数兼容性）。
• 硬件要求：由于存在噪声循环计算和高频采样模拟，建议配备 8GB 以上内存以确保大数据量 FFT 运算的流畅性。
• 依赖工具箱：主要使用 MATLAB 基础功能，无需特殊第三方库。