Аналитическое решение системы дифференциальных уравнений

Аналитическое решение системы дифференциальных уравнений.

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

Рассмотрим систему, состоящую из одного ремонтируемого элемента и имеющую два состояния:

S₀(t) –– система в момент t находится в рабочем состоянии;

S₁(t) –– система в момент t находится в ремонте ( не работает).

Система уравнений Колмогорова имеет вид:

(П1.35)
Найдем интеграл этой системы уравнений, для чего введем обозначения P₀(t)=y, P₁(t)=z.

Тогда система уравнений принимает вид:

(П1.36)
Дифференцируем первое уравнение:

(П1.37)
Из первого же уравнения находим z:

(П1.38)
Из второго уравнения системы после подстановки найденного выражения для z имеем:

(П1.39)
Подставляем полученное выражение производной (П1.39) в уравнение (П1.37), получим:

(П1.40)
Характеристическое уравнение имеет вид:

(П1.41)
Решение характеристического уравнения дает корни:

, откуда

k₁=0, k₂=-(l+m) (П1.42)
Общее решение системы однородных дифференциальных уравнений при вещественных и неравных корнях имеет вид:

, (П1.43)
(П2.10)

Учитывая начальные условия

y(0)=1 z(0)=0 (П1.44)
получим:

(П1.45)
Тогда

(П1.46)
Дифференцируем (П2.10):

(П1.47)
Подставляя в (П2.4) выражения (П2.13) и (П2.14), получим:

(П1.48)

При подстановке начального условия z(0) = 0, получим уравнение:

(П1.49)
Подставляя (П2.16) в (П2.10) и (П2.15), получим окончательно:

(П1.50)
(П1.51)

Кроме рассмотренного метода, могут использоваться и другие, например, с помощью преобразований Лапласа.

Заметим, что с пакетами Mathcad или Matlab системы дифференциальных уравнений интегрируются проще [16, 17, 18>.

Аналитическое решение системы дифференциальных уравнений

Общее решение дифференциальных уравнений.

Для нахождения аналитических решений дифференциальных уравнений в Maple применяется команда dsolve(eq,var,options), где eq – дифференциальное уравнение, var – неизвестные функции, options – параметры. Параметры могут указывать метод решения задачи, например, по умолчанию ищется аналитическое решение: type=exact . При составлении дифференциальных уравнений для обозначения производной применяется команда diff , например, дифференциальное уравнение y» + y = x записывается в виде: diff(y(x),x$2)+y(x)=x.

Общее решение дифференциального уравнения зависит от произвольных постоянных, число которых равно порядку дифференциального уравнения. В Maple такие постоянные, как правило, обозначаются как _ С1 , _ С2 , и т.д.

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

Команда dsolve выдает решение дифференциального уравнения в невычисляемом формате. Для того, чтобы с решением можно было бы работать далее (например, построить график решения) следует отделить правую часть полученного решения командой rhs(%) .

Задание 1.1.

1. Найти общее решение дифференциального уравнения y ‘+ y cos x =sin x cos x .

de : =

Итак, решение искомого уравнения есть функция 1 .

Замечание : при записи решения диффреренциального уравнения в Maple в строке вывода произвольная постоянная обозначена как _ С1 .

2. Найти общее решение дифференциального уравнения второго порядка y » — 2 y ‘+ y =sin x + e — x .

deq :=

Замечание : так как исходное уравнение было второго порядка, то полученное решение содержит две произвольные константы, которые в Maple обычно обознаются как _ С1 и _ С2 . Первые два слагаемых представляют собой общее решение соответствующего однородного дифференциального уравнения, а вторые два – частное решение неоднородного дифференциального уравнения.

3. Найти общее решение дифференциального уравнения порядка y »+ k 2 y =sin( qx ) в двух случаях: q ¹ k и q = k (резонанс).

de :=

Теперь найдем решение в случае резонанса. Для этого перед вызовом команды dsolve следует приравнять q = k .

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

Фундаментальная (базисная) система решений.

Команда dsolve представляет возможность найти фундаментальную систему решений (базисные функции) дифференциального уравнения. Для этого в параметрах команды dsolve следует указать output=basis .

Задание 1.2.

Найти фундаментальную систему решений дифференциального уравнения: y (4) +2 y »+ y =0.

de : =

> dsolve(de, y(x), output=basis);

Решение задачи Коши или краевой задачи.

Команда dsolve может найти решение задачи Коши или краевой задачи, если помимо дифференциального уравнения задать начальные или краевые условия для неизвестной функции. Для обозначения производных в начальных или краевых условиях используется дифференциальный оператор , например, условие y»(0)=2 следует записать в виде , или условие y ‘(1)=0: . Напомним, что производная n -го порядка записывается в виде .

Задание 1.3.

1. Найти решение задачи Коши: y (4) + y »=2cos x , y (0)= — 2, y ‘(0)=1, y »(0)=0, y »'(0)=0.

cond:= y(0)= — 2, D(y)(0)=1, (D (2) )(y)(0)=0, (D (3) )(y)(0)=0

y( x )= — 2cos( x ) — x sin( x )+ х

2. Найти решение краевой задачи: , , . Построить график решения.

de : =

y( x )=2 x — p + p cos( x )

Замечание : для построения графика решения предварительно следует отделить правую часть полученного выражения.

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

Команда dsolve может найти решение системы дифференциальных уравнений (или задачи Коши), если в ней указать: dsolve(,), где sys — система дифференциальных уравнений, x(t),y(t),… — набор неизвестных функций.

Задание 1.4.

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

Найдены две функции x ( t ) и y ( t ), которые зависят от двух произвольных постоянных _ С1 и _ С2 .

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

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

Чтобы найти приближенное решение дифференциального уравнения в виде степенного ряда, в команде dsolve следует после переменных указать параметр type=series (или просто series ). Для того, чтобы указать порядок разложения n , т.е. порядок степени, до которой производить разложение, следует перед командой dsolve вставить определение порядка с помощью команды Order:=n .

Если ищется общее решение дифференциального уравнения в виде разложения в степенной ряд, то коэффициенты при степенях х найденного разложения будут содержать неизвестные значения функции в нуле y(0) и ее производных D(y)(0), (D@@2)(y)(0) и т.д. Полученное в строке вывода выражение будет иметь вид, похожий на разложение искомого решения в ряд Маклорена, но с другими коэффициентами при степенях х . Для выделения частного решения следует задать начальные условия y(0)=у1, D(y)(0)=у2, (D@@2)(y)(0)=у3 и т.д., причем количество этих начальных условий должно совпадать с порядком соответствующего дифференциального уравнения.

Разложение в степенной ряд имеет тип series , поэтому для дальнейшей работы с этим рядом его следует преобразовать в полином с помощью команды convert(%,polynom) , а затем выделить правую часть полученного выражения командой rhs(%) .

Задание 1.5.

1. Найти решение задачи Коши: , в виде степенного ряда с точностью до 5-го порядка.

y(0)=0>, y(x), type=series);

В полученном решении слагаемое означает, что точность разложения была до 5-го порядка.

2. Найти общее решение дифференциального уравнения y »( х ) — y 3 ( х )= е — х cos x , в виде разложения в степенной ряд до 4-го порядка. Найти разложение при начальных условиях: y (0)=1, y ‘(0)=0.

> restart; Order:=4: de:=diff(y(x),x$2)-

Замечание : в полученном разложении запись D(y)(0) обозначает производную в нуле: y ‘(0). Для нахождения частого решения осталось задать начальные условия:

3. Найти приближенное решение в виде степенного ряда до 6-го порядка и точное решение задачи Коши: , , , . Построить на одном рисунке графики точного и приближенного решений.

de : =

cond :=y(0)=1, D(y)(0)=1, D (2) (y)(0)=1

y( x )=

Замечание : тип решения дифференциального уравнения в виде ряда есть series , поэтому для дальнейшего использования такого решения (вычислений или построения графика) его обязательно следует конвертировать в полином с помощью команды convert

На этом рисунке видно, что наилучшее приближение точного решения степенным рядом достигается примерно на интервале — 1 x

Исправляем ошибки: Нашли опечатку? Выделите ее мышкой и нажмите Ctrl+Enter

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ УРАВНЕНИЙ И ИХ СИСТЕМ В ПАКЕТАХ СИМВОЛЬНОЙ МАТЕМАТИКИ

Обращаем Ваше внимание, что в соответствии с Федеральным законом N 273-ФЗ «Об образовании в Российской Федерации» в организациях, осуществляющих образовательную деятельность, организовывается обучение и воспитание обучающихся с ОВЗ как совместно с другими обучающимися, так и в отдельных классах или группах.

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ

ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ

ПЕДАГОГИЧЕСКИЙ УНИВЕРСИТЕТ ИМЕНИ М.Е. ЕВСЕВЬЕВА»

Кафедра информатики и вычислительной техники

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ УРАВНЕНИЙ И ИХ СИСТЕМ В ПАКЕТАХ СИМВОЛЬНОЙ МАТЕМАТИКИ

Автор работы _____________________________________И. Ю. Добрынькина

Направление подготовки 44.03.05 Педагогическое образование

Профиль Информатика. Математика

Руководитель работы_______________________________ Т. В. Кормилицына

Введение

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

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

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

В начале 90-х гг. на смену универсальным языкам программирования пришли специализированные системы компьютерной математики (СКМ). Среди них наибольшую известность получили системы Eureka, Mercury, Mathcad, Derive, Mathematica 2/3/4, Maple V R3/R4/R5 и Maple 6 и др.

Научное программное обеспечение и математические пакеты играют важную роль в современном естествознании и технике. Такие пакеты как Axiom, Derive, Maсsyma, Maple, MatLab, MathCAD, Mathematica широко распространились в университетах, исследовательских центрах и компаниях развитых стран. Владение одним или несколькими математическими пакетами и регулярное использование их в работе будь то исследовательская или преподавательская задача быстро становится нормой для специалиста.

1. Mathematica . Решение простейших дифференциальных уравнений

Для решения дифференциальных уравнений в аналитической форме в пакете Mathematica используется функция DSolve, дифференциальное уравнение 29 относительно функции y(x). Функция y и все ее производные должны быть записаны с аргументом, заключенным в квадратные скобки: y[x], y’[x]

Функция DSolve стремится найти общее решение ДУ в явном виде и выдает результат в виде списка правил замены, причем каждое решение заключается в фигурные скобки. Для ДУ порядка n общее решение содержит n произвольных констант, которые обозначаются C[1], C[2],…,C[n]. Для получения частного решения необходимо в качестве первого аргумента DSolve указать список, состоящий из самого уравнения и начальных или граничных условий:

Найденные с помощью DSolve решения можно подставить в любое выражение, содержащее y(x). Однако это решение не определяет правил замены производных y’(x), y’’(x) и так далее, например:

Чтобы получить решение, не имеющее этого недостатка, нужно в качестве второго аргумента функции DSolve записать только имя искомой функции, не указывая ее аргумент. В этом случае решение представляется в виде чистой функции («purefunction»-объекта), в котором роль аргумента x, в некоторых случаях, играет символ «#1», а признаком этого объекта является символ «&». Полученное решение можно подставить в любое выражение, содержащее как функцию y(x), так и ее производные:

Для решения систем уравнений в качестве первого аргумента функции указывается список уравнений, а в качестве второго аргумента – список искомых функций:

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

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

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

1.2 Примеры из математического анализа

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

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

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

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

DSolve [Derivative [1] [у] [х] ==2*а*х^3, у[х], х]

DSolve [у» [х] — у’ [х] — 6 у [х] == 0, у [х] , х] <<У[х] ->| е-4хС[1] + С[2] -Cos[2x] -|sin[2x]>>

DSolve [у» [х] + 4 у'[х] == 10 Sin [2 х] , у [х] , х]

DSolve[y'[x] == Sin[Ex] , y[x] , x]

DSolvefz2 w»[z] +zw'[z] — (z2 + l)w[z] ==0, w[z], z]

<BesselI[l, z] C[l] +BesselK[l, z] C[2] >>

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

1. 3 Аналитическое решение дифференциальных уравнений

Общее решение дифференциальных уравнений.

Для нахождения аналитических решений дифференциальных уравнений в Maple применяется команда dsolve(eq,var,options),где eq – дифференциальное уравнение, var – неизвестные функции, options – параметры. Параметры могут указывать метод решения задачи, например, по умолчанию ищется аналитическое решение: type=exact. При составлении дифференциальных уравнений для обозначения производной применяется команда diff, например, дифференциальное уравнение y»+y=x записывается в виде: diff(y(x),x$2)+y(x)=x.

Общее решение дифференциального уравнения зависит от произвольных постоянных, число которых равно порядку дифференциального уравнения. В Maple такие постоянные, как правило, обозначаются как _С1, _С2, и т.д.

Команда dsolve выдает решение дифференциального уравнения в невычисляемом формате. Для того, чтобы с решением можно было бы работать далее (например, построить график решения) следует отделить правую часть полученного решения командой rhs(%).

eq:=<2*x-3*y+5*z+7*t=1, 4*x-6*y+2*z+3*t=2, 2*x-3*y-11*z-15*t=1>: > s:=solve(eq,); s:= < z=-11/8 t , y=y , x=3/2 y - 1/16 t+1/2 >Для нахождения частного решения следует выполнить подстановку конкретного значения одной из переменных при помощи команды subs: > subs(,s); < z=(-11)/8, x=31/16 , 1=1>«>

2. Аналитические вычисления в Mathcad

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

Команды для выполнения аналитических вычислений в основном сосредоточены в меню Символика (Symbolics) и продублированы на аналогичной панели инструментов.

Чтобы упростить выражение (или часть выражения), надо выбрать его при помощи уголкового курсора и дать команду Символика > Упростить (Symbolics > Simplify). При этом выполняются арифметические действия, сокращаются общие множители и приводятся подобные члены, применяются тригонометрические тождества, упрощаются выражения с радикалами, а также выражения, содержащие прямую и обратную функции. Некоторые действия по раскрытию скобок и упрощению сложных тригонометрических выражений требуют применения команды Символика > Раскрыть/Расширить (Symbolics > Expand).

В меню Символика (Symbolics) предусмотрен ряд операций, ориентированных на выделенную переменную, использованную в выражении. Например, команда Solve (Решить) ищет корни функции, заданной данным выражением. В примере в аналитической форме получены все корни полинома второй степени: сначала применена команда solve для решения, а затем simplify для упрощения результата:

Другие возможности использования этого меню включают:

аналитическое дифференцирование и интегрирование: Символика > Переменная > Дифференцировать (Symbolics > Variable > Differentiate) и Символика > Переменная > Интегрировать (Symbolics > Variable > Integrate);

замена переменной: Символика > Переменная > Подставить (Symbolics > Variable > Substitute) – вместо переменной подставляется содержимое буфера обмена;

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

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

Отметим, что функция Find пытается найти решение в аналитической форме. В том случае, если до блока Given задать численно значения всех параметров, входящих в уравнения, а также начальные приближения для корней, то получим решение в числовом виде.

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

Любое аналитическое вычисление можно применить с помощью ключевого слова. Cписок ключевых слов

3. Решение систем дифференциальных уравнений в символьном виде в системе MATLAB

Для решения дифференциальных уравнений в форме Коши MatLAB имеет функцию dsolve(‘eqn1’,’eqn2’, …), которая возвращает аналитическое решение системы дифференциальных уравнений с начальными условиями. Они задаются равенствами eqni(вначале задаются уравнения, затем начальные условия).

По умолчанию независимой переменной считается ‘t’ . Можно использовать и другую переменную, включив ее в конец списка параметров функции dsolve. Символ D обозначает производную по независимой переменной, то есть d/dt, при этом D2 означает d^2/dt^2 и т.д.

Начальные условия задаются в виде равенств ‘y(a) = b’ или ‘Dy(a) = b’, где y — независимая переменная, a и b – константы. Если число начальных условий меньше, чем число дифференциальных уравнений, то в решений будут присутствовать произвольные постоянные С1, С2 и т.д. Вывод осуществляется в виде массива записей.

Обратите внимание, что уравнение, которое требуется решить, задано как строка, то есть взято в одинарные кавычки. Ответ представляет собой точное (символьное) решение 1+корень(5). Для получения числовых решений введите double (ans) или vpa (ans), чтобы отобразить больше знаков. Ввод с командой solve может также быть символьным выражением, но в этом случае программа MATLAB потребует, чтобы правая часть выражения была заключена в скобки, и фактически синтаксис решения уравнения х 2 — Зх = -7 будет выглядеть так:

Ответ представляет собой точное (символьное) решение (3 + корень(19i))/2 (сложные числа, где буква i в ответе ставится для мнимой единицы V-1). Для получения числовых решений введите double (ans) или vpa (ans), чтобы отобразить больше знаков.

Заключение

В настоящее время научное программирование претерпевает серьезную трансформацию: развиваются интегрированные среды, основанные на алгоритмических языках, и растет применение универсальных математических систем (Maple, Mathematica, MATLAB, MatCad и др.). Эти системы имеют дружественный интерфейс, реализуют множество стандартных и специальных математических операций, снабжены мощными графическими средствами и обладают собственными языками программирования. Все это предоставляет широкие возможности для эффективной работы специалистов разных профилей, о чем говорит активное применение математических пакетов в научных исследованиях и в преподавании. С помощью этих пакетов проще готовить и выполнять задания, устраивать демонстрации и гораздо быстрее решать исследовательские и инженерные задачи.

Конечным продуктом исследования выступают публикации, подготовка, распространение и использование которых в настоящее время требует квалифицированного применения компьютера. Это касается редактирования текста, изготовления графических материалов, ведения библиографии, размещения электронных версий в Интернет, поиска статей и их просмотра. Де-факто сейчас стандартными системами подготовки научно-технических публикаций являются различные реализации пакета TeX и текстовый редактор Word. Кроме того, необходимы минимальные знания о стандартных форматах файлов, конверторах, программах и утилитах, используемых при подготовке публикаций.

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

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

Список литературы

1. Дьяконов В.П. Справочник по применению системы PC MATLAB. — М.: «Физматлит» , 1993. — С. 112. — ISBN 5-02-015101-7

2. Дьяконов В.П. Компьютерная математика. Теория и практика. — СПб: «Питер» , 1999, 2001. — С. 1296. — ISBN 5-89251-065-4

3. Дьяконов В.П. MATLAB 5 — система символьной математики. — М.: «Нолидж» , 1999. — С. 640. — ISBN 5-89251-069-7

4. Дьяконов В.П., Абраменкова И.В. MATLAB. Обработка сигналов и изображений. Специальный справочник. — СПб.: «Питер» , 2002. — С. 608. — ISBN 5-318-00667-608

5. Дьяконов В.П., Круглов В.В. MATLAB. Анализ, идентификация и моделирование систем. Специальный справочник. — СПб.: «Питер» , 2002. — С. 448. — ISBN 5-318-00359-1

6. Дьяконов В. П. Simulink 4. Специальный справочник. — СПб.: «Питер» , 2002. — С. 528. — ISBN 5-318-00551-9

7. Дьяконов В . П . MATLAB 6/6.1/6.5 + Simulink 4/5. Основы применения. Полное руководство пользователя. — Москва.: «СОЛОН-Пресс» , 2002. — С. 768. — ISBN 5-98003-007-7

8. Дьяконов В.П. MATLAB 6/6.1/6.5 + Simulink 4/5 в математике и моделировании. Основы применения. Полное руководство пользователя. — Москва.: «СОЛОН-Пресс» , 2003. — С. 576. — ISBN 5-93455-177-9

9. Дьяконов В . П . MATLAB 6.0/6.1/6.5/6.5+SP1 + Simulink 4/5. Обработка сигналов и изображений. Полное руководство пользователя. — Москва.: «СОЛОН-Пресс» , 2005. — С. 592. — ISBN 5-93003-158-8

10. Дьяконов В . П . MATLAB 6.5/7.0 + Simulink 5/6. Основы применения. Библиотека профессионала. — Москва.: «СОЛОН-Пресс» , 2005. — С. 800. — ISBN 5-98003-181-2

11. Дьяконов В.П. MATLAB 6.5/7.0 + Simulink 5/6 в математике и моделировании. Библиотека профессионала. — Москва.: «СОЛОН-Пресс» , 2005. — С. 576. — ISBN 5-98003-209-6

12. Дьяконов В . П . MATLAB 6.5/7.0 + Simulink 5/6. Обработка сигналов и проектирование фильтров. Библиотека профессионала. — Москва.: «СОЛОН-Пресс» , 2005. — С. 576. — ISBN 5-98003-206-1

13. Дьяконов В . П . MATLAB 6.5/7.0/7 SP1 + Simulink 5/6. Работа с изображениями и видеопотоками. Библиотека профессионала. — Москва.: «СОЛОН-Пресс» , 2005. — С. 400. — ISBN 5-98003-205-3

14. Дьяконов В . П . MATLAB 6.5/7.0/7 SP1/7 SP2 + Simulink 5/6. Инструменты искусственного интеллекта и биоинформатики. Библиотека профессионала. — Москва.: «СОЛОН-Пресс» , 2005. — С. 456. — ISBN 5-98003-255-X

15. Дьяконов В . П . MATLAB R2006/2007/2008 + Simulink 5/6/7. Основы применения. Изд-е 2-е, переработанное и дополненное. Библиотека профессионала. — Москва.: «СОЛОН-Пресс» , 2008. — С. 800. — ISBN 978-5-91359-042-8

16. Дьяконов В.П. MATLAB 7.*/R2006/2007. Самоучитель. — Москва: «ДМК-Пресс» , 2008. — С. 768. — ISBN 978-5-94074-424-5

17. Дьяконов В.П. SIMULINK 5/6/7. Самоучитель. — Москва: «ДМК-Пресс» , 2008. — С. 784. — ISBN 978-5-94074-423-8

18. Дьяконов В.П. Вейвлеты. От теории к практике. Полное руководство пользователя. Изд-е 2-е переработанное и дополненное. — Москва: «СОЛОН-Пресс» , 2004. — С. 400. — ISBN 5-98003-171-5

19. Чарльз Генри Эдвардс, Дэвид Э. Пенни Дифференциальные уравнения и проблема собственных значений: моделирование и вычисление с помощью Mathematica, Maple и MATLAB = Differential Equations and Boundary Value Problems: Computing and Modeling. — 3- е изд . — М .: « Вильямс » , 2007. — ISBN 978-5-8459-1166-7

20. Алексеев Е.Р., Чеснокова О.В MATLAB 7. Самоучитель.. — Пресс , 2005. — С. 464.

21. Курбатова Екатерина Анатольевна MATLAB 7. Самоучитель. — М.: «Диалектика» , 2005. — С. 256. — ISBN 5-8459-0904-X

22. Джон Г. Мэтьюз, Куртис Д. Финк Численные методы. Использование MATLAB = Numerical Methods: Using MATLAB. — 3- е изд . — М .: « Вильямс » , 2001. — С . 720. — ISBN 0-13-270042-5 u

источники:

http://old.exponenta.ru/EDUCAT/SYSTEMAT/SAVOTCHENKO/6_1.asp

http://infourok.ru/analiticheskoe-reshenie-uravnenij-i-ih-sistem-v-paketah-simvolnoj-matematiki-4913344.html