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

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

КУРСОВАЯ РАБОТА

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

РЕФЕРАТ

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

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

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

Работа состоит из введения, 3 разделов, 2 таблиц, заключения, библиографического списка из 4 источников, одного приложения, в котором приведена реализация примеров решения уравнений.

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

1. Построение формального решения на входном Maple-языке

2. Метод разделения переменных

3. Метод функций Грина и другие методы

ЗАКЛЮЧЕНИЕ

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

ПРИЛОЖЕНИЕ А

ВВЕДЕНИЕ

Прикладной математический пакет MAPLE обладает большим набором инструментов для работы с дифференциальными уравнениями в частных производных. Среди них: установление порядка уравнения, исследование на возможность разделения переменных, определение условий поиска решения в виде суммы или произведения функций, получение решения из функций, получаемых командой pdsolve для разделенных уравнений, выполнение замены переменных и различных подстановок и т.п.

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

1. Построение формального решения на входном Maple-языке

Проблема решения дифференциальных уравнений в частных производных средствами MAPLE представляет собой программную задачу, сочетающую использование инструментов пакета с необходимыми дополнительными алгоритмами: учет начальных и граничных условий (НУ и ГУ), сложные и, зачастую, нетривиальные преобразования промежуточных результатов (основанные, например, на исследовании асимптотического поведения функций), программное использование дополнительной и/или специальной информации (например, использование рекуррентных соотношений для некоторых специальных функций, которые пока недоступны средствами MAPLE) и т.п. Более того, при решении сложных задач требуется программирование отдельных этапов решения с последующим объединением промежуточных результатов, а также создания комплексов программ (например, при комплексном аналитическом и численном - решении уравнений и различных способах визуализации и интерпретации результатов).

Для программирования построения формального решения на входном MAPLE-языке необходим ввод необходимой начальной информации (табл. 1) с последующим выполнением определенных алгоритмических операций (табл. 2).

Таблица 1

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

Тип информации	Содержание
Основная Информация	Вызов пакетов расширения.
	Задание системы координат.
	Ввод дифференциального уравнения в частных производных.
	Ввод начальных и граничных условий.
	Ввод различных функций и операторов.
	Вызов средств аналитического или численного решения уравнений.
Дополнительная информация	Представление функции при разделении переменных.
	Выполнение замены переменных(при необходимости).
	Переопределение постоянных, которые по умолчанию присваиваются пакетом.
	Ввод математической информации, недопустимой в Maple.
	Ввод специфических данных(физические параметры, габариты и т.д.).
	Ввод и вывод информации, связанной с текущим контролем выполняемых операций(получение результата для известного частного случая, контроль другими средствами).
	Ввод информации о форме представления результата (экспоненциальная, тригонометрическая и т.п. формы решения).
	Ввод информации для исследования промежуточных и конечных результатов (о порядке разложения в ряд, асимптотике, сравнениях и т.п.).
Рабочая информация	Последовательность вывода полученных результатов.
	Форматы переменных и данных.
	Вывод промежуточных результатов.
	Типы и форматы графиков.
	Пределы изменения переменных.

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

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

Таблица 2

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

Тип операции	Содержание	Выход
1. Ввод уравнения	Программная запись уравнения на входном MAPLE-языке.	Уравнение на входном MAPLE-языке.
2. Ввод дополнительных данных	Программная запись НУ и ГУ.	НУ и ГУ на входном MAPLE-языке.
3. Использование средств исследования уравнения суммы или произведения функций.	Установление порядка ДУ.	Вывод ответов программой.
	Исследование возможности разделения переменных.
	Определение условий поиска решения в виде.
4. Использование средств преобразования уравнения.	Выполнение замены переменных.	Вывод преобразованного уравнения.
	Выполнение подстановок.
Тип операции	Содержание	Выход
5. Использование основных инструментов решения уравнения	Получение разделенных уравнений по умолчанию с применением команды «pdsolve».	Вывод разделенных уравнений.
	Получение разделенных уравнений в заданном виде с применением операторов «pdsolve» и «hint».
	Получение решения с применением команды «build» (для тех случаев, когда это возможно).	Вывод решения уравнения.
6. Использование дополнительных инструментов решения уравнения	Учет НУ и ГУ при решении уравнений с применением команды «conds» (для тех случаев,когда это возможно).	Вывод решения уравнений с (частичным) учетом НУ и ГУ.
	Проверка полученного решения с применением команды «pdetest».	Вывод результатов проверки.
7. Решение разделенных уравнений и учет НУ и ГУ на уровне разделенных уравнений	Решение задач на собственные значения и собственные функции.	Вывод решений разделенных уравнений в общем виде.
	Определение собственных значений и собственных функций.	Вывод собственных функций
	Определение коэффициентов разложения.
8. Построение частного решения	Получение частного решения исходного уравнения с учетом исходной факторизации при разделении переменных и коэффициентов разложения.	Вывод частного решения
9. Построение общего решения	Построение общего решения как суперпозиции частных решений.	Вывод общего решения
	Учет НУ и определение оставшихся коэффициентов разложения

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

2. Метод разделения переменных

Рассмотрим подробнее метод разделения переменных. Основными этапами построения решения этим методом являются:

1) ввод уравнения и разделение переменных;

2) решение разделенных уравнений;

3) построение общего решения;

4) учет начальных условий и определение коэффициентов разложения;

5) вывод общего решения в развернутом виде и его преобразование.

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

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

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

1. Ввод уравнения и разделение переменных

PDE:=diff(u(t,x),t)=a^2*diff(u(t,x),x,x);

struc:=pdsolve(PDE,HINT=T(t)*X(x));

2. Переобозначение постоянной и решение разделенных уравнений

_c[1]=-lambda^2: dsolve(diff(T(t),t)=-lambda^2*T(t)*a^2);

dsolve(diff(X(x),`$`(x,2))=-lambda^2*X(x));

3. Построение общего решения

u[lambda](t,x):=(C1(lambda)*sin(lambda*x)+

+C2(lambda)*cos(lambda*x))*exp(-lambda^2*a ^2*t);

u(t,x):=int(u[lambda](t,x), lambda=-infinity..infinity);

4. Учет начальных условий и определение коэффициентов разложения

u_0(t,x):=eval(subs(t=0, u(t,x)))=f(x);

C1(lambda):=(1/(2*Pi))*int(f(xi)*sin(lambda*xi),xi=-infinity..infinity);

C2(lambda):=(1/(2*Pi))*int(f(xi)*cos(lambda*xi), xi=-infinity..infinity);

5. Вывод общего решения в развернутом виде и его преобразование

u(t,x):=combine(int((C1(lambda)*sin(lambda*x)+

+C2(lambda)*cos(lambda*x))*exp(-lambda^2*a^2*t),

lambda=-infinity..infinity));

u(t,x):=(1/(2*a*sqrt(Pi*t)))* int(f(xi)*exp (-1/4*(x-xi)^2/ a^2/t),xi=

=-infinity..infinity);

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

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

1. Ввод уравнения и разделение переменных

PDE:=diff(u(t,r),t)=a^2*(diff(u(t,r),r,r)+(1/r)*diff(u(t,r),r));

struc:=pdsolve(PDE,HINT=T(t)*R(r));

2. Переобозначения постоянной и решение разделенных уравнений

_c[1]=-lambda^2*a^2: dsolve(diff(T(t),t)=-lambda^2*a^2*T(t));

dsolve(diff(R(r),`$`(r,2))=-lambda^2*R(r)-diff(R(r), r)/r);

3. Учет условия регулярности решения в начале координат

BesselJ(0,lambda*r)=series(BesselJ(0,lambda*r),r=0,4):

BesselY(0,lambda*r)=series(BesselY(0,lambda*r),r=0,4):_C2=0;

4. Учет граничного условия для решения на краю области: r = r0

R[n](r):=BesselJ(0,lambda[n] *r);

BesselJ(0,mu[n])=0;

mu:=BesselJZeros:mu(0,n);

lambda[n]:=mu(0,n)/r0;

5. Вывод решений радиального и временного уравнений

R[n](r):=BesselJ(0,r*lambda [n]);

T[n](t):= C[n]*exp(-lambda[n]^2*a^2*t);

6. Построение общего решения

u[n](t,r):=T[n](t)*R[n](r): u(t,r):=Sum(u[n](t,r), n=1..infinity);

7. Замена переменной

u(t,rho):=subs(r=rho*r0,u(t,r));

8. Учет начальных условий

simplify(subs(t=0,u(t,rho))= F(rho*r0));

9. Определение коэффициентов

C[n]:=2/BesselJ(1,BesselJZeros(0,n))^2*

*int(rho*BesselJ(0,BesselJZeros(0,n)*rho)*F(rho*r0),rho = 0 .. 1);

10. Вывод общего решения

(t,rho):=Sum(C[n]*exp(-BesselJZeros(0,n)^2/r0^2*a^2*t)*

*BesselJ(0,BesselJZeros(0,n)*rho),n=1..infinity);

3. Метод функций Грина и другие методы

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

Основными этапами построения решения этим методом являются:

1) ввод неоднородного уравнения;

2) ввод представления для решения уравнения в виде ряда Фурье;

3) разложение функций в ряд Фурье;

4) определение коэффициентов разложения;

5) подстановка разложений функций в исходное уравнение;

6) представление решения в виде суммы решений однородного и неоднородного уравнений;

7) учет НУ, определение коэффициентов и вывод решения однородного уравнения;

8) построение функции Грина;

9) вывод решения однородного уравнения и частного решения неоднородного уравнения с помощью функции Грина;

10) вывод решения уравнения.

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

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

1. Ввод неоднородного уравнения

PDE:=diff(u(t,x),t)=a^2*diff(u(t,x),x,x)+w(t,x);

2. Ввод представления для решения уравнения в виде ряда Фурье

u(t,x):=Sum(u[n](t)*sin(Pi*n* x/L),n=1..infinity);

3. Разложение функций в ряд Фурье

w(t,x):=Sum(w[n](t)*sin(Pi*n*x/ L),n=1..infinity);

F(x):=Sum(F[n]*sin(Pi*n*x/L),n=1..infinity);

4. Определение коэффициентов разложения

w[n](t)=(2/L)*int(w(t,xi)*sin(Pi*n*xi/l), xi=0..L);

F[n]=(2/L)*int(F(xi)*sin(Pi*n* xi/L),xi=0..L);

5. Подстановка разложений функций u(t,x) и w(t,x) в исходное уравнение PDE;

6. Представление решения в виде суммы решений однородного и неоднородного уравнений

u[n](t)=u_Un[n](t)+u_Nu[n](t): u_Un[n](t):=_C1*exp(-a^2*Pi^2*n^2/L^2*t): u_Nu[n](t):=(Int(w[n](tau)*exp(a^2*Pi^2*n^2/L^2*(tau-t)),tau)):

7. Учет начальных условий, определение коэффициентов и вывод решения однородного уравнения

u_0:=subs(t=0,u(t,x))=F(x): u[n](0)=F[n];

eval(subs(t=0,u_Un[n](t)))= F[n];

8. Построение функции Грина

G(x,xi,t,tau):=Sum(2/L*exp(-a^2*Pi^2*n^2/L^2*(ttau))*

*sin(Pi*n*xi/L)*sin(Pi*n*x/L),n=1..infinity);

9. Вывод решения однородного уравнения и частного решения неоднородного уравнения с помощью функции Грина

u_Un(t,x):=Sum(u_Un[n](t)*sin(Pi*n*x/L),n=1..infinity);

u_Nu(t,x):=int(int(G(x,xi,t,tau)*w(tau,xi),xi=0..L),tau=0..t);

10. Вывод решения исходного неоднородного уравнения

u(t,x):=u_Un(t,x)+u_Nu(t,x);

ЗАКЛЮЧЕНИЕ

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

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

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

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1. Араманович И.Г., Левин В.И.. Уравнения математической физики. - М.: Наука, 1964.

2. Голоскоков Д.П. Уравнения математической физики. Решение задач в системе Maple. -С-Пб: Питер, 2004.

3. Сдвижников О.А., Математика на компьютере: Maple8. М.:Солон-Пресс, 2003. -176 с.

4. Тихоненко А.В. Компьютерные математические пакеты в курсе «Линейные и нелинейные уравнения физики». Часть 2. Параболические уравнения в MAPLE. - Обнинск: ИАТЭ, 2005.- 80 с.

ПРИЛОЖЕНИЕ А

Листинг программы на Maple 13

> restart;

Решить неоднородное уравнение

с неоднородностью

> w(tau,xi):=mu*(x-L/2)*sin(x/7)*t*exp(-alpha*t);

и однородными начальными условиями.

Функция Грина (функция источника):

> G(x,xi,t,tau):=sum(2/L*exp(-a^2*Pi^2*n^2/L^2*(t-tau))*sin(Pi*n*xi/L)*sin(Pi*n*x/L),n = 1 .. infinity);

Решение уравнения:

> u(t,x):=simplify(sum(int(2/L*exp(-a^2*Pi^2*n^2/L^2*(t-tau))*int(mu*(x-L/2)*sin(x/7)*t*exp(-alpha*t)*sin(Pi*n*xi/L),xi = 0 .. L),tau = 0 .. t)*sin(Pi*n*x/L),n = 1 .. infinity)) assuming n::integer;



> a:=0.1;L:=100;mu:=1; alpha:=0.25;

u(t,x):=sum(-mu*sin(1/7*x)*t*L^2*(-2*x+2*x*(-1)^n+L-L*(-1)^n)*(-1+exp(a^2*Pi^2*n^2/L^2*t))/a^2/Pi^3/n^3*exp(-t*(a^2*Pi^2*n^2+alpha*L^2)/L^2)*sin(Pi*n*x/L),n = 1 .. 300):

w(t,x):=mu*(x-L/2)*sin(x/7)*t*exp(-alpha*t);

Представим полученные решения в виде двумерных анимированных графиков:

> with(plots):

animate(plot,[w(t,x),x=0..L, color=blue],t=0..40,frames=20,thickness=3);

animate(plot,[u(t,x),x=0..L],t=0..40,frames=20,thickness=3);

Warning, the name changecoords has been redefined



Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/



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

> restart;

Решить неоднородное уравнение

с неоднородностью

> w(t,x) :=x->piecewise([0, x < 995/2],[1, x < 1005/2],[0, 1000 < x])*t*exp(-.25*t);

и однородными начальными условиями.

Функция Грина (функция источника):

> G(x,xi,t,tau):=sum(2/L*exp(-a^2*Pi^2*n^2/L^2*(t-tau))*sin(Pi*n*xi/L)*sin(Pi*n*x/L),n = 1 .. infinity);

Решение уравнения:

> u(t,x):=simplify(sum(int(2/L*exp(-a^2*Pi^2*n^2/L^2*(t-tau))*int(mu*exp(-alpha*t)*sin(Pi*n*xi/L),xi = l1..l2),tau = 0 .. t)*sin(Pi*n*x/L),n = 1 .. infinity)) assuming n::integer;

> a:=0.1;L:=1000;mu:=1;l1:=L/2-L/10;l2:=L/2+L/10;alpha:=0.25;

u(t,x):=sum(2*L^2*mu*(cos(Pi*n*l1/L)-cos(Pi*n*l2/L))*(-1+exp(a^2*Pi^2*n^2/L^2*t))*exp(-t*(a^2*Pi^2*n^2+alpha*L^2)/L^2)*sin(Pi*n*x/L)/a^2/Pi^3/n^3,n = 1 .. 1000):

w(x):=piecewise(x<l1,0, x<l2,mu, x>L,0);

w(t,x):=w(x)*t*exp(-alpha*t);

Представим полученные решения в виде двумерных анимированных графиков:

> with(plots):

animate(plot,[w(t,x),x=0..L, color=blue],t=0..40,frames=20,thickness=3);

animate(plot,[u(t,x),x=0..L],t=0..40,frames=20,thickness=3);

Warning, the name changecoords has been redefined



Размещено на Allbest.ru