страница 1
|
|||||||||||||||||||||||||||||||||||||||||||
Похожие работы
|
Устойчивость задачи коши для многослойных - страница №1/1
УДК 517.55+517.929.4 УСТОЙЧИВОСТЬ ЗАДАЧИ КОШИ ДЛЯ МНОГОСЛОЙНЫХ ЛИНЕЙНЫХ НЕОДНОРОДНЫХ РАЗНОСТНЫХ СХЕМ И АМЕБЫ АЛГЕБРАИЧЕСКИХ ГИПЕРПОВЕРХНОСТЕЙ Рогозина М.С. научный руководитель д-р физ.-мат. Наук Лейнартас Е. К. Сибирский федеральный университет Разностные уравнения возникают в различных областях математики. Разностной схемой обычно называют разностное уравнение, аппроксимирующее исходное дифференциальное уравнение и дополнительные (начальные, граничные) условия. Одно из важнейших свойств разностной схемы – устойчивость. В теории Лакса [1] сходимость разностной схемы изучается в пространстве решений исходной дифференциальной задачи и теорема эквивалентности утверждает, что если исходная дифференциальная задача корректна и схема аппроксимирует эту задачу, то устойчивость необходима и достаточна для сходимости. В монографии [2] исследована устойчивость однородной двухслойной линейной разностной схемы с постоянными коэффициентами. Условие устойчивости здесь дается в терминах, связанных с понятием разностной функции Грина задачи Коши. В работе [3] к исследованию устойчивости многослойных разностных схем применяется теория амеб алгебраических гиперповерхностей и для многослойных однородных разностных схем получена формула для решения задачи Коши через ее разностную функцию Грина. Кроме того, дано необходимое условие устойчивости задачи Коши. Приведен пример, показывающий, что это условие не является достаточным. Сформулировано и доказано достаточное условие устойчивости. В данной работе исследуется случай неоднородных разностных схем. Введем необходимые обозначения и определения. Пусть Pδ,δn+1=α,β∈Aaαβδαδn+1β – полиномиальный разностный оператор, то есть A=α,β – конечное подмножество целочисленной решетки Zn+1, и δ=δ1, δ2, …, δn , α=α1, α2, …, αn и δα=δ1α1, δ2α2…, δnαn. В данной работе рассматриваются разностные операторы с постоянными коэффициентами, характеристический многочлен P(z,w) которых имеет вид Pz,w=Pmzwm+Pm-1zwm-1+…+P0z, где Pj(z) являются многочленами переменных z=(z1, z2, …, zn) и коэффициент при старшей степени Pm(z)≡1. Такие разностные уравнения возникают в теории разностных схем и ниже мы будем использовать терминологию этой теории. Рассмотрим задачу Коши для неоднородной m+1-слойной линейной разностной схемы вида δn+1m+k=1mPjδδn+1m-kfx,y=gx,y, (1) где Pj(δ) – полиномиальные разностные операторы с постоянными коэффициентами, g(x,y) – заданная функция и x=(x1, x2, …,xn)∈Zn. Задача Коши для уравнения (1) формулируется следующим образом: Определим две функции на Zn+1: φx,y=φyx, для x∈Zn и y=0, 1, …, m-1,0, для x∈Zn и y≥m; Обозначим через Π={(x,y)∈Zn+1:y≥0} полупространство в Zn+1. Теорема 1. Если f(x,y) решение задачи Коши (1) – (2), то для (x,y)∈Π справедлива формула fx,y=f0x,y+f*x,y, (3) где f0x,y=(x',y')μ(x',y')Px-x',y-y', причем суммирование проводится по всем точкам x',y'∈Zn+1, удовлетворяющим условию -m≤y'<0, а f*=(x',y')g(x',y')Px-x',y-y', где суммирование проводится по всем точкам x',y'∈Zn+1 удовлетворяющим условию 0≤y'≤m-1. При этом для любого фиксированного (x,y)∈Π число слагаемых для f0 и f* конечно. Замечание. Отметим, что f0 – решение однородной задачи Коши, f* – частное решение, с нулевыми начальными данными. Для доказательства теоремы 1 потребуется следующая лемма. Лемма. Пусть Κ – конус в R(x,y)n+1 порожденный векторами aj=(a1j, …, anj, an+1j), j=1, 2, …, m Κ=x,y∈Rn+1:x,y=j=1mλjaj, λj≥0. Если an+1j>0 для j=1, …, m, то пересечение конуса Κ с гиперплоскостью y=const является ограниченным множеством. В случае двухслойных разностных схем формула (3) была получена в монографии [2], а для задачи Коши в положительном октанте Z+n целочисленной решетки Zn в [4]. Приведем некоторые сведения из теории амеб алгебраических гиперповерхностей (см. [5]) Пусть V=λ∈Cn+1:Pz,w=0 – множество нулей многочлена P(z,w), оно называется характеристическим множеством. Определение 3. Амебой алгебраической гиперповерхности называется образ множества нулей V многочлена Pz,w=α,β∈Aaαβzαwβ при отображении Log:z,w=z1, …, zn,w→logz1, …,logzn,logw=Logz,Logw. Отметим, что множество V, а значит и LogV, замкнуто, поэтому его дополнение открыто. Оно состоит из конечного числа связных выпуклых компонент, и вершине (0,m) многогранника NP соответствует связная компонента E0.m дополнения амебы. Для произвольной функции ψ(x,y), заданной в полупространстве Π={(x,y)∈Zn+1:y≥0} определим ее норму следующим образом ψ=supψx,y. Определение 4. Назовем задачу (1) – (2) устойчивой, если существует константа L>0 такая, что при любых ограниченных начальных данных (2) и ограниченной правой части g(x,y) для соответствующего решения f выполняется неравенство f≤Lψ+g. Теорема 2. Пусть Eo,m связная компонента дополнения амебы характеристического многочлена P(z,w), соответствующая вершине (0,m) многогранника Ньютона. Задача Коши (1) – (2) устойчива тогда и только тогда, когда начало координат принадлежит Eo,m, то есть (0,0)∈Eo,m. Список литературы
|
|