Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 140
Скачиваний: 1
заменяют приближенной зависимостью, нелинейной относительно математического ожидания тх и формально линейной относительно Х°:
У = Ф о + й1Х®. |
(1.70) |
Величина ср„ представляет собой среднюю статистическую характеристику нелинейности; k l — статистический коэффициент усиления нелинейности по случайной составляющей. Для нелиней ностей, имеющих нечетные характеристики, функцию ср0 можно выразить в виде
ср0 = /е 0тх, |
(1.71) |
где /г0 — статистический коэффициент усиления по математиче скому ожиданию.
Функция ср0 и статистические коэффициенты усиления /г0, k x определяются на основании критерия вероятностной эквивалентности зависимостей (1.65) и (1.70). С точки зрения корреляционной теории таким критерием является условие равенства соответственно матема тических ожиданий п дисперсий истинной и аппроксимирующей слу чайных функций или условие минимума средней квадратической
ошибки аппроксимации функции (1.65) линеаризованной |
зависи |
мостью (1.70). |
|
При применении первого критерия имеем |
|
Шу = фо> Dy = k\Dx, |
(1. /2) |
где niy — математическое ожидание функции У; ср0— математиче ское ожидание аппроксимирующей зависимости (1.70); Dx, Dy —■ дисперсии переменных X, Y соответственно. Второе равенство (1.72)
служит для определения коэффициента
|
|
|
(1.73) |
где введены |
обозначения ау = ]/D ' |
ах = |
] /Dx— средние квад |
ратические |
отклонения функций Y |
и X. |
Знак коэффициента k x |
определяется направлением возрастания или убывания в функции ср при изменении аргумента. Если функция нечетная, то коэффициент k 0 определяется формулой
* . = * ■ • |
0 -7 4 ) |
"< Х |
|
Второй критерий состоит в выполнении условия минимума вто
рого начального момента ошибки за |
счет выбора ср0, /ех: |
|
г) = М { [ К - Фо- Л 1Х»Н, |
(1-75) |
|
где М — оператор математического ожидания. |
(1.75), |
|
Выполняя операции, указанные в правой части выражения |
||
получим |
|
|
г) = М [У2] -j- Фд -|- k f i x |
2ф0ту 2kyM [УХО]. |
(1.76) |
30
При известных М [F 2], ту, Dx, М [FX0] величина г) является функцией сро и /гх. Приравнивая нулю частные производные г) по ф0 и /гх, получим
Фо = ту- /гх = -L М [КХ°] • |
(1.77) |
и Х |
|
Для вычисления функции ср0 и коэффициентов /г0 и /гх необходимо знать одномерную плотность вероятности f (х) для каждого текущего момента времени t на входе в нелинейность. Формулы (1.73), (1.77) принимают вид
00
Фо = |
\ Ф {x)f(x)dx; |
(1.78) |
|
— СО |
|
|
СО |
|
|
j Ф2(* )/М — фо |
(1.79) |
|
— со |
|
00 |
|
|
*12)=-$■ j |
(X — тх) ф (х) f (х) dx, |
(1.80) |
■*—со
где верхние индексы у коэффициентов k[l) указывают способ линеари
зации. Расчетные формулы для /г{2) проще формул для Ар. Точность же окончательных расчетов практически одинакова.
В общем случае при вероятностном исследовании нелинейные безынерционные элементы автоматических систем управления можно представить как многомерные функциональные зависимости с неадди тивными связанными входными случайными сигналами:
Г = ф (Хх, . . ., Х п), |
(1.81) |
где ф — произвольная однозначная нелинейная функция входных не аддитивных связанных случайных переменных Х х, . . ., Х п. В част ности это относится к однозначным характеристикам реальных одно мерных нелинейных и множительных элементов. К такому же виду приводятся неоднозначные нелинейные характеристики гистерезис ного типа, если в число входных переменных включить производную, определяющую неоднозначность. Наконец, безынерционные нели нейные и линейные элементы со случайными параметрами также можно рассматривать как многомерные функциональные зависимости с неаддитивными входными случайными сигналами.
Представим каждую случайную переменную Х£(t) в виде
X l (l) = mXi(l) + X? (О,
где mx (i) — математическое ожидание; Х° (I) — центрированная
составляющая случайной функции Х(. (/). Линеаризованную зави симость для нелинейности (1.81) примем в форме выражения
>" = Фо+ |
(1.82) |
i= i
31
где ср0 и k t — неизвестные неслучайные функции вероятностных ха рактеристик входных случайных процессов.
Функция ф0 называется статистической характеристикой нелиней ности, a ki — статистическими коэффициентами усиления случай ных составляющих. Для их определения применим критерий мини мума средней квадратической ошибки аппроксимации
т) = М Ф № . ■• •, Х „ ) - Ф о - £ ktX\ |
(1.83) |
i=I |
|
где М — оператор математического ожидания. Выполнив операции, указанные в выражении (1.83), получим
Ч = М [ср2] + Ф5 -f £ |
kik/MlxUl] — |
|
L, /=I |
|
|
— 2<р0М [ср] — 2 Д |
к.М [срХ°]. |
(1.84) |
При известных М [ср2], М [ср], М [Х9Х9], М [срХ?] величина т)
является функцией параметров ср0, к... Поэтому к выражению (1.84) можно применить необходимое условие экстремума функции г) по параметрам ср0 н kt. Приравнивая нулю частные производные функции т| по ср0 и /г,., получаем уравнения
ср0 = |
М [ с р ^ , . . ., Х„)1; |
(1.85) |
£ |
т [Х?Х?] = М [срХу] |
(1.86) |
(/ = 1, • • ., п).
Разрешая систему уравнений (1.86) относительно кр получим
П
(1.87)
/=1
(i = 1, • ■ п),
где А — определитель системы уравнений (1.86); А) — алгебраиче ское дополнение элемента г'-го столбца, /-й строки определителя А.
В частном случае для несвязанных случайных переменных, вхо дящих в зависимость (1.81), формулы (1.87) принимают вид
|
М [у*?] |
ki |
( 1.88) |
м [xf] |
При вычислении статистической характеристики ср0 и коэффи циентов k( по вышеприведенным формулам необходимо знать закон распределения совокупности переменных f (хъ . . ., хп) для каждого текущего момента времени. В замкнутых автоматических системах истинный закон распределения переменных на входе в нелинейный
32
элемент неизвестен. При выполнении статистической линеаризации
ивычислении функций cp0, kt принимают некоторый нормальный эк вивалентный закон распределения, имеющий параметры — первые
ивторые вероятностные моменты, совпадающие с фактическими. Это основное допущение метода статистической линеаризации бази руется на свойстве сложных систем, включающих инерционные ли нейные цепи, нормализовывать законы распределения случайных переменных.
Втаком случае величины ф0 и /г,- оказываются функциями мате матических ожиданий тх. и корреляционных моментов связи
Qt i (t) = M [*?(/) Х° (/)].
Итак, предположим, что в общем случае сигнал на входе в нели нейный элемент имеет закон распределения, заданный в следующей форме:
/(*1. |
|
ехр |
ГД*] |
(1.89) |
|
2Д. |
|||
где А* — окаймленный определитель вида |
|
|
||
ец ... |
01* |
Xi — тх, |
|
|
Ош* ■■ |
0/т |
хп —тк |
(1.90) |
|
|
||||
пгх |
хп |
|
||
1 — |
"Хп — |
0 |
|
|
Ал |
|
|
При нормальной плотности вероятности (1.89) для коэффициентов статистической линеаризации k£ имеет место также другая формула, которая удобна для получения их в'случае сложных нелинейностей. Эта формула получается путем дифференцирования выражения ср0 частным образом по тХ£:
дфр |
J ••• (« )• • • |
J ф (хь • • х п) |
X |
|
дтх. |
||||
|
|
|
||
|
X дтхд . f (*i. ■ |
х„) dxl. . . dxn. |
(1.91) |
Имея в виду формулы (1.89) и (1.90), получим
дтх. = f(xIf |
(1.92) |
Л1 |
|
Подставляя выражение (1.92) в уравнение (1.91), получим
дфр |
(1.93) |
дтХ£ |
/= i |
|
3 в. С. Пугачев |
33 |
Сравнивая правые части формул (1.87) и (1.93), находим, что
и _ |
дфо |
(1.94) |
|
~ |
дтх.' |
||
|
(t = 1, . . ., п)
Легко показать, что необходимые условия экстремума (1.85) и (1.86) функционала (1.84) определяют в данном случае минимум выражения гр Для этого будем считать, что искомые коэффициенты в выражении (1.84) равны ср0 + к0, ki + иг Подставляя их в выраже ние (1.84), получим
|
|
|
т) = |
М |
|
|
|
32 |
|
|
|
|
|
Ф — фо — S |
|
q |
|
||||
|
|
|
|
|
|
k-iX ~Ь U -(- |
|
|||
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
+ |
М Е |
и Х |
+ |
2«0[ фо — 44ср] -\- |
|
||
|
|
|
|
L i—l |
|
|
|
|
|
|
|
|
|
+ 2 Е щ м [ Фх?] - |
Е fyM [хи° |
(1.95) |
|||||
|
|
|
|
i=l |
l |
|
/=1 |
|
|
|
Если фо и kc определяются из уравнений (1.85), |
(1.86), то выра |
|||||||||
жение |
(1.95) |
принимает вид |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
Г |
п |
Д 2 |
|
11 = М |
ф — фо — Е ь Х |
|
-|- Ко + М |
Е |
|
||||
Так |
как величина |
i—l |
|
|
|
1=1 |
|
|||
|
|
|
|
|
|
|||||
|
|
|
|
«6+ м |
S |
|
о |
|
|
|
|
|
|
|
|
|
1=1 |
|
|
|
|
при любых и 0 |
и к,-, то г) принимает минимальное значение только при |
|||||||||
и о = 0, ис = |
0, |
т. е. при фо и kt, определяемых уравнениями (1.85), |
(1.86). Рассмотрим примеры нелинейностей и их статистические линеа
ризованные эквиваленты |
[35]. |
|
Пример 1.6. Определить статистические характеристики нелинейности типа |
||
идеального реле (рис. 1.7, а) |
У = |
/ sign X. |
|
||
Линеаризованная зависимость имеет вид |
||
У = |
Фо + |
Фо = ь0тх. |
Для вычисления функции ср0 применим формулу (1.78) при нормальном законе распределения X с параметрами тх, Dx:
о
Фо |
I |
Г |
( |
(х —тх)"- \ |
||
V2nD x J |
I |
20* |
) |
|||
|
||||||
|
|
-— 00 |
|
|
|
|
|
У ш |
|
( |
(X— тх)-\ |
dx. |
|
|
~х .)ехр Г |
2D x |
|
|||
|
|
о |
|
|
|
34