В соответствии с методом статистической линеаризации заменяем НЭ статистическим эквивалентом. Пусть структура статистического эквивалента описывается уравнением (7.12), а критерий эквивалентности определен выражением (7.13). Коэффициенты статистической линеаризации K0(mx,σx) и K1(mx,σx) зависят от вида нелинейности
f(∙), а также от математического ожидания и среднеквадратического значения случайного процесса x(t) (в стационарном режиме mx=const и σx=const). Рассматриваем линеаризованную СУ с ПФ в разомкнутом состоянии (для математических ожиданий) и (для центрированных случайных составляющих). Далее можно применить известные методы анализа линейных систем с учетом следующей особенности: параметры системы зависят от mx и σx, а значения mx и σx зависят от параметров системы. Найти неизвестные значения mx и σx удобно путем совместного решения уравнений для динамической составляющей ошибки системы mx и дисперсии случайной составляющей xо(t) процесса на входе НЭ. Эти уравнения решают методом последовательных приближений или графическим способом. После нахождения mx и σx нетрудно вычислить дисперсию флюктуационной составляющей ошибки используя известные и .
В качестве примера рассмотрим задачу анализа нелинейной системы, содержащей идеальное реле и линейную часть с ПФ (заметим, что, в соответствии с результатами разд. 7.3, в такой системе отсутствуют автоколебания). Пусть шум v(t) представляет собой экспоненциально-коррелированный процесс со спектральной плотностью мощности: где и - дисперсия и интервал корреляции v(t), зависящие от характеристик радиоприемного тракта.
Линеаризованная система для математических ожиданий описывается ПФ и имеет первый порядок астатизма. Следовательно, при g(t)=const динамическая составляющая ошибки mx=0. Линеаризованная система для центрированных случайных составляющих описывается ПФ Дисперсию случайной составляющей xо(t) процесса на входе НЭ можно найти известным способом (см. разд. 3):
(7.23)
При mx=0 (см. пример разд. 7.5), поэтому откуда
Заметим, что коэффициент K характеризует степень инерционности линейной части системы и для нормализации случайного процесса y(t) требуется выполнение условия Kτ<<1. Поэтому (если это не так, метод статистической линеаризации для анализа системы неприменим).
В случае линейно изменяющегося регулярного воздействия (g(t)=Vt) динамическая составляющая ошибки не равна нулю и определяется с помощью выражения:
(7.24)
Вводим вспомогательную переменную α и представляем выражение (7.24) в виде системы из 2-х уравнений:
На рис. 7.14, а показано семейство кривых вида построенных для десяти значений σx (от σx=0,1 до σx=1,0 с шагом 0,1). Точки пересечения этих кривых и прямой (V /K=0,2) определяют зависимость (кривая 1 на рис. 7.14, б). Эта зависимость связывает значения математического ожидания и среднеквадратического значения случайного процесса x(t) на входе НЭ, удовлетворяющие уравнению (7.24). Для этих значений определяем коэффициент K1(mx,σx) (по формуле (7.16) или (7.21)) и, с помощью (7.23), определяем еще одну зависимость (кривая 2 на рис. 7.14, б построена с использованием формулы (7.16)). Эта зависимость связывает значения математического ожидания и среднеквадратического значения случайного процесса x(t) на входе НЭ, удовлетворяющие уравнению (7.23).
Точка пересечения кривых 1 и 2 на рис. 7.14, б определяет искомые значения и ( и ).
Дисперсию флюктуационной составляющей ошибки нетрудно найти известным способом (см. разд. 3):
Таким образом, в нашем примере: динамическая составляющая ошибки (при g(t)=Vt, V /K=0,2) равна 0,12, а дисперсия флюктуационной составляющей ошибки (при σv=0,5 и Kτ=0,1) равна 0,044 (размерность соответствует отслеживаемому параметру).
7.7. Применение теории марковских процесссов для анализа нелинейных систем
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.