Модель в виде фильтра Каллмана

Каллманом была доказана теорема о том, что любой динамический сигнал может быть представлен в виде:

Y_i = A ₁ · X_i + A ₂ · X_{i – 1} + … + B ₁ · Y_{i – 1} + B ₂ · Y_{i – 2} + … + C.

Рис. 6.1. Графическое представление фильтра Каллмана на схемах

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

Чем больше имеется членов ряда, то есть чем больше переменных Y учитывается в записи модели, тем глубже память системы. Заметим, что наличие члена Y_{i – 1} в модели динамической системы соответствует наличию первой производной, Y_{i – 2} — второй производной и т. д.

Допустим, известны следующие экспериментальные данные: состояния сигналов X_i и Y_i в n временных точках (табл. 6.1).

Таблица 6.1. Таблица экспериментальных данных

i Xi Yi   X 1 Y 1   X 2 Y 2 … … … n – 1 Xn – 1 Yn – 1 n Xn Yn

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

Таблица 6.2. Таблица экспериментальных данных и промежуточных расчетов

i Xi Xi – 1 … Yi Yi – 1 Yi – 2 m Xm Xm – 1 … Ym Ym – 1 Ym – 2 m + 1 Xm + 1 Xm … Ym + 1 Ym Ym – 1 m + 2 Xm + 2 Xm + 1 … Ym + 2 Ym + 1 Ym … … … … … … …

Находим ошибку между значением экспериментально снятой точки и теоретическим ее значением (гипотезой):

E_m = Y_m – A ₁ · X_m – A ₂ · X_{m – 1} – … – B ₁ · Y_{m – 1} – B ₂ · Y_{m – 2} – … – C.

Суммарная ошибка F (сумма берется по всем экспериментальным точкам) должна быть минимизирована относительно определяемых переменных A ₁, A ₂, …, B ₁, B ₂, …, C:

После взятия частных производных от F по A ₁, A ₂, …, B ₁, B ₂, …, C, приравнивания их к нулю и составления системы уравнений получается линейная множественная регрессионная модель, из которой определяются неизвестные коэффициенты A ₁, A ₂, …, B ₁, B ₂, …, C модели.

Поскольку коэффициенты модели определены, построим реализацию (см. рис. 6.2), имитирующую поведение системы, описанной фильтром Каллмана.

Рис. 6.2. Вариант технической реализации фильтра Каллмана

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

Например, в случае реализации блока задержки в среде моделирования Stratum-2000, первый способ может быть основан на перезаписи информации из одной переменной (ячейки) в другую, на что требуется один такт. Таким образом, можно организовать задержку сигнала на любое число тактов. Например, задержка сигнала X относительно Y будет составлять 3 такта, если выполнить следующую последовательность операций: A1:= X; A2:= A1; Y:= A2.

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

На рис. 6.3 приведена схема настройки (автоматического нахождения коэффициентов).

Рис. 6.3. Схема автоматической настройки коэффициентов модели «на ходу»

На рис. 6.4 приведена схема проверки фильтра Каллмана.

Рис. 6.4. Схема проверки работы модели фильтра Каллмана

Модель динамической системы в виде
Фурье представления (модель сигнала)

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

Пусть, например, в некотором сигнале присутствует сумма трех гармоник:3 · cos(t) + 2 · cos(3 t) + 0.5 · cos(5 t). Это значит, что в сигнале присутствует первая гармоника с амплитудой 3, третья гармоника с амплитудой 2, пятая гармоника с амплитудой 0.5. Сам суммарный сигнал выглядит так, как показано на рис. 7.1.

Рис. 7.1. Пример гармонического сигнала

Спектр этого сигнала показан на рис. 7.2. Ясно, что в нашем примере больший вес (амплитуду) в сигнале имеет (более других представлена) первая гармоника, наименьший вес имеет пятая гармоника.

Рис. 7.2. Пример спектра гармонического сигнала

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

Пусть задан некий сигнал X (t) (рис. 7.3).

Рис. 7.3. Временной сигнал на входе преобразования Фурье (возможный вид)

Определимся со временем рассмотрения сигнала: если сигнал периодический, то время рассмотрения равно периоду p сигнала; если сигнал непериодический, то периодом сигнала считается все время его рассмотрения.

…	…
…	…

A_i и B_i — это веса соответствующих гармоник, присутствующих в сигнале; i — номер гармоники. Формулы их расчета называются прямым преобразованием Фурье.

Значение 2π · i / p = ω_i — это частота i -ой гармоники. Отметим также, что частота i -ой гармоники связана с частотой первой гармоники простым соотношением: ω_i = i · ω ₁.

Отметим важную особенность данного способа представления: вместо всего сигнала во всех его подробностях достаточно хранить вектор чисел, представляющих весовые коэффициенты составляющих его гармоник: (A ₀, A ₁, A ₂, …, B ₁, B ₂, …). То есть эти числа полностью характеризуют исходный сигнал, так как по ним сигнал можно полностью восстановить формулой обратного преобразования Фурье:

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

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

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

Система чисел A_i и B_i является полной характеристикой сигнала. Такой же полной характеристикой сигнала является система чисел S и φ, которые также образуют спектр (рис. 7.5). S — это амплитудно-частотная характеристика (АЧХ), φ — фазо-частотная характеристика (ФЧХ).

Рис. 7.5. Сигнал, представленный в частотной области, амплитудно-частотная и фазо-частотная характеристика сигнала (возможный вид)

Системы «A и B» и «S и φ» являются полностью равнозначными. Переход из системы «A и B» в систему «S и φ» производится по следующим формулам: S_i = sqrt(A_i ² + B_i ²) — абсолютная амплитуда сигнала; φ_i = arctg(B_i / A_i) — фаза сигнала, при сложении гармоник нужно учитывать сдвиг фаз (сдвиг фаз проиллюстрирован на рис. 7.8).

В случае с системой «S и φ» обратное преобразование Фурье имеет вид:

Рис. 7.6 и рис. 7.7 разъясняют смысл коэффициентов A и B разных гармоник. Эти коэффициенты — амплитуды синусов и косинусов соответствующих частот (гармоник). Во временной области графически они соответствуют размаху гармонических колебаний (рис. 7.6 и рис. 7.7); в частотной — высоте спектральной полоски на соответствующей частоте (рис. 7.4).

Рис. 7.6. Геометрическая иллюстрация параметров А и ω для косинусной составляющей гармонического сигнала

 

Рис. 7.7. Геометрическая иллюстрация параметров В и ω для синусной составляющей гармонического сигнала

Смысл чисел S_i и φ_i разъяснен на рис. 7.8.

Рис. 7.8. Геометрическая иллюстрация параметров S и φ для составляющей гармонического сигнала

Модель динамической системы в виде
Фурье представления (модель объекта)

Пусть имеется входной динамический сигнал X (t) и объект F, преобразующий этот сигнал в выходной Y (t) (см. рис. 8.1). Если объект описывается дифференциальными уравнениями, то таким преобразованием является интегрирование входного сигнала и вычисление Y (t). Интегрирование, как было ранее показано, — операция, требующая значительных вычислительных ресурсов и имеющая значительную погрешность при реализации на цифровых машинах.

Рис. 8.1. Схема моделирования динамического объекта при переходе из временной области представления в частотную

Если перейти от описания входного сигнала во временной области к описанию в частотной области (см. рис. 8.1), а от дифференциальных уравнений перейти к частотной характеристике объекта, — то есть, фактически, заменить сигнал на частотную модель сигнала, а объект на частотную модель объекта, — то с вычислительной точки зрения процесс преобразования сигнала упростится. Конечно, полученный результат тоже будет частотной моделью выходного сигнала, которую для получения окончательного ответа придется сконвертировать во временную область Y (t). Процесс такой конвертации из частотной области во временную и обратно называется преобразованием Фурье (есть и другие преобразования). Для тех объектов, для которых известна их модель в частотной области, такой подход достаточно просто реализуется на компьютере и позволяет достичь любой наперед заданной точности.

Модель объекта в частотном виде называется передаточной функцией или АЧХ (амплитудно-частотной характеристикой). Объекты, для которых известны АЧХ, обычно называют типовыми звеньями (усилительное звено, апериодическое, колебательное и т. д.). Пусть, для примера, характеристика объекта в частотной области следующая (см. рис. 8.2).

Рис. 8.2. АЧХ (возможный вид)

Амплитудно-частотная характеристика (АЧХ) показывает, насколько пропускается объектом на выход соответствующая гармоника. Значение k_i характеризует коэффициент усилeния гармонического сигнала на определенной частоте ω_i.

Моделирование прохождения сигнала через объект в этом виде заключается в умножении коэффициента A_i гармоники с частотой ω_i входного сигнала X (t) на коэффициент усиления k_i при той же гармонике с частотой ω_i в АЧХ: A_i ^* = A_i (ω_i) · k_i (ω_i). (Для коэффициента B преобразование аналогично.) В результате получается коэффициент A_i ^* выходной гармоники данной частоты ω_i. Процедура выполняется для всех частот, представленных во входном сигнале и АЧХ. После получения спектра выходного сигнала можно восстановить сигнал как временную зависимость с помощьюформулы обратного преобразования Фурье.

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

Схема преобразования показана на рис. 8.3.

Рис. 8.3. Схема процедуры преобразования сигнала при использовании метода Фурье

Если бы временной сигнал проходил через звено, которое во временной области представлено дифференциальным уравнением, то пришлось бы его интегрировать, что, конечно, приводит к погрешностям результата. В частотной области достаточно перемножить значения коэффициентов ряда Фурье сигнала и звена при одинаковых частотах. Очевидно, что достоинством метода является замена дифференциальных уравнений модели на алгебраические. Разумеется, данный подход может быть использован только для объектов, у которых известен вид передаточной функции. (Кстати, для неизвестных случаев АЧХ может быть получена численным разложением в ряд.)

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

Пусть мы имеем радио-электронное устройство (РЭУ), состоящее из 5 блоков (см. рис. 8.4). Блоки 1, 2, 4, 5 — линейные и представлены соответствующими известными АЧХ; блок 3 — нелинейный, поэтому АЧХ для него неизвестна. Примером линейного блока может служить апериодическое звено, колебательное звено и т. д. (см. Лекцию 05. Динамические регрессионные модели, заданные в виде передаточной функции). Примером нелинейного блока может служить устройство ограничения сигнала (срез) по амплитуде.

Как видно из рис. 8.4, сначала входной сигнал X (t) прямым преобразованием Фурье переводится в частотную область и проходит в виде спектра через АЧХ 1 и 2 линейного блока, затем обратным преобразованием Фурье сигнал после 2 блока переводится во временную область. Проходим нелинейный блок 3 во временном представлении. Результат работы блока 3 снова преобразуем прямым преобразованием Фурье в частотную область и проходим через АЧХ блоков 4 и 5. В конце полученный спектр преобразуется с помощью обратного преобразования Фурье во временную область, — вид сигнала, Z (t), является результатом моделирования.

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

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

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

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

2 W _max ≤ W _дискр.,

где W _дискр. = 1/Δ t _дискр. — частота дискретизации, W _max — максимальная частота, присутствующая в сигнале.

Оценка качества модели

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

Если оценка качества применяется до исследования, то она решает задачу: есть ли связь между входом X и выходом Y и оценивает силу этой связи.

⇐ Предыдущая12345678910Следующая ⇒

Поделиться:

Воспользуйтесь поиском по сайту: