- •Конспект лекций
- •Содержание
- •Введение.
- •1. Приближенные числа и действия над ними. Оценка точности вычисления.
- •2. Методы решения нелинейных уравнений.
- •2.1. Отделение корней
- •2.2. Общие свойства алгебраических уравнений.
- •2.3. Методы уточнения корней.
- •2.4. Решение систем нелинейных уравнений.
- •3. Решение систем линейных уравнений.
- •3.1. Классификация методов
- •3.2. Точные методы.
- •3.3. Схема единственного деления.
- •3.4. Вычисление определителя и обратной матрицы.
- •3.5. Приближенные методы.
- •4. Приближение функций.
- •4.1. Интерполирование и экстраполирование функций.
- •4.2. Аппроксимация.
- •5. Численное интегрирование.
- •6. Численное решение задач оптимизации.
- •6.3. Многомерная оптимизация.
- •Литература
2.4. Решение систем нелинейных уравнений.
При решении системы уравнений очень сложно отделить корни. Для получения всех возможных решений чаще всего на практике перебирают различные начальные условия поиска.
Метод итераций аналогичен соответствующему методу при решении одного уравнения. Все уравнения системы приводятся к виду:
x1=f1(x1, x2, …, xn)
x2=f2(x1, x2, …, xn)
…
xn=fn(x1, x2, …, xn)
Произвольные начальные значения подставляются в уравнения системы. Из каждого уравнения системы находятся новые значения и т.д. При выполнении условий или
последовательность сходится к решению
ЗАДАЧА 2.5. Решить систему уравнений
sin(x1-0,6)-x2=1,6
3x1-cos x2=0,9
E – погрешность вычислений
РЕШЕНИЕ
СВЯЗЬ. Преобразуем систему к виду:
x2 =sin(x1-0,6) -1,6= f2(x1, x2)
x1=(cos x2)/3+0,3= f1(x1, x2)
Проверим условие сходимости:
и
Выберем по графику начальные значения х10=0,15 и х20=-2
x2i+1 =sin(x1i-0,6) -1,6
x1i+1=(cos x2i)/3+0,3
x1= x1i+1, x2= x2i+1 при /x1i+1 –x1i/<E , /x2i+1 –x2i/<E
3. Решение систем линейных уравнений.
3.1. Классификация методов
Методы решения
Точные
Приближенные
Метод Гаусса
Метод Крамера
Метод простой итерации
Метод Зейделя
3.2. Точные методы.
Точные методы позволяют для любых систем линейных уравнений Ах=В в принципе найти точные значения переменных. Но при использовании этих методов может появиться накапливающаяся вычислительная ошибка (не алгоритмическая!). В результате значения переменных могут быть получены с некоторой погрешностью. При высоких порядках системы погрешность может оказаться значительной. Ошибка в вычислениях неизбежно приводит к ошибке в результате. Для контроля вычислительной ошибки применяют специальные приемы.
Метод Гаусса заключается в последовательном исключении неизвестных. Метод сводится к двум этапам. 1 этап – приведение исходной системы уравнений с помощью эквивалентных преобразований к системе с верхней треугольной матрицей (к треугольному виду). 2 этап – последовательное нахождение всех переменных системы, начиная с последней..
Затем можно найти невязки δ=В-Ах и вычислить поправки ε к полученным значениям переменных, решив систему уравнений Аε= δ.
Метод Крамера заключается в нахождении значений переменных с помощью определителей
ЗАДАЧА 3.1.
Решить систему уравнений Ах=В методами Гаусса и Крамера
РЕШЕНИЕ.
1) Приведем систему к треугольному виду.
Найдем приближенные корни: x3=3; x2=2,02; x1=0,97
Вычислим невязки: δ1=3∙0,97+2∙2,02+2∙3 – 13 = 0,05; δ2=0; δ3=0,02
Решив систему Аε= δ, получим: ε1=0,03; ε2=-0,02; ε3=0,00014
Уточненные решения: x3=3; x2=2,02-0,02=2; x1=0,97+0,03=1
2) Вычислим определители системы и найдем корни по методу Крамера.
Задание 3.1.
Решить систему уравнений точными методами.
3.3. Схема единственного деления.
Для удобства вычисления по методу Гаусса производятся по схеме единственного деления. Процесс преобразования матрицы к треугольной называется прямым ходом, а вычисление значений неизвестных – обратным ходом.
Количество разделов прямого хода равно числу неизвестных системы уравнений. В раздел I схемы записываются коэффициенты при неизвестных; свободные члены, контрольные суммы и строчные суммы, равные сумме всех элементов строки. Последняя строка (b1j) получается делением первой строки раздела (a1j) на ведущий элемент (a11 –первый элемент раздела) :
Элементы раздела II вычисляются по формуле: aij(II)= aij(I)– ai1(I) · b1j . Последняя строка раздела вычисляется делением первой строки раздела на ведущий элемент.
Аналогично вычисляются элементы следующих разделов.
Обратный ход начинается с вычисления последнего неизвестного системы уравнений и заканчивается вычислением первого неизвестного, используя лишь последние строки каждого раздела.
Строчные суммы всегда равны суммам элементов своей строки (без контрольной суммы), Над контрольными суммами в каждой строке проделываются те же операции, что и над остальными элементами этой строки. При отсутствии ошибок в вычислениях контрольные суммы приближенно равны строчным.
При необходимости уточнения корней х(0) системы уравнений необходимо:
- вычислить невязки δ =В – Ах(0);
- выписать невязки в столбец ε схемы;
- считая столбец ε столбцом свободных членов, вычислить ε1, ε2, …εn как значения неизвестных
- найти уточненные значения неизвестных: х = х(0) + ε
ЗАДАЧА 3.2.
По схеме единственного деления найти корни системы уравнений
РЕШЕНИЕ
Раздел |
x1 |
x2 |
x3 |
Свободные члены |
Контрольные ∑ |
Строчные ∑ |
ε |
|
||
Прямой ход |
I
b1j |
3 2 2 |
2 3 2 |
2 2 3 |
13 14 15 |
20 21 22 |
20 21 22 |
0,05 0 0,02 |
||
1 |
0,67 |
0,67 |
4,33 |
6,67 |
6,67 |
0,017 |
||||
II
b2j |
|
1,66 0,66 |
0,66 1,66 |
5,34 6,34 |
7,66 8,66 |
7,66 8,66 |
-0,034 -0,014 |
|||
|
1 |
0,4 |
3,22 |
4,61 |
4,62 |
-0,031 |
||||
III b3j |
|
|
1,4 |
4,2 |
5,62 |
5,6 |
0,006 |
|||
|
|
1 |
3 |
4,01 |
4,01 |
0,00014 |
||||
Обратный ход |
1 |
1 |
1 |
3 2,02 0,97 |
|
|
0 -0,02 0,03 |
|
х3=3+0=3; х2=2,02–0,02=2; х1=0,97+0,03=1