Решение систем уравнений с экспонентой

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

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

Выразим у через х из (2) -го уравнения системы и подставим это значение в (1) -ое уравнение системы.

Решаем (2) -ое уравнение полученной системы:

2 х +2 x +2 =10, применяем формулу: a x + y =a x ∙a y .

2 x +2 x ∙2 2 =10, вынесем общий множитель 2 х за скобки:

2 х (1+2 2 )=10 или 2 х ∙5=10, отсюда 2 х =2.

2 х =2 1 , отсюда х=1. Возвращаемся к системе уравнений.

Ответ: (1; 2).

Представляем левую и правую части (1) -го уравнения в виде степеней с основанием 2, а правую часть (2) -го уравнения как нулевую степень числа 5.

Если равны две степени с одинаковыми основаниями, то равны и показатели этих степеней — приравниваем показатели степеней с основаниями 2 и показатели степеней с основаниями 5.

Получившуюся систему линейных уравнений с двумя переменными решаем методом сложения.

Находим х=2 и это значение подставляем вместо х во второе уравнение системы.

Находим у.

Ответ: (2; 1,5).

Если в предыдущих двух примерах мы переходили к более простой системе приравнивая показатели двух степеней с одинаковыми основаниями, то в 3-ем примере эта операция невыполнима. Такие системы удобно решать вводом новых переменных. Мы введем переменные u и v, а затем выразим переменную u через v и получим уравнение относительно переменной v.

Решаем (2) -ое уравнение системы.

v 2 +63v-64=0. Подберем корни по теореме Виета, зная, что: v₁+v₂=-63; v₁∙v₂=-64.

Получаем: v₁=-64, v₂=1. Возвращаемся к системе, находим u.

Так как значения показательной функции всегда положительны, то уравнения 4 x = -1 и 4 y = -64 решений не имеют.

Представляем 64 и 1 в виде степеней с основанием 4.

Приравниваем показатели степеней и находим х и у.

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

Системы уравнений по-шагам

Видео:Решение систем уравнений второго порядка. 8 класс.Скачать

Результат

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

• Метод Гаусса
• Метод Крамера
• Прямой метод
• Система нелинейных уравнений

Указанные выше примеры содержат также:

• квадратные корни sqrt(x),
  кубические корни cbrt(x)
• тригонометрические функции:
  синус sin(x), косинус cos(x), тангенс tan(x), котангенс ctan(x)
• показательные функции и экспоненты exp(x)
• обратные тригонометрические функции:
  арксинус asin(x), арккосинус acos(x), арктангенс atan(x), арккотангенс actan(x)
• натуральные логарифмы ln(x),
  десятичные логарифмы log(x)
• гиперболические функции:
  гиперболический синус sh(x), гиперболический косинус ch(x), гиперболический тангенс и котангенс tanh(x), ctanh(x)
• обратные гиперболические функции:
  asinh(x), acosh(x), atanh(x), actanh(x)
• число Пи pi
• комплексное число i

Правила ввода

Можно делать следующие операции

2*x — умножение 3/x — деление x^3 — возведение в степень x + 7 — сложение x — 6 — вычитание Действительные числа вводить в виде 7.5, не 7,5

Чтобы увидеть подробное решение,
помогите рассказать об этом сайте:

Видео:9 класс, 11 урок, Методы решения систем уравненийСкачать

Как не стоит вычислять матричную экспоненту

Пост написан под влиянием поста пользователя pchelintsev_an.

В данной статье я постараюсь рассказать, с какими вычислительными трудностями можно столкнуться, если пойти по «наивному» пути вычисления матричной экспоненты. Статья может быть полезна тем, кто хотел бы познакомиться с вычислительной математикой, но уже знаком с такими понятиями как система обыкновенных дифференциальных уравнений и задача Коши. Эксперименты проводились с использованием системы GNU Octave.

Что вообще такое матричная экспонента и для чего она применятеся

Матричная экспонента возникает при рассмотрении задачи Коши для линейной системы обыкновенных дифференциальных уравнений с постоянными коэффициентами:

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

Решение системы ОДУ можно выписать явно:

Здесь появляется матричная экспонента, которая определяется как сумма ряда:

Фактически, матричная экспонента получается формальной подстановкой матрицы A в ряд Тейлора для e z (Математик отметит, что ряд сходится абсолютно для любой матрицы и любого t). Удобно то, что с помощью матричной экспоненты можно вычислить решение в любой момент времени t. Например, если эту систему решать каким-нибудь методом численного интегрирования, например, Адамса или Рунге-Кутты, то затраченное время будет пропорционально времени t. Однако, решение через матричную экспоненту возможно лишь для автономных систем с постоянными коэффициентами.

В курсе дифференциальных уравнений для практических вычислений матричной экспоненты предлагается метод приведения матрицы к жордановой форме: если матрицу A представить в виде S J S -1 , где J – жоржанова форма, S – матрица перехода к жорданову базису, то матричная экспонента А выражается через матричную экспоненту J:

Матричная экспонента от жордановой матрицы J выписывается в явном виде. Однако такой способ требует вычисления системы собственных векторов, а также присоединенных к ним векторов. Эта задача значительно сложнее исходной. К тому же, даже для действительной матрицы А, ее жорданова форма J может быть комплексной.

Суммирование ряда

Итак, у нас есть бесконечный ряд и желание вычислить его за обозримое время. Очевидный путь – закончить суммирование на каком-то слагаемом, отбросив остаток. При этом суммирование ряда заканчивается, если очередное k-е слагаемое по норме мало:

Нормы матрицы

Остановимся ненадолго на норме матрицы. Различных норм матрицы существует достаточно много. Можно ввести норму, исходя из требования

В этом случае говорят о норме матрицы индуцированной (или подчиненной) некоторой векторной норме. Для наиболее употребимых векторных норм соответствующие матричные нормы имеют вид:

Первые две нормы матрицы считаются элементарно, а для последней требуется вычислить сингулярное число, что весьма непросто. Кроме индуцированных норм существует еще ряд других норм. Практически удобной является норма Фробениуса, которая, с одной стороны, легко вычисляется, а с другой – тесно связана с сингулярными числами А. Фактически, это евклидова норма вектора, полученного «вытягиванием» матрицы в одну строку:

Сюрпризы при суммировании ряда

Чтобы осознать сложности, возникающие при суммировании рядов Тейлора, не обязательно обращаться к примеру с матричной экспонентой. Попробуем просто построить график функции e -x на отрезке [0, 50], суммируя аналогичным образом ряд Тейлора для экспоненты. Суммирование завершим, если очередное слагаемое не превосходит 10 -20 . Воспользуемся следующим кодом на Octave:

Результат получается довольно странным:

При x > 30 график функции начинает вести себя хаотично и не имеет с экспонентой ничего общего.

Причина данного эффекта в накоплении погрешностей округления. Например, при вычислении e -35 было просуммировано 132 слагаемых, из которых минимум абсолютной величины имело последнее слагаемое s₁₃₂ ≈ 5.8674 · 10 -21 , а максимум абсолютной величины имело слагаемое с номером 35: s₃₅ ≈ -1.067 · 10 14 . Поскольку Octave проводит вычисления с двойной точностью (примерно 16 значащих цифр), округление 35-го слагаемого уже внесло ошибку порядка 0.01, что существенно превосходит и значение ϵ = 10 -20 и само значение суммы ряда e -35 ≈ 6.305 · 10 -16

Неужели ряд Тейлора настолько безнадежен? Нет, ошибки округления не будут давать существенного вклада, если величина суммируемых слагаемых не будет достигать таких огромных значений. Например, такой алгоритм будет надежно работать при небольших значениях x (|x| ≲ 1). Для матричной экспоненты аналогичное условие имеет вид . Кроме вычислительной устойчивости такой выбор x обеспечивает быструю сходимость, то есть достаточно взять всего несколько слагаемых из ряда для получения приемлемой точности.

Метод удвоения аргумента

Заметим, что матричная экспонента удовлетворяет следующему соотношению:

Это означает, что матричную экспоненту можно посчитать для матрицы , а затем возвести результат в степень N. Взяв в качестве N = 2 p можно воспользоваться быстрым алгоритмом возведения в степень за p матричных умножений:

1. Вычислить матричную экспоненту
2. Умножить B на себя , повторив эту операцию p раз
3. В результате,

Число p можно (и нужно) выбрать из соображений вычислительной устойчивости шага 1 алгоритма, то есть

Соответствующий Octave код:

Метод аппроксимаций Паде

Отрезок ряда Тейлора – не единственный способ приближенно вычислить значение функции. Более качественные приближения функции могут быть получены не в виде многочленов, а в виде рациональных функций. Наилучшие приближения в виде отношения двух многочленов заданных степеней называются аппроксимациями Паде. Например, аппроксимация Паде отношением двух квадратичных многочленов для экспоненты записывается в виде:

Аналогичная аппроксимация для матричной экспоненты будет иметь вид:

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

Анализ точности

Как убедиться, что матричная экспонента вычислена правильно? Достаточно ли, что или требуется более «интеллектуальная» проверка?

Вспомним, что матричная экспонента возникла не на ровном месте, а в результате решения некоторой системы обыкновенных дифференциальных уравнений. Собственные числа матрицы A характеризуют скорости процессов, протекающих в системе, описываемой данной системой дифференциальных уравнений. Например, матрица:

описывает систему с одним быстро затухающим процессом (собственное значение -1000) и двумя колебательными процессами (собственные значения ±i). Хорошим критерием правильности вычисления матричной экспоненты может являться соотношение:

связывающее собственные значения матрицы A и ее матричной экспоненты.

Вычислительный эксперимент

Для начала исследуем алгоритмы на быстродействие на случайных матрицах разного размера. В качестве методов используем метод суммирования ряда Тейлора, метод удвоения аргумента и встроенную в Octave функцию expm. Изначально хотелось сравнить также соответствие собственных значений, но от этой идеи пришлось отказаться, поскольку все три алгоритма одинаково проваливали эту проверку при размерах матрицы N ≳ 50. Возможно виновником этого была функция определения собственных значений eig, работающая недостаточно точно или плохая обусловленность самой задачи определения собственных значений для этих матриц.

Исходя из графика, метод удвоения работает примерно вдвое дольше библиотечной функции (возможно стоит ослабить завышенное требование по точности), в то время как прямое суммирование имеет даже иную асимптотику ≈O(N 3.7 ), в то время, как остальные имеют асимптотику O(N 3 ).

Далее, проверим как алгоритм работает с большими по величине собственными значениями. Рассмотрим следующую систему обыкновенных дифференциальных уравнений:

С начальным условием:

.

Её особенность в том, что одно из ее собственных значений равно -1000 (быстро затухающее решение), а два других равны ±i√2 (колебательные решения). Если в мнимые собственные значения будут внесены искажения, это будет хорошо заметно. Возьмем шаг по времени τ = 0.038, найдем матричную экспоненту разными способами. Решение системы на сетке с шагом τ получается следующим простым алгоритмом:

1. Решение в начальный момент известно из начального условия
2. Зная решение u_n в момент t_n, решение в момент t_n+1 = t_n + τ получается из него умножением на матричную экспоненту B:

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

Все траектории, кроме полученной с помощью прямого суммирования ряда Тейлора, визуально совпадают. Обратите внимание, график не выходит из начального условия x₁(0) = 1. На самом деле, график экспоненциально подходит к точке y = 1, просто характерное расстояние «поворота» графика (∼ 1 / 1000) намного меньше шага по времени τ = 0.038.

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

Дальнейшее чтение

О различных способах вычисления матричной экспоненты можно почитать в статье (англ.) [2]. По методам решения обыкновенных дифференциальных уравнений и жестких задач написан замечательный двухтомник [3,4]

🌟 Видео

Решение системы уравнений методом Крамера.Скачать

Решение систем уравнений второй степени. Алгебра, 9 классСкачать

Решение систем уравнений методом подстановкиСкачать

Решение системы уравнений методом ГауссаСкачать

Cистемы уравнений. Разбор задания 6 и 21 из ОГЭ. | МатематикаСкачать

Решение систем уравнений методом сложенияСкачать

Решение системы уравнений методом Крамера 2x2Скачать

Задание №20. Экзамен ОГЭ. Система уравнений #shortsСкачать

Матричный метод решения систем уравненийСкачать

Метод Крамера за 3 минуты. Решение системы линейных уравнений - bezbotvyСкачать

Решение систем уравнений методом сложенияСкачать

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

Графический способ решения систем уравнений. Алгебра, 9 классСкачать

Решение системы линейных уравнений графическим методом. 7 класс.Скачать

Решение системы уравнений методом Гаусса 4x4Скачать

Решение систем уравнений. Методом подстановки. Выразить YСкачать

Решение системы уравнений методом обратной матрицы - bezbotvyСкачать