§1. Учет погрешностейвычислений.
Прирешении математическихзадач могутвозникнутьпогрешностипо различнымпричинам:
При составленииматематическоймодели физическогопроцесса илиявления приходитсяприниматьусловия, упрощающиепостановкузадачи. Поэтомуматематическаямодель не отражаетреальный процесс,а дает егоидеализированнуюкартину. Погрешность,возникающаяпри этом, называетсяпогрешностьюпостановкизадачи.
Часто приходитсядля решениязадачи применятьприближенныйметод (интегралзаменяютквадратурнойсуммой, производнуюзаменяют разностью,функцию –многочленом).Погрешность,возникающаяпри этом, называетсяпогрешностьюметода.
Часто исходныеданные заданыне точно, априближенно.При выполнениивычисленийпогрешностьисходных данныхв некоторойстепени переходитв погрешностьрезультата.Такая погрешностьназываетсяпогрешностьюдействий.
Погрешность,возникающаяпри округлениибесконечныхи конечныхдесятичныхчисел, имеющихбольшее числодесятичныхзнаков, чемнадо в округлении,называетсяпогрешностьюокругления.
Определение.Пусть х – некотороечисло, числоа называетсяего приближеннымзначением,если а в определенномсмысле малоотличаетсяот х и заменяетх в вычислениях,
Определение.Погрешностью приближенногозначения ачисла хназываетсяразность
Если
Если
Определение.Границейпогрешностиприближенногозначения ачисла х называетсявсякое неотрицательноечисло
Говорят, чтоприближениеа приближаетчисло х с точностьюдо
Пример. Пустьа=0,273 – приближенноезначение х сточность до0,001. Указать границы,в которых заключаетсях.
Приокруглениичисел считают,что границыпогрешностиокругленияравна половинеединицы округляемогоразряда:
Определение.Относительнойпогрешностьюприближенногозначения ачисла х называетсяотношение
Пример.Округлитьдо десятыхчисло 27,52 и найтипогрешностьи относительнуюпогрешностьокругления:
Также каки абсолютнаяпогрешностьотносительнаяпогрешностьне всегда можетбыть вычисленаи приходитсяоценивать еемодуль. Модульотносительнойпогрешностивыражаетсяв процентах.Чем меньшемодуль относительнойпогрешности,тем выше качествоприближения.
Определение.Границейотносительнойпогрешностиприближенногозначения ачисла хназываетсявсякое неотрицательноечисло
Установимсвязь междуграницамипогрешностейабсолютнойи относительной:
§10. Вспомогательныесведения изфункциональногоанализа.
1-3– аксиомы метрики.
Говорят, чтомножествоэлементов
Последовательностьточек
Всякая сходящаясяпоследовательностьявляетсяфундаментальной,обратное верноне всегда.
Определение.Метрическоепространство,в котором всякаяфундаментальнаяпоследовательностьсходится называетсяполным.
Пример.
Зададимразличнымиспособамирасстояния:
кубическаяметрика, m-метрика
сферическаяметрика,
октаэдрическая,s-метрика
Для всехвыполняютсяаксиомы метрикии в каждой
Пусть X,Y– метрическиепространства.
Если X=Y,то
Если
Определение.Говорят,что отображение
§11. Решениеуравнений содним неизвестным.Дихотомия.
Пусть требуетсярешить уравнение
Число
Если функция
Отделитькорень уравнениязначит найтитакой интервал,внутри которогонаходится одини только одинкорень данногоуравнения.
Для отделениякорней можноприменитьследующийпризнак:
Если на отрезке
Достаточнымусловием монотонностифункции наотрезке являетсясохранениезнака производной.
Отделитькорень можнои графически:нарисоватьграфик и указатьточки пересеченияс осью Ох.
Совершенныйметод отделениякорней – методШтурма.
Дихотомия(метод деленияотрезка пополам).
Пусть
Рассмотрим
Если нужнонайти кореньс точностьюдо
Дихотомияпроста и оченьнадежна: к простомукорню она сходитсявсегда длялюбой непрерывнойфункции в томчисле и недифференцируемой,при этом онаустойчива кошибкам округления.Скорость сходимостиметода дихотомиине велика, т.е.за одну итерациюточностьувеличиваетсявдвое.
Недостатки:прежде чемприменить,необходимонайти отрезок,на концах которогофункция принимаетзначения разныхзнаков. Еслина этом отрезкенесколькокорней, то неизвестнок какому из нихсходится дихотомия.Метод не применимк корням четнойкратности.
Метод применимк корням нечетнойкратности, нохуже устойчивк ошибкам округления.Метод не применимк системамуравнений.
§12. Метод простойитерации длярешения алгебраическихи трансцендентныхуравнений.
ТЕОРЕМА1. (ПринципБанаха сжимающихсяотображений).
Пусть R– полноеметрическоепространство.Если
План доказательства.
q – коэффициентсжатия
Т.к. R– полноеметрическоепространство,то в нем всякаяфундаментальнаяпоследовательностьсходится.
ЧТД.
Метод
Еслив (*)
Методпростой итерацииимеет линейнуюскорость сходимости.
Пусть
Необходимопривести к виду
Условиесходимостидля данногометода:
ТЕОРЕМА2.
Пустьвыполняютсяусловия:
Функция
Дляначальногоприближения
Числа
Тогдауравнение
Теоремадоказываетсяаналогичнотеореме Банахас точностьюдо обозначений.
Замечание.УсловиеЛипшица применятьтрудно, вместонего применяютдругое условие:
Методитерация даетбесконечнуюпоследовательностьприближений,поэтому используютследующиеправила остановки:
задаетсяуровень останова
задаетсяуровень
Методпростой итерацииудобен в использовании,так как он легкопрограммируетсяна ЭВМ.
Недостаток:невысокаяскорость сходимости,т.е. линейная.
§13. МетодНьютона. Решениеуравнений содной переменной.
Пустьтребуетсярешить уравнение
Изэтих условийвытекает, чтона
Прежде,чем использоватьитерации, необходимо(1) привести квиду
Функция
Вкачестве
Выберемначальноеприближение
Метод, определенный(4), называетсяметодомНьютона.
Докажем, чтометод Ньютонасходится иполучим егооценку погрешности.
Еслидано, что |
|
|
|
|
Докажем,что (4) сходится.
Дляэтого покажем,что отображение
При
По непрерывностифункции
Поэтому котображению
Если выбрать
Заметим, чтометод (4) будетсходиться, еслиначальноеприближение
Докажем, чтометод Ньютонасходится.
Определимскорость сходимостиметода Ньютона.Для этого
При
Выразим
Потребуем,чтобы начальноеусловие
Тогдаиз (6) получим
МетодНьютона имеетквадратичнуюскорость сходимости.Это означает,что при переходеот одной итерациик другой количествоверных знаковудваиваетсяв последующемприближении.
Достоинство:высокая скоростьсходимости,легко программируетсяна ЭВМ.
Недостатки:узкая областьсходимости.
Если будемрешать операторноеуравнение
Геометрическийсмысл методаНьютона.
П
В точке
Если
Возьмем
Если
И так далее.Отсюда методНьютона называютметодомкасательных.
§14. Метод хорд.Метод секущих.
По прежнемурешаем уравнение
Т.е. на
Уравнение(1) запишем в виде
Тогда итерационныйметод
Докажем, чтометод хордсходится. Дляэтого необходимопоказать, что
Разложим
Рассмотримпри
Обозначимчерез
Т.е.
Следовательно,
Получим оценкупогрешностидля метода хорд
Таккак
Обозначимчерез
Сходимостьметоды хорд– линейная.
Достоинствометода хорд– легкостьпрограммированияна ЭВМ.
Общий видупростится:
При условии
При условии
Метод секущих.
Метод секущихимеет вид:
Скоростьсходимости– сверхлинейная.
Методсекущих сходитсябыстрее методахорд и методапростой итерации.
§15.Метод Гауссарешения системуравнений.
Длярешения системуравненийиспользуютметоды: точныеи приближенные.
Кточным относятся:
методГаусса;
методКрамера;
методоптимальногоисключения;
методквадратногокорня.
Кприближеннымметодам решениясистем уравненийотносятся:
методпростой итерации;
методЗейделя;
методНьютона.
МетодГаусса состоитв том, чтобыисходную системувида Ах=b(1)с произвольнойматрицей Асвести к системевида:
Процесссведения системы(1) к системе (2)называетсяпрямымходом методаГаусса.
Анахождениенеизвестных
Привычисленияхпо методу Гауссавелика вероятностьслучайныхошибок. С цельюизбежать ихвводится контрольныйстолбец:
Элементыконтрольногостолбца преобразовываютсяпо тем же формулам,что и элементыматрицы А.
Второйшаг контролясостоит в проверкеравенства суммыэлементовпреобразованнойстроки и контрольногоэлемента. Этивеличины должнысовпадать сточностью до1,2 единиц последнегоразряда.
МетодГаусса с выборомглавного элемента.
Средиуравненийвыбирают уравнение,содержащеенаибольшийпо абсолютнойвеличине коэффициент(главный элемент).
Затемуравнение делятна этот главныйэлемент и изостальныхуравненийсистемы исключаютнеизвестные,определяемыеэтим главнымэлементом.
Далее,оставляя неизменнымвыбранноеуравнение сглавным элементом,из остальныхуравненийсистемы выбираютновый главныйэлемент. Потомэто уравнениес новым главнымэлементом делятна новый главныйэлемент и исключаетнеизвестноеили определяемоеиз остальныхуравненийсистемы.
Дляудобства главныйэлемент помещаютв левый верхнийугол, переставляястроки и столбцысистемы уравнений.
Врезультатепреобразованийприходим кединичнойматрице.
Здесьпереставляютсяуравнения, чтоприводит кизменениюпорядка исключенныхнеизвестных,и во многихслучаях уменьшаютпогрешности,связанные сокруглениями.
§16. Методквадратногокорня.
Метод квадратногокорня – точныйметод решениясистем уравненийи он применяетсядля решениясистем уравнений,если матрицаА – симметричная,т.е.
где С – верхняятреугольнаяматрица;
D – диагональная,
Подставимматрицу А всистему (1)Ах=b.
Тогда
Выразимэлементы матрицыС через элементыисходной матрицыА.
Из(4) будем получатьвыражения
Пусть
Пусть
Пусть
Изформулы (*) получаем:
Получилиформулы:
§2. Оценкапогрешностейрезультатовдействий надприближеннымизначениямичисел.
(Строгийучет погрешности)
Пусть
Обозначимчерез
Утверждение1. Суммаграниц погрешностейприближенныхслагаемыхявляется границейпогрешностиих алгебраическойсуммы.
Доказательство:
ЧТД.
Утверждение2. Средиграниц относительнойпогрешностисуммы приближенныхслагаемыхсуществуеттакая, котораяне превосходитнаибольшейиз границотносительнойпогрешностислагаемых:
Утверждение3. Суммаграниц относительныхпогрешностейсомножителейявляется границейотносительнойпогрешностиих произведения:
Следствие1. При умноженииприближенныхзначений числана точный множительк, границаотносительнойпогрешностине меняется,а граница абсолютнойпогрешностиувеличиваетсяв
Следствие2. Произведениеграницы относительнойпогрешностиприближенногозначения ачисла х на
Следствие3. Частноеграницы относительнойпогрешностиприближенногозначения ачисла х и nявляется границейотносительнойпогрешностикорня n-йстепени из а:
Следствие4. Суммаграниц относительныхпогрешностейприближенныхзначений делимогои делителяявляется границейотносительнойпогрешностичастного.
§3. Приближенныевычислениябез учетапогрешностей.
Правило1. Для того,чтобы вычислитьалгебраическуюсумму приближенныхслагаемыхнужно:
среди слагаемыхвыбрать наименееточное (имеетнаименьшеечисло разрядовпосле запятой);
все остальныеслагаемыеокруглить,сохраняя одинзапасной разряд,следующий запоследнимразрядом выделенногослагаемого;
сложитьполученныепосле округлениячисла;
округлитьполученныйрезультат допредпоследнегоразряда.
Пример.S=2.737+0.77974+27.1+0.283
Определение1. Значащимицифрами вдесятичнойзаписи числаназываетсявсе его цифрыкроме нулей,записанныхслева от первойцифры не равной0.
0,00237 – 3 значащиецифры;
0,02000 – 4 значащиецифры.
Правило2. Для того,чтобы вычислитьпроизведение(деление) приближенныхчисел нужно:
выделитьсомножитель,содержащийнаименьшеечисло значащихцифр;
округлитьостальныесомножители,оставляя наодну значащуюцифру больше,чем в выделенномсомножителе;
произвестиумножение(деление);
округлитьполученныйрезультат,сохраняя столькозначащих цифр,сколько их ввыделенномсомножителе.
Пример.Р=3,34*0,7*4,748=4,7*3,3*0,7
Правило3. Привозведенииприближенногозначения вквадрат иликуб, при извлеченииквадратногоили кубическогокорня, в результатеследует оставлятьстолько значащихцифр, сколькоих имеет основание.
Правило4. Если числоявляется результатомпромежуточныхдействий, тоследует сохранитьв нем на 1-2 цифрыбольше, чемуказано в правилах1-3.
§4. Связь междучислом количестваверных цифр
и относительнойпогрешностью.
Пусть
Определение.Цифраприближенногозначения аназываетсяверной,если модульего погрешностине превосходитполовины единицыэтого разряда.
Очевидно,что все цифры,стоящие слеваот верной цифры– верные.
Пример.Пусть х=27,421, а=27,381,
Выясним, какиецифры верныев приближенииа?
4 –
8 –
3 –
3,2,7 – верныецифры.
Пусть известноколичествоnверных значащихцифр в приближенииа, тогда а запишем:
Так как цифра,стоящая в разряде-(n-1)верна, то погрешность
тогда
В качествеграницы относительнойпогрешностиможно взять
Итак, доказанатеорема 1.
Теорема1. Еслиприближениеимеет nверных значащихцифр, то число
Теоремаустанавливаетсвязь междучислами верныхзначений и егоотносительнойпогрешностью.
Замечание.Пусть приближениеимеет nверных значащихцифр и
Пример.
Итак, границаотносительнойпогрешностиприближенногозначения зависитот первой значащейцифры
Теорема2. Если границаотносительнойпогрешностиприближенияравна
Доказательство.Пусть
Из определенияследует, что–(m-1)– цифра, записаннаяв этом разрядеверная, цифры,записанныелевее тожеверные, то естьmверных цифр.
ЧТД.
Пример.Если известно,что относительнаяпогрешностьприближения
§5.Прямая задачатеории погрешностей
(функцииот приближенныхзначений аргументов).
Пустьфункция
Переменныезаданы своимиприближениями:
Известнапогрешностьэлементов
Предположим,что
Если
Пример.Вычислитьвеличину погрешностиприближенногозначения большегокорня уравнения
Вприближеннойзаписи используюттолько верныецифры, ????????????????????,обусловленныепогрешностьюприближенныхзначенийкоэффициентов.
Теперьобозначим
Рассмотрим
§6. Обратнаязадача теориипогрешностей.
Всезадачи теориипогрешностейделятся напрямые и обратные.
Прямая задача:определитьпогрешностьданной функцииот приближенныхзначений аргументов,заданных сизвестнойотносительнойпогрешностьюили с заданнойточностью.
Обратнаязадача:какими должныбыть относительнаяи абсолютнаяпогрешности,чтобы модульотносительнойили абсолютнойпогрешностизаданной функциине превышалзаданной величины.
Решениеобратной задачи.
Пусть
С какой точностью
Существуютразличныеподходы к решениютаких задач.
Принциправных влияний
заключаетсяв предположении,что погрешностивсех аргументоввносят одинаковыедоли в погрешностифункции, тоесть частныедифференциалыравны междусобой по модулю:
Предполагают,что погрешностивсех аргументовравны
Пример.С какой точностьюследует взятьдроби, чтобысумма Sмогла бытьполучена сточностью до0,001?
Решение.
Обозначим
1-й принцип
Сколько знаковпосле запятойнужно братьв дробях, чтобыполучиласьэта погрешность.Дроби необходимопредставитьв десятичномвиде та, чтобымодуль не превосходил0,00025, т.е. четырьмядесятичнымизнаками послезапятой.
§7. Метод границ.
Существуютразличныеспособы оценкиточности приближенныхвычислений:
строгий учетпогрешностей;
вычислениябез учетапогрешностей;
метод границ.
Метод границпозволяетустановитьграницы, в которыхнаходитсязначение, вычисляемоепо функции,если известныграницы, в которыезаключенызначения параметров,входящих вформулу.
х – число.
Теорема1. Суммаверхних границслагаемыхявляется верхнейграницей ихсумм. Сумманижних границслагаемыхявляется нижнейграницей ихсуммы.
Пример.
Теорема2. Разностьверхней границыуменьшаемогои нижней границывычитаемогоявляется верхнейграницей разности.Разность нижнихграниц уменьшаемогои верхней границывычитаемогоявляется нижнейграницей разности.
Доказательство.
сложим данныенеравенстваи получим результат.
ЧТД.
Пример.
Теорема3. Пустьнижняя границасомножителейнеотрицательна,то произведениенижних границсомножителейявляется нижнейграницей ихпроизведения,а произведениеверхних границсомножителейявляется верхнейграницей ихсомножителей.
Пример.
Теорема4. Если
Теорема5. Если НГделителяположительна,то частное ВГделимого и НГделителя являетсяВГ частногочисел; частноеНГ делимогои ВГ делителяявляется НГчастного
Доказательство.
Перемножими получим
ЧТД.
Пример.Вычислим значение
Действие | Содержимое | НГ | ВГ |
(1) | x | 2.57 | 2.58 |
(2) | y | 1.45 | 1.46 |
(3) | z | 8.33 | 8.34 |
(1)+(2) | x+y | 4.02 | 4.04 |
(1)-(2) | x-y | 1.11 | 1.13 |
| 9.24 | 9.43 | |
| 2.28 | 2.35 |
§8.Математическиемодели и численныеметоды.
Великароль математикив решении задачреального мира.Физиков математикаинтересуетне сама по себе,а как средстворешения физическихзадач. Один изспособов решениязадач: эксперимент.
Другойспособ: математическийанализ конструкцииили явления,однако такойанализ применяетсяне к самомуявлению, а кего математическоймодели. Математическаямодель физическогопроцесса представляетсобой совокупностьуравнений,описывающийпроцесс.
Математическаямодель должнаохватыватьважнейшиестороны явленияили процесса.Если математическаямодель выбранане точно, токакой бы мыспособ решенияне применили,результатымогут получитьсяне достаточнонадежными, аиногда и неверными.
1-ястадия работы:Постановказадачи.
2-ястадия работы:Математическоеисследование.
В зависимостиот сложностимодели применяютразличныематематическиеподходы, длянаиболее грубыхи наименеесложных моделейзачастую удаетсяполучитьаналитическоерешение (в видеформулы).
Длянаиболее точныхи сложных моделейаналитическоерешение удаетсяполучить крайнередко и тогдаприменяютчисленныеметоды решения,которые какправило требуютрасчета на ЭВМ.
3-ястадия работы.
Осмыслениематематическогорешения и егосопоставлениес даннымиэксперимента.Если решениехорошо согласуетсяс даннымиэксперимента,то такую модельможно применятьдля расчетапроцессовданного типа(модель выбранаправильно),если же решениеплохо согласуетсяс даннымиэксперимента,то такую модельнеобходимопересмотретьи уточнить.Численныеметоды являютсяодним из мощныхматематическихсредств решениязадач. Естьзадачи, где бездостаточносложных численныхметодов неудалось быполучить ответа.В современнойфизике такихзадач оченьмного, болеетого за короткоевремя нужнопровести огромноеколичествовычислений,иначе нет смысларешать задачу(суточный прогнозпогоды долженбыть просчитанза несколькочасов, а коррекциядвижения ракетыза несколькоминут). Это немыслимобез мощных ЭВМ,выполняющих1000000 операций всекунду. Современныечисленныеметоды и мощныеЭВМ позволилирешать задачи,о которых полвеканазад человекмог толькомечтать. Численныеметоды делятсяна точные иприближенные.Точные методыпозволяют законечное числоарифметическихдействий получитьрешение задачи.При этом еслиисходные данныезаданы точнои вычисленияпроизводилисьбез округления,то получаетсяточное решениезадачи.
Кточным методамотносятся:метод Гауссаи его модификации,метод Крамера,метод ортогонализациии т.д.
Приближенныеметоды (итерационные)дают бесконечнуюпоследовательностьприближений,предел которых,если он существует,является решениемзадачи. К итерационнымметодам относятсяметод Ньютонаи метод простыхитераций, методхорд и методсекущих длярешений уравнений.
§9.Понятие корректнопоставленной
инекорректнопоставленнойзадач.
Приприближенномрешении математическихили прикладныхзадач весьмасущественнымявляется вопросо том, корректноли решаемаязадача.
Большинствонекорректныхзадач записываетсяв виде уравненияпервого порядка,где по заданномунеобязательномуоператору
Определение.Задачаопределениярешения
Решениеустойчиво, еслибесконечномалым вариациямправой частисоответствуютбесконечномалые вариациих.
Определение.СледуяЖаку Адамарузадача отыскания
существуетв Х;
единственнов Х;
устойчивов Х.
Еслиже хотя бы одноиз условий 1-3не выполняется,то задачанекорректна.