Главная страница
Навигация по странице:

  • 5.3. Алгебраические критерии устойчивости. Критерий Гурвица

  • 5.4.Частотные критерии устойчивости 5.4.1 .Критерий устойчивости Михайлова

  • Формулировка критерия устойчивости Михайлова

  • _ТАУ Конспект лекций для заочников. Основные понятия теории автоматических систем


    Скачать 0.81 Mb.
    НазваниеОсновные понятия теории автоматических систем
    Дата03.03.2023
    Размер0.81 Mb.
    Формат файлаpdf
    Имя файла_ТАУ Конспект лекций для заочников.pdf
    ТипЛекция
    #967461
    страница5 из 9
    1   2   3   4   5   6   7   8   9
    5.2. Корневой метод исследования устойчивости
    Устойчивость САУ – одно из важнейших условий ее работоспособности, так как требование устойчивости включает в себя требование затухания переходного процесса во времени. Очевидно, что система с расходящимся процессом была бы неработоспособной.
    Рассмотрим комплексную плоскость распределения корней характеристического полинома
    )
    ( p
    D
    замкнутой системы управления (см. подпункт 4.4.), изображенную на рис.5.2.
    0
    Re
    Im
    1

    3
    ,
    2

    4

    6
    ,
    5

    8
    ,
    7

    Рис.5.2. Комплексная плоскость распределения корней полинома
    )
    ( p
    D
    Под устойчивостью линейной системы понимают свойство затухания переходного
    процесса с течением времени
    , то есть
    ( )
    0
    nep
    y
    t
     при


    t
    .(5.2)
    1.
    Свойство устойчивости (5.2) имеет место тогда, когда все корни характеристического полинома
    0
    )
    (

    p
    D
    n
    i
    i
    ,
    1
    ,



    обладают отрицательными вещественными частями
    n
    i
    i
    ,
    1
    ,
    0
    Re



    . Это иллюстрируется графиками на рис.5.3.
    t
    y
    0
    t
    y
    0
    Рис.5.3. Примеры графиков устойчивых переходных процессов
    2.
    Если хотя бы один из вещественных корней, или если хотя бы одна пара
    комплексных корней будет иметь положительную вещественную часть, то переходной процесс будет расходящийся, как показано на рис.5.4.
    t
    y
    0
    t
    0
    t
    y
    0
    t
    y
    0
    Рис.5.4. Примеры графиков неустойчивых переходных процессов

    3 3.
    Если в характеристическом полиноме замкнутой системы
    0
    )
    (

    p
    D
    содержится хотя бы один корень равный нулю
    0

    i

    , или хотя бы одна пара чисто мнимых корней


    j
    j
    i


    ,
    , а все остальные корни меньше нуля
    n
    i
    i
    ,
    1
    ,
    0
    Re



    , то САУ находится на границе устойчивости.
    Причем если это действительный корень, то – граница апериодической устойчивости. Например
    Tp
    p
    W
    1
    )
    (

    . А если это мнимые корни, – то граница колебательной устойчивости. В этом случае мнимая часть корня характеристического полинома численно равна частоте незатухающих колебаний. Например
    1 1
    )
    (
    2 2


    p
    T
    p
    W
    Условие
    устойчивости, следовательно, заключается в том, что
    все
    корни
    характеристического полинома
    0
    )
    (

    p
    D
    должны располагаться в левой полуплоскости комплексного переменного, а мнимая ось плоскости корней Im служит границей устойчивости.
    5.3. Алгебраические критерии устойчивости. Критерий Гурвица
    Алгебраические критерии устойчивости позволяют судить об устойчивости и неустойчивости системы непосредственно по коэффициентам характеристического полинома без вычисления его корней. В ТАУ наибольшее применение получили критерии Гурвица и Рауса.
    Рассмотрим критерий Гурвица.
    Предварительно определим необходимое условие устойчивости.
    Предположим, что характеристический полином замкнутой системы
    )
    ( p
    D
    имеет вид
    -1 2
    -1 2
    1 0
    ( )
    ( )
    ( )
    0
    n
    n
    n
    n
    D p
    A p
    kB p
    a p
    a p
    a p
    a p a




    


     .(5.3)
    Необходимым, но недостаточным условием устойчивости САУ является положительность
    всех коэффициентов характеристического полинома замкнутой САУ
    n
    i
    a
    i
    ,
    0
    ,
    0



    Доказательство. Предположим, что все корни характеристического полинома известны и имеют отрицательную вещественную часть
    n
    i
    i
    ,
    1
    ,
    0
    Re



    . Тогда (5.3) можно разложить на сомножители
    0
    )
    (
    1




    n
    i
    i
    p
    p
    1 2
    3 3
    ( )
    (
    )(
    )(
    )(
    -
    )
    (
    ) 0
    n
    n
    D p
    a p
    p
    p
    j
    p
    j
    p













    

     .(5.4)
    Произведение пары комплексных корней равно


    2 3
    2 2
    2 3
    3 3
    2
    )
    (
    )
    )(
    (


















    p
    p
    p
    j
    p
    j
    p
    После перемножения всех скобок в (5.4) получим в уравнении только положительные коэффициенты. Но так как положительные коэффициенты получаются и при положительных вещественных частях комплексных корней, то в общем случае, положительность коэффициентов уравнения (5.3) недостаточна для устойчивости системы в целом. Хотя все вещественные корни при положительности коэффициентов уравнения будут обязательно отрицательными.
    Только при
    1, 2
    n

    необходимые условия являются и достаточными условиями устойчивости. А при
    3

    n
    условие положительности коэффициентов
    n
    i
    a
    i
    ,
    0
    ,
    0



    является
    только необходимым условием устойчивости.
    Критерий устойчивости Гурвица (без доказательства).
    Для устойчивости линейных систем необходимо и достаточно, чтобы при
    положительности всех коэффициентов
    n
    i
    a
    i
    ,
    0
    ,
    0



    характеристического полинома ( ) 0
    D p

    все
    n
    главных определителей матрицы Гурвица
    n
    k
    k





    0
    ,
    0
    были положительны.

    4
    Матица Гурвица обозначена буквой «
     » и имеет вид
    Условие положительности главных определителей Гурвициана
    1 3
    5 1
    1 2
    4 2
    1 2
    3 1
    3 3
    3 2
    1 1
    4 5
    2 0
    1 0
    2 0
    0 0
    0,
    0 0
    (
    )(
    )
    (
    ) 0,
    0 0
    0 0,
    (
    )
    [(
    )(
    )
    (
    )] 0,
    0 0
    0 0.
    0 0
    0
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    n
    a
    a
    a
    a
    a
    a
    a
    a
    a
    a a
    a
    a
    a
    a
    a
    a
    a a
    a
    a
    a
    a
    a
    a
    a

















     

     


     

     
     













    (5.5)
    Рассмотрим подробней характеристический полином первого порядка,
    1

    n
    ,
    0
    )
    (
    0 1



    a
    p
    a
    p
    D
    , тогда
    0
    a


    , то есть
    0 0
    1



    a
    Для характеристического полинома второго порядка,
    2

    n
    ,
    0
    )
    (
    0 1
    2 2




    a
    p
    a
    p
    a
    p
    D
    , тогда
    0 2
    1 0
    a
    a
    a


    , то есть
    1 0
    2 0
    1 0,
    0.
    a
    a
     

       
    Очевидно, что для 1, 2
    n

    достаточно, чтобы
    0
    ,
    ,
    2 1
    0

    a
    a
    a
    Для характеристического полинома третьего порядка,
    3

    n
    ,
    3 2
    3 2
    1 0
    3
    ( )
    0,
    0
    D p
    a p
    a p
    a p a
    a






    Размерность Гурвициана
    )
    3 3
    (
    dim



    ,
    2 0
    3 1
    2 0
    0 0
    0
    a
    a
    a
    a
    a
    a
     
    ,
    1 2
    2 1 2 0 3 3
    0 2
    0 0,
    0,
    0;
    0.
    a
    a a
    a a
    a
    a
     

     


       

    Для характеристического полинома четвертого порядка,
    4

    n
    ,
    4 3
    2 4
    3 2
    1 0
    4
    ( )
    0,
    0
    D p
    a p
    a p
    a p
    a p a
    a






     . матрица Гурвица имеет вид
    3 1
    4 2
    0 3
    1 4
    2 0
    0 0
    0 0
    0 0
    a
    a
    a
    a
    a
    a
    a
    a
    a
    a
     
    ,
    1 3
    2 2
    1 1 4 2
    2 3
    1 2
    0 3 1
    2 3 1 4 0 3 4
    0 3
    0 0,
    0,
    0; (
    )
    0,
    0;
    0.
    a
    a
    a a
    a
    a a
    a a a
    a a
    a a
    a
    a
     

       

       




       

    И так далее.
    Из структуры построения определителей Гурвица следует, что
    0 1
    0
    n
    n
    a

      
     , т.е. достаточно, чтобы
    0 0
    a

    Замечание

    Система находится на границе устойчивости, если
    0 1
    0





    n
    n
    a

    Система находится на границе апериодической устойчивости, если
    0 0
    a
     .

    Система находится на границе колебательной устойчивости, если характеристический полином содержит пару чисто мнимых корней, чаще всего
    1 0
    n


     .

    Система находится на границе апериодической устойчивости, если характеристический полином
    )
    ( p
    D
    содержит бесконечный корень. Действительно, если всё уравнение (5.3) разделить на
    n

    , то получим
    -1 2
    1 0
    2 1
    1 1
    1 1
    ( )
    0
    n
    n
    n
    n
    n
    D p
    a
    a
    a
    a
    a








    


     .

    5
    Откуда видно, что при
    0
    n
    a
     имеем
    0 1 

    , а значит



    Пример.
    Передаточная функция разомкнутой части системы имеет вид
    3 2
    1 2
    1
    ,
    )
    1
    )(
    1
    (
    )
    (
    k
    k
    k
    K
    p
    T
    p
    T
    p
    K
    p
    W




    Характеристический полином замкнутой системы соответственно будет
    0
    )
    (
    )
    (
    2 2
    1 3
    2 1






    K
    p
    p
    T
    T
    p
    T
    T
    p
    D
    Коэффициенты его положительны. Условие устойчивости по критерию Гурвица (5.5) получат вид
    2 1
    2 1
    )
    (
    T
    KT
    T
    T


    или
    


    




    2 1
    1 1
    T
    T
    K
    .(5.6)
    Границы устойчивости
    1.
    0 0;
    0,
    a
    K


    2.
    2 1
    2 1
    1 0;
    ,
    ãð
    K
    T
    T
     


    3.
    3 1 2 0,
    0
    a
    T T

     .
    Эти три границы устойчивости можно изобразить графически в пространстве параметров
    2 1
    ,
    ,
    T
    T
    K
    и найти области устойчивости системы.
    Определим вначале область устойчивости системы по одному параметру K (общий коэффициент усиления разомкнутой части). Пространство параметров здесь одна прямая линия, а границы устойчивости – точки на ней:
    0

    K
    и
    гр
    K
    K

    , как представлено на рис.5.5.
    2 1
    1 1
    T
    T

    0
    гр
    K
    Рис.5.5. Область устойчивости по параметру K
    Те же границы устойчивости можно построить на плоскости двух параметров, например:
    1
    ,T
    K
    (рис.5.6).
    Первая граница
    )
    0
    (

    K
    лежит на оси
    0 1

    T
    , как показано на
    Вторая граница
    1 2
    1 1
    K
    T
    T


    имеет вид гиперболы с асимптотами,
    0 1

    T
    и
    2 1
    K
    T

    Третья граница совпадает с осью K .
    0 2
    1
    T
    1
    T
    K
    Рис.5.6. Область устойчивости по параметрам
    1
    ,T
    K

    6
    Как видно (рис.5.6), при увеличении постоянных времени
    1 2
    ,
    T T область устойчивости сужается. Отрицательно влияет на устойчивость также и увеличение коэффициента усиления K .
    При любых заданиях
    1
    T
    и
    2
    T
    существует свое граничное значение коэффициента усиления
    гр
    K , после которого система становится неустойчивой. Это является важным обстоятельством, так как для повышения точности работы системы необходимо увеличивать добротность системы, то есть K . Тут выявляется противоречие между требованиями точности и устойчивости.
    5.4.Частотные критерии устойчивости
    5.4.1
    .Критерий устойчивости Михайлова
    Рассмотрим характеристический полином замкнутой системы
    n
    -порядка (5.3) с положительными коэффициентами (необходимое условие устойчивости). Подставив в него чисто мнимое значение корня

    j
    p

    , получим
    -1 2
    -1 2
    1 0
    (
    )
    (
    )
    (
    )
    (
    )
    (
    )
    ( )
    ( )
    n
    n
    n
    n
    D j
    a j
    a
    j
    a j
    a j
    a
    U
    jV









     




    , (5.7) где
    2 0
    2 3
    1 3
    Re( )
    -
    ,
    Im( )
    -
    a
    a
    a
    a












    


    (5.8)
    Изобразим годографы этого выражения на комплексной плоскости (Re, Im) .
    0
    Re
    Im
    0
    a
    0


    2

    n
    3

    n
    4

    n
    5

    n












    )
    (

    j
    D
    Рис.5.7. Годографы Михайлова
    Анализируя годографы, можно сделать следующие выводы.

    При
    0


    имеем,
    0
    Re( ) a

     , Im( ) 0

     ,

    При



    будет Re( )

      или Re( )

      , Im( )

      или Im( )

      .
    Знаки
    
    или 
     зависят от показателя степени
    n
    . Из формулы (5.8), где все
    0
    i
    a
     , видно, что при



    для
    3

    n
    будет Re( )

      и Im( )

      , а для
    5

    n
    получим Re( )

      и
    Im( )

      и.т.д. Поэтому годографы эти имеют для различных
    n
    примерно такие формы как показано на рис.5.7. Эти годографы называются годографами Михайлова.
    Формулировка критерия устойчивости Михайлова
    . Для устойчивости линейной системы
    n
    -порядка необходимо и достаточно, чтобы изменение аргумента функции
    )
    (

    j
    D
    при
    изменении частоты




    0
    равнялось бы
    2

    n , то есть
    arg (
    )
    2
    D j
    n




    при




    0
    . (5.9)
    Другими словами, требуется, чтобы кривая Михайлова проходила последовательно
    n
    - квадрантов против часовой стрелки (окружая начало координат).

    7
    Кривые на рис.5.7. соответствуют устойчивым системам, а на рис.5.8. приведен пример годографа Михайлова, соответствующий неустойчивой системе пятого порядка.
    0
    Re
    Im
    0


    5

    n



    )
    (

    j
    D
    Рис.5.8. Годограф Михайлова, соответствующий неустойчивой системе пятого порядка
    Критерий Михайлова приведен без доказательства.
    Рассмотрим пример определения границ устойчивости системы по критерию Михайлова.
    Пусть как в предыдущем примере подразд.5.3. характеристический полином имеет вид
    0
    )
    (
    )
    (
    2 2
    1 3
    2 1






    K
    p
    p
    T
    T
    p
    T
    T
    p
    D
    Тогда, выполнив замену
    p
    j


    в ( )
    D p , получим (
    ) Re( )
    Im( )
    D j
    j





    . Действительная и мнимая части
    (
    )
    D j

    имеют вид
    2 1
    2 3
    1 2
    Re( )
    (
    )
    ,
    Im( )
    K
    T T
    T T











     
    
    Отсюда для границы устойчивости согласно выражению
    0 0
    Re( ) 0,
    Im( ) 0,







    (5.10) имеем
    2 0
    1 2
    0 2
    0 0
    1 2 0
    Re( )
    (
    )
    0,
    Im( )
    (1
    ).
    K
    T
    T
    T T














    
    Из второго уравнения получим два значения
    0 0


    и
    2 1
    2 0
    1
    T
    T


    Тогда из первого уравнения находим соответственно
    0

    K
    и
    2 1
    1 1
    T
    T
    K


    Для бесконечно удаленного корня
    )
    (
    0



    получим
    0 2
    1

    T
    T
    Эти результаты совпадают с тем, что было сделано в примере по критерию Гурвица, где были изображены области устойчивости, которые получаются такими же (см.рис.5.5, 5.6.) и по критерию Михайлова.
    Достоинством критерия Михайлова является то, что он проще и нагляднее в применении, особенно, если возможно (удастся) снять экспериментальные кривые замкнутой системы.

    1
    1   2   3   4   5   6   7   8   9


    написать администратору сайта