книги из ГПНТБ / Пулькин, С. П. Вычислительная математика пособие для учащихся 9-10 классов по факультативному курсу
.pdfхи х2, х„, только вместо чисел, стоящих в столбце «Свободные члены», используем числа, стоящие в столбце «Суммы». Получив числа хи хг, х3, хи х2, х3 и сравнив их, мы осущест вим заключительный контроль.
При такой схеме вычислений мы можем осуществить и текущий контроль: в процессе вычислений число, стоящее в последнем столбце (столбец «Суммы»), должно быть равно сумме всех остальных чисел той же строки.
П р и м е р : Решить систему уравнений:
f |
10,2л:!-}-6,07х2 — |
91,4x3-i-50,3 |
= |
0, |
||
| |
9,28хх — 79,6х2—4,92х3 |
+ |
25,8 |
= 0, |
||
{ 68,Зл^ — 2,71х2— |
8 ,14х3 |
+ |
32,6 |
= |
0. |
Вычисления будем выполнять по описанной схеме, с контролем вычислений. Результаты вычислений отра жены в таблице XVIII.
56°. Точные и приближенные методы решения систем линейных уравнений. Среди различных методов решения систем уравнений будем различать точные и приближенные методы. Метод называется точным, если при точном вы полнении всех требуемых действий мы получаем точное решение системы. При этом имеется в виду, что все коэффи циенты — точные числа. Метод называется приближенным, если при точном выполнении всех требуемых действий и при точных коэффициентах мы получаем, как правило, лишь приближенный результат.
В случае применения приближенного метода на окон чательный результат влияет всегда погрешность метода, да к тому же на практике обычно имеет место неустранимая погрешность в задании коэффициентов и погрешность округления в промежуточных действиях.
Метод решения систем уравнений, основанный на про цессе исключения неизвестных, является точным. К числу приближенных методов относится метод итераций, приме нимый к решению систем. Несмотря на то что он является приближенным, он обладает такими качествами, которые делают его выгодным с вычислительной точки зрения.
57°. Метод итераций. Дадим здесь краткое изложение метода итераций. Этот метод заслуживает внимания, так как он удобен для вычислений на ЭВМ. Сущность метода итераций изложим на примере системы трех уравнений
130
|
|
|
|
Т а б л и ц а x v r n |
|
Номе |
|
Коэффициенты при: |
Свободные |
Суммы |
|
ра |
|
Xt |
*8 |
члены |
|
строк |
|
|
|
||
1 |
10,2' |
6,07 |
—91,4 |
50,3 |
—24,83 |
2 |
9,28 |
—79,6 |
—4,92 |
25,8 |
—49,44 |
3 |
68,3 |
—2,71 |
—8,14 |
32,6 |
90,05 |
4 |
— 1 |
—0,5951 |
8,961 |
—4,931 |
2,435 |
5 |
|
—85,12 |
78,24 |
—19,96 |
—26,84 |
6 |
|
—43,36 |
603,90 |
—304,19 |
256,35 |
7 |
|
—1 |
0,9192 |
—0,2345 |
—0,3153 |
8 |
|
|
564,04 |
—294,02 |
270,02 |
9 |
|
—1 |
—1 |
0,5213 |
—0,4787 |
|
З н а ч е н и я неизвестны х |
|
|
||
10 |
—0,4052 |
0,2447 |
0,5213 |
1 |
|
|
-1,4051 |
—0,7553 |
—0,4787 |
1 |
|
с тремя неизвестными. Итак, пусть дана система трех урав нений с тремя неизвестными. Для применения метода ите раций приведем систему к виду:
( |
х — Ахх-\- B1y-sr C1z-\ |
-D 1, |
|
|
< |
у = |
A2x-j- В 2у -\-С2г + |
^ 2> |
(5.10) |
I |
2 = |
А 3х + Bsy -f- С3г -j- D3. |
|
(такую систему будем называть системой нормального вида). Пусть (дг0, у о, 2о) есть некоторое грубо приближенное решение системы. Назовем эту совокупность чисел началь ным приближением. Если неизвестно какое-либо грубое приближение, то можно взять вообще любую произвольную совокупность значений неизвестных и принять ее за на чальное приближение. Часто за начальное приближение берут совокупность свободных членов.
5' |
131 |
Подставим числа х0, у0, 20в правые части системы, в ле вых частях получим некоторые значения х, у, г; обозначим их Xi, у1, гх. Совокупность чисел хи у,, zx называется первым приближением. Подставим хи г/,, zt в правые части, получим второе приближение: х2, уг, г2. Этот процесс можно осущест влять неограниченно. При известных условиях, налагаемых на коэффициенты, &-е приближение x k, y k, zk при возрас тании k будет становиться сколь угодно близким к ис тинному решению системы; это истинное решение обозначим (х*, у*, 2*). Мы говорим в этом случае, что итерационный процесс сходится. Значит, проведя процесс подстановки некоторое, достаточно большое число раз, мы получим приближенное решение системы. Для сходимости итера ционного процесса достаточно потребовать выполнения одного из следующих условий:
а) сумма модулей коэффициентов при неизвестных в пра вой части каждого уравнения системы (5.10) меньше еди ницы;
б) сумма квадратов всех коэффициентов при неизвест ных в правых частях системы (5.10) меньше единицы.
Поясним сказанное на примере. Пусть дана система нормального вида-:
{ х = 0,12х—0,18у + 0,08г— 0,64, | y = 0,15x + 0,06y—0,11г + 0,26, \ 2= 0,04х—0 ,Юг/—0,092+1,34.
Найдем сумму модулей коэффициентов при неизвестных
вкаждом уравнении.
Впервом уравнении:
0,12+0,18+0,08=0,38.
Во втором:
0,15+0,06+0,11=0,32.
В третьем:
0,04+0,10+0,09=0,23.
Для каждого уравнения эта сумма меньше единицы, значит, итерационный процесс сходится и мы можем получить приближенное решение системы методом итераций. За начальное приближение мо^кно всегда принять х0=.0; у 0~ 0; 20= 0. Подставив эти значения в правые части уравнений, получим первое приближение:
хх= —0,64; у !=+0,26; 2l= + l,3 4 .
132
Подставив первое приближение в правые части, получим второе приближение:
х-2= —0,6540; г/2==0,0289; г2=1,1951._
Продолжая процесс, замечаем, что четвертое и пятое при ближения совпадают в пределах трех десятичных знаков. Получаем приближенное решение:
х = —0,622; у = 0,033; z—1,231.
58°. О точности решения системы линейных уравнений. Метод Гаусса является точным. Это значит, что если коэф фициенты — точные числа и все действия в процессе реше ния выполнять точно, без округлений, то мы получим точные значения неизвестных. Но в вычислительной прак тике бывает обычно иначе. Пусть коэффициенты — прибли женные числа и в процессе решения, например, по схеме Гаусса все вычисления выполняем приближенно, т. е. ок ругляем результаты промежуточных действий. Полученное решение будет приближенным. Погрешность результата складывается из неустранимой погрешности исходных данных (коэффициенты — приближенные числа) и погреш ностей округления. Так как при решении системы прихо дится выполнять очень много арифметических действий, то погрешность, возникающая в результате округлений, может быть значительной. Учесть заранее эту погрешность очень трудно. На практике достаточно принять следующие меры
предосторожности. Промежуточные |
вычисления |
выпол |
|
нять |
с 1— 2 запасными цифрами, |
полученный |
результат |
(т. е. |
готовое решение) округлять на 1— 2 цифры. |
Если применяется приближенный метод, то погрешность результата складывается из неустранимой погрешности исходных данных, погрешности округления и погрешности метода (погрешность, возникающая в силу природы данного приближенного метода).
§ 16. РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ УРАВНЕНИЙ МЕТОДОМ ОПРЕДЕЛИТЕЛЕЙ
59°. Определители второго порядка. Пусть дана квад ратная таблица из четырех чисел:
Такая таблица |
называется матрицей второго порядка, |
а составляющие |
ее числа ах, Ьх, а2, Ь2— элементами этой |
матрицы.
Определителем ( или детерминантом) второго порядка, соответствующим матрице (5.11), называют число D,
равное а1Ь2—а2Ьх. Определитель обозначается символом
|
ах Ьг |
|
|
|
й2 |
Ь2 |
|
Таким образом, |
|
|
|
D = |
«1 Ьг |
— ахЬ2—а2Ьх. |
(5.12) |
|
а2 Ь2 |
|
|
60°. Определители третьего порядка. Рассмотрим квад ратную таблицу (матрицу) из девяти чисел (элементов матрицы):
ах Ъх сх
и% Ь2 |
с2 |
(5.13) |
а3 Ья |
cs |
|
Определителем третьего порядка, соответствующим мат рице (5.13), называется число D, равное следующей алгебраической сумме произведений элементов матрицы:
D 7 ахЬ2с3 bxc2a3■]■"сха2Ь3схЬ2а3 ахс2Ь3 Ьха2с3.
Обозначим определитель третьего порядка символом
ах сх
а2 Ь2 с2 h с3 •
Таким образом, по определению
ах |
bx |
сх |
|
|
а2 |
Ь2 |
с2 |
= аА сз+ Ьхсга3+ cxa2b3—cxb2a3— |
|
а3 |
Ь3 |
с3 |
— axc2b3— bxa2c3. |
(5.14) |
|
|
|
Числа ax, bx, cx, a2, b2, c2, a3, b3, c3 называются элемен тами определителя. В записи определителя мы разли чаем три строки и три столбца, а также две диагонали. Именно, элементы ах, Ь2, с3 образуют главную диагональ, элементы а3, Ь2, сх— побочную диагональ.
134
Of |
bi |
|
\ |
A |
X |
\ |
\/V? |
\ |
|||
|
|
|
N |
|
/ |
|
> У ' |
V |
X |
?c2 |
|
4 |
V ' |
°г<^\Pi \ |
|||
// \ ^ // \ \ |
|
/ \ \ / \ |
|
||
|
V |
\ |
/ |
V |
^ 4 |
*5 |
Aj |
|
От |
Aj |
Cs |
|
o; |
|
|
6 ! |
|
|
|
|
Рис. 38 |
|
|
Несмотря на то что выражение определителя третьего порядка через его элементы кажется довольно громоздким, в действительности закон его составления весьма прост. В самом деле, в выражение (5.14) со знаком плюс входят произведения элементов, принадлежащих главной диаго нали, а также произведения элементов, образующих парал лели к этой диагонали, и третьего множителя — элемента из противоположного угла матрицы (рис. 38, а). Со знаком минус входят произведение элементов, образующих по бочную диагональ, а также произведения элементов, образующих параллели к побочной диагонали, и третьего множителя — элемента из противоположного угла мат рицы (рис. 38,6). Это правило, называемое правилом тре угольников, позволяет очень просто вычислять определители третьего порядка.
П р и м е р .
2 1 - 1
5 -- 2 3 = 2 .( —2).1 + 1-3-(—1)+5-4-(—!) — (—Г) X
—1 4 1
Х(—2)-(—1)—4-3-2—5-1-1 = — 54.
61°. Решение систем линейных уравнений при помо щи определителей. Пусть дана система двух линейных уравнений с двумя неизвестными:
( atx + Ьгу —сь
(5.15)
I azx + b2y = c2.
135
Составим определители из коэффициентов при неизвест ных и правых частей:
ах |
Ьх |
д, = |
Ci ftj |
fli Cj |
а2 |
Ь2 ) |
^1 |
ь2 |
й2 с2 |
Определитель Д называется главным определителем системы, определители Ль Д2— вспомогательными определителями.
Как известно из основного курса школы, справедливы сле
дующие свойства системы (5.15). |
определитель |
системы |
|||
Т е о р е м а |
1. |
Если главный |
|||
отличен от |
нуля, |
то система имеет решение, |
притом |
||
единственное. |
Решение системы определяется формулами: |
||||
|
|
|
х —'Ь. |
|
(5.16) |
|
|
|
А ’ |
|
|
(эти формулы называются формулами Крамера). |
системы |
||||
Т е о р е м а |
2. |
Если главный |
определитель |
равен нулю, а хотя бы один из вспомогательных опреде лителей отличен от нуля, то система не имеет решений
(система несовместна).
Рассмотрим теперь систему трех уравнений с |
тремя |
неизвестными: |
|
( a1x + b1y-\-c1z = du |
|
| а2х + b2y с2г = d2, |
(5.17) |
{ a3x + b3y + c3z = d3. |
|
Составим определитель третьего порядка из коэффициентов при неизвестных:
ах bx сх
Д= а.2 Ь2 с2 @з Ь3 с3
Этот определитель называется главным определителем системы. Составим вспомогательные определители Дц Д2, Д3, которые получаются из главного заменой того или иного столбца правыми частями уравнений:
dx |
bt |
Cj |
, Д2 |
ax |
dx Cj |
ax |
bx |
dx |
|
d2 |
b2 |
c2 |
a2 |
d2 |
c2 |
> Дз — a2 |
b2 |
d2 |
|
d b 3 |
C3 |
|
«3 d3 |
c3 |
a3 |
b3 |
d3 |
Для системы трех уравнений с тремя неизвестными справедливы теоремы, аналогичные только что сформу-
136
лированным. В случае А Ф 0 единственное решение си стемы определяется формулами (формулы Крамера):
|
х — Д ’ |
(5.19) |
Метод |
определителей можно применять к |
системе |
с любым |
числом уравнений (и соответственным |
числом |
неизвестных). Пусть дана система п линейных уравне ний с п неизвестными хъ хг, хп.
Вводится понятие определителя п-го порядка, и для системы устанавливается главный определитель А, со ставленный из коэффициентов при неизвестных, а также вспомогательные определители Alf А2, . . . , А„. Опреде литель любого порядка есть многочлен, составленный из элементов определителя по определенному закону. Для системы п уравнений оказываются справедливыми теоремы, аналогичные приведенным выше. В случае А ф О единственное решение системы определяется фор мулами (формулы Крамера)-.
Метод определителей является очень удобным в теорети ческом отношении. Применяя определители, оказывается возможным четко выделить случаи, когда система уравне ний совместна, когда несовместна, когда имеет бесконечное множество решений. В случае единственного решения (Д=£0) значения неизвестных определяются по готовым формулам (формулам Крамера).
Казалось бы, формулы Крамера дают возможность легко вычислить значения неизвестных и получить решение системы. В принципе это так, но на практике дело ослож няется тем, что для отыскания решения системы по форму лам Крамера требуется выполнить чрезмерно много вы числений. Установлено, что для отыскания решения си стемы п уравнений по формулам Крамера требуется про извести
N = (n2— \)-п \+ п |
. |
(5.21) |
одних только умножений и делений. Применяя эту формулу, можем подсчитать число арифметических действий при различных значениях п. Приведем таблицу, указывающую
137
число арифметических дейст вий при решении систем по формулам Крамера и ' спо собом Гаусса (табл. XIX).
Из этой таблицы видно, что в вычислительной прак тике применение формул Кра мера невыгодно, поэтому сле дует пользоваться другими методами, в том числе методом Гаусса.
Метод исключений по схе ме Гаусса удобен еще и тем, что он может быть применен для вычислений на ЭВМ.
Упражнения.
1. Решить систему по методу Гаусса (схема един ственного деления).
( 0,625*1— 1,27х, + 6,02*„ —6,96 = 0,
<1,43*!+0,65х2 —2,15*,+ 1,92 = 0,
{2,25*!— 1,76*2 —0,316*3 —2,37 = 0.
2.Решить систему по методу итераций:
( |
5,92*!— 1,24*2— 1,84*3 = 2,44, |
| |
2,72*i—9,71*2 + 2,43*3 = 2,40, |
{1,76*1 —3,12*2 + 9,38*3= 1,93.
3.Решить систему по методу итераций:
*i = 0,1*2—0,2*з+0,3,
+ = — 0,1*1 + 0,1*з + 0,1*4 —0,2,
*3 = 0,1*2+ 0,1*4 + 0,1,
*4= 0,1*2--0,1*3+ 0,2.
Л А Б О Р А Т О Р Н А Я Р А Б О Т А № 4
Решение системы линейных уравнений по методу Гаусса.,
З а д а н и е . Решить систему трех. уравнений с тремя неизвестными по схеме единственного деления.
138
П о р я д о к п р о в е д е н и я р а б о т ы . 1. Составить расчетный бланк по схеме, изложенной в п. 54° (табл. XVII).
2.Заполнить расчетный бланк, применяя механические правила, изложенные в п. 54°. Одновременно с заполне нием основных столбцов заполнить также столбец «Суммы».
3.Выполнить все вычисления, получить значения неизвестных хи х2, х3 (строка 10) и значения вспомога
тельных неизвестных лу, х2, х3 по правилу, |
изложенному |
|||
в п. |
55°. |
заключительный |
контроль, сравнив хх, |
|
4. |
Провести |
|||
х2, х 3 и xlt х2, |
х3 (как указано в |
п. 55°). |
|
|
5. Провести дополнительный контроль, подставив |
||||
найденные значения хг, хй, хв в |
данную |
систему урав |
||
нений. |
|
|
|