运用直接抽样法产生指数分布的随机序列

clc;clearall;closeall;

%利用直接抽样法形成指数分布、瑞利分布的随机变量序列

%%

%利用直接抽样法形成指数分布的随机变量序列

%(1)使用直接抽取的方式先形成[0,1]间的随机数列

%(2)运用累积分布函数的反函数指数分布matlab，得到的结果就是符合指数分布的序列

N=10000;%形成随机数的个数

=2;%参数λ

y=rand(1,N);%生成N个[0指数分布matlab，1]间均匀分布随机数

x=-log(1-y)/;%生成指数分布随机数

[n,xout]=hist(x,100);%分区间统计随机数出现机率

nn=n/N/mean(diff(xout));

;

bar(xout,nn,1);holdon;%绘图验证随机数是否符合机率密度函数

plot([0xout],*exp(-*([0xout])),'r');

('x');('p(x)');title('直接抽样法形成指数分布的随机变量序列');holdoff;

%%

%利用直接抽样法形成瑞利分布的随机变量序列

N=10000;%形成随机数的个数

y=rand(1,N);%生成N个[0，1]间均匀分布随机数

sigma=2;%参数sigma

z=sqrt(-2*sigma^2*log(1-y));%生成瑞利分布随机数

[m,zout]=hist(z,100);

mm=m/N/mean(diff(zout));

;

bar(zout,mm,1);

%holdon;%绘图验证随机数是否符合机率密度函数

%plot([0zout],([0zout]*(1/sigma^2))*exp(-([0

zout])^2/(2*sigma^2)),'r');

('z');('p(z)');title('直接抽样法形成瑞利分布的随机变量序列');holdoff;