 |
РУБРИКИ |
Методические указания: Численные методы |
РЕКОМЕНДУЕМ |
||
|
Методические указания: Численные методы/p>
Задание № 21. Схематически построить графики функций. 2. Найти площадь фигуры, ограниченной этими графиками (с точностью ). Расчет точек пересечения заданных функций и расчет интегралов (любыми методами) оформить отдельными программами.
9 ЧИСЛЕННОЕ ДИФФЕРЕНЦИРОВАНИЕ 9.1 Некорректность операции численного дифференцирования Задача численного дифференцирования состоит в приближенном вычислении производных функции по заданным в конечном числе точек значениям этой функции. Разобьем отрезок на частей . Пусть определены значения функции приближенного значения можно взять любое из следующих разностных отношений , называемых, соответственно, левой, правой и центральной разностными производными. Возникающая в результате такой замены погрешность, называемая погрешностью аппроксимации, характеризуется оценками где погрешность аппроксимации левой и правой разностными производными является величиной при говорят, что имеет место аппроксимация первого порядка. Центральная разностная производная аппроксимирует со вторым порядком и, следовательно, является более точным приближением к чем левая или правая разностные производные. Вторую производную в точке можно заменить второй разностной производной при этом т.е. имеет место аппроксимация второго порядка. Как правило, значения функции в точках вычисляются не точно, а с каким-то приближением. Например, элементарные трансцендентные функции вычисляются с помощью рядов, причем ряды заменяются конечными суммами. Другим источником погрешностей являются погрешности округления. Оказывается, что погрешность, возникающая при вычислении разностных отношений, намного превосходит погрешность в задании значений функции и даже может неограниченно возрастать при стремлении шага к нулю. Поэтому операцию вычисления разностных отношений называют некорректной. Поясним причину некорректности на примере вычисления разностного отношения Разностное отношение хорошо приближает только в том случае, когда шаг достаточно мал. Требование малости величины , находящейся в знаменателе разностного отношения, как раз и является причиной некорректности операции численного дифференцирования. Действительно, пусть вместо точного значения вычислены приближенные значения Тогда вместо вычислена величина Следовательно, погрешность в вычислении первой разностной производной окажется равной тогда оценка достигается при Из этой оценки видно, что вследствие малости погрешность, возникающая при вычислении первой разностной производной, значительно превосходит погрешность вычисления самой функции . Если от неограниченно возрастает при Далее погрешность такого рода будем называть погрешностью округления. Сказанное не означает, что нельзя пользоваться формулами численного дифференцирования. Чтобы не происходило существенного понижения точности, надо следить за тем, чтобы погрешность округления имела тот же порядок, что и погрешность аппроксимации. Например, погрешность аппроксимации при замене отношением превосходит величины где потребовать, чтобы и погрешность округления была бы сравнима с погрешностью аппроксимации, например, где . Это означает, что погрешность при вычислении значений функции должна быть величиной С другой стороны, это неравенство показывает, что если величина задана, и мы не можем ее менять, то вычисления надо проводить не с произвольно малым шагом шагом, удовлетворяющим условию Например, если шаг примерно равным 0,01. При этом погрешность аппроксимации и погрешность округления будут примерно равными . Вычисление производной по заданной функции также является некорректной операцией в том смысле, что для ограниченной функции производная быть сколь угодно большой. Например, для имеем при 9.2 Применение интерполирования Многие формулы численного дифференцирования можно получить как следствие интерполяционных формул. Для этого достаточно заменить функцию ее интерполяционным многочленом и вычислить производные многочлена , используя его явное представление. Рассмотрим разбиение отрезка на и обозначим через шаги этого разбиения. В качестве примера получим формулы численного дифференцирования, основанные на использовании многочлена Лагранжа построенного для функции по трем точкам Многочлен Отсюда получим Это выражение можно принять за приближенное значение в любой точке удобнее записать в виде где и если разбиение равномерное, то приходим к центральной разностной производной, При использовании интерполяционного многочлена первой степени точно таким же образом можно получить односторонние разностные производные Далее, вычисляя вторую производную многочлена получим приближенное выражение для при При равномерном разбиении это выражение совпадает со второй разностной производной приближенного вычисления дальнейших производных уже недостаточно многочлена надо привлекать многочлены более высокого порядка и тем самым увеличивать число узлов, участвующих в аппроксимации. Порядок погрешности аппроксимации зависит как от порядка интерполяционного многочлена, так и от расположения узлов интерполяции. Приведем выражение для погрешности аппроксимации, возникающей при замене выражением где Отсюда видно, что аппроксимирует вторым порядком. Хуже обстоит дело с аппроксимацией второй производной: Здесь даже при равномерном разбиении второй порядок аппроксимации имеет место лишь в точке относительно других точек (например, точек и аппроксимация только первого порядка. 9.3 Контрольные вопросы 1. Что называется левой, правой и центральной разностными производными? Какой порядок аппроксимации обеспечивают разностные производные? 2. Почему операцию вычисления разностных отношений называют некорректной? 3. Как строятся формулы численного дифференцирования, основанные на применении интерполяционного многочлена? 4. Какой порядок аппроксимации обеспечивают эти формулы численного дифференцирования? 9.4 Задание к лабораторной работе № 9 1. Составить программу и вычислить на ЭВМ производную заданной функции на отрезке точностью полученных результатов с точными значениями производной . 2.Вычислить производную заданной функции используя интерполяционный многочлен. Сравнить с точными значениями производной.
Страницы: 1, 2 |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
одинаковых
где 
. В качестве
Таким образом,
В этом случае
где 
будет
Пусть 
причем эта
не зависит
, то погрешность 
не
Естественно
не зависит от 
, а с
, то
надо брать
может
и 
частей: 
имеет вид
Его
В частности, при 
получим
Для
со
а
) выполняется
с
. Сравнить точность
© 2010 |
|