Построение вейвлет-базисов, адаптированных к дифференциальным операторам Текст научной статьи по специальности «Математика»

УДК 517.2

В.Г. Захаров, О.Н. Ипанов Институт механики сплошных сред УрО РАН

ПОСТРОЕНИЕ ВЕЙВЛЕТ-БАЗИСОВ, АДАПТИРОВАННЫХ К ДИФФЕРЕНЦИАЛЬНЫМ ОПЕРАТОРАМ

Abstract

A method of obtaining biorthogonal wavelet bases adapted to linear differential operators is considered. Examples of biorthogonal wavelet bases with compact support are obtained. Using these bases the Galerkin method of solution of differential equations dmu (x)/dxm = f (x) supplies diagonal stiffness matrix.

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

Метод Галеркина

Рассмотрим одномерное линейное дифференциальное уравнение вида

Lu( x) = f ( xX (1)

где u, f e H, H - некоторое гильбертово пространство, например L2(R), L -линейный дифференциальный оператор, такой, например, как L = dmldxm .

Сформулируем основные принципы решения уравнения (1) с помощью метода Галеркина (Петрова - Галеркина). Слабая формулировка уравнения (1) имеет вид

<Lu( x), ft) = {f, 0,), (2)

где (•, обозначает скалярное произведение в соответствующем гильбертовом пространстве, {вк} - некоторый базис пространства H. Пусть {ак} - некоторый другой базис пространства H, тогда конечномерная аппроксимация решения u(x) примет вид

u(x) = Z u]a] .

J

Подставляя последнее соотношение в (2), имеем

Т",(Ьа,, h) = {f-Pk) = fk. (3)

Наша задача найти такие базисы ак, Рк, при которых матрица жесткости А]к Ьа], в) имела бы наиболее простой вид. В случае, когда матрица жесткости имеет диагональный вид А]к = Ак5]к, уравнение (3) элементарно разрешимо,

ик = Ак!к, и (х) = £ Ак1!как.

к

Если оператор Ь представляет собой сумму дифференциальных операторов с постоянными коэффициентами Ь = £^_оакйк/йхк , то при использовании базиса Фурье

{е1кх, к є Z} матрица жесткости, очевидно, будет диагональной. Для вейвлет-базисов ситуация сложнее. Тем не менее, имея ортонормированный вейвлет-базис {/]к(х) = 2]/2/(2] х - к), ],к є Z} , можно построить биортогональные базисы [1,3,5]

= Ь -/Зк , Р]к = Ь/]к, (4)

где Ь - оператор, сопряженный оператору Ь . Тогда разлагая решение по базису а]к, а правую часть по базису / ]к

и = £{ и Рк) ак, / = £{ /, /]к} /]к,

] ,к ],к

подставляя в уравнение (1) и проектируя на базис /]к , получим

£ (u, Рк) {Ьа]к /) = £{ I /) (/к ,/ы) = (/, /ы).

],к ],к

В силу ортонормированности базиса / матрица жесткости будет диагональной:

{Ьа]к,/ы) = {ьь/]к/ы) = (/к/ы) = ёЫкт и решение уравнения принимает простой

вид

и = £{1 /к)а]к .

] ,к

Так же можно разложить решение и(х) по базису /]к, а правую часть по базису а]к, тогда решение можно записать в другой форме:

и=£{1, Рк) /к •

] ,к

В случае, когда оператор Ь самосопряженный и положительно определенный, например Ь = й2m|dx2т , схему решения можно несколько изменить [2]. Представляя оператор Ь в виде

Ь = V V,

можно определить биортогональные базисы следующим образом:

ак=V , вк=V */]к.

Разложим решение по базису а]к, а правую часть по базису Рік ,

и = £ (и, Рк) ак, 1 = £{1 ,ак) Рк,

],к ],к

подставляем в уравнение (1) и проектируем на базис а1т,

£ ( и Рк) {Ьак ,ат} = £{1, ак) (Рк ,ат) = (1, аы).

],к ],к

Очевидно, матрица жесткости имеет диагональную форму:

{а-к ,аіт) = (V'ак = ^~-/к , ^/т) = {/к /т) = 5Л5кт и решение уравнения

примет вид

u=Кf ,ak )ak.

j,k

Биортогональные вейвлет-базисы и кратномасштабный анализ

Биортогональные, в частности, ортогональные вейвлет-базисы могут быть получены с помощью так называемого кратномасштабного анализа (КМА) (multiresolution analysis). Биортогональный КМА состоит из двух последовательностей замкнутых аппроксимирующих подпространств пространства L (R) :

•••е V-2 е V-1 е V0 е V1 е V2 е•••,

•••е V-2 е V-1 е V0 е V1 е V2 е--%

таких, что

jjVj = L2(R), Wj = L2(R),

ієх ІєХ

n V = {°j, n V = !°j,

j^Z jeZ

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

f (•) є Vj о f (2;•) є Vo, f (•) є Vj о f (2;•) є Vo.

Вейвлет-пространство Wj ( Wj) представляет собой дополнение пространства V j (V j, соответственно) в пространстве V/+1 (Vj+1 соответственно):

V+i = Vj е Wj, j = V е Wj .

Условия биортогональности вейвлет-пространств разных масштабов Wj 1 W,,j^j' эквивалентны условиям

Wo 1 Vo, Wo 1 Vo. (5)

Наконец, существуют функции ф, ф - масштабирующие функции (scaling functions), чьи целочисленные сдвиги образуют базисы Риса (т.е. безусловные базисы) пространств V0 и V0 соответственно.

V0 = span {ф(- - k), k є Z}, V0 = span {ф(- - k), k є Z j.

Из вложенности пространств V0 е V^ V0 е V^ следуют так называемые масштабирующие соотношения (scaling relations)

ф( x) = л/21 \ф(2 x - k), ф( x) = лІ21 hkф(2 x - k),

k k

где

hk ={ф, ф1к) , hk = (ф, ,

или в пространстве Фурье

ф(£) = mo(^ 2)ф(£/ 2), ф(£) = mo(a 2Ж/2),

где

mo(^) = -ij I hkC-lk■ , mo(f) = - jj I V-4 .

Функции m0(^) и m0(£) называются низкочастотными фильтрами или масками и полностью определяют КМА. Для выполнения условия биортогональности (5)

достаточно, чтобы масштабирующие функции были биортогональными относительно целочисленного сдвига: ф(), Ф( - к) = 50к , что в терминах масок имеет вид

m0 (f )m0 (f) + m0(f + n)m0 (f + п) = 1, где черта сверху обозначает комплексное сопряжение. Тогда, определив вейвлет-функции следующим образом:

k(f) = m1 (f /2ф(£ / 2) = eitri m„(f2 + пф(£ / 2), (6)

/(f) = m^f/l)hl 2) = e-'f/2 mf2+T$(fl 2), (7)

получим биортогональные вейвлет-базисы всего пространства L (R). А именно,

обозначив через W0, ( W0) замыкания множеств целочисленных сдвигов вейвлет-

функций:

W0 = span {(• - к), к e Z}, W0 = span {(• - к), к e Z}, получаем, что пространство L2( R) разлагается на прямую (в общем случае неортогональную) сумму вейвлет-пространств WJ (WJ соответственно)

l2(r) =ф w, = е W,.

jeZ 1 jeZ 1

Таким образом, множества функций

/ (x) = 212 /(2]x - к^ / (x) = 21 2/(2 1 x - к), j, к e Z

образуют биортогональные базисы пространства L (R) :

(/к ,/m) = 5fi5km.

Обозначим через {gk} и {gk} коэффициенты, соответствующие высокочастотным фильтрам m1 (f ) и m1 (f ),

m1(f) = -j^ Z gke_ikf, m1(f) = ^ Z gkeiick',

gk = (/ , gk ф1^,

тогда из (6) и (7) следуют простые соотношения для коэффициентов {gk } и {gk } :

gk = (-1)kh- k+1, gk = (-1)kh-k+1.

Алгоритмы разложения и восстановления по вейвлет-базисам

Помимо того, что КМА дает методику для получения вейвлет-базисов, КМА обеспечивает простую реализацию алгоритмов разложения функции по вейвлет-базису и ее восстановления из вейвлет-коэффициентов.

Пусть нам известны коэффициенты разложения некоторой функции f по

масштабирующим функциям самого мелкого масштаба J: cJ =^ f ,фк), как правило, в

качестве коэффициентов ckJ берутся сами значения функции f из дискретной выборки

{fk} : cJ = fk . Тогда вейвлет-коэффициенты d]k={^f ,ajk^ могут быть вычислены по

следующим формулам:

dk-1 =Z gn-2kci , (8)

n

1 =Z hn-2kci , (9)

n

если сигнал дан N = 217 значениями, то мы имеем N коэффициентов разложения по вейвлет-базису а]к : с0, ё3к, 1 = О,..., J -1, к = 0,...,21 -1. Формула восстановления

имеет вид

4 =!( 4-2 „с? + ёк-2 „1 ) • (10)

п

Как мы видели, для вычисления разложения функции по вейвлет-базису достаточно знать лишь коэффициенты низкочастотных {Ик}, {Ик} и высокочастотных

{ёк}, {ёк} фильтров. Алгоритмы разложения-восстановления являются более

быстрыми, чем БПФ, и при использовании вейвлетов с компактным носителем требуют всего 0(Щ операций.

Построение биортогональных вейвлет-базисов,

1Ш I 1 Ш

адаптированных к операторам вида а /ах

Рассмотрим простейший дифференциальный оператор вида Ь = аш1ахш . Пусть мы имеем ортонормированный базис //., который может быть использован для

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

Пусть ш0 низкочастотный фильтр (маска), определяющий ортонормированный

вейвлет-базис, образованный функциями / с N нулевыми моментами, тогда

(1 + ^N

то(4) =

б(4), бп) * о,

2

где в случае использования функций с компактным носителем Q(£)-

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

масштабирующую функцию с маской ш0,

ад

ф(* = СПто (2-4),

1=1

где константа С определяется нормировкой функции ф, и классическую формулу

Псо8 (2-'4):

М £

можно показать, что маски, соответствующие вейвлетам а]к = 1^1/1к, Ркк = Ь/1к, где ь = аш/ахш , выражаются следующим образом:

=

тв *) =

1 + е

* Л

2

2е

У

•е \К

- 4 \

1 + е

N+т

б (4),

V1 + е-*,

то(4)= е

__ -іт*

С1 + е-*

б (4).

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

4 -2 0 2 4 -4 -2 О 2 4 _4 _2 о 2 4

а б с

Рис. 1. Масштабирующие функции (верхний ряд) и вейвлеты (нижний ряд), адаптированные к

оператору а2/Ох2 , построенные на основе вейвлета Добеши 5-го порядка. Вейвлет Добеши (а);

биортогональные вейвлеты а (б) и в (с)

-6 -4 -2 0 2 4 -6 -4 -2 0 2 4 -4 -2 0 2 4

а б с

Рис. 2. Масштабирующие функции (верхний ряд) и вейвлеты (нижний ряд), адаптированные к оператору О4/Ох4 , построенные на основе вейвлета Добеши 10-го порядка. Вейвлет Добеши

(а); биортогональные вейвлеты а (б) и в (с)

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

Б/ = 21/2Д2-),

тогда

ат ат

—Б1/(х) = ^Б1—/(х), ш,1е Z . дхт ахш У’

Таким образом, матрица жесткости при использовании алгоритма анализа-синтеза в

силу того, что действие алгоритма сводится к масштабированию

продифференцированных (проинтегрированных) вейвлетов, имеет диагональные

элементы вида 2-1 и решение уравнения запишется в виде

и = 221т{/,ак)ак .

1,к

Примеры решения уравнений

Приведем примеры решения уравнений 2-го и 4-го порядка с разрывными правыми частями. На рис. 3 приведены решения уравнений а2и/ах2 = /(х) и а4 и/ Ох4 = / ( х), х е [0,1] с периодическими граничные условиями

и(п)(0) = и(п)(1), п = 0,1,...,т -1,

условием совместности и условием однозначной разрешимости

1 1

| / (х)ах = 0, | и( х)ах = 0.

0 0

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

масштаба); нижний ряд: решения уравнений

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

Заключение

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

Первый автор (В.Г. Захаров) выражает благодарность за финансовую поддержку Американскому фонду гражданских исследований и развития, грант № PE-009-0.

Библиографический список

1. Dahlke S., Weinreich I. Wavelet-Galerkin methods: an adapted biorthogonal wavelet basis// Constr. Approx.-1993.-Vol.9.-P.237-266.

2. Jawerth B., Sweldens W. Wavelet multiresolution analyses adapted for the fast solution of boundary value ordinary differential equations. In N. D. Melson, T. A. Manteufel, and S. F. McCormick, editors, Sixth Copper Mountain Conference on Multigrid Methods, NASA Conference Publication 3224. - 1993. - P.259-273.

3. Lazaar S., Ponenti Pj., Liandrat J., Tchamitchan Ph. Wavelet algorithms for numerical resolution of partial differential equations // Comput. Methods Appl. Mech. Engrg. -1994.-Vol.116.-P.309-314.

4. Sweldens W. The lifting scheme: A construction of second generation wavelets. Technical Report 1995:6, Industrial Mathematics Initiative, Department of Mathematics, University of South Carolina, 1995.

5. Tchamitchian Ph. Biorthogonalite et theorie des operateurs // Rev. Math. Iberoam.-1987.-Vol.3.-P.163-189.

Получено 5.05.2003