gyzhshfb.m

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

%
%产生服从高斯分布的随机数
%
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
%
%利用高斯分布的随机数产生广义指数分布的随机数
%
s=2; %令s=2
for i=1:2048
    x(i)=u1(i)+sqrt(2*s);
    z(i)=x(i)^2+u2(i)^2;
end
HIST(z,50);
mean=mean(z);
var=var(z);