Разработка единого системного подхода к решению задачи оптимального оценивания
Дипломная работа - Компьютеры, программирование
Другие дипломы по предмету Компьютеры, программирование
частотной областях. Для функций, имеющих высокую скорость убывания в указанных областях, рассмотренный математический аппарат позволяет обеспечить требуемую точность N-кратного дифференцирования на заданном отрезке [-Т, Т], если уровни усечения не превышают некоторых заранее установленных значений. Рассмотрим возможность применения данного аппарата для вычисления производных соответствующего порядка от финитных функций имеющих плохие спектральные свойства и большие уровни усечения в пространственной (временной) области.
Первый подход к дифференцированию финитных функций, основанный на сплайн-продолжениях, состоит в следующем. Пусть - произвольная финитная функция (рис. 2.1), у которой производная непрерывна, а кусочно-непрерывна на отрезке [-Т, Т]. Используя операцию сплайн-продолжения, перейдем от функции заданной на отрезке [-Т, Т], к новой функции , заданной на всей вещественной оси (рис. 2.2):
Рисунок 2.1Рисунок 2.2
(2.44)
где ?1(t) и ?2(t) - вспомогательные функции, у которых производные и непрерывны при и соответственно, и, кроме того, выполняются равенства (условия стыковки)
(2.45)
Таким образом, из исходной финитной функции получили сплайн-продолженную функцию ?(t), заданную на интервале (-?, ?), состоящем из основного информационного отрезка [-Т, Т] и двух вспомогательных неинформационных полуоткрытых интервалов и .
Далее, вводя оператор усечения в пространственной (временной) области перейдем от ?(t) к финитной функции , заданной на отрезке (рис. 2.3). При этом потребуем выполнения следующего условия:
(2.46)
где - уровни усечения функции ?(t) и ее производных , которые выбираются с учетом полученных в предыдущих подразделах аналитических зависимостей исходя из условия минимизации результирующей погрешности N-кратного дифференцирования (N < M).
Рисунок 2.3
Очевидно, что периодически (с периодом ) продолженная на всю числовую ось функция остается непрерывной на этой оси вместе со своими производными до (М - 1)-го порядка включительно. Вспомогательные функции ?1(t) и ?2(t), удовлетворяющие условиям (2.45), (2.46), могут быть достаточно произвольного вида. Однако полагаем, что независимо от вида финитной функции в качестве ?1(t) и ?2(t) используются функции, интегрируемые в квадрате, т.е. и . Каждая из функций ?1(t) и ?2(t) должна удовлетворять 2М условиям (2.45), (2.46) и поэтому имеет 2М подлежащих определению коэффициентов.
Необходимость введения условий (2.45), (2.46) можно обосновать из следующих соображений. Если финитную функцию периодически продолжить на всю числовую ось (обозначим полученную функцию через ), то коэффициенты соответствующего ряда Фурье можно определить по формуле
(2.47)
при этом скорость убывания коэффициентов ck с возрастанием определяется дифференциальными свойствами данной периодической функции.
Близким по своей сути к подходу, основанному на сплайн-продолжении, является известный метод Малиева разложения функций в быстросходящиеся ряды Фурье. Однако данный метод предполагает лишь одностороннее продолжение исходной финитной функции с обеспечением выполнения условия непрерывности соответствующей периодической функции и ее производных различных порядков.
Построение же функции ?(t) (2.44) требует двухстороннего продолжения , выполнения указанного условия непрерывности (2.45) и обеспечения требуемых уровней усечения функции ?(t) и ее производных
в пространственной (временной) области (см. (2.46)).
Второй подход к дифференцированию финитных функций, имеющих плохие спектральные свойства, состоит в домножении исходной функции заданной на отрезке [-Т, Т], на некоторую функцию-регуляризатор (где ? - вектор постоянных параметров), заданную на том же отрезке [-Т, T] или на расширенном отрезке :
(2.48)
где является продолжением функции с отрезка [-Т, T] на отрезок с сохранением свойства непрерывности. Для функции и ее производных различных порядков на краях отрезка можно задать соответствующие уровни усечения (по аналогии с (2.46)) исходя из условия минимизации результирующей погрешности N-кратного дифференцирования.
В качестве функции-регуляризатора можно рассматривать функцию вида
(2.49)
где при .
В простейшем случае в качестве можно рассматривать гауссовскую кривую, характеризующуюся двумя параметрами. Известно [29], что такая кривая имеет максимально возможную скорость убывания в пространственной и частотной областях и реализует знак равенства в известном соотношении неопределенности. Выбор значений параметров гауссовской кривой, обеспечивающих требуемое поведение функции и ее спектра соответственно в пространственной (временной) и частотной областях, не представляет затруднений, поскольку известно, что убывает на бесконечности быстрее любой степени t при ? > 0.
Рисунок 2.4
Если найдена производная , то вычисление искомой производной осуществляется в соответствии с известной формулой [4], в которой необходимо заменить на на и на Для того чтобы избежать трудоемкой в вычислительном плане операции восстановления искомой производной в соответствии с указанной формулой, можно использовать в качестве функцию-регуляризатор срезывающего типа [4]. Данная функция является бесконечно дифференцируемой на всей вещественной оси (-?, ?) и вместе со своими производными тождественно равна нулю для всех Кроме тог