Математические методы и алгоритмы обработки информации при идентификации динамических систем тема диссертации и автореферата по ВАК РФ 05.13.01, доктор физико-математических наук Малевинский, Михаил Федорович

• Малевинский, Михаил Федорович
• доктор физико-математических наук
• 2004

Оглавление диссертации доктор физико-математических наук Малевинский, Михаил Федорович

Введение

Глава 1. Показатели идентификации динамических систем.

1.1. Основные показатели контроля текущего состояния системы. Принципы оценки.

1.2. Принцип прогнозирования.

1.3. Постановка задачи оценки показателей идентификации динамических систем.

1.4. Принципы представления результатов идентификации на дисплее ПЭВМ.

Глава 2. Метод синтеза тестовых сигналов для идентификации динамических систем.

2.1. Постановка задачи синтеза тестовых сигналов в базисе ВВСФ.

2.2. Метод вычисления ВВСФ на основе полной проблемы собственных значений для симметрической матрицы.

2.3. Вычисление ВВСФ методом коллокации с использованием сплайнов произвольной степени

2.4. Вычисление ВВСФ с использованием эрмитовых сплайнов третьей степени.

2.5. Метод вычисления ВВСФ с использованием В-сплайнов.

2.6. Метод вычисления ВВСФ на основе ряда Котельникова.

Глава 3. Алгоритмы синтеза и характеристики тестового сигнала для идентификации динамических систем.

3.1. Алгоритм синтеза тестового сигнала.

3.2. Характеристики тестового сигнала, представимого в базисе ВВСФ

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

4.1. Оценка параметров весовой функции методом наименьших квадратов.

4.2. Оценка параметров весовой функции методом минимаксных квадратов.

4.3. Оценка параметров весовой функции при априорной неопределенности относительно значений корреляционной функции ошибок измерений.

4.4. Оценка параметров весовой функции системы, представимой нелинейным разностным уравнением.

4.5. Алгоритмы прогнозирования весовой функции системы.

Глава 5. Методы оценки весовых функций при идентификации систем с расределенными параметрами.

5.1. Оценка параметров весовой функции системы двумерным оператором статистического сглаживания.

5.2. Оценка параметров весовой функции системы байесовским оператором сглаживания.

5.3. Оценка параметров весовой функции методом максимума апостериорной вероятности.

5.4. Оценка параметров весовой функции статистически линеаризованным методом максимума апостериорной вероятности.

5.5. Оценка параметров весовых функций распределенных систем методом наименьших квадратов.

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

5.7. Оценка параметров весовой функции методом максиминного правдоподобия

5.8. Оценка весовой функции, основанная на применении линейного согласованного фильтра.

5.9. Оценка весовой функции нелинейным согласованным фильтром.

5.10. Оценка весовой функции системы нелинейным оператором Гаммерштейна.

5.11. Вероятностные характеристики параметров весовой функции системы, оцениваемые двумерным оператором статистического сглаживания.

5.12. Вероятностные характеристики параметров весовой функции системы, оцениваемые байесовским оператором сглаживания.

Введение диссертации (часть автореферата) на тему «Математические методы и алгоритмы обработки информации при идентификации динамических систем»

Одним из эффективных способов под держания сложных динамических систем в требуемом состоянии является совершенствование и развитие алгоритмического диагностического контроля показателей их функционирования (состояния). Под показателями контроля состояния динамических систем понимаются отклонения их основных-доминирующих и частных характеристик, оцененных (вычисленных) в конкретных текущих условиях от требуемых их значений, установленных для нормального штатного режима. Такие показатели характеризуют устойчивость, как необходимое качество систем. Они должны не только оцениваться, но и прогнозироваться непосредственно в процессе их функционирования в различных реальных условиях: детерминированных и/или при воздействии внешних и внутренних случайных факторов в том числе и при имитации ситуаций и объектов, на которые система должна реагировать в соответствии со своим назначением.

В диссертации показатели состояния оцениваются и прогнозируются для динамических систем, аксиоматически определяемых [1] множествами моментов времени -Т, состояний — X, мгновенных входных воздействий - U, допустимых входных воздействий - С1 = {со :Т С/} * 0, ПсГ, Q- замкнутое множество, множеством мгновенных значений выходных величин -Y, множеством выходных величин Г = {у: Т У}. Принципиальным при этом является то, что перед каждой системой ставится определенная цель, то есть контролю по показателям состояния подлежат целенаправленные динамические системы. В этих системах воздействия а> = {a>(t),a(t)} обусловливаются соответственно состояниями внешней среды и внутрисистемными факторами - параметрами a(t) = {aI(0,fl2(0v.fl/(0} • Текущее состояние системы определяется по соотношению x(t) = Т(/,/0,х(/0),&>(/), д(0), где t е Т, t > t0, t0 - начальный момент времени, x(t)~ траектория поведения системы на множестве Xпри {co(t),a(t)} е Q, rj(t) - случайные возмущения, ЧР - переходная функция состояния на X. ЧР-это основная характеристика системы, но она не наблюдаема в текущих условиях, так как априори принципиально невозможно полностью описать все множество альтернативных состояний системы и связей между ними, порождаемые возможными отклонениями параметров от их номинальных значений в процессе функционирования системы. Известно только, что для режима нормального функционирования должно иметь место выполнение условия а^ < a(t) < а^, иначе будет иметь место нештатный режим: "ВНИМАНИЕ " или "БРАК". При известных же x(t), co(t,tQ) можно было бы однозначно установить текущее состояние системы в зависимости от параметров a(t) = {a,^),^/),.,^/)} и при этом спрогнозировать состояние системы на момент tk>t>t0eT, а значит и установить выходные величины y(t) = (p(t,x(t),a(t)\co{t,t0)) как на текущий, так и на прогнозируемый моменты времени; ^-выходное отображение - оператор выхода. Но поскольку и X неизвестны, то текущее состояние системы можно оценить только по результату её реакции на известные входные воздействия. Так, если входное воздействие задать в виде одиночного кратковременного (широкополосного) сигнала, то оператор выхода будет представляться интегральным оператором, например, в виде оператора Дюамеля, Гаммерштейна, Вольтерра. Ядра этих операторов описывают полностью весовую (импульсную) функцию системы при конкретных, но априори неизвестных, её параметрах.

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

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

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

Далее всюду будем исходить из учета зависимости y(t) = v(tMv0,x(t0),a(t)Mt),ri(t)), как уравнения наблюдения, которое представляется в пространстве сигналов -в пространстве входных и выходных величин при нахождении системы в состоянии x(t), характеризующимся вектором параметров a{t) ^ {ax{t),a2{t),.,a,(t)}. Из этой зависимости непосредственно следует возможность определения при контроле системы компонент вектора a(t) = {a^tXa^t),.,^^)} или, иначе, возможность определения параметров весовой функции системы и её частных характеристик при известных измеренных данных y(t) и априорных данных относительно x(tQ), co(t), ?j(t), что в целом составляет содержание задачи идентификации системы, причем идентификации как параметрической, так и непараметрической.

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

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

Очевидно при контроле - идентификации системы входные воздействия должны представлять смесь полезного сигнала &>(/) и помех rj(t), формироваться специальным образом и в том числе с модуляцией их параметров, чтобы удовлетворить условию сочленения входных воздействий [1,2], что возможно с использованием имитаторов.

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

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

Актуальность темы. К настоящему времени по проблемам теории и практики идентификации систем опубликовано достаточно большое число работ. Результаты этих работ, накопленные к 1984—1987 гг., изложены, например, в книгах Я. 3. Цыпкина, «Основы информационной теории идентификации» (М.: Наука, 1984, 320 е.), А.

A. Бессонова, Ю.В.Загашвили, А. С. Маркелова «Методы и средства идентификации динамических объектов» (JL: Энергоатомиздат, 1984,280 е.),

B.Я. Катковника «Непараметрическая идентификация» (М.: Наука, 1985, 336с), Ш.Е. Штейнберга «Идентификация в системах управления» (М.: Энергоатомиздат, 1987, 80 с.) и в «Справочнике по теории автоматического управления» под ред. А. А. Красовского (М.: Наука, 1987, 712 е.). В последующие годы была издана книга JI. Льюнга «Идентификация систем—теория для пользователя» (М.: Наука, 1991, 432 е.), перевод с английского под редакцией Я. 3. Цыпкина и опубликовано большое количество статей в различных научных журналах в нашей стране и за рубежом (их краткое содержание изложено в реферативном журнале ВИНИТИ за 1988—2003 гг. «Техническая кибернетика» в разделе «Теория кибернетических систем управления»).

В этих работах такие задачи идентификации динамических систем как формирование испытательного сигнала (ИС), оценка характеристик и параметров систем рассматриваются как независимые самостоятельные частные задачи.

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

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

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

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

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

Вторым ключевым направлением при идентификации систем является оценивание их операторов по выходным данным при подаче на вход системы испытательных сигналов. Для линейных систем (непрерывных и дискретных ) в качестве характеристики оператора, представляющего систему, принимается ее весовая функция ( импульсная переходная функция). Для нелинейных систем, представимых операторами Вольтерра или Гаммерштейна, в качестве весовой функции принимается набор их ядер.

Обычно в качестве критерия идентификации принимается минимум среднего значения квадрата (дисперсии) ошибки между выходными значениями истинного неизвестного оператора реальной (идентифицируемой) системы и искомого оператора системы. Такой критерий справедлив в тех случаях , когда случайные входные воздействия распределены по нормальному закону с известными параметрами. Однако при статистической непараметрической идентификации, когда входные случайные воздействия должны формироваться с заданными математическими ожиданиями и корреляционными функциями, то для этого потребуется большое количество реализаций. Поэтому вероятностные характеристики входного процесса реализуются с некоторыми погрешностями, границы которых можно оценить методами математической статистики. Кроме того, на практике, как правило, неизвестны значения корреляционных функций входных сигналов, а известны только границы их изменения.

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

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

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

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

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

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

Таким образом возникает необходимость комплексного подхода к решению проблемы идентификации, как проблемы синтеза оптимального испытательного сигнала, разработки методов и алгоритмов определения оптимальных операторов линейных и нелинейных динамических систем, разработки программно-реализуемых на ПЭВМ алгоритмов, обеспечивающих достоверный и своевременный контроль состояния динамических систем в текущих условиях их функционирования.

Цель работы. Разработка математических методов, вычислительных алгоритмов и программ решения комплекса задач идентификации: оптимизации испытательного сигнала и формирования случайных входных воздействий; восстановления операторов линейных и нелинейных систем, наилучшим образом по заданным критериям, аппроксимирующих истинные (реальные) неизвестные операторы идентифицируемых систем; оценки параметров весовых функций линейных систем и ядер операторов Вольтерра или Гаммерштейна для нелинейных систем (параметрическая идентификация ); оценки вероятностных характеристик систем, представимых системами обыкновенных дифференциальных уравнений, линейных и нелинейных систем, представимых операторами Вольтерра и Гаммерштейна; восстановления частотных характеристик систем (сигналов) по конечной временной выборке, содержащей случайные ошибки; создание библиотеки процедур и моделей для идентификации систем с сосредоточенными и распределенными параметрами.

На защиту выносятся следующие концепции- положения:

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

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

В основу восстановления оператора нелинейной системы принята концепция параметрического оценивания ядер полинома Вольтерра для систем с сосредоточенными параметрами и непараметрического оценивания ядер полинома Гаммерштейна для систем с распределенными параметрами, для линейных и статистически линеаризованных по методу Казакова- Бутона систем- концепция представления их байесовскими фильтрами с конечной памятью, учитывающими априорную информацию о параметрах нестационарного случайного полезного входного сигнала.

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

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

4.Концепция параметрического оценивания операторов нелинейных динамических систем, имеющих нелинейности с разрывами непрерывности производных нулевого и первого порядков, основанная на статистически линеаризованном методе максимума апостериорной вероятности.

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

6.Концепция аналитического высокоточного вычисления частотных характеристик динамических систем на основе решения интегрального уравнения Фредгольма первого рода с сильно осциллирующим ядром при их аппроксимации сплайнами Лагранжа, ВВСФ, обобщенным рядом Котельникова.

7.Концепция построения классов математических моделей в системе объектно-ориентированного программирования "Дельфи " на основе разработанной библиотеки процедур и моделей для исследования методов идентификации систем с сосредоточенными и распределенными параметрами.

Научная новизна диссертации состоит:

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

Собственно синтез осуществляется по критерию минимизации среднеквадра-тической ошибки приближения единичного спектра линейной комбинацией ВВСФ с учетом требований по заданной энергии и согласованности с полосой пропускания частот идентифицируемой динамической системы.

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

2.В синтезировании оптимальных операторов линейных динамических систем с сосредоточенными и распределенными параметрами при априорной неопределенности о пространственно-временных и частотных характеристиках внешних воздействий как условий функционирования, вызывающих дополнительные ошибки при формировании выборочных данных измерительными средствами идентифицируемой системы.

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

Это теоретическое положение охватывает имеющееся решение 1 / аналогичной задачи, для случая, когда коэффициенты разложения имеют бесконечные дисперсии и нулевые взаимные моменты. б) оптимальности метода вычисления эффективных оценок параметров весовой функции по критерию максимума правдоподобия ( при нормальном законе распределения аддитивных ошибок измерений) посредством условной максимизации квадратичной формы с положительно определенной матрицей, при условии, что ошибки измерений ограничиваются заданными пределами второй составляющей ошибок, и последующем сведением задачи максимизации к полной проблеме собственных значений, легко реализуемой на ПЭВМ. Существующие же методы решения таких задач являются весьма трудоемкими для реализации на ПЭВМ. в) необходимого условия оптимальности весовой функции линейной системы с распределенными параметрами как двумерного согласованного фильтра. Доказательство основано на принципе построения согласованного фильтра по критерию максимума отношения сигнал/шум на выходе системы, приводящего к формированию двумерного интегрального уравнения Фредгольма первого рода, решением которого является искомая оптимальная весовая функция. Установленное таким образом необходимое условие в форме интегрального уравнения является общим по отношению к существующим методам определения весовых функций линейных систем как согласованных фильтров.

3. В синтезе оптимальных операторов нелинейных динамических систем с распределенными параметрами как нелинейных согласованных фильтров-обнаружителей тестовых сигналов.

Синтез основан на реализации критерия максимального отношения сигнал/шум на выходе динамической системы, представимой оператором Гам-мерштейна n-го порядка, и приводит к необходимым условиям оптимальности в виде системы двумерных интегральных уравнений Фредгольма первого рода, которой должны удовлетворять ядра Гаммерштейна.

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

4. В синтезе оптимальных операторов нелинейных динамических систем, как сглаживающих фильтров, представимых двумерным оператором Гаммерштейна п-го порядка, по критерию минимума среднего значения квадрата ошибки воспроизведения требуемого выходного сигнала идентифицируемой системы. Необходимые условия оптимальности оператора получены в форме системы п-го порядка двумерных интегральных уравнений Фредгольма первого рода. При этом уравнение Винера- Хопфа является частным случаем полученной системы при представлении динамической системы полиномом Гаммерштейна первой степени.

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

Замена полинома Вольтерра полиномом Гаммерштейна для таких систем позволяет проводить их синтез и вероятностный анализ без необходимости вычисления многомерных интегралов повышенной кратности.

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

Выведены аналитические выражения авто и взаимных корреляционных функций ошибок фильтрации для оптимальных линейных дискретных байесовских фильтров с конечной памятью. В байесовском фильтре с нарастающей памятью (рекуррентном фильтре Калмана ) вычисляется ковариационная матрица оцениваемых параметров. Она не отражает вероятностную степень связи значений случайных ошибок в оценках параметров в различные моменты времени. В такой же мере эта вероятностная связь не учитывается и в формулах В. С. Пугачева, полученных для байесовского фильтра с конечной памятью.

Разработан сплайн-интерполяционный метод расчета вероятностных характеристик нелинейных систем, представимых системами обыкновенных дифференциальных уравнений, операторами Вольтерра и Гаммерштейна, а также метод вычисления статистических узлов на основе решения обыкновенного дифференциального уравнения первого порядка. Из сплайн-интерполяционного метода как частный случай следует интерполяционный метод В.И. Чернецкого, если интервалы существования случайных величин не разбивать на частичные подынтервалы. Применение сплайн-интерполяционного метода наиболее эффективно для расчета вероятностных характеристик динамических систем, нелинейности которых (как функции от случайных величин ) имеют разрывы непрерывности низких порядков и для несимметричных нелинейностей.

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

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

Собственно вычисление сводится к решению интегрального уравнения Фред-гольма первого рода с сильно осциллирующим ядром и аппроксимации частотной характеристики в базисе сплайнов Лагранжа, ВВСФ и обобщенным рядом Котельникова. Принятые аппроксимации приводят к аналитическому вычислению интеграла интегрального уравнения и представлению последнего в виде линейной регрессии. При этом оценивание параметров регрессии производится фильтром Калмана или байесовским фильтром с конечной памятью. Установлена оценка погрешности вычисления интеграла с сильно осциллирующим ядром, для которой выведено уточненное неравенство С.Н.Бернштейна для функций с финитным спектром и конечной энергией.

В существующих методах решения интегральных уравнений Фредгольма первого рода интеграл заменяется какой-либо квадратурной формулой. При сильно осциллирующем ядре это потребует большого числа узлов.

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

Обобщенный ряд Котельникова обеспечивает лучшую аппроксимацию функций и их спектров по сравнению с обычным рядом Котельникова.

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

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

9. В разработке метода точного вычисления распределения Стьюдента для любых значений к-степеней свободы. Новизна состоит в рекуррентном вычислении отношения гамма-функций в формуле распределения Стьюдента. Метод реализован в устройстве авторского свидетельства на изобретение. Распределение Стьюдента вычисляется через элементарные функции с помощью громоздких формул и при значениях к>20 заменяется нормальным законом, но нет оценки погрешности такой аппроксимации.

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

Практическая значимость. Разработанные алгоритмы и программы могут быть использованы для решения следующих практических задач:

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

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

3.Расчета вероятностных характеристик ( математических ожиданий, моментов высших порядков, авто и взаимных корреляционных функций, плотностей распределений) выходных процессов нелинейных стохастических систем, представимых системами обыкновенных дифференциальных уравнений, операторами Вольтерра и Гаммерштейна, сплайн-интерполяционным методом с определением статистических узлов интегрированием обыкновенного дифференциального уравнения первого порядка.

4,Оценки параметров траекторий движения объектов максиминным методом правдоподобия при их аппроксимации уравнениями линейных регрессий и статистически линеаризованным методом максимума апостериорной вероятности при аппроксимации траекторий уравнениями нелинейных регрессий. Прогнозирования траекторий движения объектов и расчета вероятностных ошибок прогнозирования по аналитическим выражениям.

5.Оценки текущего и прогнозируемого состояний систем передачи информации с применением синтезированного в базисе ВВСФ испытательного сигнала в аппаратно-программных комплексах контроля функционирования этих систем.

6.Вычисления частотных характеристик линейных стационарных динамических систем по весовым функциям, реставрации изображений, аналитического продолжения спектра изображения для его восстановления на основе решения интегрального уравнения Фредгольма первого рода с аппроксимацией искомых функций сплайнами Лагранжа, ВВСФ, обобщенным рядом Котельникова и аналитическим вычислением интеграла этого интегрального уравнения. Данный метод вычисления интегрального уравнения применим также для синтеза антенн (акустических, радиотелескопов и др.), для повышения разрешающей способности антенн, оценки энергетического спектра по функции автокорреляции, восстановлении сигналов и их спектров в оптике, в том числе и с аналитическим продолжением спектров.

7. Оценки доверительных интервалов для оцениваемых параметров весовых функций динамических систем с вычислением точных значений квантилей ( без привлечения таблиц и интерполяции ) интегрированием дифференциального уравнения первого порядка. Алгоритм вычисления квантилей для больших выборок (обратной функции Лапласа) реализован в устройстве, защищенным авторским свидетельством на изобретение, алгоритм вычисления доверительных интервалов для математических ожиданий при малых выборках реализован в устройстве, защищенным в другом авторском свидетельстве на изобретение.

Результаты диссертации используются при чтении лекций и проведении практических занятий со студентами по специальному курсу: "Методы и алгоритмы оценки параметров случайных процессов".

Апробация работы. Результаты работы докладывались и обсуждались: на XIV Московской городской НТК, посвященной Дню Радио (Москва, 1988); на двух межотраслевых научно-производственных конференциях " Развитие и совершенствование телевизионной техники " (НИИТТ "Электрон" г. Львов, 1990-1991г.г.); на 1-12-ой ежегодных Всероссийских научно-технических конференциях "Современное телевидение" ( Москва; 1992-2004 г. г.); на НТК специалистов и молодых ученых "Развитие и совершенствование телевизионной техники" (НИИТТ " Электрон ", г. Львов, 1993) ; на двух Международных конференциях "Информационные технологии в проектировании", "Восток-Запад" (Москва, 1994,1996 г.г.); на Международной конференции " 100-летие начала использования ЭМ волн для передачи сообщений и зарождения РТ " (Москва, 1995); на 2-ой Международной конференции " Спутниковая связь " (Москва, 1996); на LII научной сессии, посвященной Дню Радио(РНТО РЭС им. Попова А.С., Москва, 1997); на 22nd Europen Meeting of Statistical, 7nd Vilnus Conference of Probability Theory and Mathematical Statistics (Вильнюс, 1998); на 3-м сибирском конгрессе по прикладной и индустриальной математике, посвященном памяти С.Л.Соболева (1908-1989) ( Новосибирск, 1998); в Российском научном обществе исследования операций, ВЦ РАН (Москва, ^ 2001); на 4-ой Московской Международной конференции по исследованию операций (Москва,2004 г.).

Публикации. Основные результаты диссертации опубликованы в двух монографиях издательства "Радио и связь", статья в иностранном журнале, в центральных научных журналах, в ведомственных научных журналах и научно-методических сборниках, в Трудах Межведомственных, Всероссийских и Меж-дународых конференций; в сборниках трудов Тверского государственного университета, в виде авторских свидетельств на изобретения. Список работ приведен в конце автореферата. Достоверность результатов основана:

-на корректности постановок задач, адекватно описывающих изучаемые физические процессы;

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

-на реализации требований несмещенности и состоятельности оценок параметров;

-на доказательстве, что сплайн-интерполяционный метод дает лучшие по точности оценки математического ожидания и дисперсии для несимметричных и негладких (в том числе разрывных) нелинейностей по сравнению с интерполяционным методом Чернецкого;

-на оценке положительного вклада от использования априорной информации в точность оценок параметров весовых функций в зависимости от объема выборки;

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

-на оценивании параметров весовых функций динамических систем с доверительными вероятностями в диапазоне 0,95 <Р<0,999.

В основу постановки и решения проблемы идентификации приняты следующие исходные положения-данные:

1. Когда система находится в нормальном режиме функционирования, то считается, что её амплитудно-частотная характеристика аппроксимируется константой в пределах рабочей полосы пропускания частот и равна нулю—вне ее, а фазочастотная характеристика—линейной функцией. Иначе говоря, это приводит к следующим известным фактам: входной сигнал передается и обрабатывается без каких-либо искажений, справедлив принцип суперпозиции и выходной сигнал записывается во временной области сверткой вида I y(t)= jh(T)x(t-T)dr, t-T где х(т) - входной сигнал; hit) - весовая функция; Т—память системы.

Для пространственно-временного входного сигнала (двумерного поля, изображения) f(t,u,v), где и, v - пространственные координаты, значение выходного t* сигнала в момент времени t в точке {a, ff) через весовую функцию двумерной системы h(u,v) записывается в виде а Р z(t,a,p)= j \h(u,v)f(t,a-u,p-v)dudv, а-Т, Р-Ту где ТХ,ТУ- память двумерной системы.

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

Пользуясь теперь условием равенства для неравенства Буняковского-Коши, имеем h(t) = kx(t-t0), а в частотной области Sy(ja>) = K(ja>)Sx(ja>)- для одномерной системы и S2Uax,jay) = K{j(ox, jcoy)Sf (jcox, j<oy)- для двумерной системы, где Sy(ja>), Sx(jco), Sz(ja>x,jcoy), Sf{jcox,jcoy)- спектральные плотности сигналов на выходе и входе системы;

K(jco) = ke~ia*° - частотная характеристика одномерной системы с АЧХ, равной к и ФЧХ - (-ja> ta), а для двумерной системы K{jcox,jcoy) = ке~1а'хе~]ЮуУ и (-jcoxx - j(oyy) соответственно.

Весовая функция одномерной системы к Sin2nF(t-t0) h{t) = л (t-t.) Л к Sin27rF {и -и ) Sm2nF (у- v ) а двумерной - h(u,v) =--jV oJ-—-2-, п (и-и0) Ф-VJ где F, Fx, Fy—верхние граничные частоты, t0, Uq, Vq - характеристики запаздывания систем.

Отсюда следует, что амплитудно-частотный спектр входного сигнала совпадает с АЧХ одномерной или двумерной системы. Это положение принимается в основу формирования одного из требований к испытательному сигналу в частотной области.

2. Если в процессе функционирования параметры и характеристики системы изменяют свои свойства и возникает рассогласование текущего состояния системы с входным испытательным сигналом, то сигнал на выходе естественным образом должен быть согласован с текущим, но априори неизвестным состоянием системы и нести информацию об этом. Установление и идентификация текущего состояния по выходному сигналу составляет обратную задачу относительно задачи формирования входного сигнала, согласованного с системой. Она заключается в восстановлении в ортогональном базисе и оценке параметров весовой функции Щ) по выражению t y(t) — \k(T)x(t-T)dT + 0(t), k(t) * h(t), t-T когда имеют место мешающие воздействия 0{t) и линейные искажения в системе, или в восстановлении ядер отрезка ряда Вольтерра совместно с коэффициентами нестационарности по выражению [3] t y(t) = bD(t) + kit) J*,(r,)*(f-г,)</г, + • • • + t-T t K(t)\--- \кп{тх,тг,.,тп)х(t - rt) • • • x(t - TjdTy •■■drn+ G{t), t-T t-T когда в системе возникают нелинейные искажения и её функционированию сопутствуют мешающие воздействия 0{t), которые считаем заранее неизвестными.

Теперь воспользуемся еще одним известным обстоятельством: динамическая система в исходном нормальном состоянии представляется идеальным фильтром низких частот с весовой функцией h(t), а в других состояниях—реальным фильтром низких частот с неизвестной весовой функцией k(t) или с неизвестной совокупностью ядер отрезка ряда Вольтерра kx(j\k2(jxj2),.,kn(tx,t2,.,tn) и коэффициентов нестационарности ba(t),bx(t),.,b„(t). Очевидно, что наибольшую трудность представляет решение второй из отмеченных обратных задач. Эта задача из-за необходимости вычисления многомерных интегралов и синтезирования входного сигнала, наилучшим образом удовлетворяющего требованиям сосредоточения энергии на заданном отрезке времени Т, является очень сложной. Сложность задачи усиливается при программной ее реализации на ПЭВМ в реальном масштабе времени. 3. Для облегчения решения можно воспользоваться математической моделью системы в виде оператора Гаммерштейна. Переход к оператору Гаммерштейна обоснован, так как при отмеченном в п.2 представлении системы сигнал на её выходе полностью совпадает с сигналом на входе, когда эта система представляется в рассматриваемой исходной обратной задаче отрезком ряда Вольтерра. Отсюда следует, что ядра Вольтерра в ортогональном базисе функций пространства L2 (обозначим их i=0, N), должны быть представлены не многомерными суммами, а одномерными, т. е. N 0 N

К (j„T2) = Yu (Г1Ж (Т2 ) > 0 N кп(тх,т2,.тп) = • • • у/Хт„),

1=0 что значительно понижает порядок системы алгебраических уравнений для определения неизвестных коэффициентов Фурье ;. Точность решения исходной задачи в целом не ухудшается.

4. Наличие мешающих воздействий 0(t) обусловливает статистический характер обратной задачи, поэтому метод ее решения существенно зависит от полноты описания характера и частотно-временных параметров воздействий. В диссертации рассматриваются аддитивные флуктуационные, импульсные широкополосные и сосредоточенные узкополосные помехи. Статистические характеристики и основные параметры помех в процессе исследований были приняты по данным [4-10]. Согласно этим данным оказалось возможным представить помехи как аддитивные гауссовы и помехи, для которых могут быть установлены только диапазоны изменения энергетического параметра в полосе частот динамической системы.

5. Построение метода решения задачи идентификации должно осуществляться с учетом требований по частотно-временным характеристикам к входному испытательному сигналу, по достоверности, своевременности и полноте оценок характеристик текущего состояния контролируемой системы. Это, в свою очередь, вызывает необходимость восстановления соответствующего принципа оптимальности. Для его восстановления общая задача идентификации декомпозируется на две: задачу оптимального синтеза широкополосного импульсного входного испытательного сигнала и задачу оценивания параметров ядер Воль-терра как функций, представимых в ортогональном базисе пространства L2. При этом вводится допущение о том, что на интервале времени (t-T\ t) получения выборки измерений состояние системы не изменяется, изменение может иметь место при сдвиге интервала Г на временной оси. В результате искомый принцип оптимальности должен объединить отмеченные здесь требования, т. е. должен представляться объединением принципа неопределенности, принципа формирования общего вида функционала в Ьг на всевозможных линейных комбинациях базисных функций и принципа оптимальности Лагранжа.

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

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

- для уменьшения динамической ошибки аппроксимации необходимо увеличивать степень полинома, а это приводит к увеличению дисперсии случайной ошибки аппроксимации, к трудности реализации на ЭВМ вычислений самих полиномов и коэффициентов аппроксимации;

- на апертуре фильтра во входном сигнале трудно учесть разрывы производных низких порядков.

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

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

- обладают свойством локальности, то есть наличие особенности в какой-либо точке входного сигнала не влияет на поведение аппроксимирующего сплайна в целом;

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

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

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

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

В первой главе формулируются принципы и показатели контроля состояния систем; приведены методы их идентификации.

Во второй главе сформулирована постановка задачи синтеза тестовых сигналов в базисе ВВСФ для контроля технического состояния систем. Показаны преимущества предложенного метода синтеза тестового сигнала перед существующими. Разработаны методы вычисления ВВСФ на основе ряда Котельнико-ва, сплайнов Лагранжа, В-сплайнов, эрмитовых сплайнов третьей степени и полной проблемы собственных значений для симметрической матрицы.

В третьей главе изложен вычислительный алгоритм синтеза тестовых сигналов по заданным частотно-временным параметрам. Приведены характеристики тестового сигнала, представимого линейной комбинацией ВВСФ.

Четвертая глава содержит методы оценки весовых функций линейных и нелинейных систем с сосредоточенными параметрами при различных условиях априорной неопределенности об ошибках измерений.

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

В шестой главе разработаны методы оценки временных и частотных характеристик контролируемых систем: переходной характеристики, амплитудно-частотной и фазо-частотной характеристик, группового времени запаздывания.

В седьмой главе предложен метод расчета вероятностных характеристик: математического ожидания, центральных моментов высших порядков, корреляционных функций, интегрального закона распределения выходных процессов нелинейных динамических систем, представимых нелинейными дифференциальными уравнениями. Эти вероятностные характеристики рассчитываются с использованием сплайнов Лагранжа, что позволяет учитывать априорную информацию о свойствах контролируемой системы и входного процесса для оптимального выбора статистических узлов. При таком выборе узлов сокращается объем вычислений и повышается точность расчета вероятностных характеристик. В главе предложен также численный метод вычисления распределения Стьюдента, устойчивый при реализации на ПЭВМ к большим значениям степеней свободы. Разработанный метод вычисления статистических узлов на основе решения обыкновенного дифференциального уравнения первого порядка может использоваться для произвольных законов распределения с конечными и бесконечными диапазонами изменения аргумента а также для формирования случайных чисел на ПЭВМ.

В восьмой главе изложены статистические методы восстановления спектров тестовых сигналов, частотных характеристик систем по конечной выборке во временной области. С использованием сплайнов Лагранжа, ряда Котельникова и ВВСФ задача восстановления частотных характеристик, спектров сигналов сведена к оценке параметров уравнений линейной регрессии. Оценка параметров производится фильтром Калмана и байесовским оператором сглаживания. В главе предложен обобщенный ряд Котельникова, обладающий лучшими ап-проксимационными свойствами по сравнению с обычным рядом Котельникова. В приложении к главе получена оценка погрешности интегрирования сильно осциллирующих функций с финитным спектром на основе сплайнов Лагранжа и Бернштейна. При этом уточнено фундаментальное неравенство Бернштейна для сигналов с конечной энергией.

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

В главе рассматриваются также вопросы определения геометрических характеристик объектов на тестовом двумерном поле. Приведена вычислительная схема расчета контуров объектов. Для аппроксимации контуров используются параметрические сплайны различных типов: кубические дефекта 1, кусочно-линейные, рациональные. С помощью этих сплайнов рассчитываются геометрические характеристики объекта: длина периметра, кривизна контура, радиус кривизны контура, площадь объекта, координаты геометрического центра объекта.

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

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

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

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

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

1. Показатели контроля состояния динамических систем

Заключение диссертации по теме «Системный анализ, управление и обработка информации (по отраслям)», Малевинский, Михаил Федорович

Заключение

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

1 .Разработаны теоретические основы по обоснованию: -необходимого условия оптимальности весовой функции двумерного линейного согласованного фильтра обработки информации с целью обнаружения входных возмущений по критерию максимума отношения сигнал/ шум. Необходимое условие выведено в форме двумерного интегрального уравнения Фредгольма первого рода с учетом априорной информации о математических ожиданиях и центральных корреляционных моментах второго порядка коэффициентов обобщенного полинома входного сигнала.

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

-необходимых условий оптимальности нелинейного согласованного двумерного фильтра обработки информации, представимого оператором Гаммерштейна n-го порядка, при обнаружении негауссовых случайных входных воздействий, характеризуемых одно и двумерными моментами до 2п-го порядков.

Собственно необходимые условия выведены по критерию максимума отношения сигнал/ шум на выходе фильтра и сформулированы для ядер Гаммерштейна в виде системы n-го порядка двумерных интегральных уравнений Фредгольма первого рода.

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

-необходимого условия оптимальности весовой функции двумерного дискретного байесовского фильтра обработки информации при идентификации динамических систем по критерию максимина среднего квадрата ошибки оценки параметров восстанавливаемой функции с учетом априорной информации о входном воздействии.

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

-необходимых условий оптимальности весовых функций нелинейных фильтров обработки информации, представимых ядрами оператора Гаммерштейна, для идентификации нелинейных динамических систем с учетом одно и двумерных моментов до 2п-го порядков негауссовых тестовых входных воздействий. Необходимые условия оптимальности ядер Гаммерштейна выведены по критерию минимума среднего значения квадрата ошибки идентификации в виде системы n-го порядка двумерных интегральных уравнений Фредгольма первого рода. Они охватывают теоретические положения линейных стационарных фильтров обработки информации с конечной памятью а при степени равной единице оператора Гаммерштейна - уравнения Винера-Хопфа.

-высокоэффективных методов восстановления спектров тестовых сигналов, частотных характеристик идентифицируемых динамических систем по конечной выборке отсчетов на их выходе.

Восстановление спектров осуществляется как решение интегральных уравнений Фредгольма первого рода их аппроксимацией линейными комбинациями сплайнов Лагранжа, ВВСФ, обобщенным рядом Котельникова с неизвестными параметрами с последующим аналитическим вычислением интеграла с сильно осциллирующим ядром и определением параметров линейной комбинации в зависимости от условий получения выборочных данных методом коллокации, фильтром Калмана или байесовским фильтром с конечной памятью.

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

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

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

Соответственно проблема идентификации декомпозируется на следующие задачи :

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

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

-разработки методов оценки вероятностных характеристик выходных процессов линейных оптимальных двумерных фильтров, операторов Вольтерра и Гаммерштейна обработки информации, вероятностных характеристик фазовых координат динамических систем, представимых нелинейными обыкновенными дифференциальными уравнениями (методов оценки математических ожиданий, центральных моментов высших порядков, корреляционных функций, интегральных законов распределений); -разработку программно-реализуемых методов отображения результатов идентификации на дисплеях ПЭВМ для принятия соответствующих решений.

3. В качестве критериев оценки параметров в зависимости от объема сведений о полезном сигнале (параметрах весовой функции системы) и помехах принимаются критерии в виде максимума отношения сигнал/ шум, минимума, максимина, минимакса среднего значения квадрата ошибки идентификации, максимума апостериорной плотности вероятностей, максимина функции правдоподобия. По приведенным критериям оптимизации рассмотренные задачи идентификации динамических систем относятся к классу задач оптимального управления при случайных факторах, так как в этих задачах определяемая весовая функция по существу представляет собой функцию управления.

Синтезирование испытательного сигнала осуществляется по критерию минимума среднеквадратической ошибки аппроксимации равномерного спектра в заданной полосе частот и при заданной величине энергии сигнала на конечном временном интервале; сигнал ищется в классе линейных комбинаций ВВСФ.

4. Для решения поставленных задач обработки информации в диссертации разработана соответствующая совокупность методов и алгоритмов. Эта совокупность включает:

-алгоритм синтеза оптимального испытательного сигнала в базисе ВВСФ, включающий два метода вычисления ВВСФ;

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

- максиминный метод максимума правдоподобия при априорной неопределенности относительно ошибок выборочных данных;

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

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

- сплайн-интерполяционный метод расчета вероятностных характеристик динамических систем, представимых системами нелинейных обыкновенных дифференциальных уравнений, операторами Вольтерра и Гаммерштейна;

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

Результаты моделирования на ПЭВМ показали высокую эффективность разработанных методов и алгоритмов для идентификации систем. Так, для сложных условий функционирования динамической системы типа канала передачи информации, когда имеют место изменения его частотных и временных характеристик из-за воздействия различных факторов (помех, сбоев, неисправностей и др.) установлены следующие значения показателей качества методов: -точность оценок весовой функции-99,5%; -точность оценок переходной характеристики-(99-98,5); -точность оценок амплитудно-частотной характеристики-97%; -точность оценок фазо-частотной характеристики-97%; -точность оценок группового времени запаздывания-97%. 5. Для отображения на экране монитора ПЭВМ результатов идентификации динамической системы разработаны алгоритмы формирования и выдачи данных по оцененным весовой функции, переходной характеристике, АЧХ, ФЧХ, групповому времени запаздывания.

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

Обобщенная информационная модель функционирует на этапе установки (смены, снятия) режимов контроля системы и отображает: -дату и время начала работы оператора с моделью идентификации системы;

-установленный режим контроля (оперативный, функциональный, периодический);

-таблицу контролируемых характеристик динамической системы и их краткие условные обозначения; -контролируемые характеристики;

-меню-перечень и назначение клавиш, с помощью которых оператор работает с моделью.

Информационная модель о текущем состоянии функционирует на этапе контроля текущего состояния системы по установленному режиму и отображает:

-дату и время проведения контроля; -режим контроля;

-таблицу контролируемых характеристик и их условные обозначения; -оцениваемые характеристики;

-обобщенное состояние: " ИСПРАВНА", "ВНИМАНИЕ", "БРАК"; -меню- перечень и назначение клавиш, с помощью которых оператор работает с моделью.

Информационная модель детальной информации по установленной характеристике отображает: -дату и время проведения контроля;

-характеристика системы, по которой проводится просмотр детальной информации;

-состояние системы по данной характеристике :"ИСПРАВНА", "ВНИМАНИЕ", "БРАК";

-таблицу уровней значений оцениваемой характеристики в контрольных точках во временной или частотной областях;

-вычисленные значения установленной характеристики в контрольных точках;

-меню-перечень и назначение клавиш, с помощью которых оператор работает с моделью.

Информационная модель графического изображения оцененной характеристики отображает: -установленную характеристику;

-график установленной характеристики во временной или частотной области;

-области значений нормального, допустимого и состояния брак функционирования системы по установленной характеристике; -состояние системы по установленной характеристике: "ИСПРАВНА", "ВНИМАНИЕ", "БРАК";

-меню-перечень и назначение клавиш, с помощью которых оператор работает с моделью.

Разработанные алгоритмы и информационные модели использовались в аппаратно-программном комплексе контроля текущего состояния динамической системы - телевизионного канала передачи информации.

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

1. Калман Р., Фалб П., Арбиб М. Очерки по математической теории систем, М.: Мир, 1971.

2. Вязгин В.А., Федоров В.В. Математические методы автоматизированного проектирования. М.: Высшая школа, 1989.

3. Ван-Трис Г. Синтез оптимальных нелинейных систем управления. Пер. с англ., под ред. А.Ю. Ишлинского. М.: Мир, 1964.

4. Тихонов В. И. Нелинейные преобразования случайных процессов. М.: Радио и связь, 1986.

5. Шерайзин С. М. Адаптивная коррекция и фильтрация телевизионного сигнала. М.: Радио и связь, 1987.

6. Обнаружение радиосигналов. Под ред. Колосова А. А. М.: Радио и связь, 1989.

7. Варакин JI. Е. Теория систем сигналов. М.: Сов. радио, 1978.

8. Правила эксплуатации магистральной и внутризоновой первичных сетей ЕАСС. Часть 4.—М.: Радио и связь, 1987.

9. Системы подвижной радиосвязи. Под ред. И. М. Пышкина. М.: Радио и связь, 1986.

10. Дворецкий И. М., Дриацкий И. Н. Цифровая передача сигналов звукового вещания. М.: Радио исвязь, 1987.

11. Чернецкий В.И. Анализ точности нелинейных систем управления . М.: Машиностроение, 1968.

12. Абрамов О.В., Розенбаум А.Н. Прогнозирование состояния технических систем. М.: Наука, 1990.

13. Френке JI. Теория сигналов. М.: Сов. радио, 1974.

14. Ямпольский Э. М. Вариационные принципы согласования сигналов каналом связи. М.: Радио и связь, 1978.

15. Дворкович В. П. Оптимизация измерительных сигналов для оценки характеристик телевизионного канала// Радиотехника, 1988. № 2.

16. Лившиц К.М., Терпугов А.Ф. О выборе сигналов для идентификации линейных систем по методу наименьших квадратов // Известия АН СССР, Техническая кибернетика, Мм, 1974, №5.

17. Лившиц К.М., Терпугов А.Ф. Идентификация линейных систем био-ртогональными функциями // Автоматика и телемеханика, 1974, №9.

18. ГришинВ.Н., Дятлов В.А., Милов Л.Т. Модели, алгоритмы и устройства идентификации сложных систем. Л.: Энергоатомиздат, 1985.

19. Бессонов А.А., Загашвили Ю.В., Маркелов А.С. Методы и средства идентификации динамических объектов. Л.: Энергоатомиздат, 1989.

20. Хургин Я. И., Яковлев В. Н. Финитные функции в физике и технике. М. :Наука, 1971.

21. Виленчик Л. С., Катулев А. Н., Михно В. Н., Михно Г. А. Алгоритмические измерения в телевидении и радиовещании. М: Радио и связь, 1995.

22. Левин Б. Р. Теоретические основы статистической радиотехники. М.: Радио и связь, 1989.

23. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1970.

24. Уилкинсон, Райнш. Справочник алгоритмов на языке АЛГОЛ. Линейная алгебра. М.: Машиностроение, 1976.

25. Марчук Г.И. Методы вычислительной математики. М.: Наука, 1980.

26. Льюнг Л. Идентификация систем теория для пользователя. М.: Наука, 1991.

27. Катулев А.Н., Малевинский М.Ф., Соломаха Г.М. Двумерный фильтр с конечной памятью и его вероятностные характеристики // Вестник РУДН. №1(11), 2002.

28. Лэннинг Д.Ж., Бэттин Р.Г. Случайные процессы в задачах автоматического управления. М.: ИЛ, 1958.

29. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1976.

30. Ланкастер П.Л. Теория матриц. М.: Наука, 1982.

31. Казаков И.Е., Мальчиков С.В. Анализ стохастических систем в пространстве состояний. М.: Наука, 1983.

32. Малевинский М.Ф., Виленчик Л.С., Катулев А.Н. Минимаксный метод обработки изображений // Известия РАН Теория и системы управления , 2000.№ 2. 0,4 пл.

33. Малевинский М.Ф., Катулев А.Н., Кузнецов В.Н.,Соломаха Г.М. Двумерный полиномиальный фильтр // РАН Автоматика и телемеханика ,№ 9,2003. с.77-88.

34. Левин Б.Р. Теоретические основы статистической радиотехники. М.: Радио и связь, 1989.

35. Основы автоматического управления. Под ред. B.C. Пугачева. М.: Наука, 1968.

36. Малевинский М.Ф, Виленчик Л.С., Катулев А.Н. Ковариационная функция оценок фильтров Семенова и Калмана // Радиотехника, № ЗД991.с.З-5.

37. Давенпорт В.Б., Рут B.JI. Введение в теорию случайных сигналов и шумов. М.: ИЛ, 1960.

38. Крутько П.Д. Статистическая динамика импульсных систем. М.: Сов. Радио, 1963.

39. Михлин С.Г. Вариационные методы в математической физике. М.: Гостехиздат, 1957.

40. Пупков К. А., Капалин В. И., Ющенко А. С. Функциональные ряды в теории нелинейных систем. М.: Наука, 1976.

41. Верлань А. Ф., Сизиков В. С. Интегральные уравнения: методы алгоритмы программы справочное пособие. Киев, Наукова думка, 1986.

42. Поцелуев А.В. Статистический анализ и синтез сложных динамических систем. М.: Машиностроение, 1984.

43. Малевинский М.Ф, Катулев А.Н., Соломаха Г.М. Оценка точностных характеристик динамических систем на основе сплайнов Лагранжа // РАН.Известия Академии наук.Теория и системы упрвле-ния.№3,2002.с. 19-28.

44. Фаронов В.В. Программирование на персональных ЭВМ в среде ТУРБО-ПАСКАЛЬ. М.: МГТУ, 1991.

45. Малевинский М.Ф.,Виленчик Л.С.,Катулев А.Н. Метод восстановления спектра функции по конечной выборке // РАН РФ Автоматика и телемеханика, том LVIII, №6, 1997.0,7 пл.

46. Богуславский И.А. Прикладные задачи фильтрации и управления. М.:Наука,1977.

47. Березин И.С., Жидков Н. П. Методы вычислений .Т. 1, 2 . М.-.Наука, 1966.

48. Анисимов Б. В., Курганов В. Д., Злобин В. К. Распознавание и цифровая обработка изображений. М.: Высшая школа, 1983.

49. Яншин В.В. Анализ и обработка изображений : принципы и алгоритмы. М.: Машиностроение , 1995.

50. Методы компьютерной обработки изображений. Под ред. В.А. Сой-фера. М.: Физматлит, 2003.

51. Воробьев В.И. Трибун В.Г. Теория и практика вейвлет-преобразования. ВУС, С-Пб, 1999.

52. Добеши И. Десять лекций по вейвлетам. М.: Регулярная и хаотическая динамика, Ижевск, 2001.

53. Дремин И.М., Иванов О.В., Нечитайло В.А. Вейвлеты и их использование. // Успехи физических наук, Т. 171, №5, М.: Наука, 2001.

54. Бакут Б.А., Колмогоров Г.С., Ворновицкий И.Э. Сегментация изображений: методы пороговой обработки. // Зарубежная радиоэлектроника. 1987., № 10.

55. Денисов Д.А., Низовкин В.А. Сегментация изображений на ЭВМ //

56. Зарубежная радиоэлектроника, 1985, №10.

57. Дидэ Э. И др. Методы анализа данных. — М.: Финансы и статистика, 1985.

58. Виленчик JI.C., Катулев А.Н., Богданчук В.З. Методы оптимальной обработки информации в информационно-измерительных системах М.: Радио и связь, 1991.

59. Козлов М.В. Прохоров А.В. Введение в математическую статистику. М.: МГУ, 1987.

60. Миддлтон Д. Введение в статистическую теорию связи. Перевод с английского. Под редакцией Б.Р.Левина.М.: Сов.радио.1961.

61. Малевинский М.Ф., Виленчик Л.С., Катулев А.Н., Колганов С.К. Идентификация каналов передачи информации: новые методы и алгоритмы. М.: Радио и связь, 1993.81 с.

62. Малевинский М.Ф., Виленчик Л.С., Катулев А.Н., Колганов С.К.

63. Идентификация каналов передачи информации, М.: Радио и связь, 1996. 100с.ч