rlfb.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
%
%利用高斯分布的随机数产生瑞利分布的随机数
%
for i=1:2048
    r(i)=sqrt(x(i).^2+y(i).^2);
end
mean=mean(r);
var=var(r);
hist(r,50);