1、clc;clear all;close all;%运用直接抽样法产生指数分布、瑞利分布的随机变量序列%运用直接抽样法产生指数分布的随机变量序列% (1)使用直接抽取的方法 先产生0,1间的随机数列% (2)利用累积分布函数的反函数,得到的结果就是符合指数分布的序列 N=10000; %产生随机数的个数lambda=2; %参数入y=rand (1ZN) ; %生成N个0, 1间均匀分布随机数x=-log (1-y) /lambda; %生成指数分布随机数n, xout =hist (xf 100) ;%分区间统计随机数出现概率nn=n/N/mean(diff(xout);figure;bar
2、xout,nn,l) ;hold on; %画图验证随机数是否符合概率密度函数plot(0 xout,lambda*exp(-lambda*(0 xout), 1r1 );xlabel ( * x 1) ; ylabel ( fp (x) ) ; title (1直接抽样法产生指数分布的随机变量序列1 ); hold off;%运用直接抽样法产生瑞利分布的随机变量序列 N=10000; %产生随机数的个数y=rand (1,N) ; %生成N个0, 1间均匀分布随机数 sigma=2; %参数 sigmaz=sqrt (-2夫sigma八2夫log (1-y) ) ; %生成瑞利分布随机数m,zout=hist(z,100);mm=m/N/mean(diff(zout);figure;bar(zout,mm,1);% hold on; %画图验证随机数是否符合概率密度函数% plot ( 0 zout, (0 zout*(1/sigmaA2)*exp(-(0zout)人2/(2*sigmaA2), 1r f);xlabel ( 1 z 1) ; ylabel (p (z) 1 ; title (!直接抽样法产生瑞利分布的随机变量序列); hold off;