gyrlfb.m

来自「产生了十种随机数字」· M 代码 · 共 21 行

%
%产生服从高斯分布的随机数
%
clear;
clc;
u1=rand(1,2048);  %标准正态分布的随机数
u2=rand(1,2048);  %
for i=1:2048
    x(i)=sqrt((-2)*log2(u1(i)))*cos(2*pi*u2(i));
    y(i)=sqrt((-2)*log2(u1(i)))*sin(2*pi*u2(i));  %xi,yi相互独立且正交
end
%
%利用高斯分布的随机数产生广义瑞利分布的随机数
%
a=1;   %令a=1
for i=1:2048
    g(i)=sqrt((x(i)+a)^2+y(i).^2);
end
hist(g,50);
mean=mean(g);
var=var(g);