Идентификация объектов управления. Семенов А.Д - 89 стр.

UptoLike

закон распределения в диапазоне [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