2

我想生成包含随机生成的粉红噪声的任意大小的 2D 图像。 Wikipedia 建议粉红噪声的 2D 泛化将具有以 1/f^2 下降的能量。我在 MATLAB File Exchange 上找到了一些计算一维粉红噪声向量的代码。但我不知道如何正确地将其推广到二维——我对 fft 不是很熟悉,当我计算 ifft 时,我在下面的天真尝试会产生复杂的向量。

function pink = pinkNoiseImage(nrow,ncol)

rnrow = 2.^(ceil(log2(nrow)));
rncol = 2.^(ceil(log2(ncol)));
r = randn(rnrow,rncol);
rf = fft(r);
rnup = rnrow/2+1;
cnup = rncol/2+1;
frf = kron(1./sqrt(1:cnup),1./sqrt(1:rnup)');

rf(1:rnup,1:cnup) = rf(1:rnup,1:cnup).*frf;
rf(rnup+1:rnrow,1:cnup) = real(frf(rnrow/2:-1:2,1:cnup))-1i*imag(frf(rnrow/2:-1:2,1:cnup));
rf(1:rnup,cnup+1:rncol) = real(frf(1:rnup,rncol/2:-1:2))-1i*imag(frf(1:rnup,rncol/2:-1:2));
rf(rnup+1:rnrow,cnup+1:rncol) = real(frf(rnrow/2:-1:2,rncol/2:-1:2))-1i*imag(frf(rnrow/2:-1:2,rncol/2:-1:2));

pink = ifft(rf);

如何生成包含粉红噪声的二维矩阵?

4

1 回答 1

3

首先,不要总是相信维基百科告诉你的。或者,仔细阅读,因为粉红噪声的定义对于 2D 不是一对一的。其次,您可以使用以下FEX 文件生成1/f^beta具有正态误差分布的空间噪声。在该文件的文档中阅读更多内容。

于 2013-10-14T21:45:25.240 回答