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

ЛИЧНЫЙ КАБИНЕТ 

 

Здравствуйте гость!

 

Логин:

Пароль:

 

Запомнить

 

 

Забыли пароль? Регистрация

Повышение уникальности

Предлагаем нашим посетителям воспользоваться бесплатным программным обеспечением «StudentHelp», которое позволит вам всего за несколько минут, выполнить повышение уникальности любого файла в формате MS Word. После такого повышения уникальности, ваша работа легко пройдете проверку в системах антиплагиат вуз, antiplagiat.ru, etxt.ru или advego.ru. Программа «StudentHelp» работает по уникальной технологии и при повышении уникальности не вставляет в текст скрытых символов, и даже если препод скопирует текст в блокнот – не увидит ни каких отличий от текста в Word файле.

Результат поиска


Наименование:


Реферат Предмет и методы изучения дифференциальной векторно-матричной алгебры, ее структура. Векторное решение однородных и неоднородных дифференциальных уравнений. Численное решение векторно-матричных уравнений. Формулы построения вычислительных процедур.

Информация:

Тип работы: Реферат. Предмет: Математика. Добавлен: 15.08.2009. Сдан: 2009. Уникальность по antiplagiat.ru: --.

Описание (план):


4
Реферат
на тему:
"Решение систем дифференциальных уравнений"

1. Дифференциальная линейная алгебра

С собственными значениями и векторами матрицы приходится иметь дело в задачах, связанных с решением систем линейных дифференциальных уравнений и исследованием устойчивости этих решений. Дифференциальная векторно-матричная алгебра включает в себя операции интегрирования и дифференцирования, которые во множестве случаев в своей нотации напоминают соответствующие операции обычного дифференциального исчисления. Производная по скалярной переменной и интеграл от вектора и матрицы в заданных пределах изменения скалярной переменной определены так:
Производные от векторных и векторно-матричных выражений определяются следующими правилами:
,
,
,
,
.
2. Векторное решение однородного уравнения
Пусть система линейных однородных дифференциальных уравнений задана в векторной форме:
Если уравнение записано в форме однородного дифференциального уравнения n-го порядка и его характеристический многочлен имеет различные корни, то общее решение представляется суммой n частных решений с экспоненциальными базовыми функциями:
,
где - константы, определяемые начальными условиями.
Можно предположить, что векторное уравнение, представляющее общее решение, имеет аналогичную форму
.
Для выяснения вопроса, что есть в таком представлении и , подставим частное решение в уравнение:
Отсюда видно, что будет частным решением, если будут собственным значением и собственным вектором матрицы A.
Таким образом, если матрица A имеет собственные значения и векторы , k=1,2,, n, то общее решение однородного векторного уравнения после ряда эквивалентных преобразований предстанет в следующем виде:
.
Используя значение решения при t=0, находим . Таким образом, общее решение однородного векторного уравнения имеет следующий вид: .
Матричная экспонента выражается через проекторы и собственные значения матрицы по формулам спектрального разложения:
.
После подстановки X в решение вместо экспоненты получим:
.
В случаях, когда собственные значения и векторы найти не удается, матричную функцию можно разложить в ряд по степеням матрицы:
,
что позволяет численно получать многомерный переходной процесс, если ряд сходится.
Матричный ряд сходится, если существует предел последовательности частичных сумм. Достаточным условием является сходимость ряда из норм членов степенного матричного ряда. Используя, например, признак сходимости Даламбера ряд, представляющий матричную экспоненту, сходится, если существует и меньше единицы предел отношения
,
где R - радиус сходимости.
Объем вычислительной работы при оцифровке многомерного переходного процесса существенно зависит от числа членов в матричном ряде. Для повышения скорости сходимости применяют различные аппроксимации этого ряда. В частности, для экспоненты широко используются аппроксимации отрезков ряда дробно-рациональными функциями Падэ вида:
.
Так, матричная экспонента для трех и четырех членов имеет вид:
В свете приведенных разложений матричной экспоненты общее решение линейного векторно-матричного дифференциального уравнения приближенно можно вычислить по формуле:
.
3. Решение неоднородных дифференциальных уравнений
Познакомившись с общим подходом к построению решений линейных векторных дифференциальных уравнений, покажем теперь, как получаются решения неоднородных уравнений.
Представим исходное уравнение с неоднородностью, локализованной в правой части уравнения, и умножим обе части уравнения на матричную экспоненту :
.
Обращаясь к правилам дифференцирования векторно-матричных выражений, приведенных выше, несложно заметить, что слева от знака равенства находится производная от произведения матричной экспоненты на вектор y:
.
Сделаем соответствующую замену и проинтегрируем левую и правую части по независимой переменной t:
.
Умножая слева обе части равенства на матрицу , получим общее решение неоднородного дифференциального уравнения:
.
Формула общего решения в своей нотации точно соответствует случаю скалярного уравнения. При невозможности аналитического решения переходный процесс можно вычислить по точкам, заменив непрерывное время дискретным с шагом , где R - радиус сходимости степенного матричного ряда с матрицей :
.
В интеграле можно заменить независимую переменную на дискретную с тем же шагом, что и при разложении экспоненты: , тогда, применяя метод интегрирования по правилу прямоугольников и обозначая матричную экспоненту на k-том шаге через , получим
.
Удобно из формулы вычисления дискретных значений векторного переходного процесса получить рекуррентную формулу. Этого можно добиться, если найти в выражении для часть, которую можно заменить значением :
Повышения точности вычисления переходного процесса достигают за счет замены интеграла квадратурами более высокого порядка, например, первого - формула трапеций, или второго - формула парабол (Симпсона).
Использование формулы трапеций приводит после соответствующих преобразований к следующей рекуррентной формуле:
Если использовать формулу Симпсона, то рекуррентная формула для расчета переходного процесса от точки к точке будет такой:
В приведенных рекуррентных формулах матричные экспоненты имеют следующий вид:
.
4. Примеры численного решения векторно-матричных уравнений
В качестве примера построим переходный процесс для системы уравнений:
.
Эта система может быть представлена дифференциальным уравнением второго порядка относительно переменной :
,
или относительно переменной :
.
Характеристическое уравнение имеет два комплексных корня: . Общее решение этих уравнений будет:
,
где - постоянные, которые вычисляются по заданным начальным условиям путем решения системы уравнений:
Несложные преобразования приводят к следующим точным решениям этого уравнения для двух различных наборов начальных условий:
Получим такое же аналитическое решение векторного переходного процесса в форме экспоненциальной функции, используя спектральное разложение матрицы по собственным значениям.
Характеристический полином заданной матрицы имеет вид:
.
Собственные значения матрицы (корни характеристического уравнения) и собственные векторы равны:
Проекторы находим матричным произведением левых и правых собственных векторов. Для этого обратим матрицу и в качестве левых собственных векторов возьмем ее строки:
Векторное аналитическое решение имеет вид:
Решение совпадает с точным решением уравнений второго порядка.
Для численного построения векторного переходного процесса по заданному векторно-матричному уравнению с использованием Падэ-аппроксимации матричной экспоненты дробно-рациональными выражениями первого, второго и третьего порядков, вычислим сначала эти аппроксимирующие матрицы:
Вектор приближенного решения вычислим по рекуррентной формуле, в которую, для демонстрации влияния на точность результата, поочередно подставим каждое из трех приведенных выше приближений к матричной экспоненте:
:

В таблице помещены численные значения переходных процессов, полученные для трех названных случаев аппроксимации матричной экспоненты вместе с точным аналитическим решением.
t
Аналитическое
решение
Аппроксимация
Падэ порядка 1
Аппроксимация
Падэ порядка 2
Аппроксимация
Па и т.д.................


Перейти к полному тексту работы



Смотреть похожие работы


* Примечание. Уникальность работы указана на дату публикации, текущее значение может отличаться от указанного.