Главная » Просмотр файлов » Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu)

Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu) (1160088), страница 82

Файл №1160088 Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu) (Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы) 82 страницаН.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu) (1160088) страница 822019-09-19СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Вск>ду в двльнгйшелс прелпалагается, что ршп ыатрнцы В равен 1 — г, в ранг матрицы В ранен г. Прежде чен отыскивать практически пригаднлле ыоп>ды решения задачи (1], абсучилс вопрос о чужтвитольнасчи решений краевых звдшс к разного рада вазыущсннял>. В качестве надели возьмем краовук> задачу у' — Ау =. О, А = савей Ву(О) = Ъ, Ву(Х) = с1. (2) Рассмотрим случай, когда все юбственнья:иишенив на>рицы А различны; в абсцеы слу иис ход расгуждеиий изменяется несущжтвенпо. Пусть Л = и .ь 1В> — собственные значения ьшт)ищы А, .упорядоченные в парнике возраслапня аь а е - соответствующие собственные векюры, причем ((е>((.= 1. Общее решение систаыы у' — Ау = О запишется в виде у(г) = 2с е ехр (Л .с) Собсп>енпые зна >сипя Л.

разде>ппя на три груш>ы, присваивая собственныы значениялс кажс(ой грузны сао>впгствующий верхний индекс. Сабственньсе значения Л>, для которых величина ехр((а>(Х) очень болыпая, абозначаеы как Л+, если а > О, н как Л,, если о < О. Остальные собг ' > > огненные значения, *г.е. те, для которых величина ехр((а (Х) не очень вешска, обазпачаеы как Ло. Соответственно снабдим верхнилш индексаыи >' собствевные векторы е>, суыыироваспсе по лщдексаы )Л сосггвегствующиы о этилс группаы, будем обозначать ~, ), ~ . Пусть 1л, 1", 1 число собственных значений в сопгветствусощих группах. срорыа зшшси решения (3) ставит в неодинаковое положение концы отрезка интегрирования> все функции ехр(Л в) в п>чке к = О имеют порядок 1, в то иреыя как 4 5.

Рбсужденае носгановок краевых зад ш Выпишем систему уравнений Ву(й) = Ь. Ру(Х) = с1, из которой следует определить постоянпыс с, гоотвеп:твуюшие отыскиваемому решению: о Ву(0) = ~ с, Ве. + ~ ссВео+ ) г~Ве ехр(-Л4Х) =Ь, а Ру(Х) = ~ с"Ре ехр(Л., Х) -1-~ сбРеоехр(Л~~Х)-~-) сГРе+ = с1. (4) Систему уравнений (4) можно зя~игать в г)юрис Сс= и; с = (гг',..,, с!~)т, й =-(б!,..., Ь! „, г!и..., й! ) здесь матрица С имеет клеточный вид: 1о в котором клетки записываются следующим образом: ((Вег)ч], 1<1<В', !<с!<1 — ~", ((Веге)ч], Г < ! <1 +1в, 1 < г! <! — г, НВ г)ч ч1( — ЛгХН, ! +1'<)<1, 1<4<1 — г, ((Ре )еехр(Л ХЦ, 1<!<1, 1<4<с, ((Ре.) ехр(ЛоХ)], Г < ! < Г 4 ус, 1 < 4 < с, ](Рее) ] 1-4!о <! <! ! « , С, = е Сг = Предполсоким, что !!г = г1егС ф О, е предсгавиы решение системы (4) н вьгле с = С 'й.

Согласно известным формулам злементы обратной ьягтрицы имеют вид (б) в точке г = Х одви из них очень Геггыпнс, другие очень малые. Удобнее другая форма записи общего решения: О У(х) = ) сг ег ехР(Л:г) +~ гоеосхР(Лет) + ~ г+егсхР(Лг(х — Х)) Глава 9. Чисвенные методы решения краевых задач 44й где Сц миноры матрицы С. При гфорыулированных ранее предположениях величины ехр( — Л~Х), ехр(Л Х) ни'пожно малы, п<нтол<у ни. чтожно малы элементы матриц С< и Ся н определ<пель близок к определителю / с< с" ,о о с с:,' Рас<ь<отрил< отдельно случаи: а) определи<ель <Зе ве мвл; б) определитель йо мвл, в чвсщнхзи равен нулях Поскольку среди эпеь<ецчов матрицы С нет больших, то, вследствие формулы (5), я случае а) элементы матрицы С < обычно ве очень большие.

Предпол<нкнм, что правые чанги тряпичных углоний Ь и б <юдержат некоторые <югрешиогти бЬ и бнй Пусть бй = (бЬ, бб)т, тогд«погрешность вектора с равна С <бй. Если взел<виты матрицы С < не очень велики, то нлияние погрешности бй на ком)и)<нцнгнти с будет не очень большим. Решение задачи янлиет«:я линейной комбннацие»1 с к<нффнциентами с. <шагаеыых е. ехр(Л.

х), ез охр(Л<х), е«'ехр(Л«(х — Х)). П<нтоыу погрсппюсть приближенного решения, явлиющая<я следствием погрешностей бЬ и бб, закже будет прис»о<смой. Заметим, что п«п<и вы<шпыввния я<к ят давольно неопределенный характер: «мазый», «очень малый», «небозыпой», «очень большой»; при ш<елнзе конкретной задачи исслепователь Лолжен сзм решать, нм:колько приемлем дзя него тот или иной поряжж рассматриваемых иванчин. В часгнссги, »с ш решается <мс<ема болыпо<о» порядка, ю при «умеренных звачсниях коэффициентов сисшмы и «не очень»<ало»<» опрелелителе <Ге возм<экно, что миноры, состоящие нз сумм прож<венею<в бовьшого числа эле»<ентое, окзжушя «недопустимо бозыпими».

Если определитель <Зо очень мал или равен нулю, то, вследствие равенства йейсб<дс = Ь среди элементов матрицы С < встретятся большие. ТЬгда малые возмущения правых частей граничных условий могут приводить к большим в<пмущениям коэффициентов сч а следовательно, и решения задачи. З 5. Обсуждевие постановок краевых заЛа < 447 Зная элементы матрицы С и габствснные зн<»<ения л<втряцы А, из -< полученных выше соотношений л<ожно палучитл довольно точную информацию о возмущении решения диф<реренцнальной задачи. Однако получение этой информации «амо по <тбе требует больш<»го обьема вычислег<ий; перенос этих построений на случай переменной матрицы А(х) потребует еще большего об'ьема вычи<шений.

Попьппемся поэтому получить критерии устойчивости решения к нозмущепиям б!» и б<) качествснногп харакзпра, требующие меныпей информации о зэ»<а<с, хотя, пожег быть, и иескш<ъко менее надежные. Таким критерием могут служить со<аношения ыежду чи<шами 1, 1, 1, ! — г н г. Среди элементов первых е 1 — г строк матрицы С ненулевые »пол<виты могут находиться в пер,е вых 1 41а столбцах, соответству<ощих матрицаы С<, Сс<. Если 1л > г, тг» Р +!" < ! — г, и тогда асс миноры порядка (1 — <) х (1- г), ле»кащис в первых 1 — г строках, обращшотся в нуль. Раскрываи апредалитпаь <Ле по парным 1-.< строкам, иолучаел«»э=-О.

Точно так же, если 1 >1 — <., то все л<нноры порядка < х г, лежащие и последних г строках, обрап<аются в нуль, поэтому Аю = О. Если 1» «., а ! < 1 — г, тп опрсдел»ггель <Ле окажет<я линейной комбинацией произведений элементов матриц В и В и к<юрдинаг собственных векторов е, причем коэффициентами при этих пр<»ж<ведениях будут произведения чи<пл ех!»(Лс<Х), ве очень Ухл»ыпих и не оченл ыаленьких по ыодулю. Можно принять г»п»сапзу, что этот определитель оказывается лпшыл< числом довольно редко.

Тогда решение задачи (2] будет мало чуиствительш» к веал<ущепияы 6Ь и бд правых частей граничных условий. Мы можал< гформулировать полученные иь<ноды в качесгне ске!<ук»г»<его прсдложшшя. Если 1< >г или ! >1 — г, то ре<аекисдифф<рекциалъкой зада <и сильно чувства<велько к возмущениям <Опек:с частей греки"<киы условий< если)ь < г, а1 <1 — т, пю, лпк и!<<зало, ре<аение задачи (2) буди« мало «!агтеителъно к измекекилм прн»ьп; частей гроте<кис <!слоеии.

Первую часть этого утверждения можно переформулировать еще в такой <)н»рл<е< длл .малой <рве<пег»гаельнасга<г .задали к еозл<ри1<»<г<лм .<ракички»: условий необходима, шлаби чис»ю независимых часпшмх ре<иекий е» ехр(Л<х), сильно !<<со<ба<их иа (О, Л) с растая< х, кс пр< восходило <исаа грал<никит условий ка щпеом конце, а числа частных ушлет<й е ехр(Л,х), сильно убъ<ва<ап!нх ка (О, Х) с росным х, ке превосходило <исаа граничных, условий <т левом конце. Эта формулировка ири определенных уточнениях может быть перенесена и на случай задачи (1) с переменной матрицей А(х). Строгая п< реформулировка этоге утнерждепия будет доиольно граыоздкой; однако если элементы матрицы А(:е] атноситп»<ьно плавно ысняютея па (О, Х), то при первоначальяом исследовании устойчивости задачи к возл<ущениям граничных условий зачастую можно огравичиться подсчшпм числа собственных значений матрицы А(х) с большим положительным и бальшим отрицательным значениями величины Х Ке Л» (х).

Глава 9. Чишегнные леетодье решенюе краевых задач 446 Краевую чэдешу называют хоро>но обуслоелеииои (хора шо ею<ее>ос>>ее>ее>>>1) шлв малые возмущения келчффвциснтои и правых чшггсй уравнении и грани еных условий >>рива>щт к столь же малыы па порядку изыенениям решения зады>ее. Более аккуратное определение хорошей обусловленности можно дать следующим образо>>.

Наряду с краевой ведшей (1) расс>>ее грим краевые задачи у' — (А(х) 4 6А(г))у = Г(ее) Ь ЬГ(г), ()) -У 6В)у(0) = Ь-У 60< ((3 У Ы3)у(Х) .= 4 + бб с пе очсеп, большой мерой вепмущения г = шах ()[ЗА(х)[[ -1 [[ЬГ(т)[[) т [[БВ[[+ [[Ы)[) + [[Ь(>[[ -:- [[бг([[. е< 'х Если для агах рошоний таких краевых задач выполняютгя исравснсгва пеах [[у(:г) — у(т)[[ < ВХс э«.х (6) с не очень бгшьшим значением постоянной М, го исходную задачу на>ываэуг хороию обуглоелснноц, в противном сяучве задачу нвзыввн>т плохо абйслоелсниаб. Минимальное значение 31(ге>), при котором неравенство (6) выполняется ири всех с < со (сэ > 0 фикгиранано)е иногда называют жсроб ебуглаолеюеосиш данной зв>иши (относите>п,но возмущений с нормой, не большей со).

Обусловленность задачи характеризует устойчивость решения к возмущениям исходных данных, напримор к нвгочпгхли задания ковффициеегп>в уравпепвя. 11оскальку пагрсшнагти от округления при вычислениях эквивалентны возмущениям коэффициентов исходного уравнения, то мера обусловленности характеризует и устаэчивгють чиешсинога решения к возьшжньем округлепияы при численном решении. Если известь>а ориентировочная оценка с возмущения кеоффие>иве>тон задачи и погрешность порядка МЯс допустима, та имеет смысл несо. средствепиае чиглепное решение задачи. Рагсматрим в «ачостве примера задачу Каши для сис>емы у, '= Э>, у> = у> при Эе(0) =.

1, э>(О) .= 1 на отршкг [О, 30). Собсгвшшые значения матрипы А=( ) равны Л1. Величина ехр( — 30) очоиь малая, а ехр(30) — очень Гюлылвя, 1 1, 1э = О, 1е = 1, е = О и 1 = 1 > г = О, Поэтому малые эшиушеиия граничного условия должны приводить к очень Гх>пылим изменениям решения. В данно г еучае проволивп имя «еми в эео пара>рафе построения не имеют стсбаго гмыг еа, песка.еьку и бгз них ягно, что, вследогвне сшп ного роста решений потопного линейного уравнения, погршпн<хть решения растет очень быстра. Пусть лля той же системы рассматривается задача с краевыми условееями Ьещ(0) + Ь>у>(0) = Ь, фв>(30) + д>Ш(30) = ф собственные веки>ры, <хютветсгвуюшие собственным значениям Л, = -1 и Ля = 1, равны соотвепхвенна 449 Ь б.

Алгоритггы регпеиии краевых звлач (!, — 1) г и (1, 1)ч, !г =, ! = ! — г. Восбвце говоря, слсцует ожидать, что задача будет ус юйчивв к возмугценивм Ьь и Ьг). Решение отмсквваеп л в виде с, ( ) ехр( — х) бг", ~ ) ехр(х — ЗО). Уравнения (4) имеют вил (Ьг — Ь)1+(Ь +Ь);+, гг( — ЗО)=Ь, (г1т — г!г)с, гхр( — ЗО)+ (ф ьг) )гз = И. Очггош б(г), .1- 4 ) — 4(Ь1 + и) ехр( — ЗО) 1 гь 4(ьг — Ь!) — 1~0! — гь!) ехр( — ЗО) г и тле гь = (Ь, — Ьз)(ф чцг) — (ф — !г)(ь~ -~-Ь~) схр( — 60) близко к бе = (Ьг — Ьз)(цг + Из).

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

Тип файла
DJVU-файл
Размер
6,05 Mb
Тип материала
Высшее учебное заведение

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

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