Методы и анализ нелинейного режима работы системы ЧАП. Метод фазовой плоскости

БЕЛОРУССКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИНФОРМАТИКИ И РАДИОЭЛЕКТРОНИКИ

Кафедра РТС

РЕФЕРАТ

На тему:

"Методы и анализ нелинейного режима работы системы ЧАП. Метод фазовой плоскости"

МИНСК, 2008

К нелинейным относят системы, описываемые нелинейными дифференциальными уравнениями.

Система является нелинейной вследствие наличия в ее составе звеньев, описываемых нелинейными дифференциальными уравнениями, или имеющих нелинейную статическую характеристику (например, дискриминационную).

Нелинейный режим работы имеет место в системе при выходе ошибки слежения за пределы линейного участка (переходной режим, срыв слежения, большой уровень помех и т.д.).

Методы анализа нелинейных систем:

Метод кусочно-линейной аппроксимации. Нелинейная характеристика разбивается на ряд линейных участков, в пределах каждого из которых система описывается линейным дифференциальным уравнением. Далее на каждом из этих участков система исследуется линейными методами; находятся решения, описывающие работу системы, которые затем "сшиваются". Метод удобен при небольшом числе участков разбиения. Недостаток метода в громоздкости вычислений при увеличении количества участков.

Метод гармонической линеаризации. Нелинейный элемент (НЭ) заменяется его линейным эквивалентом. Критерий эквивалентности состоит в равенстве первой гармоники напряжения на выходе НЭ и его линейного эквивалента по амплитуде и фазе при подаче на входы НЭ и его эквивалента гармонического сигнала. Метод эффективен, когда все высшие гармоники подавляются последующими цепями.

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

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

Моделирование на аналоговых и цифровых вычислительных машинах. Не имеет ограничений на количество и вид нелинейностей, порядок дифференциального уравнения, позволяет исследовать поведение системы при детерминированных и случайных воздействиях.

Отсутствие возможностей найти аналитические зависимости для исследуемых явлений является недостатком метода.

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

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

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

Анализ нелинейного режима работы системы ЧАП

Для определения некоторых характеристик системы, произведем качественный анализ системы ЧАП (рис.1)

Рис.1. Структурная схема нелинейной системы.

Исходные данные:

─ крутизна регулировочной характеристики генератора;

─ дискриминационная характеристика;

─ нестабильность частоты генератора;

─ флюктуационная составляющая;

─ отклонение от частоты от номинального значения. .

─ постоянная времени фильтра.

Составим ДУ описывающее поведение системы:

(1)

; (2)

Подставив (8.2) в (8.1), получим

;

. (3)

В установившемся режиме ; , следовательно,

. (4);

Решение уравнения (4) может быть найдено графическим способом (рис.2).

Рис. 2.

- прямая проходящая через точку , с наклоном .

Абсциссы точек и есть решение этого ДУ.

Исследуем на устойчивость в "малом" систему в точках .

С этой целью линеаризируем дискриминационную характеристику в окрестности точек равновесия системы и представим ее зависимостью

; (5)

где - крутизна дискриминационной характеристики;

.

Подставим (5) в (3) и введем новую переменную ; в результате получим дифференциальное уравнение следующего вида:

. (6)

Уравнение (6) описывает поведение системы в окрестности точек равновесия системы. Определим исходя из алгебраического критерия условия устойчивости системы:

; .

В точке, соответствующей решению , следовательно,

Таким образом соответствует устойчивому состоянию равновесия.

В точке, соответствующей ,, но , поэтому соответствует устойчивому состоянию равновесия.

В точке, соответствующей , и , здесь условие устойчивости не выполняется.

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

Для разомкнутой системы эта зависимость линейна.

Рис.3. Зависимость частотной ошибки от первоначальной частотной расстройки.

Для замкнутой системы при увеличении увеличивается и , и в точке Б система скачком переходит в точку В: происходит срыв слежения. При дальнейшем увеличении система будет вести себя как и разомкнутая. При уменьшении система войдет в режим синхронизма в точке Г, ошибка скачком уменьшится, при этом будет меньше, чем при срыве слежения.

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

Участок В– Г соответствует решению типа 3 (устойчивому состоянию).

Участок Б – Г соответствует решению типа 2 (неустойчивому состоянию).

Участок Б – Б соответствует решению типа 1(устойчивому состоянию).

Аналогичную зависимость можно получить для системы ФАПЧ (рис.4),

Где - расстройка между частотой входного сигнала и частотой собственных колебаний опорного генератора;

- ошибка слежения по частоте.

Не для всех систем . Это определяется типом фильтра и дискриминатора. Для цифровых следящих систем и называется полосой синхронизации.

Рис.4. Зависимость частотной ошибки от первоначальной частотной расстройки.

Метод фазовой плоскости

Предположим, что поведение следящей системы описывается нелинейным дифференциальным уравнением второго порядка

. (7)

Обозначим

х = х1;

;

. (8)

Состояние системы, описываемой уравнениями (8), определяется в каждый момент времени величинами и т.е. величиной координаты и скоростью ее изменения. Это состояние системы можно отобразить точкой на плоскости с координатами , называемой фазовой плоскостью. При изменении состояния системы изображающая точка перемещается на фазовой плоскости по кривым, которые называют фазовыми траекториями. Совокупность фазовых траекторий для различных начальных условий называют фазовым портретом.

Чтобы получить уравнение фазовых траекторий, исключим из (2) время, поделив для этого второе из них на первое:

. (9)

Его решение . Каждой комбинации начальных условий соответствует свое решение уравнения (3) и своя фазовая траектория.

В качестве примера рассмотрим затухающий колебательный процесс, показанный на рис.5.

Рис.5. Затухающий колебательный процесс.

Цифрами отметим характерные точки кривой и сопоставим их с фазовой траекторией. В точке 1 х(0) 0, х(0) =0, поэтому фазовая траектория начинается на положительной полуоси абсцисс (Рис.6). В точке 2 х=0, х0, поэтому эта точка расположена на отрицательной полуоси абсцисс. В точке 3 х0, х(0) =0, и на фазовой плоскости она расположена на отрицательной части горизонтальной оси и т.д. В результате для затухающего колебательного процесса фазовая траектория имеет вид сходящейся спирали.

Рис.6. Фазовая траектория затухающего колебательного процесса.

Для затухающего монотонного процесса (рис.7а) фазовая траектория приведена на рис.7б.

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

Построение фазовых траекторий позволяет судить о свойствах нелинейных систем по переходному процессу.

Рис.7. Апериодический процесс и его фазовая траектория.

Построение фазового портрета системы обычно начинают с определения его характера вблизи точек равновесия системы, в которых производные . Координаты точек равновесия определяются, как следует из (8), равенствами , . Точки равновесия при построении фазового портрета системы называют особыми.

Поведение фазовых траекторий вблизи особых точек зависит от характера корней соответствующего характеристического уравнения

,

где

, ;

- отклонение от состояния равновесия.

Если и , то процесс является затухающим гармоническим колебанием

, (10)

где и - амплитуда и начальная фаза колебания; - его частота, равная

.

Продифференцировав выражение (10) для по времени, получим

. (11)

Фазовая траектория, построенная по приведённым выражениям для процессов и , имеет вид скручивающейся спирали (см. рис.8), получившей название – устойчивый фокус.

При и процесс является гармоническим колебанием с нарастающей амплитудой. Особая точка соответствует при этом неустойчивому состоянию равновесия и называется неустойчивым фокусом (см. рис.9).

При выполнении условия корни действительные и имеют одинаковый знак. Если они отрицательны, то особая точка является устойчивым узлом (см. рис.10). Положительным корням соответствует особая точка типа неустойчивого узла (см. рис.11). При корни действительные и имеют разные знаки. Особая точка называется седлом (см. рис.12).

Рис.8. Устойчивый фокус.

Рис.9. Неустойчивый фокус.

Рис.10. Устойчивый узел.

Рис.11. Неустойчивый фокус

.

Рис.12. Особая точка типа седла.

Для построения фазового портрета необходимо определить изоклины. Изоклиной называют геометрическое место точек в котором касательные к фазовым траекториям имеют постоянный наклон.

Уравнение изоклины:

.

Для горизонтальных касательных уравнение изоклины:

;

для вертикальных:

.

Ось абсцисс является изоклиной вертикальных касательных. Для особых точек типа узла и седла существуют изоклины, совпадающие с фазовыми траекториями: (). Они называются сепаратрисcами.

Рассмотрим пример.

Определим условия вхождения в синхронизм системы, представленной структурной схемой (рис.13), если задающее воздействие изменяется по линейному закону (t) = at и в момент включения системы при t = 0 начальная ошибка имеет конечное значение х(0) = х.

Рис.14. Дискриминационная характеристика (а) и фазовый портрет (б)

Обозначим ошибку слежения.

х(t) = х = (t) – y(t).

Тогда производная этой функции:

= = a – .

Так как в качестве фильтра системы используется интегрирующее звено, то

y(t) = k F(x) /p.

В результате уравнение ошибки примет вид

= а – k F(x).

Обозначим

= х

и, пользуясь уравнением

х= а – k F(x),

построим фазовый портрет системы в координатах (x, х) для различных значений скорости изменения задающего воздействия а.

При различных значениях а кривая х=f(x) перемещается параллель - но самой себе. На рис.14 изображено семейство кривых для положительной скорости а. Обозначим максимальное значение функции F(x) = F. Направление движения изображающей точки обозначим в соответствии с правилами: в верхней полуплоскости слева направо; в нижней – справа налево. Проанализируем фазовый портрет.

При а=0 ошибка слежения х 0 при начальных значениях | х(0) | , что следует из направления движений на фазовой траектории. Если 0 а k F, то x стремится к устойчивой точке 1, если начальное рассогласование х(0) меньше величины, соответствующей точке 2. Когда х(0) , захвата не происходит, так как x неограниченно растет. Если скорость /а/ k F, то захвата не будет ни при каких начальных условиях, поскольку нет устойчивых точек на фазовой траектории. Таким образом, условия захвата сигнала, изменяющегося с постоянной скоростью а, состоят в выполнении неравенства kFа. При этом область захвата х(0) . Величина находится из уравнения а – kF() =0. Первый корень этого уравнения соответствует точке 1 устойчивого равновесия, а второй корень, соответствующий точке 2, является искомой величиной .

ЛИТЕРАТУРА

1. Коновалов. Г.Ф. Радиоавтоматика: Учебник для вузов. – М.: Высш. шк., 2000.

2. Радиоавтоматика: Учеб. пособие для вузов. / Под ред. В.А. Бесекерского. - М.: Высш. шк., 2005.

3. . Первачев С.В. Радиоавтоматика: Учебник для вузов. - М.: Радио и связь, 2002.

4. Цифровые системы фазовой синхронизации / Под ред. М.И. Жодзишского – М.: Радио, 2000