Решение систем линейных уравнений maple (7 видео)

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

Решение обыкновенных уравнений.

Для решения уравнений в Maple существует универсальная команда solve(eq,x) , где eq – уравнение, x – переменная, относительно которой уравнение надо разрешить. В результате выполнения этой команды в строке вывода появится выражение, которое является решением данного уравнения. Например:

Если уравнение имеет несколько решений, которые вам понадобятся для дальнейших расчетов, то команде solve следует присвоить какое-нибудь имя name . Обращение к какому-либо k –ому решению данного уравнения производится указанием его имени с номером решения k в квадратных скобках: name[k] . Например:

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

Системы уравнений решаются с помощью такой же команды solve(,) , только теперь в параметрах команды следует указывать в первых фигурных скобках через запятую уравнения, а во вторых фигурных скобках перечисляются через запятую переменные, относительно которых требуется решить систему. Если вам будет необходимо для дальнейших вычислений использовать полученные решения уравнений, то команде solve следует присвоить какое-нибудь имя name . Затем выполняется присвоения команда assign(name) . После этого над решениями можно будет производить математические операции. Например:

Численное решение уравнений.

Для численного решения уравнений, в тех случаях, когда трансцендентные уравнения не имеют аналитических решений, используется специальная команда fsolve(eq,x) , параметры которой такие же, как и команды solve . Например:

Решение рекуррентных и функциональных уравнений.

Команда rsolve(eq,f) позволяет решить рекуррентное уравнение eq для целой функции f . Можно задать некоторое начальное условие для функции f(n) , тогда получиться частное решение данного рекуррентного уравнения. Например:

Универсальная команда solve позволяет решать функциональные уравнения, например:

F := proc ( x ) RootOf(_ Z ^2 — 3*_ Z + 2* x ) end

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

Решение тригонометрических уравнений.

Команда solve , примененная для решения тригонометрического уравнения, выдает только главные решения, то есть решения в интервале [0,2 p ]. Для того, чтобы получить все решения, следует предварительно ввести дополнительную команду _EnvAllSolutions:=true. Например:

В Maple символ _ Z

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

Решение трансцендентных уравнений.

При решении трансцендентных уравнений для получения решения в явном виде перед командой solve следует ввести дополнительную команду _EnvExplicit:=true . Пример решения сложной системы трансцендентных уравнений и упрощения вида решений:

Содержание

Решение уравнений
Системы уравнений
Maple. Решение алгебраических задач. Решения уравнений, систем уравнений и неравенств в Maple
Страницы работы
Содержание работы
💥 Видео

Видео:Графики, функции, решение системы линейных уравнений в MapleСкачать

Решение уравнений

Методы нахождения корней полиномов, решения уравнений, содержащих элементарные и специальные функции и систем сложных уравнений

Видео:Решение системы линейных уравнений в MapleСкачать

Системы уравнений

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

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

В первом примере используем solve для простой задачи линейной алгебры: John вдвое старше Kimberly. Возраст Kimberly плюс возраст John равен 27. Найти возраст каждого. Если использовать пакет LinearAlgebra, то придется рассматривать матрицу, но можно применить solve (и fsolve ), которые могут непосредственно работать с уравнениями:

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

solve и fsolve можно применять для решения нелинейных систем, т. е. систем уравнений, в которых переменные – квадраты, кубы, синусы, экспоненты и т. п. Например, вот система двух нелинейных уравнений:

Вначале попробуем применить команду solve :

Maple сделал по-умному: чтобы получить уравнение для х , он исключил у из Е2 с помощью Е1 . Затем он факторизовал это квадратное уравнение, выдал ответ (x,y)=(3,4) и свел оставшуюся часть задачи к кубической. Если завершить задачу командой evalf , получим:

Но если нарисовать кубическую часть в RootOf (для оценки положения корней), то увидите, что есть еще два решения. Где же они? Примените fsolve и получите:

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

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

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

В Maple есть другой полезный инструмент для случая двух нелинейных уравнений с двумя неизвестными. Точно так же, как в задачах с одной переменной, полезно сначала строить график, чтобы увидеть, где есть решения. В данном случае двумерные графики помогут искать корни двух неизвестных величин. Примените команду построения графиков implicitplot , которая берет уравнение с двумя переменными, вроде x 2 – y = 5, и строит определяющую их кривую в плоскости xy (для добавления команды графики надо вначале загрузить пакет графики with(plots) ):

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

Из картинки ясно, что две параболы пересекаются в четырех местах, поэтому должно быть четыре решения. Окно графика должно быть достаточным, чтобы увидеть всю картинку. А если оно мало, то получится вот что:

Если у вас есть три нелинейных уравнения с тремя переменными, implicitplot3d может сделать нечто подобное (см. Maple help).

Найдите все решения (Re и Im) системы уравнений

Чтобы определить количество искомых корней, сначала постройте график с помощью implicitplot .

Вот еще нелинейная система:

Вначале попробуем solve :

(Maple на мгновение задумается, но ничего не произойдет.) Теперь попробуем fsolve с диапазонами для каждой переменной:

Похоже, что (x, y, z) = (1, 1, 3) достаточно близко к решению. Предупреждение: в трех и более измерениях Maple может ошибиться и работать, несмотря на то, что:

(a) известно, что здесь есть решение и

(b) указано, где искать приближенное значение корня.

Для лучшего понимания, где следует искать решение, можно применить implicitplot3d :

Щелкните на рисунке и покрутите его, чтобы разглядеть подробнее. После этого перерисуйте график так, чтобы он был вблизи известного решения: [x, y, z]=[1, 1, 3]:

В середине графика все три поверхности – E1 , E2 и E3 – пересекаются в точке. Это и есть то, что искали с помощью implicitplot3d , но в целом рассматривать трехмерные задачи сложно.

Пусть гладкая функция y(x) представлена тремя точками (x 1 , y 1 ), (x 2 , y 2 ), (x 3 , y 3 ). Предположим, что эти точки определяют параболу вида y(x) = a + bx + cx 2 и что есть три аппроксимирующих уравнения для нахождения коэффициентов параболы a, b, c. Постройте эти три уравнения и используйте solve , чтобы найти формулы для a, b, c, а затем определите выражение Maple для параболы.

Задание выглядит вполне приемлемо для равноотстоящих точек. Пусть x 2 = x 1 + h и x 3 = x 1 + 2h . Упростим выражения для a, b, c. В результате получили приближенную форму функции, с которой можно работать:

(a) Оцените площадь под кривой между x 1 и x 3 путем интегрирования параболы между этими двумя пределами. Получится правило Симпсона. Чтобы посмотреть, хорош ли этот приближенный интеграл, задайте x 1 = 0, x 2 = 0.5, x 3 = 1.0 и y 1 = cos(x 1 ), y 2 = cos(x 2 ), y 3 = cos(x 3 ) , при этом приближенное значение площади будет близко к Дает ли формула приближенного интегрирования хороший результат?

(b) Оцените производную функции у(х) в средней точке х 2 путем дифференцирования параболы и вычисления значения производной в x = x 2 . Она называется формулой центральных разностей для первой производной. Проверьте ее точность при x 2 = 0.5, используя x 1 = 0.4 и x 3 = 0.6 с функцией

(c) Повторно дифференцируя формулу параболы, оцените вторую производную функции у(х) в x 2 . Это центральная вторая производная для равноотстоящих точек. Проверьте ее точность, как в части (b).

Эти формулы дифференцирования и интегрирования понадобятся в курсе физики.

Видео:Решение систем линейных уравнений в MapleСкачать

Maple. Решение алгебраических задач. Решения уравнений, систем уравнений и неравенств в Maple

Страницы работы

Содержание работы

Решение алгебраических задач

Решение уравнений (часть 1)

Для решения уравнений, систем уравнений и неравенств в Maple используется команда (оператор) «solve». Например, для решения уравнения х2 – 6х + 5 = 0 набираем: > solve(x^2–6*x+5=0); 1, 5 Обратите внимание, что если аргумент solve не является уравнением (или неравенством), то Maple трактует его так, как если бы это выражение было приравнено к 0. Можно было бы написать «> solve(x^2–6*x+5);». При решении алгебраических уравнений Maple приводит все корни, включая комплексные: > solve(x^4=1); 1, –1, I, –I

Решение уравнений (часть 2)

С помощью команды «solve» можно решать не только алгебраические уравнения. Например, решим тригоно-метрическое уравнение tg x – 2 sin x = 0: > solve(tan(x)–2*sin(x)); Обратите внимание, что Maple привёл решения, лежащие в пределах одного промежутка периодичности (от –π до π). Для вывода всех решений необходимо присвоить зарезервированной переменной _EnvAllSolutions значение true: > _EnvAllSolutions := true; > solve(tan(x)–2*sin(x)); где _Z

обозначает любое целое число.

Решение уравнений (часть 3)

Приведём примеры применения функции «solve» для решения уравнений с несколькими переменными. Решим, например, уравнение xy + x – 1 = 0 относительно x: > solve(x*y+x–1,x); относительно y: > solve(x*y+x–1,у); В общем виде Maple решает это уравнение так: > solve(x*y+x–1); Видно, что форма ответа определяется вторым параметром (или его отсутствием) команды «solve», указывающим, относительно какой переменной решать уравнение.

Неравенства решаются тем же оператором «solve». Например, решим неравенство x2(x – 1) solve(x^2*(x–1) solve(x^2*(x–1)>=0); 0, RealRange(1, ∞) В переводе на математический язык ответ: U[1; ∞). Открытый интервал (или луч) задаётся в Maple с помощью функции «Open», применяемой к концам интервала, задаваемого функцией «RealRange».

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

Все уравнения системы записываются в фигурных скобках через запятую. Решим например систему > solve(); , Решим систему с параметром > solve(,);