- •Указания по выполнению и защите практических заданий
- •Листинг программы
- •Интерфейс программы
- •Ввод
- •Вывод
- •Построение графиков
- •Нелинейные уравнения и системы
- •Уравнения с одним неизвестным
- •Системы нелинейных уравнений
- •Минимизация функций
- •Минимум функции 1 переменной.
- •Многомерная минимизация.
- •Вычислительные задачи линейной алгебры
- •Прямые методы для задач линейной алгебры
- •Итерационные методы решения СЛАУ
- •Алгебраическая проблема собственных значений
- •Приближение функций
- •Функции одной переменной
- •Обратная интерполяция (ИМН, ИМЛ)
- •Наилучшее среднеквадратическое приближение
- •Функции многих переменных
- •Численное дифференцирование
- •Численное интегрирование
- •Приложения
- •Правила оформления листинга программы
- •Нормы векторов и матриц
- •Нормы векторов
- •Нормы матриц
Результат. Оценить точность приближения, вычислив
v |
10·Nx 10·Ny |
f(a + i |
|
hˆx,c + j |
|
hˆy) |
|
P(a + i |
|
hˆx,c + j |
|
hˆy) 2 |
|||
u |
|
|
|
|
|
· |
|
· |
|
− |
|
· |
|
· |
|
u i=0 j=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
u |
X X |
|
|
|
|
|
|
|
|
|
|
|
|
||
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где hˆx = |
b−a |
, hˆy = |
d−c |
. |
|
|
|
|
|
|
|
|
|
|
|
|
10·Nx |
|
10·Ny |
|
|
|
|
|
|
|
|
|
|
Методы решения:
1.Треугольная интерполяция (многочлен Ньютона)
2.Последовательная интерполяция (многочлен Лагранжа)
3.Последовательная интерполяция (Сплайн)
4.Среднеквадратичное приближение
Варианты функций. Необходимо выполнить задание на интерполяцию для обеих функций (f(x,y) и g(x,y)) из указанной преподавателем системы нелинейных уравнений (см. стр. 11).
7Численное дифференцирование
Практические задачи на численное дифференцирование можно условно подразделить на две категории:
1.Аппроксимация производной на отрезке.
2.Аппроксимация производной в точке.
Впервом случае результатом будет являться некоторая непрерывная функция, а во втором случае — одно-единственное значение в заданной точке или же полная таблица для всех узлов исходной дискретной величины.
23
Этапы выполнения:
Для методов 1-й категории
Задание. Для заданных метода и функции f(x), построить графики: самой функции f(x), аналитически вычисленной производной, а также все графики интерполяционных многочленов.
Оформление. Требования к оформлению графиков см. в разделе «Приближение функций», подпункт «Функции одной переменной»
Результат. Поварьировать расположение отрезка интерполяции [a,b] и его длину. Сделать выводы.
Для методов 2-й категории
Задание. Для заданных метода и функции f(x), построить графики: самой функции f(x) и аналитически вычисленной производной. По формулам 2-го и 4-го порядков точности вычислить значение производной в точке x˜. Точность результата оценить и уточнить по формулам Рунге.
Оформление. Требования к оформлению графиков см. в разделе «Приближение функций», подпункт «Функции одной переменной». Значения малых величин (оценки погрешности) необходимо выводить на экран без каких-либо форматных ограничений21
Результат. Поварьировать расположение точки x˜ и длину шага между узлами разностной формулы h. Сделать выводы.
Методы:
1.Аппроксимация f0(x) на основе ИМН степени 1, 2 и 3
2.Аппроксимация f00(x) на основе ИМН степени 1 и 2
3.Аппроксимация f00(x) на основе ИМЭ степени 3 и 5
4.Вычисление f0(x˜) в заданной точке x˜
5.Вычисление f00(x˜) в заданной точке x˜
Варианты функций см. на стр.9.
21Т.е., например, оставить величину в виде 9.56761818095963E-0002, а не преобразовывать к виду 0.09568
24
8Численное интегрирование
«Численное интегрирование» — завершающая практическая работа семестрового курса «Вычислительная математика» и первого семестра годового курса «Численные методы». Выполнив предыдущие работы, студент должен был усвоить все основные требования к оформлению работы.
Не будем углубляться в подробности оформления и постараемся кратко сформулировать требования к выполнению данной работы:
1.Реализовать программно требуемый численный метод.
2.Проиллюстрировать в достаточной степени процесс решения.
3.Получить конечный результат.
4.Оценить, наглядно продемонстрировать точность полученного решения.
5.Сделать выводы.
Перечисленные требования оставляют довольно большую свободу выбора студенту. Однако, при написании и защите программы следует учесть, что от работы ожидается максимальная эффективность. Если по какому-либо из перечисленных пунктов работа выполнена недостаточно убедительно, то она будет возвращена на доработку с соответствующими замечаниями.
Методы решения:
b
R
1. Метод трапеций. Вычислить f(x)dx, оценить эффективный порядок
a
точности по Эйткену.
|
b |
|
|
|
|
|
|
|
2. |
Метод Симпсона. Найти Ra |
f(x)dx с заданной точностью. |
||||||
|
|
|
b |
|
|
|
|
|
3. |
Метод прямоугольников. Вычислить |
Ra |
|
( |
) |
|
, оценить эффективный |
|
порядок точности по Эйткену. |
f |
|
x |
dx |
|
b
R
4. Метод Грегори. Найти f(x)dx и уточнить по формуле Рунге.
a
5.Сравнение методов трапеций и Симпсона. Реализовать программно оба метода, оценить точность по Рунге и сравнить.
6.Метод статистических испытаний, 1-й вариант.
25
7. Метод статистических испытаний, 2-й вариант.
b d
RR
8.Метод ячеек. Вычислить f(x,y)dxdy, оценить порядок точности
a c
по Рунге.
RR
9. Метод ячеек. Вычислить f(x,y)dxdy по области D сложной формы.
D
10. Последовательное интегрирование. Вычислить |
f(x,y)dxdy по об- |
|||
|
ласти D сложной формы. |
|
|
D |
|
|
|
RR |
|
11. Метод статистических испытаний, |
2-й |
вариант. Вычислить |
||
|
b d |
|
|
|
|
R R f(x,y)dxdy. |
|
|
|
|
a c |
|
|
|
на |
Варианты функций: для заданий на |
f(x)dx см. стр.9; для заданий |
||
RR |
системы нелинейных уравнений |
|||
на |
f(x,y)dxdy см. функции, входящие в |
R |
|
|
стр.11.
y |
|
|
|
|
|
3 |
|
|
|
|
|
2 |
d |
|
|
|
|
1 |
|
D1 |
|
|
|
|
|
|
|
||
0 |
c |
|
|
|
|
|
a |
|
|
b |
|
-1 |
|
|
|
|
x |
-1 |
0 |
1 |
2 |
3 |
|
|
(a) Область D1 |
|
|
||
y |
|
|
|
|
|
3 |
|
R=(b-a)/2 |
|
|
|
2 |
d |
|
|
|
|
1 |
|
R |
|
|
|
|
D3 |
|
|
|
|
|
|
|
|
|
|
0 |
c |
|
|
|
|
|
a |
|
b |
|
|
-1 |
|
|
|
|
x |
-1 |
0 |
1 |
2 |
3 |
|
|
(c) Область D3 |
|
|
y |
|
|
|
|
|
3 |
d |
|
D2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
0 |
c |
|
|
|
|
|
a |
|
|
b |
|
-1 |
|
|
|
|
x |
-1 |
0 |
1 |
2 |
3 |
|
|
(b) Область D2 |
|
|
||
y |
|
|
|
|
|
3 |
d |
|
|
|
|
2 |
|
|
|
|
|
1 |
|
|
D4 |
|
|
0 |
c |
|
|
|
|
|
a |
|
|
b |
|
-1 |
|
|
|
|
x |
-1 |
0 |
1 |
2 |
3 |
|
|
(d) Область D4 |
|
|
RR
Рис. 1. Варианты сложных областей D для заданий на f(x,y)dxdy.
26