ОПЕРАЦИИ С МАТРИЦАМИ

Вычисление определителя матрицы

Пусть необходимо с использованием встроенных операторов пакета

⎛11⋅a

5 1 ⎞

MathCADвычислить определитель матрицы ⎜

7 11⋅b

4 ⎟и полу-

⎜ ⎟

⎜ 2 6 11⋅c⎟

⎝ ⎠

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

при a=2, b=2, c=2.

Соответствующий фрагмент документа MathCAD представлен на рис. 2.18.

Рис. 2.18. Вычисление определителя матрицы с использованием па-

кета MathCAD

Транспонирование матрицы

Пусть необходимо с использованием встроенных операторов пакета

MathCADтранспонировать матрицу (1 b c

2).

Соответствующий фрагмент документа MathCAD представлен на рис. 2.19.

Рис. 2.19. Транспонирование матрицы с использованием пакета Math- CAD

Нахождение обратной матрицы

Пусть необходимо с использованием встроенных операторов пакета

c−1

b −1 ⎞

⎝ ⎠

Соответствующий фрагмент документа MathCAD представлен на рис. 2.20.

Рис. 2.20. Нахождение обратной матрицы с использованием пакета

MathCAD

ГЛАВА 5.

ЗАДАЧИ ОПТИМИЗАЦИИ

Рассмотрим решение задач оптимизации в MathCADна

примере

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

ки. Данное уравнение можно легко решить аналитически и численно с использованием программы MathCAD.

Пусть дана цепь постоянного тока, состоящая из источника (гене-

ратор постоянного тока независимого возбуждения) и приёмника (рис.

2.21). Источник ЭДС

Eист

имеет внутреннее сопротивление

Rист

и ра-

ботает на омическую нагрузку

Rнагр . В цепи протекает ток

Iнагр .

Rист

Eист

Rнагр

Iнагр

Рис. 2.21

Необходимо определить такую величину сопротивления нагрузки

Rнагр , когда мощность нагрузки

будет максимальной

Pнагр = Iнагр2 ⋅ Rнагр

Pнагр = max.

(2.1)

Выведем зависимость мощности нагрузки как функцию от сопро-

тивления нагрузки.

Согласно закону Ома ток нагрузки определяется как

Iнагр =

E ⋅

Rист + Rнагр

(2.2)

Напряжение на нагрузке определяется как

U = I

⋅ R =

E ⋅ R

(2.3)

нагр нагр нагр Rист + Rнагр нагр

Мощность нагрузки есть функция сопротивления нагрузки

⎛ ⎞2

(2.4)

Pнагр (Rнагр )= Iнагр ⋅Uнагр = ⎜

E

⎟ ⋅ Rнагр .

Найдём производную функции

Pнагр (Rнагр ) по

Rнагр

dP (R ) 2

⎛ ⎞2

(2.5)

нагр нагр

dRнагр

= −2 ⋅

E

(R R

⋅ Rнагр + ⎜ E ⎟ .

)3 ⎜ Rист + Rнагр ⎟

ист +

нагр ⎝ ⎠

ние

Приравняв выражение (2.5) к нулю получим нелинейное уравне-

E2 ⎛

E ⎞2

(2.6)

−2 ⋅

3 ⋅ Rнагр + ⎜ R R

⎟ = 0 ,

⎜ + ⎟

(Rист + Rнагр )

⎝ ист нагр ⎠

где

Rнагр

- переменная, а

Rист и

E - константы.

После несложных математических выкладок, которые здесь не

приводятся, можно доказать, что

Rнагр = Rист

является корнем уравне-

ния, а, следовательно, и искомым решением задачи.

Данную задачу можно решить как численно, так и аналитически в MathCAD. Аналитическое решение в MathCAD по этапам показано да- лее.

Объявление функции мощности

нагрузки Pот переменных E–

ЭДС генератора, r – сопротивление генератора, R – сопротивление на-

грузки согласно (2.4) показано на рис. 2.22.

Рис. 2.22

Вызов процедуры вычисления производной возможно с помощью вкладки «Исчисление» математической панели или сочетания клавиш

<Shift> + </> (рис. 2.23).

Рис. 2.23

Определение производной функции P(E,r,R) по R производится после нажатия значка → на символьной панели (рис. 2.24).

Рис. 2.24

Для решения нелинейного уравнения воспользуемся процедурой

solve, доступной на символьной панели (рис. 2.25). Слева от ключевого слова solve записывается выражение, а справа – переменная.

Рис. 2.25

В выражении для нелинейного уравнения необходимо применить

знак «логическое равно», взятое с логической панели помощью сочетание клавиш <Ctrl> + <=>.

(рис. 2.26) или с

Рис. 2.26

Окончательно аналитическое решение выглядит так, как показано на рис. 2.27, то есть при равенстве сопротивления нагрузки и сопротив- ления генератора на нагрузке выделяется наибольшая мощность.

Рис. 2.27

Рассмотрим численное решение данной задачи. Пусть ЭДС гене-

ратора

составляет 10 В, а сопротивление генератора – 400 Ом. После

объявления функции P(E,r,R) (рис. 2.22) можно построить её график (рис. 2.28), из которого можно определить, что при R=r=400 Ом функ- ция имеет максимум, а это есть графическое решение задачи.

Рис. 2.28

Для численного решения задачи встроенной процедурой Minerr новой переменной R2 присвоим начальное приближение 300 Ом (рис. 2.29). Начальное приближение взято из графика (рис. 2.28) и мо- жет быть как больше, так и меньше 400 – приближённого графического решения. Затем пишется ключевое слово Given, после которого следует

нелинейное уравнение с правой частью, отличной от нуля и выражение для искомой величины Ropt.

Рис. 2.29

В данной главе на простом примере показано, что MathCAD мо-

жет решать задачи оптимизации, как в аналитическом, так и в числен-

ном виде.

ГЛАВА 6.