Подякувати Студентському архіву довільною сумою
Admin
26.02.2023 12:38
Дякуємо, що користуєтесь нашим архівом!
Інформація про навчальний заклад
ВУЗ:
Національний університет Львівська політехніка
Інститут:
О
Факультет:
РТ
Кафедра:
Не вказано
Інформація про роботу
Рік:
2015
Тип роботи:
Звіт до лабораторної роботи
Предмет:
Інші
Група:
УІ 21
Частина тексту файла (без зображень, графіків і формул):
МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА»
ЗВІТ
З лабораторної роботи №2
на тему: «МЕТОД ГАУССА ДЛЯ РОЗВ’ЯЗУВАННЯ СИСТЕМ
ЛІНІЙНИХ АЛГЕБРАЇЧНИХ РІВНЯНЬ
»
Львів-2015
Мета роботи – ознайомлення з прямими методами розв’язування систем лінійних алгебраїчних рівнянь.
Короткі теоретичні відомості:
Прямі методи розв’язування систем лінійних алгебраїчних рівнянь
Класичний метод Гаусса.
Розглянемо систему рівнянь четвертого порядку:
(1)
Зауважимо, що елементи вектора-стовпчика вільних членів занесені в матрицю коефіцієнтів А.
Будемо вважати, що . З першого рівняння знаходимо х1:
, (2)
де , .
З допомогою рівняння (2) можна виключити з решти рівнянь, для чого достатньо підставити (2) для в друге, третє і четверте рівняння системи. Це і є першим кроком – кроком виключення невідомого .
,
Перехід від початкової системи
до новоствореної
відбувається за такою формулою:
Другий крок – виключення невідомого відбувається аналогічно:
Третій крок – виключення невідомого
,
;
Останнє рівняння можна переписати у вигляді:
або .
Отже, в результаті прямого ходу одержимо систему рівнянь:
Розглянемо систему лінійних рівнянь
(1)
з невиродженою матрицею А розміру.
Більшу частину всього обчислювального процесу поглинає зведення матриці А до трикутного вигляду.
Можливі дві форми розкладу (зведення) матриці А до трикутного вигляду – рядкова або стовпцева. Як ми вже знаємо, стовпцева форма розкладу зображується наступною обчислювальною схемою:
Для до
Для до
(2)
Для до
Права частина системи (1) також може оброблятися в ході зведення матриці А до трикутного вигляду. Тому можна приєднати до і-го рядка (член ) – як ми й робили раніше (в такому випадку в циклі "k" верхня межа зростає до ). Можна залишити на місці, не вносячи в масив А. В цьому випадку в результаті виконання прямого ходу методу Гаусса одержується система рівнянь:
(3)
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..
.
Обернений хід при стовпчиковій формі розкладу описується загальною формулою:
, (4)
Розглянемо тепер рядкову форму розкладу матриці А. Вона базується на зведенні системи лінійних рівнянь до трикутного вигляду. Для цього спочатку нормують перше рівняння, ділячи його на а11(0), тобто роблять коефіцієнт при х1 рівним 1. Потім це перше рівняння домножують відповідно на коефіцієнт аі,1(0) при х1 всіх інших рівнянь і послідовно віднімають від усієї решти рівнянь. В результаті х1 буде виключене із всіх рівнянь, крім першого. На другому кроці виключають х2 з третього, четвертого, ..., п –го рівнянь. Цю процедуру повторюють до тих пір, доки вся система не буде зведена до такого трикутного вигляду
(5)
Рядкова форма зображається наступною обчислювальною схемою:
Для до
Для до
(6)
Для до
Тобто, на відміну від стовпчикової форми, обчислення коефіцієнтів нової матриці відбувається по рядках. Результат же одержується той самий.
При (обертанні) обчисленні оберненої матриці доцільно використовувати розклад матриці А до трикутного вигляду за рядковою формою.
Завдання
Розв’язати систему лінійних алгебраїчних рівнянь методом Гаусса.
Список ідентифікаторів констант, змінних, процедур і функцій.
Main – головна функція
for – оператор циклу;
ar[] – матриця коефіцієнтів
x[] –матриця розв’язків
if, else – умовний оператор
БЛОК_СХЕМА
початок
1
f(x)
x1:=f(x)
ні
k<3 k:=k+1
так
i<k+1
i:=i+1
ні
так
j>=k ні
j:=j-1
так
a[i,j]:= f(x)
f(x)
i>=0
i:=i-1
ні
j>i
j:=j-1
x[i]:= f(x)
x[i]:=f(x)
так
i<k+1
i:=i+1
так
x1,x2,x3,x4
КОД
using System;
class program
{
public static void Main()
{
double[] x = new double[4];
int i, j, k;
double[,] ar = { {8.3, 3.02, 4.1, 1.9, -9.55},
{3.92, 8.45, 7.38, 2.46, -12.21},
{3.77, 7.61, 8.04, 2.28, -14.45},
{2.21, 3.25, 1.69, 6.99, -8.35}};
Console.WriteLine("\n");
for (i = 0; i < 4; i++)
{
for (j = 0; j < 5; j++)
{
if (j == 4)
Console.Write(" = " + ar[i, j]);
else
{
Console.Write("{0} x{1}", ar[i, j], j + 1);
Console.Write(" + ");
}
}
Console.WriteLine("\n");
}
for (k = 0; k < 3; k++)
{
for (i = k + 1; i < 4; i++)
{
for (j = 4; j >= k; j--)
{
ar[i, j] = ar[i, j] - (ar[k, j] / ar[k, k]) * ar[i, k];
}
}
}
Console.WriteLine("\n\n\n result:\n");
for (i = 3; i >= 0; i--)
{
for (j = 3; j > i; j--)
x[i] = x[i] + ar[i, j] * x[j];
x[i] = (1.0 / ar[i, i]) * (ar[i, 4] - x[i]);
}
for (i = 0; i < 4; i++)
Console.WriteLine("x{0}={1}", i + 1, x[i]);
Console.ReadKey();
}
}
Результат
Висновок
При виконанні лабораторної роботи я ознайомилася з прямим методом розв’язування систем лінійних алгебраїчних рівнянь. Написала програму для розв’язання системи лінійних алгебраїчних рівнянь класичним методом Гаусса.
26.12.2015 19:12-
Ділись своїми роботами та отримуй миттєві бонуси!
0%
Оголошення від адміністратора