- •Міністерство освіти і науки України Національний технічний університет України "Київський політехнічний інститут"
- •Розрахунково-графічна робота з курсу «Обчислювальна математика»
- •Содержание
- •1. Вычисления с учетом погрешностей Необходимые сведения из теории
- •1.1 Абсолютная и относительная погрешности приближенных
- •Чисел и правило их записи.
- •1.2. Верные значащие цифры приближенных чисел.
- •1.3. Нахождение абсолютной погрешности по верным цифрам
- •1.4. Правило округления чисел
- •1.5. Правило записи приближенных чисел
- •1.6. Оценка влияния погрешностей аргументов на значение функции
- •1.7. Оценка погрешностей арифметических действий
- •Задание
- •2. Метод половинного деления Необходимые сведения из теории
- •2.1. Этапы приближенного решения уравнений с одним неизвестным
- •2.2. Отделение корней. Графическое отделение корней
- •2.3. Условия применения метода половинного деления
- •2.4. Алгоритм метода половинного деления
- •2.5. Условие окончания процесса деления при заданной допустимой погрешности
- •Задание
- •3. Комбинированный метод хорд и касательных
- •3.5. Условие окончания процесса вычислений при заданной допустимой погрешности
- •Задание
- •4. Уточнение корней уравнений методом простой итерации Необходимые сведения из теории
- •4.1. Методы отделения корней уравнений
- •4.2. Алгоритм построения итерационной последовательности, порождаемой уравнением )
- •4.3. Достаточное условие сходимости итерационной последовательности
- •4.4. Оценка погрешности n-го приближения к корню
- •4.5. Условие окончания итерационного процесса при заданной допустимой погрешности
- •Задание
- •5. Метод простой итерации приближенного решения систем линейных алгебраических уравнений Необходимые сведения из теории
- •5.1. Способы определения расстояния в пространстве Rn
- •5.2. Абсолютная погрешность числового вектора и его координат
- •5.3. Сходимость последовательности векторов в Rn
- •5.4. Приведенная система уравнений. Cпособы преобразования систем к приведенному виду
- •5.5. Построение итерационной последовательности
- •6. Интерполирование математических таблиц Необходимые сведения из теории
- •6.1. Табличная функция
- •6.2. Задача интерполирования табличной функции
- •6.3. Теорема о единственности задачи полиномиального интерполирования
- •6.4. Конечные разности таблиц
- •6.5. Первый и второй интерполяционные многочлены Ньютона. Оценка погрешностей интерполяционных формул Ньютона
- •6.6. Формула линейного интерполирования и способы оценки ее погрешности
- •6.7. Обратное линейное интерполирование
- •Задание
- •Квадратичное приближение табличных функций по методу наименьших квадратов Необходимые сведения из теории
- •Задача аналитического приближения табличных функций
- •Задача приближения по методу наименьших квадратов
- •7.3. Алгоритм построения наилучшего многочлена по данному методу
- •7.4. Уклонения, среднеквадратичное уклонение
- •Задание
- •Задание по вариантам
- •Порядок выполнения работы
- •8. Приближенное вычисление определенных интегралов Необходимые сведения из теории
- •8.1. Численный метод приближенного вычисления определенных интегралов
- •8.2. Квадратурные формулы трапеций и Симпсона
- •8.3. Строгая оценка погрешностей этих формул
- •8.4. Оценка погрешностей методом двойного пересчета
- •8.5. Определение шага разбиения отрезка интегрирования, при котором квадратурная формула обеспечивает заданную точность
- •8.6. Вопросы оценки точности приближенного интеграла с учетом вычислительных погрешностей
- •Задание
- •Численное решение обыкновенных дифференциальных уравнений методом эйлера-коши Необходимые сведения из теории
- •Общее и частное решения обыкновенного дифференциального уравнения 1-го порядка. Интегральные кривые
- •Задача Коши. Теорема Пикара
- •Геометрический смысл правой части дифференциального уравнения, разрешенного относительно производной
- •Понятие численного решения. Ломаная Эйлера
- •Метод Эйлера- Коши, его геометрический смысл
- •Оценка погрешности численного решения методом двойного пересчета
- •Задание
- •Порядок выполнения работы
- •Использованная литература:
3. Комбинированный метод хорд и касательных
Необходимые сведения из теории
3.1. Отделение корней уравнений аналитическим способом
Пусть
на
определена функция
f,
причем:
1) f непрерывна на ;
2)
на
cуществует
производная
,
сохраняющая на нем знак;
3)
.
Тогда
найдется точка
такая, что
3.2. Условия, при которых для уточнения корней применяются методы хорд и касательных
Т.к.
метод хорд дает недостаток, а метод
касательных – избыток, применяют данный
метод,
для которого
значения
и
функции f
на концах отрезка имеют разные знаки
и, кроме того, функция f
имеет непрерывные производные
и
с отличными от нуля и сохраняющими
постоянный знак значениями при всех
.
3.3. Правила выбора начальных приближений для методов хорд и касательных
а) для метода касательных выбираем за
тот из концов отрезка
,
в котором выполняется условие
, то есть
и
одного знака;
б) начальные приближения корней:
по методу касательных
;по методу хорд
.
первое приближение корня
.
3.4. Алгоритм уточнения корней комбинированным методом хорд н касательных
Вычислить значения f(a) и f(b);
Проверить выполнение условия
если условие не выполняется, неверно
выбран отрезок
Найти первую и вторую производные f;
Проверить постоянство знака производных на отрезке
если
условие не выполняется, он выбран
неверно;Выполнить действия, указанные в п.3.3.;
Проверить условие окончания процесса, если задана погрешность; если оно не выполняется, то отрезок изоляции корня сужается и имеет вид
.
Приближенные корни находят по формулам
3.5. Условие окончания процесса вычислений при заданной допустимой погрешности
Для того чтобы найти приближенное значение корня с точностью до ε > 0, необходимо остановить процесс половинного деления па таком шаге n, на котором отрезок [xn1 ; xn2], будет иметь длину
и
вычислить
Тогда
можно взять
,
причем
.
Проверяется
условие
,
где ε – заданная точность. Вычисления
продолжаются до тех пор, пока не будет
найдено такое значение
,
при котором
и
не совпадут с точностью ε.
Задание
Отделите аналитически один из корней данного уравнения и определите его с точностью до ε = 0,5 · 10-5 комбинированным методом хорд и касательных.
Вариант |
Уравнение |
10 |
|
1. Отделим корни уравнения аналитически и выберем один из отрезков изоляции, на котором выполняются условия применимости метода.
Нули
:
Тогда при
функция возрастает, а при
– убывает. Примем за
.
Так
как
на
отрезке
существует один корень.
отрезок
выбран верно.
2. Возьмем соответствующие начальные приближения и найдем вручную первые приближения. Проверим условие окончания процесса вычислений.
Проверим
условие
:
3. Составим программу уточнения корня с точностью до ε, которая выводила бы приближения к корню, найденные методами хорд и касательных, и расстояния между ними, в таблицу.
Текст программы:
function f(var x:real):real;
begin
f:= 2*X*X*X+9*X*X-21;
end;
function f1(var x:real):real;
begin
f1:= 6*X*X+18*X;
end;
function f2(var x:real):real;
begin
f2:= 12*X+18;
end;
var
a,b,e,e1,x0,x11,x12:real;
n:integer;
begin
a:=-4.5;
b:=-3.5;
e:=0.000005;
x0:=a;
n:=1;
if f(a)*f(b)>0 then writeln('otrezok vibran neverno')
else
begin
x11:=x0-f(x0)/f1(x0);
x12:=a-((b-a)*f(a)/(f(b)-f(a)));
e1:=(x11+x12)/2;
writeln ('x11=',x11:1:6,' x12=',x12:1:6,' iteraciya #',n, ' x=',e1:1:6);
while abs(x11-x12)>=2*e do
begin
a:=x11;
b:=x12;
x11:= a-F(a)/f1(a);
x12:= a-((b-a)*f(a)/(f(b)-f(a)));
e1:=(x11+x12)/2;
n:=n+1;
writeln ('x11=',x11:1:6,' x12=',x12:1:6,' iteraciya #',n, ' x=',e1:1:6)
end;
writeln ('x=',e1:1:6)
end;
readln
end.
Результат, выведенный на экран:
x11=-3.981481 x12=-3.642857 iteraciya #1 x=-3.812169
x11=-3.786967 x12=-3.736732 iteraciya #2 x=-3.761849
x11=-3.756286 x12=-3.755065 iteraciya #3 x=-3.755676
x11=-3.755532 x12=-3.755531 iteraciya #4 x=-3.755531
x=-3.755531
4. Найдите приближенный корень и выпишите его с верными значащими цифрами.
Ответ: на отрезке существует только один корень данного уравнения, который был вычислен с точностью до ε = 0,5 · 10-5:
.
