matlab生成随机矩阵

原文地址：matlab生成随机矩阵 作者：LFF_Coder

一. matlab里和随机数有关的函数：

（1） rand：产生均值为0.5、幅度在0~1之间的伪随机数

（2） randn：产生均值为0、方差为1的高斯白噪声

（3） randperm(n):产生1到n的均匀分布随机序列

（4） normrnd(a,b,c,d)：产生均值为a、方差为b大小为cXd的随机矩阵

rand

rand(n)：生成0到1之间的n阶随机数方阵

rand(m,n)：生成0到1之间的m×n的随机数矩阵

randn

randn（）命令是产生白噪声的，白噪声应该是0均值，方差为1的一组数，同rand有randn(n)，randn(m,n)

rand是0-1的均匀分布，randn是均值为0方差为1的正态分布

二. 功能：生成服从正态分布的随机数

R＝normrnd(MU,SIGMA)

R＝normrnd(MU,SIGMA,m)

R＝normrnd(MU,SIGMA,m,n)

说明:

R＝normrnd(MU,SIGMA)：生成服从正态分布(MU参数代表均值，DELTA参数代表标准差)的随机数。输入的向量或矩阵MU和SIGMA必须形式相同，输出R也和它们形式相同。标量输入将被扩展成和其它输入具有 相同维数的矩阵。

R＝norrmrnd(MU,SIGMA,m)：生成服从正态分布(MU参数代表均值，DELTA参数代表标准差)的 随机数矩阵，矩阵的形式由m定义。m是一个1×2向量，其中的两个元素分别代表返回值R中行与列的维数。

R = normrnd(MU,SIGMA,m,n)

其中MU为均值，SIGMA为标准方差，m、n为矩阵大小；

三、求矩阵中所有元素的均值和方差

x是一个m*n的矩阵

均值：sum(x(:))/(m*n)

方差：var(x(:))

四、已知直方图，随机生成一个图像，使得新图与原图的直方图相同

img=imread('raw.bmp');

h=imhist(img);

h=h/sum(h);

x=1:1:256;

figure,plot(x,h);

cp=cumsum(h);

figure,plot(x,cp);

dst=uint8(zeros(size(img,1),size(img,2)));

for row=1:size(img,1)

   for col=1:size(img,2)

       dst(row,col)=x(sum(cp<=rand(1))+1); 

   end

end

figure,imshow(dst);