Анализ нелинейных систем. Метод гармонической линеаризации. Применение теории марковских процессов для анализа нелинейных систем, страница 8

          В соответствии с методом статистической линеаризации заменяем НЭ статистическим эквивалентом. Пусть структура статистического эквивалента описывается уравнением (7.12), а критерий эквивалентности определен выражением (7.13). Коэффициенты статистической линеаризации K0(mxx) и K1(mxx) зависят от вида нелинейности

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(mxx) (по формуле (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. Применение теории марковских процесссов для анализа нелинейных систем