Question

怎样用MATLAB语句产生均值为零服从瑞利分布的随机数

Answer 1

Matlab中提供了直接的函数，生成瑞利分布随机数raylrnd

下面我提供一个直接的函数吧

function x = RelayDist(x0,sigma,n)
x = zeros(n,1);
for i=1:n
r = MixMOD(x0,10,1);
k = 0;
while r(10) == 0
k = k + 1;
r(10) = power(2,k);
r = MixMOD(r(10),2,1);
end
y = -2*log(r(10));
x(i) = sigma*sqrt(y);
x0 = x(i);
end

function r = MixMOD(x0,n,type)
format long;
M1 = power(2,31);
M2 = power(2,35);
a1 = 314159269;
a2 = power(5,15);
c1 = 453806245;
c2 = 1;
r = zeros(n,1);
x = zeros(n+1,1);
x(1) = x0;

if type == 1
for i=2:n+1
y = a1*x(i-1)+c1;
x(i) = mod(y, M1);
r(i-1) = x(i)/M1;
end
else
for i=2:n+1
y = a2*x(i-1)+c2;
x(i) = mod(y, M2);
r(i-1) = x(i)/M2;
end
end
format short;