Материал: 6386

Глава 7. Системы линейных алгебраических уравнений

В приведенном выше примере функция lsolve возвращает символьные зависимости, показывающие, что решение может быть выражено через параметр t_. А функция rref возвращает матрицу, эквивалентную А, последняя строка которой состоит из нулей. Запишем, какой системе уравнений соответствует матрица rref(AB):

Листинг 42. Бесконечное множество решений СЛАУ

7.4. Практическая работа №6. Системы линейных алгебраических уравнений

Продолжительность – 4 часа.

Цель работы

Закрепление умений, навыков и компетенций работы с матрицами и векторами в MathCAD. Закрепление знаний, полученных в курсе «Высшая математика» в разделах «Линейная алгебра», посвященных поиску решений СЛАУ. Развитие творческого, осознанного подхода к решению систем уравнений. Получение навыков рассуждения и анализа полученных результатов.

Задание на практическую работу

1.В соответствии с индивидуальным заданием (Таблица 11, Приложение 6) определить существование и единственность решения заданной системы.

2.Создать программу в среде MathCAD и найти решение заданной СЛАУ четырьмя способами:

1.1.по методу Гаусса при помощи встроенной функции rref( );

1.2.по правилу Крамера;

1.3.непосредственно умножением на обратную матрицу;

1.4.при помощи встроенных функций MathCAD.

3.Сравнить решения, найденные разными способами.

4.Для индивидуальных заданий (Таблица 12, Приложение 6) найти определитель, ранг основной и расширенной матриц и их размерности и объяснить условия существования решения.

5.В случае, если решение существует, обосновать разделение вектора неизвестных на зависимые и свободные переменные.

6.Выписать общее решение вида (20) и построить не менее трех частных решений, при фиксированных свободных переменных.

7.Составить отчет, в котором отразить листинг программного кода с комментариями.

Контрольные вопросы

56

Профессиональные математические пакеты

1.Как определить, разрешима ли система алгебраических уравнений?

2.Если система имеет решения, то как определить, сколько их?

3.Как отделить зависимые переменные от свободных? Сколько их?

4.Что такое расширенная матрица системы?

5.Что такое ранг СЛАУ? Чем он отличается от ранга матрицы СЛАУ?

6.В чем смысл метода Гаусса решения СЛАУ?

7.Как решить СЛАУ по методу Крамера? Можно ли решить вырожденную систему?

8.Какие методы поиска решений СЛАУ вам известны?

9.Как решить СЛАУ, если определитель матрицы равен 0?

7.5.Электрические цепи постоянного тока. Закон Ома. Законы Кирхгофа

Умение решать системы линейных алгебраических уравнений будет востребовано в процессе поиска решения задач из теории электрических цепей, теоретических основ электротехники (ТОЭ) [30, 31], физики, теории автоматического управления (ТАУ) и методов анализа и расчета электронных схем (МАРЭС).

Рассмотрим, в качестве примера, простейшую задачу из теории цепей постоянного тока.

Рис. 22 – Задача на баланс мощностей в цепи постоянного тока

На рисунке изображена электрическая цепь, включающая идеальные источники ЭДС Е1, Е2 и резисторы R1-R6. В правой части рисунка приведены номиналы этих элементов. Ставится задача при помощи закона Ома и законов Кирхгофа определить токи, протекающие по всем элементам цепи и вычислить генерируемую и потребляемую схемой мощность. Вспомним эти законы:

Иными словами, при полном обходе контура потенциал, изменяясь, возвращается к исходному значению. Это правило вытекает из 3-го уравнения Максвелла для стационарного магнитного поля.

Частным случаем второго закона Кирхгофа для цепи, состоящей из одного контура, является закон Ома для этой цепи.

57

Рассмотрим применение законов (21)-(23) на примере расчета токов в цепи (Рис. 22).

При построении уравнений и проведении расчетов необходимо учитывать тот факт, что схема содержит идеальные источники напряжения Е1 и Е2, т.е. вводится допущение, что внутреннее сопротивление этих источников равно нулю.

Алгоритмом расчета цепи постоянного тока [30, 31]

1)Произвольно выбрать положительные направления токов и обозначить их на схеме.

2)Составить уравнения по первому закону Кирхгофа (21), причем на одно уравнение меньше числа узлов (т.к. для последнего узла уравнение будет линейно зависимым от предыдущих уравнений).

3)Выбрать независимые (главные) контуры и направление их обхода. Удобно для всех контуров выбрать одинаковое направление обхода.

4)Записать уравнения по второму закону Кирхгофа (22) для выбранных контуров.

5)Решить полученную систему уравнений – определяют искомые токи.

Рис. 23 – а) условные направления токов в ветвях схемы; б) направления обхода контуров

Решение:

1) В цепи имеется 4 узла, обозначенных на Рис. 23, а цифрами 1-4, значит, потребуется составить 3 уравнения по первому закону Кирхгофа. Стрелками обозначены произвольно выбранные направления протекания токов i1 i6 .

2) Составим 3 уравнения токов в узлах по первому закону Кирхгофа (21):

3)На Рис. 23, б предложены выбранные произвольно независимые контуры I – III и направление их обхода.

4)Для выбранных контуров запишем уравнения по второму закону Кирхгофа (22):

5) В полученных уравнениях перейдем от напряжений к токам, пользуясь законом Ома (23) и учитывая направления токов (если ток направлен против обхода контура, он отрицательный)

контур I: i1 R1 i3 R3 i4 R4 E1 (ток i4 направлен против обхода контура I);

58

Профессиональные математические пакеты

контур II: i2 R2 i3 R3 i5 R5 E2 (токи i2 , i3 и i5 – против обхода контура II); контур III: i4 R4 i5 R5 i6 R6 E1 (все токи по направлению обхода контура III).

6) Из первых трех и последних трех уравнений составим матрицы СЛАУ и разрешим её любым известным способом (см. Практическая работа №6. Системы линейных

образом направления протекания токов в цепях схемы, мы угадали направление токов i2, i4, i5 и i6, а токи i1 и i3 имеют отрицательное значение, следовательно в реальной схеме они направлены в противоположную сторону.

Листинг 43. Решение полученной системы уравнений

7)Проверить полученные правильность вычисления результатов можно подставив значения найденных токов в уравнения, построенные по 1 закону Кирхгофа.

8)Зная сопротивления элементов цепи [R] и токи [A], через них протекающие легко построить значения падений напряжений [B] на каждом из резисторов по закону Ома. После этого можно произвести проверку по уравнениям, построенным по 2 закону Кирхгофа.

Листинг 44. Расчет падения напряжений в цепи постоянного тока

9) Проверить полученные результаты можно опираясь на закон сохранения энергии. То есть, вся суммарная мощность, выработана источниками ЭДС нашей схемы, должна быть потреблена нагрузкой.

Из школьного курса физики известно, что мощность источника ЭДС может быть вычислена как произведение напряжения на ток: Pист E i [Вт]. Таким образом, все источники схемы вырабатывают суммарную мощность, равную

n

P ист Ek ik [Вт]. (24)

k 1

С другой стороны, мощность, потребляемая нагрузкой может быть подсчитана Pпотр R i2 [Вт], отсюда можно вычислить суммарную мощность, рассеянную на сопротивлениях R1-R6:

m

P потр Rk ik2 [Вт]. (25)

k 1

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

59

Глава 7. Системы линейных алгебраических уравнений

Листинг 45. Баланс мощностей в цепи постоянного тока

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

7.6. Лабораторная работа №7. Расчет цепи постоянного тока

Продолжительность – 4 часа.

Цель работы

Закрепление умений и навыков по поиску решений СЛАУ. Актуализация математических компетенций в предметной области направления 210100.62. Развитие творческого, осознанного подхода к решению систем уравнений. Получение навыков рассуждения и анализа полученных результатов.

Восстановление знаний школьного курса физики, разделы «Электричество и магнетизм», «Электростатика». Практическое осознание закона сохранения энергии. Вырабатывание стратегического, аналитического взгляда на прикладные задачи и практического подхода к теоретическим знаниям.

Задание на лабораторную работу

1.Рассчитать токи заданной в индивидуальном задании цепи.

2.Топологический вариант схемы выдается преподавателем в соответствии с индивидуальным заданием (Таблица 14, Приложение 7), номиналы элементов схемы – из таблицы (Таблица 13, Приложение 7), лишние параметры не учитывать.

3.Отразить в отчете следующую последовательность действий:

3.1.Выбрать положительные направления протекания токов, обозначить их на схеме. Пронумеровать узлы схемы и составить n-1 уравнение по первому закону Кирхгофа.

3.2.Выбрать контуры, задать направление обхода контуров и составить уравнения по второму закону Кирхгофа. Перейти от напряжений к токам по закону Ома.

3.3.Составить систему линейных алгебраических уравнений относительно токов рассматриваемой схемы и решить ее в MathCAD.

3.4.Определить токи в цепях схемы. Рассчитать падение напряжения на каждом из резисторов.

3.5.Рассчитать баланс мощностей рассматриваемой схемы.

4.Составить отчет, в котором отразить последовательность построения СЛАУ относительно токов рассматриваемой схемы, программу с комментариями и привести скриншоты с результатами работы программ.

5.Все физические единицы в отчете обязательно(!) должны указываться с единицами измерения.

Контрольные вопросы

1.Сформулируйте первый закон Кирхгофа.

2.Как связано напряжение, сопротивление и ток в резисторе?

3.Сформулируйте закон сохранения энергии в приложении к балансу мощностей.

4.Мощность. По каким двум формулам можно вычислить мощность?

5.Сформулируйте закон Ома.

6.Могут ли токи иметь отрицательное значение? Почему?

7.Назовите единицы измерения тока, напряжения, сопротивления и мощности.

8.Сформулируйте второй закон Кирхгофа.

60