Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Boyarshinov_ChM_T1

.pdf
Скачиваний:
90
Добавлен:
13.03.2016
Размер:
1.32 Mб
Скачать

погрешности. Прецизионные измерения обеспечивают погрешность в пределах 10-8 - 10-10 %.

Как и в предыдущем случае, погрешности измерения вносят неустранимые искажения в результаты решения задач. В связи с этим имеет смысл накладывать ограничения и на точность выполнения математических вычислений: погрешность вычислений должна быть меньше погрешности измерения примерно в 10 раз. Более высокая точность вычислений является бессмысленной из-за наличия погрешностей представления исходных данных.

Погрешность численного метода

Погрешность метода решения задачи на вычислительной машине определяется неточностью замены алгебраического, дифференциального или интегрального оператора в исходном уравнении поставленной задачи приближенным (линейным или разностным).

Например, при замене интеграла конечной суммой:

b

N

 

h b a

 

f(x)dx f(a k h) h,

,

a

k 1

 

N

 

где N - число разбиений отрезка [a, b],

h - шаг интегрирования, точное значение площади

под графиком заменяется суммой площадей соответствующих прямоугольников (рис 1.2).

Очевидно, что при приближенном определении значения интеграла появляется погрешность, определяемая величиной “лишних” частей взятых прямоугольников. Чем меньше шаг интегрирования h, тем выше точность вычисления значения интеграла.

a

b

Рис. 1.2. Схема численного интегрирования

Аналогичная ситуация имеет место при замене производной разностным аналогом

(рис. 1.3):

df f(x h) f(x) . dx h

Понятно, что в этом случае погрешность замены производной разностным аналогом также уменьшается с уменьшением величины шага h. Иначе говоря, при использовании

11

численных методов погрешность решения является регулируемой: при корректном построении разностной аппроксимации исходного уравнения всегда имеется некоторый параметр, варьированием которого можно регулировать величину погрешности получаемого результата.

x x+h

Рис. 1.3. Схема численного дифференцирования

Вместе с тем следует иметь в виду, что повышение точности решения модели приводит к ощутимому повышению затрат ресурсов ЭВМ (времени проведения вычислений и оперативной памяти). Поэтому необходимо придерживаться “золотой середины”:

достижение приемлемых затрат ресурсов при получении удовлетворительной точности.

Погрешность проведения расчетов на вычислительных машинах

Оценка погрешности результата вычисления при заданной ошибке в представлении данных может быть проведена следующим образом. Пусть x - точное значение, а ~x -

приближенное значение той же переменной. Обозначим абсолютную погрешность d в

определении переменной разностью

~

x x.

Поскольку точное значение x неизвестно, введем “верхнюю” оценку D для величины погрешности: D . Определим также величину относительной погрешности

 

 

 

~x x

1

x

~

~

~ .

 

x

 

x

 

x

Абсолютная погрешность делится на приближенное значение переменной, поскольку ее точное значение неизвестно.

Погрешности округления чисел в ЭВМ

12

Округлением будем называть операцию замены заданного числа другим числом,

первые S значащих цифр которого совпадают с соответствующими цифрами исходного

числа, а начиная с S+1 разряда содержат нули.

Многие практические способы округления чисел выполняют отбрасывание “лишних” разрядов, хотя возможны варианты, при которых в “младший” разряд округленного числа, в

зависимости от ситуации, может добавляться единица.

Пусть, например, x=123456789, тогда при S=7 округленное число принимает значение

~ ~

x 123456700. В этом случае погрешность округления равна x x 89 100 D , то есть

не превышает единицы (с соответствующим порядком) в младшем разряде округленного числа.

Всякое вещественное число в компьютере представляется в нормализованном виде x =

a pb, где p - основание, b - показатель степени (целые числа), a - мантисса (вещественное число). Для определенности и однозначности будем считать p=10, 0,1 a 1,0. Ошибки

округления появляются при хранении именно мантиссы вещественного числа. В

представлении чисел на персональных компьютерах IBM достоверными могут быть 7

значащих цифр (для хранения числа отводится 4 байта оперативной памяти), 15 цифр (8

байт) или 19 (10 байт).

Рассмотрим оценки погрешности округления при S=7. Округленное число представляется в виде 0,XXXXXXX 10b , где под символом X может пониматься любая

цифра от 0 до 9. Очевидно, что абсолютная погрешность определяется значением

 

 

~x x D 1 10 7 10b 10b S,

S 7.

Модуль относительной погрешности

 

 

 

 

 

 

 

~

 

 

b S

 

 

b S

 

 

 

x x

 

10

 

 

10

 

b 101 S .

 

 

~

 

 

b

 

 

 

x

0,XXXXXXX 10

 

01, 10

 

Для некоторых

частных

случаев

погрешность представления вещественных чисел

оценивается:

 

 

 

 

 

 

 

 

 

S 7,

10 6

00001,%,

 

 

 

 

 

 

S 15,

10 14

0000000000001,%,

 

 

 

 

 

S 19,

10 18

00000000000000001,%.

 

 

 

 

Для всех чисел, представимых в электронно-вычислительной машине, относительная погрешность одна и та же, что очень существенно при получении оценок погрешностей математических моделей.

Погрешность результатов вычисления арифметических операций

13

Оценим погрешность результата сложения двух чисел, заданных с ошибкой:

~

~

x1

~

~

~

1

~

2 1 2,

D1 D2 .

x1

x2

x2 x1

x2

x1

x2

Аналогично определяется погрешность результата вычитания:

 

~

~

x1

~

~

~

1

~

2 1 2,

D1 D2 .

x1

x2

x2 x1

x2

x1

x2

Для определения абсолютных погрешностей операций умножения и деления двух

чисел проведем соответствующие выкладки:

~ ~

x1

x

~ ~

~

 

~

 

~ ~ ~ ~

 

~

 

~

1 2

 

x1 x2

2 x1 x2

x1

1 x2

2 x1 x2 x1 x2

1x2

2x1

~~

1x2 2 x1 1 2;

 

~

x1

 

~

 

~

1

 

~

~

 

~

~

1)

 

~

~

2

 

 

x1

 

x1

 

x1

 

x1

(x2

2 ) x2

(x1

 

x2 1 x1

.

~

 

~

~

2

 

 

~ ~

2 )

 

~ ~

 

2 )

 

x2

x2

 

x2

 

x2

 

 

 

x2 (x2

 

 

x2 (x2

 

Оценим относительные погрешности результатов умножения и деления:

~

~

 

 

 

 

 

 

 

 

~

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

x2 x1 x2

 

x

2

x

2

 

1

 

2

 

 

1

2

 

 

 

 

 

 

~ ~

 

 

 

 

 

1

 

2

 

1

 

 

1

~

 

~

 

 

 

 

1 2 1 2 1 2 ,

 

 

 

 

 

 

 

~

~

 

 

 

 

 

 

~ ~

 

 

x1

x2

 

 

 

 

 

x1

 

x2

 

 

 

 

x1

 

x2

 

 

x1

x2

 

 

 

 

 

 

 

~

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

1

 

 

~

 

2

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

~

 

 

 

~

 

 

 

 

 

~

 

 

x2

 

x1

 

 

 

 

 

 

 

 

 

~

 

 

x2

 

1

 

2

 

 

 

~

 

~

~

 

~

 

 

 

1 2

 

 

 

 

 

 

x2

~

 

 

x2

x1

 

x2

 

 

x1

x2

 

 

x1

x2

 

 

 

2 .

 

 

 

 

 

~ ~

 

 

 

 

 

 

~

 

 

 

 

2

 

 

 

 

1

 

1

 

 

 

 

 

x1

 

x (x

2

) x

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

В последних выражениях учитывается, что величины 1, 2 1.

Таким образом, при выполнении арифметических операций сложения и вычитания складываются (вычитаются) абсолютные погрешности, а при умножении и делении -

относительные погрешности.

“Потеря порядка” и “переполнение” при проведении вычислений на ЭВМ

Пусть w - наименьшее положительное число, представимое в ЭВМ; w 11754943, 10 38

при использовании 4 байт, и w 22250738585072014, 10 308 при 8-байтовом представлении

вещественной переменной. Это означает, что

все действительные числа

из интервала

w x w нельзя представить в ЭВМ в

нормализованном виде. В

этом случае

приближение истинного числа в компьютере равно ~x 0. Но это означает, что имеет место катастрофическая потеря точности арифметических вычислений:

1

x

,

x 0.

~

 

x

 

 

С другой стороны, пусть W - наибольшее положительное число, представимое в ЭВМ.

Все действительные числа |x| > W нельзя представить в нормализованной форме.

14

Приближение числа x определяется как

~

 

x W. Очевидно, что и в этом случае имеет место

x

,

~

 

x .

потеря точности: 1 ~

x W const,

x

 

 

 

 

Машинная реализация вычислений

 

 

 

 

 

 

1000

Рассмотрим задачу

об оценке погрешности вычисления произведения ai . В

 

 

 

 

i 1

соответствии с определением относительной погрешности приближенное значение

переменной можно представить в виде ~x x . Приближенный (округленный) результат

1

произведения двух первых чисел в указанном выражении может быть оценен формулой

a1 a2 . Очевидно,

что

для

получения

значения всего

выражения

следует многократно

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

повторить эту операцию:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

a

 

 

a

 

 

a

 

a

 

 

1000

 

999

 

 

1

 

2

 

 

3

 

 

4

1000

ai 1

.

 

 

 

 

 

 

1 1

1

i 1

 

 

 

1

 

 

Для 10 6

получаем

 

1 999

1000999499, ,

что соответствует величине

относительной погрешности вычисления всего результата 0,0998501 % . При перемножении

1 000 000 сомножителей погрешность может достигнуть 63,212 % . На практике такого катастрофического нарастания погрешности, как правило, не наблюдается в силу того, что погрешности округления могут иметь разные знаки и зачастую компенсируют друг друга.

Контрольные вопросы и задания

Дайте определение понятиям “математическая модель”, “вычислительный эксперимент”, “алгоритм”, “компьютерное моделирование”.

Опишите этапы выполнения вычислительного эксперимента.

Укажите причины погрешностей математической модели. Какие из них являются неустранимыми, а какие - регулируемыми?

Объясните причины погрешности исходных данных.

Что понимается под погрешностью численного метода?

Как оценить величину погрешности вычислений на ЭВМ?

Какую погрешность, относительную или абсолютную, целесообразно оценивать и почему?

15

2 . С И С Т Е М Ы Л И Н Е Й Н Ы Х А Л Г Е Б Р А И Ч Е С К И Х У Р А В Н Е Н И Й

Система m линейных алгебраических уравнений представляется в виде

 

 

Ax = f,

 

 

(2.1)

где A - квадратная матрица ранга m;

 

 

 

f1

 

 

x1

 

 

 

 

 

 

 

 

f2

 

 

x2

 

 

 

 

- правая часть системы уравнений;

 

 

- искомый вектор.

f f3

 

x x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fm

 

xm

 

Система уравнений (2.1) имеет единственное решение, если определитель det(A)

отличен от нуля. В развернутой (компонентной) записи эта система уравнений имеет вид

a11 x1 a12 x2 a13 x3 a1m xm f1,

 

 

 

x1 a2 2 x2 a2 3 x3 a2 m xm

f2,

 

 

a21

 

 

 

x1

a32

x2

a33 x3

a3 m xm

f3,

 

(2.2)

a31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

a

 

 

x

 

a

 

x

 

a

 

x

 

f

 

.

a

1

m 2

2

m 3

3

m m

m

m

m1

 

 

 

 

 

 

 

 

 

Прямые методы решения

Прямыми называют методы решения систем линейных алгебраических уравнений, для которых результат получается за конечное, заранее определенное, число арифметических операций. Наиболее популярным среди них является метод Гаусса5.

Метод Гаусса

Рассмотрим процедуру решения системы линейных алгебраических уравнений методом Гаусса на следующем примере:

5 Гаусс Карл Фридрих [30.4.1777 - 23.2.1855]. С 1795

по 1798 годы учился в Геттингенском университете. В

1799 году получил доцентуру в Брауншвейге, а

с 1807 года - кафедру математики и астрономии

Геттингенского университета, а также должность директора Геттингенской астрономической обсерватории.

С 1802 года являлся иностранным (с 1824 - иностранным почетным) членом Петербургской академии наук.

17

2x1 4x2 6x3 8,1x1 4x2 1x3 12,2x1 6x2 0x3 14.

Главный определитель такой системы

 

 

 

2

4

6

 

det A 1

4

1

8,

2

6

0

 

что гарантирует единственность решения.

1 шаг. Первая строка системы уравнений делится на первый коэффициент:

1x1 2x2 3x3 4,1x1 4x2 1x3 12,2x1 6x2 0x3 14.

2 шаг. Первая строка вычитается из второго уравнения:

1x1 2x2 3x3 4,0x1 2x2 4x3 16,2x1 6x2 0x3 14.

3 шаг. Из третьего уравнения вычитается первая строка, умноженная на 2:

1x1 2x2 3x3 4,0x1 2x2 4x3 16,0x1 2x2 6x3 22.

4 шаг. Второе уравнение делится на 2:

1x1 2x2 3x3 4,0x1 1x2 2x3 8,0x1 2x2 6x3 22.

5 шаг. Из третьего уравнения вычитается второе уравнение, умноженное на 2:

1x1 2x2 3x3 4,0x1 1x2 2x3 8,0x1 0x2 2x3 6.

6 шаг. Определяются искомые величины:

 

 

1x1 2x2 3x3 4,

1x1 2x2 3x3 4,

x1 4 2x2 3x3 1,

 

0x1 1x2 2x3 8,

 

x2 8 2x3 2,

 

x2 8 2x3 2,

 

 

 

 

x3 3.

 

x3 3.

 

x3 3.

 

 

 

18

Таким образом, получено решение исходной системы уравнений.

Теперь рассмотрим процедуру получения решения методом Гаусса в более общем случае. Пусть a11 0. Тогда первое уравнение системы (2.2) можно поделить на этот коэффициент:

 

 

 

1 x c

 

x

 

 

c

 

x

 

c

 

 

x

 

 

 

 

y f1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

12

 

 

 

 

2

 

 

 

 

 

13

 

 

 

3

 

 

 

1m

 

 

 

 

m

 

 

 

1

 

 

a11

 

 

 

 

 

 

 

 

 

С помощью этого уравнения можно преобразовать систему уравнений (2.2) к виду

1 x1 c12 x2 c13 x3 c1m xm y1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

(1)

 

x2 a

(1)

x3

a

(1)

xm

 

 

f

(1)

f2 a21 y1,

 

 

 

0 x1

2 2

 

2 3

2 m

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

(1)

 

x

 

 

a

(1)

x

 

 

a

(1)

x

 

 

f

 

(1)

f

 

a

 

 

y

,

 

 

 

 

0 x

1

 

2

3

m

 

 

3

31

 

 

 

 

 

 

 

32

 

 

 

 

 

 

 

33

 

 

 

 

 

 

 

 

 

 

3m

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

(1)

 

x

 

 

a(1)

 

x

 

a(1)

 

x

 

 

 

f

(1) f

 

 

a

 

 

y

 

.

0 x

1

 

2

 

 

3

 

m

m

m 1

1

 

 

 

m 2

 

 

 

 

 

 

 

m 3

 

 

 

 

 

 

 

 

 

 

m m

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

Здесь обозначено ai(1j) ai j

ai1

c1 j,

 

i, j 2,m . В полученной системе можно выделить

подсистему (m-1) линейных уравнений с (m-1) неизвестными величинами:

 

 

 

 

 

 

 

a(21)2

x2

 

a(21)3

x3

 

a(21)m

 

xm

 

f2(1),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

(1)

 

x

2

 

a

(1)

x

3

 

a

(1)

 

x

m

 

 

f(1),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32

 

 

 

 

 

 

 

 

33

 

 

 

 

 

 

3m

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

a(1)

x

 

 

a

(1)

 

 

 

x

 

 

 

f(1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a(1)

 

2

3

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 2

 

 

 

 

 

 

 

m 3

 

 

 

 

 

 

m m

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть теперь a22(1)

0. Поделим второе уравнение системы на этот коэффициент:

0 x1 1 x2 c23 x3 c2m xm y2

 

 

 

f

 

(1)

 

f2

a21

y1

.

 

 

 

 

 

 

2(1)

 

 

 

a

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a22

 

 

 

 

 

22

 

 

 

 

 

 

 

 

С помощью этого соотношения уравнения системы преобразуются к виду

1 x1 c12 x2 c13 x3 c1m xm y1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 x1 1 x2 c2 3 x3 c2 m xm y2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a(332)

x3

 

 

 

a(32m) xm f3(2)

f3(1)

 

a(321)

y2,

 

 

 

0 x1 0 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 x

 

 

a

(2)

 

x

 

 

 

a(2)

x

 

 

f

 

(2)

f

1

 

a

(1)

y

 

.

 

0 x

1

2

 

 

3

 

m

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

m 3

 

 

 

 

 

 

 

 

 

 

 

m m

 

 

 

 

 

 

 

m

 

 

 

m

 

 

 

 

m 2

 

 

 

 

 

 

Здесь обозначено a(i

2j)

a(i1j)

a(i12) c2 j,

 

 

i,j 3,m .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате преобразований получена подсистема (m-2) уравнений с (m-2)

неизвестными:

19

 

 

 

 

a(332)

x3

a(32m)

xm

f3(2),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

(2)

 

(2)

 

 

 

 

 

 

 

 

 

am 3

x3 am m

xm fm .

 

 

 

 

 

В предположении, что a33(2)

0, делим третье уравнение системы на этот коэффициент:

0 x 0 x

 

1 x

 

c

 

x

 

y

 

f

(2)

 

f3(1)

a32(1) y2

 

f3

a31

y1 a32(1) y2

.

 

 

 

 

 

 

3

 

a33(2)

 

 

a33(2)

1

2

 

3

 

3m

 

m

 

 

3

a33(2)

 

 

 

 

 

 

Снова выполняется следующий шаг по понижению порядка системы алгебраических уравнений, и так далее, до тех пор, пока вся система уравнений не будет преобразована к виду

1 x1 c12 x2 c13 x3 c1m xm y1,

 

1 x2 c2 3 x3 c2 m xm y2,

 

 

 

 

1 x3 с3 m xm y3,

 

 

 

 

 

 

1 xm ym.

 

В результате всех произведенных выкладок матрица коэффициентов А системы алгебраических уравнений приведена к виду

1 c12 c13

c1m

 

1

c23

 

 

0

c2m

U 0

0

1

c3m

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

0

 

- “верхняя” треугольная матрица6, укоторой равны нулю все элементы, расположенные под главной диагональю. Процедура получения такой матрицы носит название “прямого хода” метода Гаусса. Очевидным условием для успешного выполнения “прямого хода” является a(jjj 1) 0, j 1,m .

“Обратный ход” метода позволяет определить искомые величины:

6 Обозначение матрицы U принято по первой букве английского слова upper - “верхний”.

20

Тут вы можете оставить комментарий к выбранному абзацу или сообщить об ошибке.

Оставленные комментарии видны всем.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]