LOFAR & DEMON

LOFAR (Low frequency analysis and recording)谱可反映信号的非平稳特性，进而可提取信号中的宽带线谱分布特征。但是轴频及其倍频却因为基本上都淹没在低频宽带噪声中而无法直接获取。

而舰船宽带噪声高频段存在调制现象，DEMON (Detection of Envelope Modulation On Noise)分析通过对接收的宽带信号进行解调以获得低频的包络谱，从而获得了诸如目标轴频、叶频等低频段较强的物理特征。

可以说，由LOFAR分析可获得反应舰船目标各部件结构特征的宽频带特征，而DEMON分析则获得较低频段的强制线谱特征，弥补了LOFAR分析在低频段的不足，共同描绘了目标信号的谱系特征。

LOFAR谱图分析

LOFAR谱图分析是较具代表性的被动声纳信号处理方法之一。该方法通过对采样数据作短时傅立叶变换而构成信号表达的三维图（声谱图）。

当噪声信号随时间发生比较明显的变化，即信号出现非平稳特性时，建立在平稳性假设基础上的传统的傅立叶表示方法已不适合应用。而信号的LOFAR谱图从时间和频率两个角度对信号进行联合域分析，适宜于水声目标辐射噪声的局部平稳特性。

LOFAR谱图的具体实现步骤如下：
（1）将原始信号的采样序列分成连续的若干段，每段N个采样点。根据具体情况，段间可有部分重叠。
（2）对每段信号采样样本 $L(n)$ 作归一化和中心化处理，归一化处理的目的是使接收信号的幅度（或方差）在时间上均匀；中心化是为了使样本的均值为零。

u (n) = \frac{L (n)}{m a x_{1 \leq i \leq N} [L (i)]}

$u(n)=\frac{L(n)}{max_{1\leq i\leq N}[L(i)]}$

x (n) = u (n) - \frac{1}{N} \sum_{i = 1}^{N} u (i)

$x(n)=u(n)-\frac{1}{N}\sum_{i=1}^{N}u(i)$

(3)对信号 $x(n)$ 作短时傅立叶变换得到LOFAR谱图。

DEMON谱图分析

DEMON谱是从接收信号的高频调制谱中解调出的线谱信息，是一种常用的获取轴频与倍频的方法。

在舰船辐射噪声中，螺旋桨节拍对舰船的辐射噪声存在着明显的振幅调制，调制频率等于螺旋桨转速或螺旋桨叶片频率（轴频乘以叶片数目）。因此舰船辐射噪声具有鲜明的节奏感。这种幅度调制信号即包络信号是一种慢变化信号，它反映了能量随时间的变化，由周期性成分和非周期性成分组成。通常采用宽带噪声包络调制分析技术测量其基频和各次谐波，提取螺旋桨的轴频和叶片频，这种方法称为DEMON。

在DEMON分析中，一般选取调制最强的时变谱频段进行包络检测，然后进行频谱分析得到DEMON谱。根据信号调制的原理，舰船噪声以周期性局部平稳过程为模型可表示为：

s (t) = [1 + f (t)] n (t)

$s(t)=[1+f(t)]n(t)$

$n(t)$ 是宽带平稳白色高斯随机过程；
$f(t)$ 是调制函数，它是慢变化的周期函数，它所在的频域比 $n(t)$ 所在频域的频率要低得多。
利用平方解调方法即可获得调制信号 $f(t)$ 。
再对获得的解调信号进行FFT处理即可得到解调信号的DEMON谱。
这可以称为基本的DEMON分析过程，流程图：-->带通滤波-->平方检波-->低通滤波-->FFT

仅仅做基本的DEMON分析是不够的，不能获得理想的DEMON线谱。舰船噪声调制信号对不同频带噪声信号的调制程度不一样，对某些频带调制强，对其它频带则相对调制弱些。另外，各个频带的DEMON谱线也不可能完全相同，不同频段存在不同的谱线缺失现象，同时还可能存在“虚假谱线”。采用分频段处理结果的方法可以弥补单一频段解调DEMON谱处理的不足。
由于所获得的DEMON谱中存在连续谱而使线谱成分不明显，还存在谱线靠的太近而无法分辨的情况，使DEMON谱看起来比较“模糊”。这就要求将分频段处理再融合获得的DEMON谱中的连续谱去除，并进一步进行线谱净化处理，使谱线更加清晰。

应用 $\alpha$ 双向滤波器获得DEMON谱的趋势项，再以趋势项适当的倍数为门限，保留原始DEMON谱中过此门限的谱线即可做到线谱与连续谱的分离，获得DEMON谱线。

$\alpha$ 双向滤波器是一阶递归滤波器，设某采样信号 $s(n)$ 长度为 $N$ ，当 $n=1,2,...,N-1$ 、 $m=1,2,...,N$ 时， $\alpha$ 双向滤波器的基本原理可表示为

s_{1} (1) = s (1)

$s_{1}(1)=s(1)$

s_{1} (n + 1) = s_{1} (n) + [s (n + 1) - s_{1} (n)] / Q

$s_{1}(n+1)=s_{1}(n)+[s(n+1)-s_{1}(n)]/Q$

s_{2} (n) = s_{2} (n + 1) + [s (n) - s_{1} (n + 1)] / Q

$s_{2}(n)=s_{2}(n+1)+[s(n)-s_{1}(n+1)]/Q$

s_{2} (N) = s (N)

$s_{2}(N)=s(N)$

s_{m} (k) = \frac{s_{1} (k) + s_{2} (k)}{2}

$s_{m}(k)=\frac{s_{1}(k)+s_{2}(k)}{2}$

$s_{m}(k)$ 为经过 $\alpha$ 双向滤波器后获得的信号， $Q$ 为递归系数，递归系数的值越大，跟踪能力越好，滤波效果越差。当递归系数取2的整数幂的形式时，用移位就可以实现整个滤波，并且不用担心溢出问题。

综上所述，将DEMON分析方法的具体实现步骤归纳整理如下：
（1）选取可能具有较强振幅调制的目标噪声信号的频带，一般选择大于1000Hz某频带范围（例如3000-5000Hz）的信号，将此频带内的信号按照一定的规律（平均或按比例或按频程）分为N个频带进行时域滤波，对每个频带的信号分别进行平方检波以得到其包络信号，对包络信号进行FFT运算获得到该频带的DEMON谱，假设为 $P_{i},i=1,2,...,N$ 。
（2）对每个频带的解调信号的DEMON谱 $P_{i}$ 进行加权处理。将第 $i(1\leq i\leq N)$ 个频段的DEMON谱 $P_{i}$ 经过 $\alpha$ 双向滤波器，适当选取递归系数 $Q$ 的值，可得到趋势项即连续谱 $\bar{P_{i}}$ 。
（3）将以上获得的N个归一化处理的DEMON线谱和与其相应的加权系数相乘，再进行求和就得到了较完善的DEMON谱 $P$ 。
（4）DEMON线谱 $P$ 的净化。经过 $\alpha$ 双向滤波器，适当选取递归系数Q的值，得到趋势项即连续谱 $\bar{P}$ 。...
[1]朱世才. 目标通过特性的LOFAR及DEMON分析[D].哈尔滨工程大学,2011.