§13. Разложение выражения относительно переменной X на эле- ментарные дроби

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

Mathcad разложить относительно переменной x на элементарные дроби

_{3 4}

_{выражение}

^{1 − x − x − x} _.

x⁴ −1

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

Рис. 2.16. Разложение выражения относительно переменной x на элементарные дроби с использованием пакета MathCAD

§14. Предел функции

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

₂

_{Mathcad найти предел функции}

_lim ^{4 ⋅ x}

^{− 2} _.

_x→ ¹

2

2 ⋅ x − 2

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

Рис. 2.17. Предел функции с использованием пакета MathCAD

Глава 4. Операции с матрицами

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

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

_{⎛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

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

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

^{MathCAD транспонировать матрицу} ₍^{1 b c}

²₎ ^.

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

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

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

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

⎜

_{MathCAD найти обратную матрицу от матрицы} ^{⎛ a −1}

_{c −1}

b −1 _⎞

.

_{d −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нагр .}

⎝ ⎠

^{⎜ R}ист ^{+ R}нагр ^⎟

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

^P_{нагр (} ^R_{нагр )} ^по

^Rнагр

_{dP ( R ) 2}

_{⎛ ⎞2}

_{( 2.5)}

нагр нагр

_dRнагр

= −2 ⋅

E

_{( R R}

_{⋅ Rнагр + ⎜} ^E _{⎟ .}

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

ист ⁺

нагр ^{⎝ ⎠}

ние

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

_E ² _⎛

_{E ⎞}²

₍ 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 мо-

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

ном виде.