Система алгебраических уравнений в маткаде

Решение нелинейных уравнений и систем уравнений в пакете MathCAD

Решение нелинейных уравнений

Вычисление корней численными методами включает два основных этапа:

· уточнение корней до заданной точности.

Рассмотрим эти два этапа подробно.

Отделение корней нелинейного уравнения

Учитывая легкость построения графиков функций в MathCAD , в дальнейшем будет использоваться графический метод отделения корней.

Пример. Дано алгебраическое уравнение

.

Определить интервалы локализации корней этого уравнения.

Пример. Дано алгебраическое уравнение

.

Определить интервалы локализации корней этого уравнения.

На рисунке приведен график функции , построенный в MathCAD . Видно, что в качестве интервала изоляции можно принять интервал . Однако уравнение имеет три корня. Следовательно, можно сделать вывод о наличии еще двух комплексных корней. ¨

Уточнение корней нелинейного уравнения

Для уточнения корня используются специальные вычислительные методы такие, как метод деления отрезка пополам, метод хорд, метод касательных (метод Ньютона) и многие другие.

Функция root . В MathCAD для уточнения корней любого нелинейного уравнения (не обязательно только алгебраического) введена функция root , которая может иметь два или четыре аргумента, т.е. или , где – имя функции или арифметическое выражение, соответствующее решаемому нелинейному уравнению, – скалярная переменная, относительно которой решается уравнение, – границы интервала локализации корня.

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

Заметим, что для вычисления всех трех корней использовался прием понижения порядка алгебраического уравнения, рассмотренный в п. 8.1.1. ¨

Функция root с двумя аргументами требует задания (до обращения к функции) переменной начального значения корня из интервала локализации.

Пример 8.1.5. Используя функцию root , вычислить изменения корня нелинейного уравнения при изменении коэффициента а от 1 до 10 с шагом 1.

Функция polyroots . Для вычисления всех корней алгебраического уравнения порядка (не выше 5) рекомендуется использовать функцию polyroots . Обращение к этой функции имеет вид polyroots (v) , где v – вектор, состоящий из n +1 проекций, равных коэффициентам алгебраического уравнения, т.е. . Эта функция не требует проведения процедуры локализации корней.

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

Блок Given . При уточнении корня нелинейного уравнения можно использовать специальный вычислительный блок Given , имеющий следующую структуру:

Решаемое уравнение задается в виде равенства, в котором используется «жирный» знак равно, вводимый с палитры Логичес­кий .

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

Функция Find уточняет корень уравнения, вызов этой функции имеет вид Find ( x ), где x – переменная, по которой уточняется корень. Если корня уравнения на заданном интервале не существует, то следует вызвать функцию Minerr ( x ), которая возвращает приближенное значение корня.

Для выбора алгоритма уточнения корня необходимо щелкнуть правой кнопкой мыши на имени функции Find ( x ) и в появившемся контекстном меню (см. рисунок) выбрать подходящий алгоритм.

Аналогично можно задать алгоритм решения и для функции Minerr ( x ).

Использование численных методов в функциях Find ( x ), Minerr ( x ) требует перед блоком Given задать начальные значения переменным, по которым осуществляется поиск корней уравнения.

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

Решение систем уравнений

В зависимости от того, какие функции входят в систему уравнений, можно выделить два класса систем:

· алгебраические системы уравнений;

· трансцендентные системы уравнений.

Среди алгебраических систем уравнений особое место занимают системы линейных алгебраических уравнений (СЛАУ).

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

Системой линейных алгебраических уравнений (СЛАУ) называется система вида:

В матричном виде систему можно записать как

,

где – матрица размерности , – вектор с проекциями.

Для вычисления решения СЛАУ следует использовать функцию lsolve , обращение к которой имеет вид: lsolve (А, b ), где А – матрица системы, – вектор правой части.

Решение систем нелинейных уравнений

MathCAD дает возможность находить решение системы уравнений численными методами, при этом максимальное число уравнений в MathCAD 2001 i доведено до 200.

Для решения системы уравнений необходимо выполнить следующие этапы.

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

Пример. Дана система уравнений:

Определить начальные приближения для решений этой системы.

Видно, что система имеет два решения: для первого решения в качестве начального приближения может быть принята точка (-2, 2), а для второго решения – точка (5, 20). ¨

Вычисление решения системы уравнений с заданной точностью . Для этого используется уже известный вычислительный блок Given .

Функция Find вычисляет решение системы уравнений с заданной точностью, и вызов этой функции имеет вид Find ( x ), где x – список переменных, по которым ищется решение. Начальные значения этим переменным задаются в блоке . Число аргументов функции должно быть равно числу неизвестных.

Следующие выражения недопустимы внутри блока решения:

· ограничения со знаком ¹ ;

· дискретная переменная или выражения, содержащие дискретную переменную в любой форме;

· блоки решения уравнений не могут быть вложены друг в друга, каждый блок может иметь только одно ключевое слово Given и имя функции Find (или Minerr ).

Пример. Используя блок Given , вычислить все решения системы предыдущего примера. Выполнить проверку найденных решений.

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

Система алгебраических уравнений в маткаде

РЕШЕНИЕ УРАВНЕНИЙ И СИСТЕМ УРАВНЕНИЙ

4 Решение уравнений и систем средствами Mathcad

Система Mathcad обладает широкими возможностями численного решения уравнений и систем уравнений.

Функция root, блоки Given…Find, Given…Minerr

В ходе численного решения обычно выделяют два этапа:

• отделение корней – определение интервала нахождения каждого корня или определение приблизительного значения корня. В системе Mathcad наиболее наглядным будет отделение корней уравнения графическим способом;
• уточнение корней – нахождение численного значения корня с указанной точностью.

Точность нахождения корня устанавливается с помощью системной переменной TOL (Convergence Tolerance – Допуск сходимости), которая по умолчанию равна 10 -3 . Чем меньше значение TOL, тем точнее, вообще говоря, находится корень уравнения. Однако оптимальным является TOL = 10 -5 . Переопределить значение TOL можно в окне математических свойств документа Math Options на вкладке Build-In Variables (Встроенные переменные) или присваиванием, например, TOL:=0.0001.

Для решения одного уравнения с одной неизвестной предназначена встроенная функция root, которая в общем виде задается

root(f(x), x, [a, b])

и возвращает значение переменной x, при котором функция f(x) обращается в ноль. Аргументы функции root:

• f(x) – функция левой части уравнения f(x) = 0;
• x – переменная, относительно которой требуется решить уравнение;
• a, b (необязательные) – действительные числа, такие что a -1 слева: A -1 Ax=A -1 b. Учитывая, что A -1 A, вектор-столбец решений системы можно искать в виде

Этот прием используется в Mathcad так:

1. задается матрица коэффициентов при неизвестных системы A;
2. задается столбец свободных членов b;
3. вводится формула для нахождения решения системы X:=A -1 b;
4. выводится вектор решений системы X=.

Кроме того, пакет Mathcad имеет встроенную функцию

lsolve(A, b),

возвращающую вектор-столбец решений системы линейных алгебраических уравнений. Аргументами функции lsolve являются матрица коэффициентов при неизвестных системы и столбец свободных членов. Порядок решения аналогичен рассмотренному, но вместо формулы X:=A -1 b используется X:=lsolve(A, b).

Реализовать широко известный метод Гаусса решения систем линейных уравнений позволяет встроенная функция rref(M), возвращающая ступенчатый вид матрицы M. Если в качестве аргумента взять расширенную матрицу системы, то в результате применения rref получится матрица, на диагонали которой – единицы, а последний столбец представляет собой столбец решений системы.

Решение системы линейных уравнений можно осуществить с помощью блоков Given…Find, Given…Minerr. При этом неизвестным системы задается произвольное начальное приближение, а проверка необязательна.

Порядок выполнения лабораторной работы

1. Загрузить Mathcad Start / All Programs / Mathsoft Apps / Mathcad (Пуск / Все программы / Mathsoft Apps / Mathcad).
2. Сохранить в личной папке на диске z:\ новый документ с именем ФИО1, лучше использовать латинские буквы. Производить сохранение регулярно в процессе работы (Ctrl + S).
3. Вставить текстовую область Insert / Text Region (Вставка / Область текста) и ввести в поле документа текст:

Лабораторная работа № 4
Решение уравнений и систем в Mathcad.

1. В новой текстовой области ввести фамилию, имя, отчество, учебный шифр и номер варианта.
2. Выполнить задание 1.

Задание 1. Решить уравнение .

Решение.

Решение данного уравнения будем проводить в два этапа: отделение корней уравнения графически, уточнение корней уравнения.

Определим функцию f(x), равную левой части данного уравнения, когда правая равна нулю:

Зададим ранжированную переменную x на некотором диапазоне с мелким шагом, например:

Вставим в документ графическую область. Для этого выберем дважды пиктограмму с изображением графика сначала на панели Math (Математика), затем на палитре графиков Graph или выполним из главного меню последовательность команд Insert / Graph / X-Y Plot (Вставка / График / X-Y Зависимость).

Снизу по оси абсцисс наберем x, а сбоку по оси ординат введем f(x).

Для появления графика щелкнем левой клавишей мыши вне графической области.

Отформатируем график функции f(x). Для этого щелкнем правой клавишей мыши в области графика и выберем в контекстном меню команду Format (Формат). Установим пересечение осей графика (Crossed – Только оси), добавим вспомогательные линии по координатным осям (Grid Lines – Вспомогательные линии). Отменим при этом автосетку (Autogrid – Автосетка) и установим количество линий сетки, равное 10.

Для подтверждения внесенных изменений нажмем последовательно кнопки Apply (Применить) и ОК.

После указанных преобразований график функции f(x) будет выглядеть следующим образом:

Из графика функции f(x) видно, что уравнение имеет три корня, которые приблизительно равны: x1 ≈ -1; x2 ≈ 1; x3 ≈ 2,5.

Этап отделения корней завершен.

Уточним теперь корни уравнения с помощью функции root.

Присвоим начальное приближение переменной x и укажем точность поиска корня:

Уточним заданное приближение к значению корня с помощью функции root:

Выполним проверку, подтверждающую, что первый корень найден с заявленной точностью:

Начальное приближение можно не задавать при использовании в качестве аргументов root границ отрезка нахождения корня, например, второй корень можно уточнить:

Задание 2. Решить уравнение .

Решение.

Напечатаем левую часть уравнения, не приравнивая выражение к 0, и выделим синим курсором переменную x:

Выберем из главного меню Symbolics / Polynomial Coefficients (Символика / Коэффициенты полинома). Появившийся вектор коэффициентов полинома выделим целиком синим курсором и вырежем в буфер обмена, используя кнопку Вырезать на панели инструментов Formatting (Форматирование) или комбинацию клавиш Ctrl + X.

Напечатаем v := и вставим вектор из буфера обмена, используя кнопку Вставить на панели инструментов или комбинацию клавиш Ctrl + V.

Для получения результата напечатаем polyroots(v) =:

Задание 3. Решить систему линейных уравнений Сделать проверку.

Решение.

1-й способ. Использование блока Given … Find.

Зададим всем неизвестным, входящим в систему уравнений, произвольные начальные приближения, например:

Напечатаем слово Given. Установим визир ниже и наберем уравнения системы, каждое в своем блоке. Используем при этом логический знак равенства (Ctrl + =).

После ввода уравнений системы напечатаем X := Find(x, y, z) и получим решение системы в виде вектора, состоящего из трех элементов:

Сделаем проверку, подставив полученные значения неизвестных в уравнения системы, например, следующим образом

После набора знака «=» в каждой строке должен быть получен результат, равный или приблизительно равный правой части системы. В данном примере системная переменная ORIGIN = 1.

2-й способ. Использование блока Given…Minerr.

Порядок решения системы этим способом аналогичен порядку использования блока Given … Find и представлен ниже вместе с проверкой:

3-й способ. Решение системы линейных уравнений матричным способом.

Создадим матрицу А, состоящую из коэффициентов при неизвестных системы. Для этого напечатаем A := , вызовем окно создания массивов (Ctrl + M). Число строк (Rows) и столбцов (Columns) матрицы данной системы равно 3. Заполним пустые места шаблона матрицы коэффициентами при неизвестных системы, как показано ниже:

Зададим вектор b свободных членов системы. Сначала напечатаем b :=, затем вставим шаблон матрицы(Ctrl + M), где количество строк (Rows) равно 3, а количество столбцов (Columns) равно 1. Заполним его:

Решим систему матричным способом по формуле

Решим систему с помощью функции lsolve:

Для проверки правильности решения системы, полученного матричным способом, достаточно вычислить произведение A·X, которое должно совпасть с вектором-столбцом свободных членов b:

Система алгебраических уравнений в маткаде

Электронный курс по MathCAD

Лекция 5.
Решение уравнений и систем.

5.1 Решение алгебраических (и других) уравнений и систем.

5.2 Решение дифференциальных уравнений и систем (задача Коши и граничные задачи).

5.3 Задание.

5.1 Решение алгебраических (и других) уравнений и систем.

Линейные алгебраические уравнения.

Определение: Уравнение вида ax+b=0 с заданным базовым множеством G_x, a из G_a , b из G_b называется линейным уравнением.

Этапы решения при помощи Mathcad:

1. Ввести уравнение (знак «=» вводится при помощи комбинации [Ctrl++]).
2. Выделить курсором переменную, относительно которой должно быть решено уравнение.
3. Выбрать команду Solve (Вычислить) подменю Variable (Переменные) меню Symbolics (Символы).

При решении линейных уравнений (без параметров) или дробных уравнений, которые сводятся к линейным, MathCAD находит все существующие решения. Однако при этом следует правильно интерпретировать сообщения, выдаваемые системой.

Нормальный случай.

В качестве решения MathCAD выдает число — это означает,

что уравнение однозначно разрешимо (однозначное решение линейного уравнения над множеством действительных чисел, которое одновременно является областью определения этого уравнения).

Рассмотрим другой пример:

После выполнения описанных выше действий для нахождения решения Mathcad выдает сообщение о том, что решение не найдено.

Проанализировав данное уравнение приходим к выводу, что выданное Mathcad сообщение означает, что решений нет L=<>.

MathCAD выдает сообщение «Решение не найдено», даже если уравнение имеет «формальное решение», которое не принадлежит области определения (смотри примеры ниже).

Многозначность.
Если в качестве решения MathCAD выдает имя переменной, это означает, что множество решений уравнения совпадает с областью определения. Однако, такие понятия, как множество решений уравнения и область определения, отсутствуют в MAthCAD и он не выписывает оболасть определения. Вы можете найти область определения, решая с помощью Mathcad систему неравенств или уравнений

Такой результат, выданный Mathcad после выполнения действий по решению уравнения, означает, что любое значение x из базового множества удовлетворяет этому уравнению, т. е. L=R.

Дробные уравнения

Команда Solve (Вычислить) из подменю Variable (Переменные) меню Symbolics (Символы)выдает множество решений: L = .

Решение 6 копируем в буфер, а затем выделяем маркером переменную x и активизируем команду Substitute (Замена) подменю Variable (Переменные) меню Symbolics (Символы) для замены переменной значением 6.

Рассмотрим другой пример:

Последнее уравнение (рисунок справа) условно эквивалентно уравнению:2x=4. Решение уравнения Mathcad: 2. Формальное решение x = 2 не входит в область допустимых значений. Mathcad выдает правильное сообщение!

Здесь также правильное решение: множество решений совпадает с областью допустимых значений L = D. Только следует учесть, что D=>.

Квадратные уравнения и алгебраические уравнения высших порядков.

Определение: Уравнение P(x)=0 называется алгебраическим уравнением n-го порядка, если P(x) представляет собой полином степени n, при n=2 данное уравнение называется квадратным уравнением.

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

Квадратное уравнение.

Команда Solve (Вычислить) подменю Variable (Переменные) меню Symbolics (Символы) дает решение в виде вектора: L= .

Иррациональное уравнения (уравнения с радикалами).

Корни (радикалы) могут вычисляться в MathCAD либо при помощи знака корня (клавиши [Ctrl+\]), либо как степени (клавиша [^] с дробными показателями. Знак квадратного корня вводится нажатием клавиши [\]. Знак корня и квадратный корень можно найти на панели Calculator (Калькулятор). Последовательность действий при решении уравнений с радикалами та же, что и при решении рассмотренных ранее уравнений.

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

• Радикалы определены не везде в действительной области. Это обстоятельство приводит к необходимости находить область определени, прежде чем решать само уравнение. Данная проблема справедливо игнорируется MathCAD, поскольку он не может знать, во множестве каких чисел (действительных или комплексных) вы намерены решать уравнение. Выход: вы можете самостоятельно найти область определения, воспользовавшись при этом возможностями MathCAD, связанными с решениями неравенств.
• Вторая проблема, возникающая при решении уравнений с радикалами, имеет принципиальный характер. Функция x 2 (как и любая другая функция с четным показателем) на является инъективной (проблема главных значений). В связи с этим возведение в квадрат обеих частей уравнения, содержащего квадратные корни, не является эквивалентным преобразованием. Как всегда, при применении к обеим частям уравнения не инъективного преобразования увеличивается множество решений. В результате в него могут войти «фиктивные» решения. Как ни удивительно, MAthCAD сам производит проверку решений на «фиктивность».

Классический случай решения уравнения с радикалами.

Mathcad распознает «фиктивные» решения (которые могут возникнуть в результате неэквивалентного преобразования «возведение в квадрат») и выдает верное сообщение: Решение не найдено. L =

В приведенных примерах демонстрируется способность MathCAD находить область определения иррационального уравнения путем решения неравенств.

Уравнения с радикалами третьей степени, как и уравнения с комплексными коэффициентами, не представляют для MathCAD никакой сложности.

Уравнения с параметрами.
При решении уравнений с параметрами MathCAD ведет себя по-разному, в зависимости от того, каким образом производятся символьные вычисления — с помощью символьного знака равенства или команд меню Symbolics.

В данном примере использование палитры символьных преобразований позволяет решить уравнение (solve) и упрстить результат (simplify)