![](/user_photo/_userpic.png)
книги / Теория автоматического управления техническими системами
..pdfВектор помехи п (см. рис. 11.8) представляет собой реали зацию многомерного случайного процесса с конечными момен тами второго порядка. Поскольку z*=z-J-n, то z* независимо от характера вектора управления интерпретируется как случайный вектор. Будем считать, что р также имеет моменты второго по рядка, поэтому при идентификации можно использовать свой ства оценок параметров вектора b вместо его плотности распре деления вероятности. Такими свойствами могут быть:
1) несмещенность, когда для каждого k математическое ожидание вектора р совпадает с Ь, т. е.
М[р]=Ь;
2) состоятельность, если с ростом длины выборки k для лю бого 6>0
lim р(1р—Ь|>е}=0, Ач-оо
где р —вероятность события; 3) эффективность оценки (достижение наилучшей точности
или минимальной дисперсии результирующей оценки для всех несмещенных оценок).
Когда 1-е и 3-е из перечисленных свойств выполняются только при k-+oot то их называют асимптотическими.
Метод байесовских оценок. При использовании этого метода для решения задач идентификации объекта в распоряжении
разработчика имеется наибольшая априорная информация (т. е. полученная до проведения измерений), включающая следующее:
1) плотность распределения вероятности шума п. По ней можно определить условную плотность распределения вероят ности измерений z*, которая зависит от параметров объекта b
иуправления и и обозначается через p(z*/b);
2)плотность распределения вероятности параметров Ь, обо
значаемая через q{b);
3) потери, связанные с численной оценкой {}при истинном значении параметра Ь. Эта функция штрафа, или потерь с(р, Ь),
имеет минимум при |3=Ь.
При использовании байесовских оценок исходим из формулы
Байеса: |
|
p(b/z*)=piz* |
(П.45) |
где р(b/z*)—апостериорная3 плотность распределения вероят
ности |
вектора b при заданных результатах |
измерения z*; |
p(z*) |
—плотность распределения вероятности |
z*. На основе |
информации о p(b/z*) требуется определить вектор Ь. Опреде лим условный риск выбора p(z*) оценки при истинном значе нии вектора параметров b как математическое ожидание функ-
3 Апостериорной называют информацию, полученнуюпосле выполнения измерений.
ции потерь С(Р, Ь) по наблюдениям z?, а средний |
риск /?(Р) |
||
как |
математическое ожидание условного |
риска по распре |
|
делению значений параметра объекта Ь: |
|
|
|
%{$)=Мь\Мг*/Ъ\с®, Ь)] = |
|
|
|
|
с(р, b)p(z*/b)q{b)dkz*dm"b, |
|
|
где |
\ означает (w + 1)-кратный интеграл |
S - S ' |
|
|
ш+1 |
m+1 |
|
dm+\b—db0dbJ ... dbm. |
|
||
|
|
||
Оценка, минимизирующая выражение #(р), |
называется |
оценкой минимального риска. Поскольку p(z*)^0, то средний риск можно минимизировать, сделав внутренний интеграл как
можно меньше при z*=c. Необходимое условие |
минимума |
#((}) имеет вид: |
|
ц $ с(р.Ъ) р(b/с)d^'b |м =0. |
(11.46) |
т+\ |
|
Получаемая из выражений (11.45) и (11.46) оценка р назы вается байесовской.
Метод максимального правдоподобия. Допустим, что функ ция штрафа с(р, Ь), условная плотность распределения вероят ности p(z*/b) и априорная плотность распределения вероят ности q{b) параметров b системы не известны. В этом случае <7(b)=Â=const в области возможных значений параметра Ь. Для произвольного z* имеем
maxp(b/z*)=^max/>(z*;b), |
(11.47) |
причем шах берется по bGR”*+1, где R,n+1—(т+1)-мерное ве
щественное арифметическое пространство. Так как параметр b не является случайной величиной, а представляет неизвестный постоянный параметр, то в выражении p(z*/b) правой части уравнения (11.47) вместо косой черты стоит точка с запятой.
Априорная |
плотность |
распределения |
вероятности |
имеет |
вид p{z*(l), |
..., z*(k), b). Апостериори |
становится известной |
||
реализация выборочных |
значений z*(l)=Cf, ...; z*{k)—ch |
|||
или z*=c. |
плотность |
распределения |
вероятности, |
связы |
Обозначим |
вающую реализацию выборочных значений с вектором р оцен ки параметров системы, через
Цси .... Ск; р}.
Это выражение называется функцией правдоподобия. Выберем такое значение рт, которое максимизирует L. Необходимое условие этого максимума:
P)lp=t—0. |
(11.48) |
|
В выражении .(11.48) вместо L удобнее |
рассматривать |
InL, |
так как из-за монотонности логарифма |
максимумы L |
и InL |
достигаются при одном и том же значении р, определяемом из
£taL(C5» 1,^-0. |
(11.49) |
Равенство (11.49) называют уравнением правдоподобия. Отыс кивая решение этого уравнения, обеспечивающее наибольшее
значение L(c; р) или InL(с; р), определим оценку максималь
ного правдоподобия р.
Перечислим некоторые асимптотические свойства этой оценки:
1) нормальность, т. е. р(р/Ь) при k-+co приближается к нормальному распределению;
2)несмещенность, т. е. Af[p]-*b при £-»-оо;
3)эффективность, т. е. стремление к большей точности или
минимуму дисперсии; 4) состоятельность.
Метод наименьших квадратов. Не требует никакой априор ной информации о статистических характеристиках рассматри
ваемых процессов. Результаты наблюдений в отличие от преды дущих представляют векторами-столбцами. Оценку вектора параметров ищут в классе линейных несмещенных оценок. Не обходимое и достаточное условие минимума функции от ошиб ки /(р) определяют системой нормальных уравнений. Оценку,
минимизирующую /(р) по |
p6Rm'H при R#E (где |
Е —единич |
||
ная матрица порядка k)t |
называют оценкой метода наимень |
|||
ших квадратов с весовой |
матрицей |
R. Полагая R=E и мини |
||
мизируя /(р) по p6Rm+I, |
получают |
классический |
метод наи |
меньших квадратов.
Метод марковских оценок, или обобщенный метод наимень ших квадратов. Предположим, что известна ковариационная матрица аддитивного шума, воздействующего на САУ:
МЩ\)п(\)]... М[л(1)й(Л)П М[га(*)«(!))... М[п(к)п(к\\\
Можно доказать, что при N#1 наилучшую линейную оценку получают минимизацией выражения
/=eTN_1e,
ззз
что приводит к следующему правилу оценивания: |
|
^[UTNMlJl-'UTN^z, |
|
где U —неизвестная kX(m+\)-матрица, которая строится |
с |
использованием элементов вектора и. |
эф |
Марковская оценка является линейной, несмещенной, |
фективной. В качестве критерия для оценки параметров, ко торая зависит от р, примем /(р). Полагая N=1 и минимизируя /(р) по р, получаем классический метод наименьших квадра
тов:
UW=Uz
или p=[UTL)]_1UTz.
Метод наименьших квадратов в детерминированной и ста
тистической интерпретации |
широко |
|
используют |
для оценки |
|
параметров в |
дифференциальных и |
в разностных |
уравнениях, |
||
а также для |
определения |
значения |
импульсной;' |
переходной |
функции.
Все рассмотренные методы идентификации объектов управ ления и САУ имеют соответствующие алгоритмическое и про граммное обеспечение.
Контрольные вопросы
1.Что такое идентификация?
2.Перечислите методы идентификации САУ.
3.Как производится идентификация методом частотных ха рактеристик?
4.В чем заключаются преимущества идентификации мето дом корреляционных функций?
5.В чем состоит ортогональный метод моментов?
6.Сформулируйте сущность метода дифференциальной аппроксимации?
7.В каких случаях целесообразно применять байесовские оценки?
8.В чем состоят оценки методом максимального правдопо добия?
12. ОСНОВНЫЕ ПРИНЦИПЫ ПОСТРОЕНИЯ АДАПТИВНЫХ СИСТЕМ
Адаптивные (приспосабливающиеся) системы автомати ческого управления (АСАУ) могут рассматриваться как си стемы с элементами искусственного интеллекта. Назначение
АСАУ состоит в том, чтобы заменить человека-оператора при принятии решений об улучшении характеристик системы в
процессе ее нормальной эксплуатации (такая необходимость может возникнуть, например, тогда, когда условия функциони рования системы не допускают участия оператора).
Первые работы по самоприспосабливающимся системам, относящимся к классу самонастраивающихся, появились в се редине 60-х годов1. АСАУ стали внедрять в различные отрас ли промышленности во второй половине 70-х годов, чему спо собствовало быстрое развитие вычислительной техники.
12.1. Функциональные особенности и классификации адаптивных САУ
Системы управления с неизменяемой в процессе эксплуата ции структурой, а также параметрами настроек, полученны ми на этапе проектирования, часто не могут обеспечить ка чественного и даже устойчивого управления объектами. Это связано с тем, что: во-первых, математическая модель САУ обычно не может быть точно задана на этапе разработки си стемы; во-вторых, САУ всегда действует в случайной среде, ее реальные характеристики априорно можно учесть лишь статис тически. Эти два фактора приводят к тому, что в процессе нор мальной эксплуатации параметры системы изменяются случай ным образом, а степень их неопределенности может быть раз личной. Случайные изменения параметров САУ в процессе ее функционирования и сами эти параметры называют неконтро лируемыми. Причины неконтролируемости изменений парамет ров САУ различны (например, из-за нестабильности источни ков питания, температуры и давления окружающей среды, естественного старения и пр.). Что касается контролируемых или желаемых изменений, то под ними понимают изменения статических и динамических свойств АСАУ, производимые ре гулятором на основе данных текущей информации о внешних и внутренних условиях работы системы. Контролируемые па раметры регулятора адаптивной САУ называют настраивае мыми.
Если диапазон изменения неконтролируемых параметров невелик, их разброс может быть уменьшен в процессе работы благодаря обратным связям в системе. Но если этот диапазон или степень неопределенности параметров значительны, то удовлетворительные технические характеристики могут быть получены лишь в результате самоприспособления системы к изменяющимся условиям в ходе ее нормальной эксплуатации. Примерами таких систем являются: АСАУ сверхзвукового са молета, характеристики которого существенным образом меня ются в зависимости от скоростного напора; система моделнро-1
1См. например: Солодовников В. В. Некоторые принципыпостроения н вопросытеории самонастраивающихся систем автоматического управления // Автоматическое регулирование и управление. М.: Изд-во АНСССР, 1956. С. 253—259. (Труды сессии АН СССР).
вания объекта высокой размерности с помощью модели невы
сокой размерности и др.
Приведем краткую классификацию приспосабливающихся систем, объединив их общим термином —«адаптивные систе мы». Итак, адаптивные САУ можно подразделить на (рис.
12.1):1
Рис. 12.1. Классификация адаптивных САУ
1) самонастраивающиеся (или системы с параметрической адаптацией). Используя текущую информацию о характерис тиках внешних воздействий или/и о динамических характерис тиках объекта, они осуществляют контролируемые изменения параметров регулятора, что обеспечивает улучшение качества функционирования системы;
2) самоорганизующиеся (или функционально адаптивные). Уменьшение априорной неопределенности, приводящее к улуч шению качества управления технологическим процессом или техническим объектом, достигается в них путем использования информации, получаемой в ходе непрерывного измерения дос тупных входных и выходных сигналов;
3) самообучающиеся. Информация о неизвестных характе ристиках процесса и окружающей среды вырабатывается ими в ходе обучения. Эта информация используется в дальнейшем для оценивания, классификации и принятия соответствующих
решений с целью существенного улучшения качества работы системы.
Теперь можно дать следующее общее определение приспо сабливающейся системы: система управления называется при
спосабливающейся, или адаптивной, если она обладает способ ностью уменьшать априорную неопределенность и повышать эффективность управления, используя информацию, получае муюв процессе нормальной (штатной) эксплуатации.
Необходимость применения адаптивных систем возникает в случае неполной информации о проектируемой системе, т. е. когда не известны: параметры объекта и точные характеристи ки их изменения во времени; внешние возмущения; задающие или управляющие воздействия.
В адаптивной САУ в процессе ее работы могут осущест вляться изменения контролируемых (настраиваемых) парамет ров, структуры системы, цели и критерия управления, что поз
воляет компенсировать в широком диапазоне изменение не контролируемых параметров и поддерживать технические ха
рактеристики системы на требуемом уровне.
12.2. Пассивные и активные, разомкнутые и замкнутые адаптивные системы
В зависимости от способа, которым осуществляют желае мые изменения, АСАУ можно подразделить на (см. рис. 12.1):
1) пассивные. В них контролируемые изменения програм
мно зависят от имеющейся на стадии проектирования АСАУ априорной информации о внешних и внутренних условиях ра боты системы;
2) активные. В «их осуществляются контролируемые из менения собственных характеристик в зависимости от текущей, а не только от априорной информации об условиях работы си
стемы управления.
Кроме того, АСАУ могутбыть (см. рис. 12.1): 1) с разомкнутым контуром самонастройки; 2) . с замкнутым контуром самонастройки.
Адаптивными системами с разомкнутым контуром (циклом) самонастройки, или, просто, разомкнутыми АСАУ, называют системы, в которых программа настройки параметров опреде ляется заранее и таким образом, чтобы система управления находилась в расчетном квазиоптимальном режиме при неко торых типовых или наиболее вероятных внешних и внутренних условиях. Для настройки системы во .всем диапазоне изменения этих условий необходим больший объем априорной информа
ции.
На рис. 12.2 приведена схема обычной, неприспособавливаю-
щейся системы стабилизации углового положения ЛА. При изменении условий полета меняется передаточная функция
![](/html/65386/197/html_n9WY3FfOLU.fTwN/htmlconvd-rfS78B338x1.jpg)
m s,(t) fn(t)
Рис. 12.2. ^приспосабливающаяся система |
стабилиза |
ции ЛА |
|
№о(5) ЛА, а следовательно, и динамическая |
характеристика |
всей системы стабилизации: |
|
■У(д) . Wo (s) W,(s) |
(12.1) |
U77Й“(s) 1+W,(s)W,(s)’ |
где U(s) и V(s)—лапласовские изображения входного и вы
ходного сигналов. |
..., fn(f), при- |
Возмущения со стороны внешней среды |
водящие к неконтролируемым изменениям параметров системы,
приложены к различным точкам регулятора и объекта.
На рис. 12.2 возмущающее воздействие f(t), приложенное непосредственно к входу объекта управления, в отличие от
•fi(0* •••* Ш) не меняет его параметров. Поэтому в процессе работы системы измеряют только f\(t)t .... hit) (см. рис. 12.3).
В соответствии с принципом обратной связи и выражением (12.1) неконтролируемые изменения характеристики U?0(5) из-
за возмущений и помех вызывают сравнительно небольшие изменения параметров Ф($). Если поставить задачу более пол ной компенсации неконтролируемых изменений, чтобы переда точная функция <D(s) системы стабилизации ЛА оставалась практически неизменной, то следует надлежащим образом из менить характеристику регулятора ^i(s). Это и осуществляется
в приспосабливающейся САУ, выполненной по схеме, приведен ной на рис. 12.3. Параметры внешней среды, характеризуемые
сигналами fi(t), ...» ЫО. например давление скоростного на пора рн(0* тмпература окружающего воздуха T°(t) и скорость
полета о (0, непрерывно измеряются датчиками Д1—ДЗ, и те кущие значения параметров поступают в вычислительные
устройства В1—ВЗ, вырабатывающие сигналы, с помощью ко торых подстраивается характеристика №i(s), чтобы компенси ровать изменения характеристики W0(s). Однако в АСАУ дан ного типа (с разомкнутым циклом настройки) отсутствует са моанализ эффективности осуществляемых ею контролируемых изменений.
Адаптивные системы с замкнутым циклом самонастройки,
или, просто, |
замкнутые АСАУ, обладают способностью осу |
ществлять такой самоанализ. С точки зрения структуры замк |
|
нутые АСАУ |
отличаются от разомкнутых наличием дополни |
тельного |
замкнутого контура самонастройки, |
служащего для |
анализа эффекта изменения качества системы |
управления в |
|
процессе |
адаптации. |
|
Приведем несколько примеров простых адаптивных систем.
ВАСАУс разомкнутым циклом (рис. 12.4) самонастройка проходит с
Рис. 12.4. АСАУс разомкнутым циклом на стройки
учетомизменений характеристик системыи входного сигнала. Текущая ин |
|||
формация об этих изменениях получается на основе анализа ошибки e(f). |
|||
Самонастройка |
выполняется дополнительным |
корректирующим устройством |
|
№„($), которое |
включается, наприМер, когда |
сигнал .ошибки контура управ |
|
ления превышает заданные пороговые значения. Эта система, следовательно, |
|||
обнаруживает изменения: или входного сигнала при наличии шумов; или па |
|||
раметров объекта; или те и другие одновременно. |
|||
К замкнутой АСАУследует отнести автоколебательнуюсистему с нели |
|||
нейным элементом (НЭ) в прямой цепи |
(рис. 12.5). (На этом рисунке и далее |
||
Ucc —сигнал самонастройки.) Система |
обладает способностьюизменять пе- |
редаточный коэффициент прямой цепи таким образом, чтобы обеспечить |
|
оптимальный процесс управления, мало зависящий от динамических свойств |
|
объекта. Замкнутый контур системыпрактически можно выполнить эквива |
|
лентным усилительному звену с передаточным коэффициентом, равным еди |
|
нице, а следовательно, инвариантным к изменениюдинамических характе |
|
ристик объекта; причем значение коэффициента усиления в прямой цепи это |
|
го контура должно быть наибольшим. Желаемая динамика системы управ |
|
ления обеспечивается характеристикой модели фильтра Ф.Эта модель вклю |
|
чается последовательно с замкнутым контуром управления. |
|
Чем больше значение передаточного коэффициента прямой цепи, тем в |
|
меньшей степени характер переходных процессов зависит от динамических |
|
свойств объекта. Однако увеличение значения передаточного коэффициента |
|
ограничено из-за автоколебаний или неустойчивого функционирования сис |
|
темы управления,в связи с чем возникает проблема контроля автоколебания. |
|
Возможны два пути управления автоколебаниями: 1) |
настройка передаточ |
ного коэффициента с цельюограничения амплитуды |
колебаний; 2) поддер |
Рис. 12.5. АСАУс замкнутый циклом самонастройки |
жание коэффициента усиления на уровне, соответствующем границе появле ния автоколебаний.
Структурная схема АСАУс настройкой коэффициента усиления для огра ничения амплитуды автоколебаний показана на рис. 12.6. К основному кон-
Рис. 12.6. АСАУс настройкой коэффициента усиления