Файл: Шахнович, А. Р. Математические методы в исследовании биологических систем регулирования.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 90
Скачиваний: 0
V
оптимальной задачи определяются оптимальное управление, при надлежащее ограниченному пространству, и оптимальная траек тория.
Пространство Q может быть разделено гиперповерхностью на два подпространства (Е' и Е"), каждому из которых принадлежат точки (траектории), обладающие определенными свойствами:
Рассмотренные выше математические модели динамических си стем принадлежат к детерминированному классу, процессы в ко тором характерны тем, что знание их в определенном интервале позволяет определить поведение этих процессов вне этого интер вала.
К классу стохастических принадлежат такие модели динамиче ских систем, процессы в которых характерны тем, что знание их на некотором интервале позволяет определить вероятностные ха
рактеристики поведения этих процессов вне этого |
интервала. |
||||||||||
Эти модели связаны со стохастическими процессами, |
определяе |
||||||||||
мыми случайной |
функцией, |
значение |
которой в каждый |
момент |
|||||||
времени является случайной |
величиной. |
|
|
|
|
||||||
Случайное событие |
а |
характеризуется |
вероятностью |
р (а) |
|||||||
(О ^ |
р (а) |
1). |
Случайная |
величина |
характеризуется |
множе |
|||||
ством ее возможных значений хг ... хп |
п их вероятностями, -т. е. |
||||||||||
задаются п вероятностей рк |
— р (arj), |
где рі — вероятность |
слу |
||||||||
чайного события, |
заключающегося в появлении значения хі |
слу |
|||||||||
чайной величины £. |
|
|
|
|
|
|
|
|
|||
Другой характеристикой случайной величины является функ |
|||||||||||
ция |
распределения F (х), |
т. е. вероятность |
случайного |
события |
|||||||
£ < |
X, заключающегося в том, что величина \ меньше |
некоторого |
|||||||||
фиксированного |
уровня: |
|
|
|
|
|
|
|
|||
|
|
|
F |
(х) |
= |
р (£ <х). |
|
|
(1-0-14) |
Если F (х) — непрерывная и дифференцируемая функция на ин тервале — оо <^ X <[ со, то I — непрерывная случайная вели чина и
Р{х) = ^ - |
(1-0-15) |
— плотность вероятности.
Отметим из характеристик случайной величины еще математи ческое ожидание (среднее значение):
|
+00 |
|
М{1}= |
jj xP(x)dx. |
(1-0-16) |
|
—оо |
|
Введенное выше понятие вероятности является так называемой априорной вероятностью, так как определяет вероятность события до проведения опыта. Апостериорная вероятность (т. е. вероят ность, определяемая по результатам опыта), вообще говоря, иѳ равна априорной.
9
Апостериорная плотность вероятности стремится к истинному значению параметра а. Апостериорная плотность вероятности, являясь информационным показателем опыта, представляет собой условную плотность вероятности Р (АIX) при условии, что на блюдаемый вектор X задан. Определяется апостериорная плот ность вероятности с помощью формулы Бэйеса:
|
|
|
Р {АІХ) |
= |
РіЛ)Р{ХІЛ) |
_ |
|
|
|
|
|
|
|
|
^ |
P{A)P(X/A)dQ(A) |
|
|
|
|
|
|
|
|
П(А) |
|
|
|
|
Р (А) |
— априорная |
плотность |
распределения |
вектора |
А. |
||||
Если |
Р |
(A) |
dQ (А) — безусловная вероятность нахождения |
век |
|||||
тора |
А |
в объеме dQ (А), то Р |
{XIА) — условная |
вероятность на |
|||||
хождения |
вектора |
X |
при фиксированном |
А. |
|
|
|||
Функция Р (Х/А) |
|
называется функцией |
правдоподобия. Фор |
мула Бэйеса дает возможность определить апостериорные веро ятности значений А, если известны априорные вероятности и век тор X, определяемый в результате опыта.
Отдельные наблюдения над случайным процессом £ (t), про текающим в контролируемых условиях опыта, дают каждый раз различные функции х (t) — реализации случайного процесса.
Одним из распространенных является марковский случайный процесс, для которого поведение (вероятностные характеристики для будущих моментов времени) определяется состоянием в дан ный момент и не зависит от «предыстории».
Динамические задачи для стохастических моделей могут ре шаться с помощью условных марковских процессов, причем диф ференциальное уравнение записывается, например, в виде
dxt = butdt |
+ dr\u |
(1-0-18) |
где X, — координата пространства |
состояний; |
ut — управление; |
т}< — шум с определенными статистическими |
свойствами. |
|
Рассмотренные выше классы |
дифференциальных уравнений |
не представляют всего многообразия моделей управляемых систем. Многие классы моделей управляемых систем имеют дискретновероятностные представления (связанные с нейронными сетями, конечными автоматами и т. д.), которые рассмотрены ниже в со ответствующих главах.
Л И Т Е Р А Т У Р А
Беллмап Р. Введение в теорию матриц. М., «Наука», 1969.
Беллман Р., Кук К. Дпфференцпально-разностные уравнения. М., «Мир», 1967.
Буткоеский А. Г. Теория оптимального управления системами с распреде ленными параметрами. М., «Наука», 1965.
Заде Л., Дезоер Ч. Теория линейных систем. М., «Наука», 1970.
10
Петровский |
И. Г. Лекции по теории дифференциальных уравнений в частных |
производных. М., 1947. |
|
Понтрягин |
Л. С. Обыкновенные дифференциальные уравнения. М., 1961. |
Стратонович |
Р. Л. Условные марковские процессы и их применение к тео |
рии оптимального управления. М., 1966. |
Фелъдбаум А. А. Основы теории оптимальных автоматических систем. М., 1963.
Ципкин Я. 3. Адаптация и обучение в автоматических системах. М., «Нау ка», 1969.
Глава 1-1
ИДЕНТИФИКАЦИЯ
Одной из основных проблем, возникающих при исследовании управляемых систем, является определение характеристик объек та и приложенных к нему возмущений — проблема идентификации.
В этой проблеме возникают две задачи: определение структуры и параметров объекта и определение параметров объекта при изве
стной |
структуре. |
|
|
В первой из этих задач исследуется «черный ящик», и она имеет |
|||
более |
теоретическое значение. |
|
|
Вторая задача имеет значительные практические применения. |
|||
Методы, применяемые при идентификации, |
весьма разнообразны |
||
и определяются |
классом объектов. Методам, применяемым при |
||
идентификации, |
посвящена значительная |
литература. Одну из |
задач идентификации для линейного объекта рассмотрел Браверман (1966). Более общие проблемы рассмотрели Цыпкин (1968), Райбман (1967) и др.
Одним из методов определения характеристик объекта являет ся частотный метод. Как известно, при этом на вход объекта подается гармоническое возмущение
X (t) — X sin at,
где А (г- амплитуда; ш = 2л/ Г — частота; Т — период колебания. Выходная переменная объекта представляется также в виде
гармонических колебаний |
у (і) |
= |
Y siu(wt + <р), где Y |
— ампли |
|
туда; ф — фазовый сдвиг. |
|
|
|
|
|
Этот метод основан на двойственном представлении |
линейного |
||||
объекта в виде дифференциального |
уравнения |
|
|||
n |
• |
|
m |
. |
|
• d |
У(1) |
—Wh |
dx (x) |
|
|
2 я :•' |
dt- |
- |
^Di~~dt~ |
|
|
3=0 |
•> |
|
1=0 |
|
|
11
или в виде передаточной |
функции |
|
|
|
m |
|
|
2 W " |
W (ja) |
= |
i=o |
Ä (/od) |
о
которые связаны между собой с помощью прямого и обратного пре образования Лапласа.
Наиболее простой путь определения частотной характеристики физического или биологического объекта состоит в измерении установившейся реакции на синусоидальный входной сигнал для различных частот. При этом могут быть получены амплитудночастотная характеристика YIX (со) и фазочастотная характери стика ф (а>), характеризующие линейный объект.
Этот метод имеет ограниченную область применения (особенно для действующих систем).
В общем случае удобнее находить частотную характеристику путем анализа неустановившейся реакции объекта на произволь ное возмущение, которое может иметь форму единичной ступени или импульса.
Можно считать, что неустановившееся колебание представляет собой сумму синусоидальных колебаний с различными амплиту дами, перекрывающих весь частотный диапазон. Поэтому реакция линейной системы на неустановившийся входной сигнал может
рассматриваться как реакция этой системы на сумму |
синусоидаль |
|
ных колебаний. |
|
|
Метод преобразования данных переходного процесса объекта |
||
в данные, выра?кенпые в функции частоты, |
основан на использо |
|
вании интеграла Фурье |
|
|
оо |
|
|
/ (t) = О, |
при t > |
0. |
о |
|
|
Этот интеграл должен вычисляться для каждой частоты в пре делах времени от нуля до бесконечности. Интегрирование может
быть выполнено |
только при |
условии, что известно |
поведение |
/ (t) на бесконечном интервале |
времени после возмущения. Это |
||
означает, что сист |
ма должна достичь установившегося |
состояния |
в течение некоторого конечного периода времени Т так, чтобы функция / (t) могла быть выражена аналитически на интервале времени от Т до бесконечности.
Таким образом, реакция объекта на входное возмущение мо жет быть представлена в виде
тоо
От
12