Файл: Шахнович, А. Р. Математические методы в исследовании биологических систем регулирования.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