Решение систем линейных алгебраических уравнений методом Гаусса

Работа выполняется с использованием палитры программирования системы автоматизации математических вычислений Mathcad

Задание на работу:

1. Разработать программу для решения системы линейных алгебраических уравнений методом Гаусса с выбором главного элемента по столбцам.

2. Решить систему уравнений, определенную вариантом задания.

3. Произвести проверку решения.

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

Варианты заданий .

(Система уравненийв матричной записи имеет видAx = b)

1. A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 2. A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

3. A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 4. A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

5. A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 6.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

7. A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 8.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

9.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 10.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

11.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 12.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

13.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b = Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 14.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

15.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b = Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 16.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

17.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b = Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 18.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

19.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b = Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru 20.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

21.A= Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru , b = Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

Требования к оформлению отчета

Отчет должен содержать:

- Название и цель работы

- Задание на работу

- Текст программы на Mathcadе

- Результаты работы программы

- Проверку решения

Вспомогательные материалы

Пример разработки программы для решения системы линейных агебраических уравнений в системе Mathcad

Этап 1. Объявление переменных

Определим переменные и зададим матрицу системы линейных уравнений:

Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru Создадим новый лист в MathCAD и для начала напишем оператор ORIGIN := 1 – это для того чтобы MathCAD начинал нумерацию в матрицах с 1, а не с нуля.

Теперь определим матрицу А.

Обратите внимание, что в MathCAD строчные и прописные буквы различны, следите за регистром переменных.

Наберите A:=нажмите меню Insert è Matrixзадайте размерность матрицы 3 на 3 и нажмите OK. Заполните матрицу как показано на рисунке. Тем же способом задайте матрицу-столбец B.

Объединим все в одну расширенную матрицу matr с помощью оператора augment.

Присвоим переменной nзначение 3, она будет определять количество корней системы.

Чтобы посмотреть значение любой переменной, достаточно написать ее имя и поставить знак равенства. Выведем значение переменной matr.

Еще одной важной особенностью MathCAD является последовательность операторов,например, если мы попробуем написатьвыражениеmatr:=augment(A,B) перед определением матриц A и B, то MathCAD выдаст ошибку.

Переменные определены, теперь перейдем к написанию функции проверки матрицы на линейно-зависимые строки.

Этап 2. Проверка на линейно-зависимые строки

На языке Pascal эта процедура будет выглядеть так:

Function Proverka(): boolean;

var d1: single;

da: word;

lin: boolean;

Begin

d1:=0;

da:=0;

lin:=False;

For i:=0 to n-2 do

Forj:=i+1 to n-1 do

Begin

For k:=0 to n-1 do

If matr[j,k]<>0 then d1:=d1+matr[i,k]/matr[j,k];

d1:=d1/n;

For k:=0 to n-1 do

If matr[j,k]<>0 then If d1=matr[i,k]/matr[j,k] then da:=da+1;

If da=n then lin:=True;

d1:=0;

da:=0;

end;

Iflin then Proverka:=True Else Proverka:=False;

Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru end;

Решение систем линейных алгебраических уравнений методом Гаусса - student2.ru

Наши рекомендации