生成带对数正态分布值的蓝/紫噪声：PSD形状丢失问题及测试

问题背景 #

需要生成同时满足以下两个条件的随机信号：

• 幅值服从对数正态分布（或其他下有界、非零众数的长尾分布）
• 功率谱密度（PSD）呈现低频功率低、随频率升高递增（比如蓝噪声、紫噪声的特性）

现有方法存在的问题：

• 直接对蓝/紫噪声（正态分布）做指数变换得到对数正态信号后，PSD会丢失原有递增特性，趋近于白噪声
• 某改进方案虽能匹配PSD，但对数正态分布的众数偏离预期（本该在1附近，实际接近0）

核心解决思路 #

1. 修正对数变换参数约束，减少PSD失真 #

对数正态分布由正态变量 $X \sim \mathcal{N}(\mu, \sigma^2)$ 经 $Y = e^X$ 生成，其众数为 $e^{\mu - \sigma^2}$。要让众数落在1附近，需满足 $\mu = \sigma^2$。

直接对蓝/紫噪声做指数变换会破坏PSD的原因是：指数变换的非线性会过度放大高频成分，掩盖原有PSD的递增趋势。解决步骤：

• 先对目标PSD的正态噪声做标准化处理，确保其均值 $\mu$ 等于方差 $\sigma^2$
• 进行指数变换后，在频域对信号做小幅修正：将信号的傅里叶幅值强制匹配目标PSD的平方根，保留原相位，再转回到时域

2. 迭代优化：平衡分布与PSD的匹配度 #

采用循环调整策略：

1. 生成符合目标PSD的正态噪声，按 $\mu = \sigma^2$ 调整参数后做指数变换，得到初始对数正态信号
2. 用Welch法计算当前信号的PSD，与目标PSD对比误差
3. 在频域修正：替换傅里叶幅值为目标PSD的平方根，保留原相位
4. 逆傅里叶变换得到时域信号，调整幅值使其符合对数正态分布（比如通过分位数映射）
5. 重复步骤2-4，直到PSD误差和分布偏差都在可接受范围内

3. 替代分布方案：换用其他长尾分布 #

如果对数正态的约束难以平衡，可尝试伽马分布或韦伯分布：

• 这类分布天生下有界且有非零众数，调整参数就能让众数落在1附近
• 生成目标PSD的正态噪声后，用分位数映射将正态分布转换为目标长尾分布——单调变换对PSD的破坏更小，能更好保留频谱特性

验证方法 #

• 幅值分布：用直方图、核密度估计验证是否符合目标分布，检查众数是否在1附近
• PSD特性：用Welch法计算PSD，确认低频到高频的递增趋势是否保留

内容的提问来源于stack exchange，提问作者Paul Brodersen