Главная » Просмотр файлов » Тихонов, Свешников - Дифференциальные уравнения

Тихонов, Свешников - Дифференциальные уравнения (947323), страница 33

Файл №947323 Тихонов, Свешников - Дифференциальные уравнения (Тихонов, Свешников - Дифференциальные уравнения) 33 страницаТихонов, Свешников - Дифференциальные уравнения (947323) страница 332013-09-15СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 33)

Точка лз может быть либо граничной, либо внутренней точкой 10, Л. Пусть ха — граничная точка. Но тогда М=О, поскольку уэ(0) уз(0=0. Допу-' стим теперь, что лз — внутренняя точка, тогда в этой точке име-. ется максимум и у"(хз) ~О. Но нз самого уравнения (см. (6.89) прп 1(х) = О) следует О ~~ У > (и ) = Ч ( га) Уа (за) = Ч (ха) М ' з О, (6. 92) а это может выполняться ли>пь при М = О.

Итак, М= О. Аналогичвым образом легко показать, что и минимальное и> значение функции уо(л) на замкнутом отрезке 10, Л равно нулто: п>=0. Отсюда следует, что уз(х) =О, х>в10, Л, что и доказывает утверждение об единственности решения краевой задачи (6.89), (6.90). Перейдем к построению разностной схемы, аш>роксимирующей задачу (6.89), (6.90). Введем на отрезке 10, Л сетку юз а рассмотрим сеточную функцию '"'у, определенную в узлах сетки.

Вторую производную у" (х) аппроксимируем разностным отно- шением КРАЕВЫЕ ЗАДАЧИ 473 родной системы достаточно установить отсутствие нетривиааьных решений у соответствующей однородной системы. Предположим, что однородная система имеет нетривиальное решение '"'уа О, 'а'уп .„., '"'у„=0. Наибольшее из этих и+ $ чисел обозначим через М: 44> -М ~ч (6.95) Так как '"'уа='"'у„=О, то М~О. Будем теперь рассукдать, как при доказательстве отсутствия нетривиального решения у однородного уравнения (6.89).

Пусть М достигается в одном из граничных узлов. Тогда М= О. Если же М достигается во внутреннем узле 4= й, то в силу однородного уравнения (6.94) имеет место равенство (2+а~-'"'Ч)'"'у =(2+Йа и'Ч)М='"'у +'"'у 4~2М (6.96) Если М) О, то отсюда получим противоречие, поскольку в силу условия дЫ )О имеем (2+на ° 'муа)М) 2М. Поэтому и в этом случае М = О. Итак, М= О.

Аналогичным образом легко показать, что минимальное т значение п+1 чисел '"'уа, 'а'уь ..., '"'у„ также равно нулю: т = О. Следовательно„все значения 'муа, '"'уп ..., '"'у„равны нулю, что означает отсутствие нетривиаль- ных решений у однородной системы линейных алгебраических уравнений (6.94). Итак, разностная схема (6.94) однозначно раз- решима для любой правой части 'а'~ при любом и. Перейдем теперь к изучению аппроксимирующих свойств схе- мы (6.94). Пусть функции д(х) н г(л) имеют вторые непрерывные производные на отрезке (О, Л. Тогда из (6.89) следует, что ре- шение краевой задачи (6.89), (6.90) обладает на отрезке [О, Л непрерывными и ограниченными производйыми до четвертого порядка.

Разлагая точное решение краевой задачи в строку Тей- лора и подставляя эти разложения в левую часть уравнения (6.94): ь' ьа — а~У(Х;) + ЬУ'(Л4) + — )/'(Ха) + — У'э(яа) + Аа 4 Ь + —., уню (ха + 64Ь) — 2у (ла) + у (ха) — йу' (х4) + — у" (ла)— Аа Ьа С у (ла) + 4 у'~т~ (ла — 6 44)( — у (ла) у (З4) (6.97) (О:а 6 ~ ~4, 0(6,(4), получим ь' у" (*4) — у(,.)у(,)+ — „(уг >( 4+О,Ь)+у< >(,— О,й)) = = ~ (ха), (6.98) откуда следует, что порядок аппроксимации схемы (6.94) равен 2. числвнныв мктоды [Гл.

е 174 Рассмотркм вопрос об устойчивости разностной схемы (6.94) по относпению к возмущениям граничных условий и правой части уравнения. Возмущенная задача имеет вид <лг <л>- <юус+> 2 [>с+ у<-> <л> <лг а' — Чс Ус <л> Уе <лг Ул <">)с+ б[">(с о е» Ел у= (6.99) Обозначив ошибку решения б [му [л>у [л>у и вычитая (6.99) из (6.94), получим в силу линейности йз [л> )б <л> 6 см [ б [и [ ьеб [л>л (6ЛОО) (6:101) 6' ус=ась б' у< еь Выберем то значение >'= й 0= 1, ..., п — 1), при котором абсолютная величина ошибки наибольшая: (6'л'у ! ~ В'л'у.! ([=1, ..., п — 1).

(6Л02) Из уравнения (6.101) имеем (2+ й~ ' с~>дл) ! б с~>ул! ~ (б с~>у ! .(- ! б < >у ! + йз(б с~>Д (6 103) Неравенство (6.103) только усилится, если заменить !б"'у,+с! и !б'л'Ул с! на максимальное значение В'"'Ул!, а Цб[~>)л! на максимальное значение 16 [">)[!. Тогда (6.103) дает лз .

<л> !б <л> ! .- лз[[б [л>~1 (6.104) Отсюда получим ц б <л>) ц [о. Й а так как В'л'ул! ( шах(!ес(, (ес!) (й= О, и), то окончательно Цз["')Ц Цб<">у)~шах$ ., Це»$, Це>! (СЦб<л>[7Ц„(6.105) ! [е.с> что н доказывает устойчивость схемы (6.94) по граничным условиям и правой части уравнения. Установлепные свойства аппроксимации и устойчивости схемы (6.94) позволяют на основании, теоремы ОЛ утверждать справедливость следующей теоремы. Теорема 6.3[ Если функции д(к) ~ 0 и )(х) дважды непреръгвно дифференцируемы на отрезке (О, Й, то семейство [л>у решений ктлввыв злдлчи разностной схемы (6.94) сходится при >»- 0 к точному решени»о задачи (6.89), (6.90), причем порядок сходимости равен 2.

Остановимся теперь на вопросах реализации схемы (6.94). В отличие от рассмотренных выше разностных схем решения начальной задачи (схемы Эйлера, Рунге — Кутта), данная схема не дает явного алгоритма последовательного вычисления значения сеточной функции в узлах сетки, а представляет собой систему линейных алгебраических уравнений, в которую входят неизвестные значения сеточной функции во всех узлах сетки.

Поэтому для численного решения задачи (6.94) можно воспользоваться общими методами решения линейных алгебраических систем. Однако в случае достаточно болыпого порядка и эти методы оказываются весьма трудоемкими. Естественно воспользоваться специальными свойствами матрицы полученной алгебраической системы. В данном случае матр»ща система трехдиагональная— лишь па главной диагонали н двух побочных диагоналях элементы матрицы отличны от яуля. Это позволяет предложить для решения системы (6.94) весьма эффективный специальный метод, к изложению которого мы и перейдем. 3.

Метод алгебраической прогонки решении линейных алгебраических систем с трехднагональной матрицей. Мы рассмотрим этот метод для следующей системы, являющейся некоторым обобщением схемы (6.94): А»у»» — С»у»+В,у,„,=р» В=1, ..., и — 1), (6.106) ус=ау»+р» у =(у -»+б. (6.107) Пусть коэффициенты, входящие в уравнения (6.106) и граничные соотношения (6ЛО7), удовлетворяют условиям Аь Вь С,)0; С»)А+Вб 0<а<1; О<у<1.

(6.108) Очевидно, в рассмотренном выше случае схемы (6.94) условия (6.108) выполнены. Так же, как и в случае схемы (6.94), нетрудно показать, что при выполнейин условий (6Л08) задача разрешима н притом единственным образом. Будем искать такие коэффициенты а», 6», чтобы длн всех значений индекса 1= 1, 2, ..., и имело место соотношение у» — =ау+1». (6.109) Подставив нскомын вид решения (6Л09) в уравнения (6Л06), получим (А»а» — С»)у»+ В;у»ь»+ (А»р»+ Г») =О. (6.110) Выражая у» через у,+» по формуле (6Л09), перепишем (6.110) в ваде НА»ов — С»)а,+, +В)у, »+ ПА໠— С»)р; >+А»()»+Р»! =О. (6111) !гл. а численные метОды 176 Для того чтобы (6ЛИ) удовлетворялось тождественно при 1 1.

2, ..., и — 1, достаточно потребовать, чтобы каждая иэ квадратны скобок обращалась в нуль при 1= 1, 2, ..., и — 1. Это требование дает рекуррентные соотношения для последовательного определения по заданным вначенинм а=а~ и р=р1 «прогоночных» коэффициентов со+1 и ()ы1 («примак прогонка»): В~ А,.р +г «+х С. — А.а." с1+х С вЂ” А а. $ $ 3 1 3 Легко показать, что при выполнении условий (6.108) знаменатели в формулах (6.И2) отличны от нуля и, более того, О<сц~1. Действительно, перепишем второе условие (6ЛО8) в виде С,=А,+)),+Пь В,~О.

(6.ИЗ) Тогда первую формулу (6.И2) можно эаписать в виде а. В '+г В«+ А (х — а,.)+Р ' (6. И4) Так как 0 а~ ( 1, А, В,) 0 и Р, ~ О, отсюда следует, что и для всех аь вычисляемых по формуле (6.И4), вьшолняетсн условие О--. а1( 1.

Итак, все знаменатели в формуле (6.И2) строго больше нуля, что гюэволяет найти все прогоночные коэффициенты. Вычислив а„и р„вз (6Л09), получим р~-1 а р +~п (6.И5) С другой стороны, в силу граничного условия (6Л07) р„7р„1+6 7(а р +() )+б, (6.И6) откуда 7)) +б и — 1 (6. И7) Выше мы покаэалн, что 0 а с1. Поэтому в силу условия (6ЛО8) (0(7 с 1) знаменатель и этого выражения отличен от нуля. Тем самым р„определено.

Теперь по формуле (6Л09) можем последовательно вычислить все неизвестные р;-1 П = и, и — 1, ..., 1) («обратная прогонка»). Замечания. 1. В силу устало»ленного соотношения 0«а; ( «нн прн првмой, нн прн обратнон прогонке но происходит накопления ошибок вычисления. 2. Мы рассмотрели схему метода прогонки, в которой сначала определяются прог«нотные ко»ффвлксвты прн переносе левого граннчвого условия, а »атом восстанавлвваотсв решенно по формуле (б.109). Очевидно, аналогвчным обравом может быть рассмотрена схема, в которой прогоночньте ко»ф. фнпнонты онрсдоляютсн прогонкой справа колено правого граннчного условна, а решение восстанавливается прогонкой слева направо. ГЛАВЛ 7 АСИМПТОТИКА РЕШЕНИЙ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ПО МАЛОМУ ПАРАМЕТРУ Необходимость развития прнблия~епньгх методов решения дифференциальных уравнений отмечалась выше.

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

Точность такой формулы имеет естественное ограничение (подробнее об этом см. ниже). Заметим сразу же, что численные н асимптотнческие методы не исключают, а взаимно дополншот друг друга. й $. Регулярные возмущения $. Понятие асимптотического представлении. В $5 гл. 2 подробно исследовался вопрос о зависимости решения нвчвяьнсй задачи от входящих в уравнение параметров. В наствищем параграфе мы положим (с = О, чего всегда можно добитьсн ааменой независимого переменного, н ограничимся скалярным случаем (у — скаляр, р — скаляр), что не принципиально, а делается лишь в целях краткости изложения, Итак, рассмотрим начальную задачу (7.1) Прн этом будем считать, что р изменяется в некоторой окрестности значения р = О, чего тоже можно добиться соответствующим выбором начала отсчета по оси р. 12 л.

н. тихонов в ас. 178 Асиьп1тотикА Рюкении по мАлому пАРАметРу 1гл. т В гл. 2 (теорема 2.7) было докавано, что при соответствующих условиях на правую часть (7 1) решение у(1, р) задачи (7 1) существует и является непрерывной функцией 1 и (1 на множестве 1ы(0, Т), !)1! Сс. Доказанная теорема заключает в себе следующу1о возможность , построения приближенного решения задачи (7.1). Рассмотрим задачу, которая получается из (7.1), если в ней формально поло7кить р = О: '" = 1(у, 1, О), у (О) = у . (7.2) Задача (7.2), вообще говоря, проще исходной задачи (7,1), и ее решение, которое обозначим у(1), исследовать проще, а возможно, даже удастся эффективно построить. Из теоремы 2.8 следует, что на некотором сегменте (О, Т) у(1, (1) =у(1)+ з(1, 11), .(7.3) где е(1, (1) =. О при (1- О.

Характеристики

Тип файла
DJVU-файл
Размер
2,41 Mb
Тип материала
Учебное заведение
Неизвестно

Список файлов книги

Свежие статьи
Популярно сейчас
Почему делать на заказ в разы дороже, чем купить готовую учебную работу на СтудИзбе? Наши учебные работы продаются каждый год, тогда как большинство заказов выполняются с нуля. Найдите подходящий учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6418
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее