Алгебраическая многосетка: почему продукт интерполяции и ограничения не приводит к чему-либо с нормой 1?

12

В настоящее время я работаю с «Многосеточным учебником» Бриггса и др., Глава 8.

Конструкция оператора интерполяции имеет вид: введите описание изображения здесь

Тогда конструкция оператора ограничения и оператора точной сетки задаются как:

введите описание изображения здесь

Давайте предположим, что у нас есть три точки сетки x0, x1, x2 со средней, x1 в порядке, а остальные грубые. Средний интерполируется x1 = x0*w0 + x2*w2. Следовательно, оператор интерполяции (в Matlab):

I = [1, 0, 0; w0, 0, w2; 0, 0, 1]

I =

[  1, 0,  0]
[ w0, 0, w2]
[  0, 0,  1]

Оператор ограничения тогда:

transpose(I)

ans =

[ 1, w0, 0]
[ 0,  0, 0]
[ 0, w2, 1]

Теперь давайте посмотрим, что произойдет, если кто-то ограничит, а затем интерполирует напрямую, что приведет к умножению Iи transpose(I):

I*transpose(I)

ans =

[  1,          w0,  0]
[ w0, w0^2 + w2^2, w2]
[  0,          w2,  1]

Я ожидал бы, что эта матрица является чем-то вроде единичной матрицы или, по крайней мере, имеет норму 1 или что-то в этом роде. Но если бы мы применили x = [1, 1, 1] для, скажем, w0 = w2 = 0,5, мы получили бы [1,5 1,5 1,5]. Я бы предположил, что неоднократно примененные операции ограничения-интерполяции хотя бы сходятся к чему-то. Но нет, в этом случае все компоненты вектора умножаются на 1,5 при каждой ограничительной интерполяции. Это кажется мне очень странным.

Кто-нибудь может объяснить, что происходит?

Майкл
источник
1
I3×32×3

Ответы:

8

x

Джесси Чан
источник
Я это понимаю. Но я бы, по крайней мере, предположил, что неоднократное применение ограничений и интерполяции с чем-то сходится. Но нет - в приведенном выше случае все элементы вектора будут умножаться на 1,5 для каждой ограничения-интерполяции. Это кажется странным для меня.
Майкл
Конечно - несколько коротких ответов. (1) Сглаживание и нормализация не принимаются во внимание, которые обычно применяются в сочетании с интерполяцией / ограничением. (2) Отчасти это может быть выбор весов. Они часто соответствуют различным вариантам операторов интерполяции / ограничения, некоторые из которых приводят к доказуемо лучшему поведению, чем другие. (3) Существуют другие операторы interp / restrict, для которых interp + restrict является проекцией. Например, вы можете делать глобальные проекции от грубой до точной сетки, но это дорого и не стоит для решателя.
Джесси Чан
4

ITI

EEIEI(EI)(EI)=EI

Вольфганг Бангерт
источник