MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn moreOpportunities for recent engineering grads.

Apply Today
Asked by vipul utsav on 15 Feb 2013

clc close all; clear all; n=0+(10*randn(400)); nn=imcrop(n,[0 0 16 16]); std(nn(:))

I generate noise '0+10*randn' function,then i found std(:) of noisy image but problem is that when i run matlab code using F5 in matlab,every time std(:) is different(above code) 9.5163,10.5658,9.5044,9.40,10.41

so my question is that which value should taken of std(nn(:)) of above code.

*No products are associated with this question.*

Answer by Vishal Rane on 15 Feb 2013

The command

randn(400)

generates a random 400*400 matrix every time. Hence the random result. Are you not expecting this ?

Walter Roberson on 15 Feb 2013

Right. With larger matrices, the std() will vary *less* between runs, but they are going to be different.

randn() approaches a standard deviation of 1 as the number of trials approaches infinity.

Answer by Jan Simon on 15 Feb 2013

Edited by Jan Simon on 15 Feb 2013

It is not clear, what "one particular value" could be, when you add random noise. When you want to initialize the random number generator, such that RAND replies the same values in each run:

s = RandStream.create('mt19937ar', 'seed', 5489); RandStream.setDefaultStream(s);

This is explained in the documentation, see `doc rand`.

vipul utsav on 15 Feb 2013

in above code, i have put value for each run. and it are different.

so my confusion is that which value should be taken of std(nn(:)).

Jan Simon on 15 Feb 2013

I cannot imagine what "I put value" means. I still do not understand, what you are looking for: You state, that the standard deviation of a random matrix varies from call to call, because the values replied by the random number generator vary - according to the definition of "random numbers". I've suggested to use a fixed seed, such that the random number generator is initialized by the same value in each run. Did you try to insert the suggested lines in your code and is your problem solved then? If not, please explain the difference between the results and your expectations.

Please, vipul, try to explain, what you actually want. Letting us guess it, wastes your and our time.

Image Analyst on 15 Feb 2013

I think for some reason he's thinking he should get exactly 10 out since he specified 10 as the stddev. However that doesn't agree with his response to Walter, where he indicated that he realized that the stddev will vary slightly. So, like you I have no idea what he wants and we need to see if he'll clarify.

## 1 Comment

Direct link to this comment:http://www.mathworks.com/matlabcentral/answers/63560#comment_129743

Please take a look on your code again. Half of the lines perform a brute clearing of Matlab's currently parsed functions, the command window, the variables and open windows. This is neither useful nor does it help to solve the problem. See e.g. Answers: 22301