ВУЗ:
Составители:
закон распределения в диапазоне [0,1] с последующем пересчетом искомых
чисел по формуле:
arabR
nn
+
−
=
)( , (3.33)
Генерация случайных чисел с нормальным законом распределения,
имеющим нулевое математическое ожидание (
m=0) и единичную дисперсию
(
σ
2
=1) проводится по формуле
()
1
2cos
1
ln2
−
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
=
n
n
n
r
r
n
π
. (3.34)
Для получения случайных чисел N
n
с другими параметрами закона
распределения используют формулу:
mnN
nn
+
=
σ
. (3.35)
Генерация случайных чисел с экспоненциальным законом распределения
проводится методом обратной функции по формуле
λ
)1ln(
n
n
r
E
−
−=
. (3.36)
Программа и результаты расчетов приводятся ниже.
r(1)=.1234567; %
Начальное значение случайного числа
l=.01; % Параметр экспоненциального закона распределения
for i=2:1000
x=37*r(i-1);
z=floor(x); %
Выделение целой части числа х
r(i)=x-z; %
Вычисление случайных чисел с равномерным законом
n(i)=sqrt(2*log(1/r(i)))*cos(2*pi*r(i-1)); % Вычисление случайных чисел с
нормальным законом
end
rr=rand(1000,1); %
Вычисление случайных чисел с равномерным законом в
MATLAB
hist(r) % Построение гистограммы для r
pause
hist(rr) %
Построение гистограммы для rr
pause
nr=randn(1000,1); %
Вычисление случайных чисел с нормальным законом в
MATLAB
hist(n)
pause
hist(nr)
pause
закон распределения в диапазоне [0,1] с последующем пересчетом искомых
чисел по формуле:
R n = (b − a )rn + a , (3.33)
Генерация случайных чисел с нормальным законом распределения,
имеющим нулевое математическое ожидание (m=0) и единичную дисперсию
(σ2 =1) проводится по формуле
⎛1⎞
nn = 2 ln⎜⎜ ⎟⎟ cos(2πrn−1 ). (3.34)
⎝ rn ⎠
Для получения случайных чисел Nn с другими параметрами закона
распределения используют формулу:
N n = σn n + m . (3.35)
Генерация случайных чисел с экспоненциальным законом распределения
проводится методом обратной функции по формуле
ln(1 − rn )
En = − . (3.36)
λ
Программа и результаты расчетов приводятся ниже.
r(1)=.1234567; % Начальное значение случайного числа
l=.01; % Параметр экспоненциального закона распределения
for i=2:1000
x=37*r(i-1);
z=floor(x); % Выделение целой части числа х
r(i)=x-z; % Вычисление случайных чисел с равномерным законом
n(i)=sqrt(2*log(1/r(i)))*cos(2*pi*r(i-1)); % Вычисление случайных чисел с
нормальным законом
end
rr=rand(1000,1); % Вычисление случайных чисел с равномерным законом в
MATLAB
hist(r) % Построение гистограммы для r
pause
hist(rr) % Построение гистограммы для rr
pause
nr=randn(1000,1); % Вычисление случайных чисел с нормальным законом в
MATLAB
hist(n)
pause
hist(nr)
pause
Страницы
- « первая
- ‹ предыдущая
- …
- 87
- 88
- 89
- 90
- 91
- …
- следующая ›
- последняя »
