我有一个随机行为但具有非随机(非高斯?)偏度和峰度值分别为 ~0.9 和 ~7.4 的信号的时域数据。
信号的 FFT 表明它不是宽带随机的——它通常在 100 Hz 以下是平坦的,并且在 100 Hz 以上时以大约 -3 dB/oct 逐渐变小,直到 400 Hz 的奈奎斯特。
我想创建一个与我的参考信号在统计上等效(或大致等效)的合成信号,几乎相同
- 频率内容
- 有效值
- 偏度和
- 峰度。
此外,我已经推断出高达 10 kHz 的 PSD 幅度,我想将其添加到混音中并反转到时域。
过去,我通过使用例如 MATLAB 的randn()
函数来创建随机白噪声(即创建随机相位角),对白噪声应用 FFT,缩放其频率分量,仅满足 (1) 和 (2)匹配所需的 FFT 频谱,并执行 IFFT 以转换为时域。该rand()
函数分别接近 0.0 和 3.0 的偏度和峰度,因此我最初的想法是创建具有指定偏度和峰度的白噪声。
我已经使用 MATLAB 的pearsnd()
函数尝试了相同的过程,该函数允许指定偏度和峰度值,但上述过程似乎没有保留缩放时域信号的偏度和峰度。我相信IFFT假设一个高斯分布。有解决办法吗?
谢谢!