船舶辐射噪声连续谱仿真过程中，噪声成型滤波器的设计

船舶辐射噪声时域信号主要包括三部分：调制谱，连续谱与离散线谱信号。调制谱与离散线谱比较简单，此处不再过多说明。
空化状态下，船舶连续谱噪声由于在其在高频部分呈现-6dB/oct下降的趋势，需要构造特殊频谱响应的FIR滤波器以保证此频谱特征。连续谱噪声如下图所示：

图 船舶辐射噪声连续谱

生成连续谱噪声的流程为：

关键在于FIR滤波器的构造，也就是所谓的噪声成型滤波器（自己取的名字233），其频谱形状与上图右部大致相同。
**

特定频响FIR设计：

**
根据傅里叶变换的性质，我们知道：时域实偶对称信号，其频谱响应为实偶信号。这里我们已知频谱信号，故想要得到时域实信号，需要将频谱信号补全为偶信号。操作步骤为：将已知频谱响应截取0-fs/2部分，然后反转对称，补全频响曲线至0-fs区间。

图 FIR滤波器完整频谱响应

然后将此FIR响应取傅里叶逆变换（iFFT），然后取实部即为我们所需要的时域FIR滤波器系数。得到的FIR系数可能阶数很高，此时可以根据需要，对FIR频谱进行抽样，以减少时域滤波器长度。

验证：
将高斯白噪声通过该FIR滤波器（卷积）后，将频谱曲线平滑之后，可以得到改噪声的连续谱曲线，可以看到效果还是很不错的：

图 仿真噪声连续谱

2020.05.27 修改内容

生成的FIR滤波器系数为两头大，中间小，呈现“凹状”，这样会导致连续谱FFT结果震荡较大。将FIR滤波器系数进行循环（首尾相接），取中间部分，变为两头小，中间大的”凸状“，会使连续谱FFT结果更加稳定。对比结果如下：