49461 (ЭВМ с использованием математического пакета MathCad в среде Windows 98 для решения системы дифференциальных уравнений), страница 2

Поиск обратной матрицы возможен, если матрица квадратная и ее определитель не равен нулю. Произведение исходной матрицы на обратную по определению является единичной матрицей. Для ввода оператора поиска обратной матрицы нажмите кнопку Inverse (Обратная матрица) на панели инструментов Matrix (Матрица).

Начальные условия:

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

Ф=P*Q*P^-1

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

Тогда получим 4 решения:

Рисунок 1.1. Графики изменения переменных состояния системы при заданных начальных условиях и отсутствии внешнего воздействия.

2.3.2 Определение аналитических зависимостей изменения переменных состояния системы при заданных начальных условиях и отсутствии внешнего воздействия с использованием функции Mathcad

СПРАВКА: В Mathcad 11 имеются три встроенные функции, которые позволяют решать поставленную в форме (2—3) задачу Коши различными численными методами.

• rkfixed(y0, t0, t1, M, D) — метод Рунге-Кутты с фиксированным шагом,

• Rkadapt(y0, t0, t1, M, D) — метод Рунге-Кутты с переменным шагом;

• Buistoer(y0, t0, t1, M, D) — метод Булирша-Штера;

  • у0 — вектор начальных значений в точке to размера NXI;

  • t0 — начальная точка расчета,

  • t1 — конечная точка расчета,

  • M — число шагов, на которых численный метод находит решение;

  • D — векторная функция размера NXI двух аргументов — скалярного t и векторного у При этом у — искомая векторная функция аргумента t того же размера NXI.

Таким образом, воспользуемся функцией Rkadapt (y0, t0, t1, M, D) -получим матрицу решения системы дифференциальных уравнений численным методом Рунге-Кута на интервале от t0 до t1 при M шагах решения и правыми частями уравнений, записанными в D. Тогда решение уравнения динамики электротехнической системы с помощью встроенной функции Rkadapt выглядит так:

Зададим интервал интегрирования t0 - t1, количество шагов интегрирования М, вектор заданных начальных условий X0 и правую часть дифференциального уравнения y(t):

Сформируем матрицу системы дифференциальных уравнений:

Применим функцию:

-Интервал времени.

-Значение искомой координаты.

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

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

2.3.3 Определение аналитических зависимостей изменения переменных состояния системы при заданных начальных условиях и отсутствии внешнего воздействия с использованием преобразования Лапласа

Заданную систему уравнений преобразуем по Лапласу и найдем переходную матрицу и изображение по Лапласу переменной состояния системы:

На основании переходной матрицы определим изображение и оригинал переменных состояния систем:

Графики изменения переменных состояния во временной области при отсутствии внешних возмущений и заданных начальных условиях, полученные с помощью преобразования Лапласа представлены на рисунке 7.1.

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

Как видно рисунок 1.3. совпадает с рисунком 1.1, где неизвестные получены с помощью характеристического уравнения системы и рисунком 1.2.- численный метод с использованием функции MATHCAD.

2.4Частное решение неоднородной системы дифференциальных уравнений при заданном внешнем воздействии и нулевых начальных условиях

2.4.1 Решение с применением функций MATHCAD

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

2.4.2 Решение с применением преобразования Лапласа

Преобразуем по Лапласу заданную систему уравнений и найдем переходную матрицу и изображение переменной состояния системы.

B(s) – преобразованный по Лапласу вектор-столбец внешних возмущений.

Переходная матрица и изображение переменных состояния системы:

На основание матрицы определим изображение и оригинал переменных состояния системы:

Аналогично вычисляем остальные значения x(t)

Также применим обратное проеобразование Лапласа , нажав ключевое слово invlaplace на панели Символика.

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

Как видно графики совпадают.

2.5 Частное решение неоднородной системы дифференциальных уравнений при заданном внешнем воздействии y=cos(2t) и нулевых начальных условиях

2.5.1 Решение с помощью переходной матрицы

В качестве примера рассмотрим случай, если на систему действует воздействие одного вида, например y=cos(2t) .

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

пусть

Рисунок 3.1. Графики изменения переменных состояния системы при при y(t)=cos(2t) и нулевых начальных условиях, полученные способом решения с использованием переходной матрицы.

2.5.2 Численный метод решения системы дифференциальных уравнений при нулевых начальных условиях и заданном внешнем воздействии y=cos(2t) c помощью MATHCAD

Рисунок 3.2. Графики изменения переменных состояния системы при нулевых начальных условиях и воздействии y=cos(2t)

Как видно из графиков решения совпадают.

2.5.3 Решение системы дифференциальных уравнений при нулевых начальных условиях и заданном внешнем воздействии y=cos(2t) c помощью преобразования Лапласа

Применив обратное преобразование Лапласа (invlaplace) получим значения x(t), графическое изображение которых на рисунке 3.3. Рисунок совпадает с двумя полученными ранее.

Рисунок 3.3. Графики изменения переменных состояния системы при при y(t)=cos(2t) и нулевых начальных условиях, полученные с помощью преобразования Лапласа.

2.6 Решение неоднородной системы дифференциальных уравнений при заданном внешнем воздействии и начальных условиях

2.6.1 Решение с помощью функции MATHCAD

Рисунок 4.1. Графики изменения переменных состояния системы при заданных начальных условиях и воздействии, полученных с помощью функции MATHCAD.

2.6.2 Решение с помощью преобразования Лапласа

Аналогично получаем значения Х2общ, Х3общ, Х4общ и строим график изменения переменных системы при заданном внешнем воздействии и начальных условиях, полученный с помощью преобразования Лапласа.

Рисунок 4.2. Графики изменения переменных состояния системы при заданных начальных условиях и воздействии, полученных с помощью преобразования Лапласа..

2.6.3 Решение с помощью преобразования Лапласа (способ второй)

Применив обратное преобразование Лапласа, получим изменение параметров системы в зависимости от времени.

Рисунок 4.3. . Графики изменения переменных состояния системы при заданных начальных условиях и воздействии, полученных с помощью преобразования Лапласа.

Выводы по работе №4

В данной работе были изучены возможности математического пакета MathCad в среде Windows для решения системы дифференциальных уравнений, что часто используется в инженерных расчетах электротехнических систем. Были выполнены решения системы дифференциальных уравнений численным методом и с использование преобразования Лапласа, используя математический пакет MathCad.

Классическим методом расчёта является метод расчёта с использованием переходной матрицы.

Решение уравнения динамики электротехнической системы с помощью встроенной функции Rkadapt очень наглядно и быстро. Воспользовавшись функцией Rkadapt (y0, t0, t1, M, D)-получим матрицу решения системы дифференциальных уравнений численным методом Рунге-Кута на интервале от t0 до t1 при M шагах решения и правыми частями уравнений, записанными в D.

Преобразование Лапласа позволяет преобразовывать дифференциальные уравнения по t в линейные уравнения по S. Переменные вещественного аргумента t меняется на переменные комплексного аргумента s. Дифференцирование заменяется умножением на s, повторное- на s в квадрате и т.д.С помощью laplace находим изображения функций, описывающих внешние воздействия на систему.

Решение с помощью преобразования Лапласа требует создания переходной матрицы. Этот вопрос решается очень легко, используя функцию identity из mathcad. Для нахождения изменения параметров системы в зависимости от времени используется обратное преобразование Лапласа (invlaplace).

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

35