Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Математическое моделирование радиотехнических устройств и систем

..pdf
Скачиваний:
139
Добавлен:
05.02.2023
Размер:
2.45 Mб
Скачать

 

 

1

 

 

2

1

2

 

 

 

 

 

Y

n

7

1

8 Y

n

7

 

 

7

T

5G3nT,Y1n24 7 G3

 

 

 

1

 

0,1,2,

(2.42)

2

3n 7 14T,Y1n 7 1246, n 8

 

 

 

 

 

 

 

 

 

 

 

где 1 1n 5 12 6 Y1n2 5 TG3nT,Y1n24 . Анализ (2.42) показывает, что тан

Y

генс угла наклона элемента ломаной Эйлера на текущем шаге опреде ляется как среднее арифметическое тангенсов угла наклона ломаной Эйлера на текущем и следующем за ним шагах обыкновенного метода Эйлера. Такая процедура позволяет уменьшить ошибку вычислений и сделать ее пропорциональной T3.

1

 

 

11

 

 

 

 

12

 

 

2

1

3

23

Рис. 2.12

Метод Рунге–Кутта. Данный метод является наиболее точным и используемым на практике среди одношаговых методов численного ре шения обыкновенных дифференциальных уравнений. Правило вычис ления значения очередного отсчета выходного процесса для метода Рун ге–Кутта выглядит следующим образом:

Y1n 512 6 Y1n2 5 T 3K0 5 2K1 5 2K2 5 K34, n 6 0,1,1,

(2.43)

 

6

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

K0 5 G3nT,Y1n24;

 

 

 

 

 

K1

6

6

1 7

 

8

1

7

 

5 G9

9 n 8

T,Y1n2

2

K0 ;

 

 

 

 

2

 

 

 

 

K2

6

6

1 7

 

8

1

7

 

5 G9

9 n 8

T,Y1n2

2

K1 ;

 

 

 

 

2

 

 

 

 

K3 5 G3

3n 8 14T,Y

1n2 8 K2

4.

 

Используя этот метод, оказывается возможно сделать ошибку вычис лений пропорциональной T5.

61

3. ОБРАБОТКА РЕЗУЛЬТАТОВ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

3.1. Оценка закона распределения вероятностей

Как было отмечено во введении, особенностью радиотехнических систем является постоянное воздействие на них случайных факторов. Следовательно, результаты моделирования будут также носить случай ный характер. Поэтому, принимая во внимание, что время эксперимен та и объем полученных данных ограничены, необходимо так обрабаты вать результаты, чтобы получаемые оценки наилучшим образом дава ли представление о свойствах и параметрах моделируемых устройств.

Поскольку, вследствие ограниченности полученных в ходе модели рования данных, можно говорить лишь об оценивании тех или иных свойств и параметров, необходимо определить какими качествами дол жны обладать эти оценки. Из теории статистического оценивания из вестно, что качество оценок определяется следующими показателями:

1. Состоятельность. Оценка 1 называется состоятельной, если ве x

роятность ее отклонения от истинного значения оцениваемого пара метра x при увеличении размера выборки данных N стремится к нулю, т. е. выполняется следующее условие:

lim Pr

1

 

1

 

4 5

6 0,

 

 

 

 

 

x 3 x

 

(3.1)

N12

 

 

 

2

 

 

 

 

 

 

где 1 – произвольное, положительное, сколь угодно малое число. Смысл условия (3.1) состоит в том, что для состоятельной оценки при увели

чении объема наблюдений значения оценки все ближе и ближе группи

руются вокруг истинного значения параметра.

 

1

1

2.

Смещенность. Смещением bx оценки x называется разность

 

1

(3.2)

 

bx 3 E1x24 x,

где 112 – математическое ожидание оценки. Оценка называется не

x

E

смещенной, если b1 1 0 . Это означает, что центр группирования значе

x

ний несмещенной оценки совпадает с истинным значением оцениваемо го параметра.

3. Эффективность. При достаточно общих предположениях оказы

1

 

вается, что дисперсия оценки x при фиксированном значении размера

выборки N не может быть меньше, чем некоторая величина D0 , т. е.

1

(3.3)

D1x23 D0.

62

Условие (3.3) называется неравенством Крамера–Рао, а величина D0 – соответственно границей Крамера–Рао. Граница D0 может быть вычислена достаточно просто

 

 

 

 

 

 

 

4

 

 

 

3b1

52

 

 

 

 

 

 

 

 

 

71 6

 

 

x

8

 

 

 

 

 

 

 

 

 

 

3x

 

 

 

 

 

 

 

 

D0 6

9

 

 

 

 

,

 

 

 

 

 

 

 

E

3

2 L1r x2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

(3.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где L1r

 

x2

 

 

 

3x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– логарифм функции правдоподобия наблюдаемой выбор

ки

r 3 1r1,1,rN 2

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

; E 3

– оператор вычисления математического ожида

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

ния. За меру эффективности оценки x принимается величина

 

 

 

 

 

 

 

 

3 4

 

D0

.

 

 

 

 

 

 

 

 

 

 

 

D1x2

 

 

 

(3.5)

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

В силу неравенства Крамера–Рао 0 1 2 11 . Оценка называется эф фективной, если 1 2 1 , т. е. ее дисперсия минимальна и равна границе Крамера – Рао.

При обработке результатов измерений следует стремиться к тому, чтобы использовать состоятельные, несмещенные и эффективные оцен ки. Это удается далеко не всегда. Поэтому иногда вышеназванные тре бования смягчают, заменяя свойства несмещенности и эффективности соответствующими асимптотическими свойствами:

4. Асимптотическая несмещенность:

lim b1 1 0.

(3.6)

N12 x

5. Асимптотическая эффективность:

lim 1 2 1.

(3.7)

N12

 

Рассмотрим теперь тот перечень задач, с которыми чаще всего встре чается исследователь, осуществляющий обработку результатов изме рений.

Обычно при обработке результатов машинного моделирования воз никают следующие проблемы:

определение эмпирического закона распределения вероятности оцениваемых в ходе эксперимента параметров;

проверка совпадения эмпирического закона распределения веро ятности с модельным распределением;

оценка параметров распределения оцениваемого параметра. Рассмотрим, как решаются вышеназванные задачи.

Пусть в результате эксперимента для оцениваемого параметра x

получен ряд оценок 1 ,1 ,2,1 . Этот ряд называется выборкой объема x1 x2 xN

63

1 ,1 ,2,1 x1 x2 x
K 112 3,2 lg N.
1 1 2
f x

(размера) N. Существуют два способа оценки эмпирического распреде ления вероятности: оценка плотности распределения и оценка функ ции распределения.

При оценивании плотности распределения (дифференциальной фун

кции) весь интервал возможных значений оценки 1 разбивают на K x

интервалов. Длины интервалов не обязательно должны быть равными.

1

1

1

подсчитывают количество по

После получения выборки x1,x2,2,xN

паданий выборочных значений в каждый из K интервалов – mk, k 11,K и определяют частоты

1

 

2

mk

, k 2

 

.

 

k

1,K

(3.8)

 

 

 

N

 

 

 

 

Если k й интервал имел длину 1xk , то за оценку значения плотности распределения f(x) на этом интервале берут

1

3k

 

mk

 

 

 

f1x2 4

 

4

 

, xk 1

5 x 6 xk,

(3.9)

7xk

N 7xk

 

 

 

 

 

где xk 1 и xk – границы k го интервала. Таким образом, оценка

имеет вид ступенчатой функции или гистограммы. Поэтому данный метод в литературе получил названия метода гистограмм. Можно по казать, что полученная оценка является состоятельной, несмещенной и эффективной. Непростым вопросом при использовании метода гис тограмм является выбор количества интервалов K. Считается, что для получения приемлемых результатов необходимо, чтобы в каждый ин

тервал попало не менее восьми значений случайной величины 1 . Одна x

ко это требование сложно использовать для определения K. Существу ет эмпирическое правило Штюргеса, которое утверждает, что число ин тервалов гистограммы K и объем выборки N связаны следующим соот ношением:

(3.10)

Так, при N=1000 получаем, что K 1 10 .

Обратимся теперь к оценке интегральной функции распределения. Расположим результаты эксперимента N в порядке возраста ния. Получим в результате так называемый вариационный ряд

x

1 x

121 x

.

(3.11)

1 (1)

1 (2)

1 (N)

 

Пусть x – неслучайная величина. Количество членов вариационно

го ряда (3.11) меньших x называется эмпирической частотой

 

31x2

1 (m)

5 x, m 4

1,2,M.

(3.12)

4 M, x

 

64

Эмпирической функцией распределения называется функция

 

 

 

 

 

 

3

2

(1)

 

 

 

 

 

6

1x2

 

50,

x 4 x

,

 

 

1

 

 

 

5

 

2

(m)

2

(m11)

F1x2

7

 

 

7

9M N,

x

 

8 x 4 x

.

 

N

 

 

 

 

 

 

5

2(N)

8 x.

 

(3.13)

 

 

 

 

 

 

51,

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

имеет ступенчатый вид и изменяется от 0 до 1. Скач

Функция F1x2

ки функции приходятся на выборочные значения и их величина рав

на1 N . В математической статистике доказано, что

1

F1x2 – состоя

тельная и несмещенная оценка интегральной функции распределения.

3.2. Проверка соответствия выбранной модели распределения данным эксперимента

Во многих случаях исследователь, проводивший математический

эксперимент, с той или иной степенью уверенности может предполо

1

 

 

жить, что наблюдаемая оценка x имеет некоторое модельное распреде

1

1

1

ление вероятностей. В этом случае выборочные данные x1,x2,2,xN мож

но использовать для того, чтобы либо принять гипотезу о справедливо сти сделанного предположения, либо отвергнуть ее. Задача проверки соответствия выбранной модели распределения и данных эксперимен та решается с помощью так называемых критериев согласия. Алгоритм проверки для всех критериев согласия следующий.

Пусть H0

– гипотеза о том, что полученное эмпирическое распреде

 

1

1

1

согласуется с некоторым модельным. Для

ление выборки x1,x2,2,xN

подтверждения или отрицания гипотезы H0 выбирают некоторую фун кцию U, заданную на множестве эмпирических и модельных распреде лений. Величина U, вычисленная на полученном эмпирическом и выб ранном модельном распределении, является случайной и называется решающей статистикой. В качестве решающей статистики нельзя брать любую функцию, поскольку в общем случае ее распределение ве роятности зависит и от фактического, и от модельного распределений. Однако при разработке методов согласия были найдены такие функции U, статистика которых практически не зависит от условий эксперимен та и может быть определена заранее. При использовании такой решаю щей статистики процедура принятия гипотезы H0 состоит в проверке условия

U 1 U0, (3.14) где U0 – порог сравнения, который вычисляется на основании решения одного из следующих уравнений:

65

Pr 1U 3 U02 4 5 или Pr 1U 3 U02 415 6.

(3.15)

В уравнениях (3.14) – (3.15) соответствующие вероятности вычис ляются на основании известного распределения вероятностей решаю щей статистики U. Постоянная 1 , которая называется уровнем значи мости, выбирается заранее и устанавливает допустимую вероятность ошибки при проверке гипотезы H0. Обычно 1 2 0,0510,1 . В случае, если неравенство (3.14) не выполняется, гипотеза H0 отвергается, и иссле дователь должен предложить другую гипотезу о виде модельного рас пределения. Все существующие критерии согласия различаются по виду решающей статистики. Рассмотрим наиболее часто используемые кри терии согласия.

3.2.1. Критерий Пирсона

Разобьем область определения модельной плотности распределения вероятностей F(x) на K интервалов и обозначим вероятности попада

1

p

3 Pr

1

1

k2. Согласно крите

x45

ния оценки x в k й интервал 1k через

k

 

 

рию Пирсона, решающей статистикой является следующая мера рас хождения модельного и эмпирического распределений:

K

 

 

2

K

 

2

 

1mk

3 Npk 2

 

1mk

N 3 pk 2 ,

 

U 4 5

 

4 N5

(3.16)

k 1

 

Npk

 

k 1

 

pk

 

где mk – количество значений оценки, попавших в k й интервал 1k

K

( 2mk 1 N ).

k 1

Пирсон доказал следующую теорему. Если проверяемая гипотеза H0 об истинности модельного распределения F(x) верна, то при объеме выборки N 12 закон распределения решающей статистики U зависит только от числа интервалов K и приближается асимптотически к за кону распределения 12 (хи квадрат) с s = K – J –1 степенями свободы, где J – число неизвестных параметров модельного распределения.

Выдвинув гипотезу H0 об истинности модельного распределения, тем самым устанавливаем, от какого числа параметров оно зависит. Если значения всех или части параметров не известны, то они заменя ются оценками. После определения оценок параметров модельной фун

кции распределения вычисляются вероятности

1

 

pk 3 Pr 1x45k2

 

pk 3 F1xk 1 2 4 F1xk 2 , 5k 3 xk 1 4 xk

(3.17)

попадания оценки в k й интервал.

 

 

66

N объема N
борки 1 ,1 ,2,1 x1 x2 x

Из теоремы Пирсона следует, что какова бы ни была гипотетическая функция распределения вероятностей F(x) случайная величина U при N 1 2 имеет распределение 12

 

 

 

Pr

1

U 6 u 7

53s2,u

24

,

 

(3.18)

 

 

 

53s24

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

где 3 1 s

 

22

u 2

 

 

 

 

 

22

1

2,u

4 5 xs 2 1e xdx – неполная, а 3 1 s

4 5xs 221e2xdx – пол

 

 

 

0

 

 

 

 

 

 

0

ная гамма функции, которые могут быть вычислены либо численно, либо по таблицам.

Таким образом, процедура проверки гипотезы H0, согласно крите рию Пирсона, выглядит следующим образом:

1.Задаем значение доверительной вероятности Pд 1 12 3 или уровня значимости 1 . Обычно 1 2 0,0510,1 .

2.По значению 1 и числу степеней свободы s на основании (3.18) находим величину порога сравнения U0.

3.На основании (3.16) и (3.17) вычисляем значение решающей ста

тистики Uи сравниваем его с порогом U0. Если значение U меньше порога U0, то гипотеза H0 принимается. Если U 1 U0 , то гипотеза отвергается.

Критерий Пирсона является одним из наиболее широко используе мых на практике и дает хорошие результаты при объеме выборки N порядка 100 и выше. Недостатками критерия являются:

1.Необходимость иметь сравнительно большую выборку.

2.Произвольное разбиения области определения модельной плот ности распределения вероятностей F(x) на K интервалов, которое не учитывает особенностей функции F(x).

3.2.2. Критерий Колмогорова

Согласно этому критерию количественной мерой соответствия мо

дельного F(x) и эмпирического 1 распределений вероятностей для вы

F

служит максимум их модуля разности

U 3 max

 

F1x2

1

 

.

(3.19)

 

 

 

4 F1x2

 

x

 

 

 

 

 

 

 

 

 

 

 

Колмогоров доказал, что если проверяемая гипотеза H0 верна, то при N 1 2 и дополнительном предположении о непрерывности F(x) функция распределения величины NU асимптотически стремится к

функции Колмогорова

67

Pr 1

NU 5 z26 K3z4

1

 

 

6 1 7 283714k e22k2z2 .

(3.20)

 

 

k

3

1

 

Представление (3.20) удобно использовать при z 11 , когда ряд схо дится быстро. При z 11 , когда ряд в (3.20) сходится медленно, удобнее пользоваться другим представлением

 

23 3

2

12

12k2122

 

 

 

 

8z2

 

 

K1z2 4

 

5e

 

 

 

.

(3.21)

z

 

 

 

 

 

k41

 

 

 

 

 

Таким образом, правило проверки гипотезы H0 , согласно критерию Колмогорова, таково:

1.Задаем значение доверительной вероятности Pд 1 1 2 3 или уровня значимости 1 . Обычно 1 2 0,0510,1 .

2.По значению 1 на основании (3.20) или (3.21) находим величину порога сравнения U0 .

3.На основании (3.19) вычисляем значение решающей статистики

U и сравниваем его с порогом U0. Если значение U меньше порога U0, то гипотеза H0 принимается. Если U > U0, то гипотеза отвергается.

Как и критерий Пирсона, критерий Колмогорова используется при достаточно больших объемах выборки (N = 50...80). Однако при ис пользовании этого критерия не требуется дополнительного разбиения области определения F(x) на интервалы.

3.2.3. Критерий Крамера–Мизеса

Согласно этому критерию, количественной мерой соответствия для выборки объема N служит значение среднего значения квадрата от клонения модельного распределения от эмпирического

1

 

1

42

 

U 5 9

3

(3.22)

7F1x2

6 F1x2

8 f 1x2 dx,

21

где f 1x2 – плотность распределения f1x2 3 dF1x2 . Подстановка (3.13) dx

в (3.22) и интегрирование позволяет получить другое выражение для решающей статистики

 

1

 

 

1

N

2

4

1

(n) 5

 

2n 1132

 

U 6

 

 

7

 

 

F8 x

9

1

 

.

(3.23)

12N

2

 

 

 

 

 

N n11

 

 

 

 

 

2N

Для величины NU в математических таблицах можно найти пре дельное при N 1 2 распределение вероятностей, на основании которо

68

го по заданному уровню значимости 1 можно определить величину по рога сравнения U0 . Дальнейшая процедура использования критерия ничем не отличается от алгоритма использования критерия Колмого рова.

Критерий Крамера–Мизеса может использоваться при малых объе мах выборки ( N 1 50 ).

3.3. Оценка моментов распределения

Оценка эмпирического распределения и проверка его соответствия модельному распределению не являются единственными проблемами, возникающими при обработке результатов эксперимента. Обычно при

обработке оценивают моменты распределения оценки 1 . Как правило, x

при этом ограничиваются лишь первыми начальными или центриро ванными моментами. Оценка моментов распределения случайной вели

чины 1 позволяет не только определить центр группирования резуль x

татов измерений и меру их разброса, но и судить о качественном харак тере распределения вероятностей.

Для случайных величин в теории вероятностей вводятся начальные моменты

 

1

 

 

Mp 3 4 xpf1x2 dx, p 3 1,2,1

(3.24)

 

21

 

 

и центральные моменты

 

 

1

1x 4 M1

2p f1x2 dx, p 3 2,3,1,

 

mp 3 5

(3.25)

21

 

 

 

где f1x2 – плотность распределения оценки 1 . Между семействами на x

чальных и центрированных моментов существует взаимно однозначное соответствие.

Наибольшее применение нашли величины, связанные с первыми четырьмя моментами распределения.

1. Математическое ожидание:

1 2 M1

(3.26)

представляет собой центр группирования результатов измерений. 2. Дисперсия:

D 1 m2

(3.27)

представляет собой меру разброса случайной величины. Для случай ных процессов дисперсия равна средней мощности процесса.

69

3. Коэффициент асимметрии:

31 4

m3

(3.28)

1m2

2

 

 

3/2 .

Унимодальное распределение с 11 2 0 имеет левую асимметрию, с

11 2 0 – правую (рис. 3.1). Если 11 2 0 , распределение симметрично.

1122

1122

 

1122

1 < 0

 

1 = 0

1 > 0

2

 

2

2

 

Рис. 3.1

 

4. Коэффициент эксцесса:

 

 

 

3 4

m4

5 3

 

2

1m2 2

 

(3.29)

2

 

характеризует остроту вершины плотности распределения. За ноль по шкале 12 принят эксцесс гауссовского распределения (в этом слу чае m4 1m2 22 3 3 ). Плотности с 12 2 0 имеют более плоскую вершину, чем гауссовская плотность, а плотности с 12 2 0 – более острую (рис. 3.2).

1122

1

> 0

 

1 = 01 < 0

2

Рис. 3.2

Начальные моменты распределения оцениваются по выборке

1

1

1

 

 

 

 

x1,x2,2,xN следующим образом:

 

 

 

 

 

1

1

N

2 p

 

 

 

M p 1

 

2xn.

(3.30)

 

 

 

 

 

 

N n11

 

 

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]