水声信号处理中DEMON谱分析的原理、实现与改进

一、DEMON谱分析的基本概念

DEMON谱（Detection of Envelope Modulation on Noise，噪声包络调制检测谱）是水声信号处理中提取舰船目标低频特征的关键技术，专门用于从舰船辐射噪声的高频调制信号中解调出轴频（螺旋桨转动频率）、叶频（轴频与桨叶数的乘积）等强制线谱特征。这些特征是舰船识别的重要依据（如螺旋桨转速、桨叶数等），弥补了LOFAR谱在低频段（<1000Hz）特征不明显的不足。

二、DEMON谱分析的核心原理

舰船辐射噪声中的螺旋桨节拍会对高频宽带噪声产生振幅调制，调制频率等于轴频或叶频。DEMON谱分析的本质是解调这种包络调制信号，具体流程如下：

带通滤波：选取调制最强的高频频带（如3000-5000Hz），过滤无关噪声；
平方检波：将高频信号转换为包络信号（平方运算增强调制特征）；
低通滤波：去除高频成分，保留慢变化的包络；
FFT变换：对包络信号进行频谱分析，得到DEMON谱（低频包络谱）。

三、DEMON谱分析的MATLAB实现

MATLAB是DEMON谱分析的主流工具，以舰船辐射噪声为例：

% 1. 信号加载与预处理
load('ship_noise.mat'); % 加载舰船辐射噪声信号
fs = 12000; % 采样频率
t = 0:1/fs:length(noise)-1/fs; % 时间轴

% 2. 带通滤波（选取调制最强频带）
band = [3000, 5000]; % 带通滤波区间
b = fir1(100, band/(fs/2), 'bandpass'); % 100阶FIR滤波器
filtered_noise = filter(b, 1, noise);

% 3. 平方检波（提取包络）
envelope = abs(filtered_noise).^2;

% 4. 低通滤波（保留慢变化包络）
cutoff = 100; % 低通截止频率（轴频范围）
d = fir1(50, cutoff/(fs/2), 'low');
lowpass_envelope = filter(d, 1, envelope);

% 5. FFT变换（生成DEMON谱）
N = length(lowpass_envelope);
f = (0:N-1)*(fs/N); % 频率轴
demod_spectrum = abs(fft(lowpass_envelope, N));
demod_spectrum = demod_spectrum(1:floor(N/2)+1); % 取单边谱
f = f(1:floor(N/2)+1);

% 6. 结果可视化
figure;
plot(f, 20*log10(demod_spectrum/max(demod_spectrum))); % 归一化对数谱
xlabel('频率（Hz）');
ylabel('相对幅度（dB）');
title('DEMON谱（轴频/叶频特征）');
grid on;

说明：

带通滤波选取3000-5000Hz（舰船辐射噪声调制最强的频带）；
平方检波将高频信号的振幅调制转换为包络信号；
低通滤波（截止频率100Hz）保留轴频（通常<100Hz）；
FFT变换生成DEMON谱，显示轴频（如5Hz）、叶频（如25Hz，桨叶数5）等特征。

四、DEMON谱分析的改进方法

传统DEMON谱分析存在连续谱干扰、虚假谱线、不稳定时刻线谱不足等问题，以下是常见的改进策略：

1. 背景均衡处理（排序截断平均法）

目的：去除DEMON谱中的连续谱，增强线谱对比度；
实现：对DEMON谱数据进行排序，去除极端值（如α=0.05，去除5%的最大值和最小值），计算截断平均值作为背景，用原始数据减去背景得到均衡后的谱。

2. 阵元域融合（大孔径声呐）

目的：降低声传播引起的包络畸变，提高DEMON谱估计质效；
实现：（1）对阵列各阵元的DEMON谱进行相关系数分析，选取相关系数高的阵元（如Cth=0.9）组成最大相关子阵；（2）计算子阵DEMON谱的统计估值（如均值），用连续谱作为检测门限，筛选信号成分。

3. 线谱净化（α双向滤波器）

目的：去除虚假谱线，保留真实的轴频/叶频；
实现：用α双向滤波器（一阶递归滤波器）提取DEMON谱的趋势项（连续谱），以趋势项的1.5-2倍为门限，保留超过门限的谱线（线谱）。

4. 改进的最大公约数法（轴频提取）

目的：精确提取轴频（线谱差频的最可能值）；
实现：

（1）对每组DEMON谱提取强度最高的I根线谱（如I=10）；

（2）计算线谱间的差频，统计各差频的权重（如差频fij的权重为wi+wj，wi为线谱fi的归一化强度）；

（3）用线频差频比β=fi/fij更新权重（β为整数时权重加1），选取权重最大的差频作为轴频。

参考代码水声信号处理 demon谱分析 www.youwenfan.com/contentcso/45020.html

五、DEMON谱分析的应用场景

DEMON谱分析主要用于舰船目标识别，具体应用包括：

轴频/叶频提取：通过DEMON谱中的低频线谱，计算螺旋桨转速（轴频）和桨叶数（叶频=轴频×桨叶数）；
目标分类：结合轴频、叶频等特征，区分不同类型舰船（如商船、军舰）；
状态监测：通过轴频变化，监测螺旋桨磨损或故障（如桨叶断裂会导致叶频消失）。

六、总结

DEMON谱分析是水声信号处理中提取舰船低频特征的核心技术，通过带通滤波、平方检波、低通滤波和FFT等步骤，从高频调制信号中解调出轴频、叶频等关键特征。MATLAB是实现DEMON谱分析的主流工具，结合背景均衡、阵元域融合、线谱净化等改进方法，可有效提高分析精度，为舰船目标识别提供可靠的特征依据。