how to generate gaussian noise with certain covariance and zero mean ?

162 次查看（过去 30 天）
Tejas Appaji
Tejas Appaji 2017-9-25
回答： getahun kibret 2023-2-22
i have a signal and i want to add gaussian noise to it with zero mean and 0.1 covariance.
variance = 0.1
W = sqrt(variance).*randn(1,size(Xmodt,2)); %Gaussian white noise W
mysignal = mysignal + W; %Add the noise
this code lets me define variance. but i need an algorithm or code to generate gaussian noise with specific covariance and zero mean.
thanks in advance

请先登录，再回答此问题。

采纳的回答

David Ding
David Ding 2017-9-27
Hi Tejas,
If I understand your question correctly, you wish to generate AWGN with certain co-variance. In this case, you would have a vector of zero-mean Gaussian noises that are statistically dependent. In order to model this in MATLAB, your workflow would be to generate an n x 1 noise vector and then pre-multiply that by the co-variance matrix.
For example:
% Generate a 2 x 1 Gaussian noise vector with covariance
noiseVec = randn(2, 1);
var1 = 0.1;
var2 = 0.2;
covar = 0.05;
cMatrix = [var1, covar; covar, var2];
noiseVec = cMatrix * noiseVec;
  4 个评论
显示 2更早的评论
Shah Mahdi Hasan
Shah Mahdi Hasan 2020-8-14
I agree with JUNHO. There should have been sqrt(). Think about the scalar analogy. If you have a standard normal random variable x~N(0,1) and want to have a certain variance sigma, then you would multiply the following:
y ~ N(0,sigma^2) = sigma*x
Brendan Nichols
Brendan Nichols 2020-11-26
I agree with Junho as well. Test out a variation of David's answer:
noiseVec = randn(2, 1e6);
var1 = 0.1;
var2 = 0.2;
covar = 0.05;
cMatrix = [var1, covar; covar, var2];
noiseVec = cMatrix * noiseVec;
cov(noiseVec(1, :), noiseVec(2, :))
Note the covariance is not equal to cMatrix. Try the following instead:
noiseVec = randn(2, 1e6);
var1 = 0.1;
var2 = 0.2;
covar = 0.05;
cMatrix = [var1, covar; covar, var2];
noiseVec = chol(cMatrix, 'lower') * noiseVec;
cov(noiseVec(1, :), noiseVec(2, :))
Note the addition of the Cholesky decomposition to get the covariance of the noise to match.
Read more here

请先登录，再进行评论。

更多回答（1 个）

getahun kibret
getahun kibret 2023-2-22
i want to get a matlab code that adds AWGN with zero mean and variance 0.85

类别

在 Help Center 和 File Exchange 中查找有关 Get Started with Curve Fitting Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by