Курс лекций по основам алгоритмизации и программирования задач машиностроения. Кравченко Д.В. - 50 стр.

UptoLike

ВУЗ: 

УлГТУ | Ульяновск

Составители: 

Рубрика: 

48
Таким образом, искомое положение точки М (см. рис. 3.2) будет задано
координатами: х
м
= х
1
= 36,102 мм, а у
м
= х
2
= 63,141 мм.
Осуществим проверку правильности решения исходной системы (см. п.
3.2, пример 9):
+⋅−=
+⋅=
950,79х466,0у
50х364,0у
мм
мм
.
Подставим х
м
и у
м
в эту систему и проверим выполнение тождеств:
.
141,63141,63
141,63141,63
;
950,79102,36466,0141,63
50102,36364,0141,63
=
=
+⋅−=
+
⋅
=
Тождества обеспечиваются с погрешностью |E
1
| = 0 для первого уравнения
и |Е
2
| = 0 для второго уравнения в системе. В сравнении с методом Крамера (см.
п. 3.2, пример 9) метод Гаусса обеспечил более высокую точность решения сис-
темы.
Разработаем блок-схему алгоритма решения поставленной задачи методом
Гаусса (рис. 3.3) и программу SLU_GAUSS на языке TURBO PASCAL.
ПРОГРАММА SLU_GAUSS
(постановка задачи см. пример 9, алгоритм см. рис. 3.3)
Program SLU_GAUSS;
Uses crt;
Label 1;
Var i, j: Byte;
D, a11z, a12z, b1z, a22z, b2z, x2, x1, xm, ym, E1, E2: Real;
A,B: Array[1..5] Of Real;
Begin
Clrscr;
1: Writeln (‘Введите матрицу А совокупности коэффициентов а
ij
.’);
For i:=1 To 2 Do Begin
For i:=1 To 2 Do Begin
Read (’ ’, a[i,j],’ ’);{a
11
= -0,364, a
12
= 1, a
21
= 0,466, a
22
= 1}
End; End;
Writeln (’Введите вектор-столбец В совокупности коэффициентов b
i
.’);
For i:=1 To 2 Do Begin
Read (’ ’,b[i],’ ’); {b
1
= 50, b
2
= 79,950}
End;
D:=(a[1,1]*a[2,2])-(a[2,1]*a[1,2]);
If D=0 Then Begin

Страницы