8.5. РасходимостьВ некоторых приложениях можно обнаружить, что истинные ошибки оценивания во много раз превышают теоретически предсказанные ошибки, характеризуемые дисперсией Может показаться, что утверждение о расходимости резко противоречит результатам по устойчивости, полученным в § 7.5. Однако на самом деле это не так, поскольку расходимость вызывается такими ошибками, которые учтены в этом параграфе не были. Одними из главных причин расходимости являются: неточность задания процесса, моделирующего сообщение и наблюдение, линеаризация уравнений, отсутствие полной информации о реальной физической задаче, а также всякого рода упрощающие предположения, позволяющие дать математическое описание задачи. Кроме того, неустойчивость могут вызывать ошибки, связанные с моделированием вероятностных характеристик шума и неизвестных входных сигналов. Неустойчивость могут вызвать также ошибки округления, неизбежно возникающие при цифровом моделировании алгоритмов фильтрации и которые могут привести к потере положительной определенности и симметрии ковариационной матрицы ошибок. Для того чтобы выяснить, каким образом возникает расходимость и как ее можно проконтролировать с помощью простой модификации алгоритма фильтрации, рассмотрим простой пример, в котором расходимость возникает в результате неточного задания модели сообщения. Допустим, что необходимо оценить состояние дискретного процесса
где
Модель наблюдения
где Используя алгоритмы табл. 7.2, легко показать, что оценка после
а дисперсия ошибки, согласно фильтрации по Калману, равна
Следовательно, при увеличении Действительное состояние после
Поэтому действительная ошибка оценивания равна разности между (8.271) и (8.269):
Следовательно, действительная среднеквадратическая ошибка равна
откуда (видно, что при увеличении Может показаться, что подобная задача тривиальна, поскольку мы просто игнорируем константу, которую должны были бы учесть. Однако даже если и включить константу в модель сообщения, то оценка все равно могла бы разойтись, до тех пор пока два входных сигнала не станут абсолютно равными, что, к сожалению, практически трудно сделать. Если, например, модель сообщения задать в виде
то очень просто показать, что действительная среднеквадратическая ошибка равна
и оценка снова будет расходиться до тех пор, пока Основная причина, вызывающая расходимость в фильтре Калмана, состоит в том, что коэффициент передачи (усиления) Эта трудность может быть преодолена одной из простых модификаций алгоритма фильтрации. Можно, например, просто ограничить снизу коэффициент передачи
Так как, с другой стороны, мы ограничили
Тогда при Другой подход к устранению расходимости в указанной выше задаче возникает в случае, когда в рассматриваемой модели a есть неизвестное состояние некоторой динамической системы, управляемой белым шумом. Хотя эта процедура часто дает хорошие результаты, она имеет неприятную особенность, связанную с повышением порядка рассматриваемой задачи. Так как порядок фильтра и алгоритма вычисления дисперсии для Выводы, полученные при рассмотрении этой простой задачи, имеют фундаментальное значение и в общем случае. Расходимость может возникнуть из-за слишком малого веса новой информации или, наоборот, из-за слишком большого веса прошлых наблюдений. Поэтому расходимость (насыщение информацией) — явление противоречивое. Как уже известно из гл. 7, в том случае, когда входной шум мал по сравнению с наблюдениями, дисперсия ошибки калмановской фильтрации и, следовательно, коэффициент усиления имеют тенденцию быстро уменьшаться при увеличении времени. В том случае, когда модель сообщения вообще не содержит управляющего шума, коэффициент усиления асимптотически стремится к нулю. Эта особенность, которая в общем не является неожиданной, поскольку (каждый отсчет (в среднем) содержит гораздо больше нежелаемого шума наблюдения, чем информации об управляющем шуме, имеет тенденцию «отключать» фильтр от последовательности наблюдений и таким образом вызывать неустойчивость. Итак, мы приходим к тому, что ситуации, в которых управляющий шум слабо сравним с шумом наблюдения, приводят к проблемам неустойчивости. Насколько сильно проявляется расходимость и в каких случаях она возникнет вообще, все это зависит от длины интервала наблюдения и от точности моделирования реального процесса. Основываясь на предыдущих рассуждениях, можно предложить определенные модификации фильтра Калмана, устраняющие расходимость. В работе [221] рассматриваются несколько простых методов устранения неустойчивости в специфических задачах определения орбит. Основная идея этих подходов связана с ограничением коэффициента усиления для того, чтобы избежать «нечувствительности» к наблюдениям. Эти методы можно разбить на три широких класса: - уменьшение матричных коэффициентов усиления; - ограничение дисперсии ошибок; - искусственное увеличение мощности управляющего шума сигнала. Рассмотрим каждый из этих методов. Существует несколько способов непосредственного уменьшения матричных коэффициентов усиления. Так, можно использовать процедуру, развитую в рассмотренном выше примере, и ограничить элементы матрицы снизу некоторой заранее выбранной величиной. Можно просто добавить к элементам матрицы некоторую величину или использовать менее очевидные способы для увеличения усиления. Например, в работе [222] матрица коэффициентов усиления представлена в виде
Подобно тому, как можно ограничивать снизу элементы матрицы коэффициентов усиления, аналогичная методика может быть применена для ограничения элементов матрицы дисперсий ошибок. Легко видеть, что из-за неизбежной неопределенности в формулировке начальных условий физически невозможно иметь ошибку, равную нулю. Однако точную нижнюю границу дисперсии ошибки также трудно найти, если не использовать эксперимент. Можно показать, что, имея минимальное значение дисперсии ошибок, можно избавиться от расходимости. Поскольку основная причина, из-за которой матричные коэффициенты становятся малыми и возникает неустойчивость, состоит в том, что управляющий шум модели значительно слабее, чем шум наблюдения, то очевидное решение задачи состоит в искусственном увеличении дисперсии управляющего шума. Несмотря на то, что можно вычислить величину необходимой добавки к дисперсии управляющего шума, точная ее величина должна находиться экспериментально. Нетрудно видеть, что при использовании перечисленных выше методов дисперсия ошибки будет большей, если модель сообщения на самом деле точная. Другими словами, контроль расходимости может быть осуществлен только за счет увеличения дисперсии ошибки свыше «теоретического» минимума. Конечно, если оптимальный фильтр неустойчив, то модифицированный неоптимальный фильтр будет иметь дисперсию ошибки, меньшую, чем действительная дисперсия ошибки «оптимальной» фильтрации. Другой подход, позволяющий контролировать расходимость, состоит в вычислении квадратного корня из матрицы ковариаций ошибок оценивания. Подобный способ был развит в работах [21], [23], [5]. Эта процедура полезна в случаях, когда неустойчивость вызывают ошибки округления. Тогда лучше вычислять квадратный корень из дисперсии, чем саму дисперсию, потому что число необходимых операций при счете сокращается в два раза. Рассматривая этот способ, предположим, что шум модели равен нулю; мы уже видели, что неустойчивость возникает чаще всего при малом шуме модели, поэтому наше предположение не является ограничительным. В случае ненулевого шума можно воспользоваться алгоритмом, разработанным в [58]. В дискретном времени уравнения для ковариационной матрицы ошибок оценивания при
Легко показать прямой подстановкой, что корень из матрицы ковариаций равен
где используется нижняя треугольная матрица корней из дисперсий. Аналогичные выражения для случая непрерывного времени были получены в [5]. Другой подход [96] состоит в том, что оценку получают только по прошлым Для устранения неустойчивости можно воспользоваться адаптивными алгоритмами фильтрации [97], [208], [243]. Применение этих алгоритмов может оказаться более целесообразным, поскольку при этом ошибка из-за расходимости становится большой, а значит, и большее количество информации можно использовать для адаптации.
|