
- •1. Уравнение колебания струны
- •2. Уравнение теплопроводности
- •1.Волновое уравнение
- •2Уравнение теплопроводности.
- •Уравнение Пуассон
- •В общем виде Теплопроводность
- •Вторая краевая задача
- •В общем виде Волновое уравнение
- •Первая краевая задача
- •Вторая краевая задача
- •Постановка краевых задач
- •1.Уравнения описывающие стационарные процессы
- •Вторая краевая задача
- •2.Не стационарное уравнение.
- •Первая краевая задача
- •Вторая краевая задача
- •Корректность по Адамару
- •Пространство Соболева
- •След функции
- •1. Уравнение теплопроводности
- •Теорема 3
- •Задача Коши
- •Теорема 3
1. Уравнение теплопроводности
u/t=L(u(t,x))+f(t,x) (1), L(u)=i,j=1naij(t,x)uxixj(t,x)+i=1nbi(t,x)uxi+c(t,x)u
i,j=1naij(t,x)ij>0, (t,x)П[0,T]={(t,x)|0tT,xEn}, f,aij,bi,c-заданы в П[0,T], u(0,x)=(x),xEn (2).Задача Коши (1),(2): Найти функцию uC1,2(П[0,T])C(П[0,T]) удовлетворяющую (1) в П[0,T] и совпадающую с (х) при t=0(выполнение условия (2)). П[0,T]=(0,Т]En, C1,2(G) линейное пространство функций на G имеющие непрерывные производные ut,uxi,u,uxixj.
2.Уравнение колебания
2u/t2=L(u(t,x))+f(t,x) (3),u(0,x)=u0(x),xEn(4), ut(0,x)=u1(x),xEn(5), (3)-(5) задача Коши для (3): Найти u(t,x)C2(П(0,T]C1(П[0,T]) удовлетворяющуу в П(0,T] уравнению (3) и условию (4)(5).
Уравнение колебания струны. Формула Даламбера.
utt=a2uxx (1), u(0,x)=u0(x) (2), ut(0,x)=u1(x), xE1(3). Введем следующую замену переменных: =x+at,=x-at подставляя в (1) получим u(,)=0 отсюда получим решение u()=g()+f() произвольные дважды дифференцируемые, подставляя из замены получаем u(t,x)=v(x+at,x-at)=f(x+at)+g(x-at) (4)-есть решение (1). Нужно найти f,g на основании (2),(3). При t=0 u(0,x)=f(x)+g(x)=u0(x) (5), ut(0,x)=af`(x)-ag`(x)=u1(x) (6).f,g- функции одного аргумента. Проинтегрируем (6) по х, разделив его на а. f(x)-g(x)=1/a0xu1()d+c (7) сложим (5) и (7) и разделим на 2. f(x)=u0(x)/2+1/2a0xu1()d+c/2 (8). Из (5) вычтем (7) находим g(x)= u0(x)/2-1/2a0xu1()d-c/2, u(t,x)=f(x+at)+g(x-at)= u0(x+at)/2+1/2a0x+atu1()d+c/2+ u0(x-at)/2-1/2a0x-atu1()d-c/2, u(t,x)= u0(x+at)/2+ u0(x-at)/2+1/2ax-atx+atu1()d..Функция v(t,x)=(x+-at) где функция одного аргумента, а константа называется простой волной. Покажем что решение задачи (1)-(3) есть сумма двух простых волн u(t,x)=(x+at)+(x-at) где (x+at)=u0(х+at)/2+1/2a 0x+atu1()d,(x-at)=u0(х-at)/2+1/2a x-at0u1()d,x+at=C1,x-at=C2, (x+-at)=(const), (x+at),(x-at)-простые волны. Пусть u1=0u(t,x)=u0(x+at)/2+u0(x-at)/2 ,X1=x0-at0, x2=x0+at0, u(t0,x0)=[u0(x1)+ u0(x2)]/2, u(t0,x0)=[u0(x1)+u0(x2)]/2+1/2a[x1,x2]u1()d. Пусть u1(x)0 тогда u(t,x)=u0(x+at)/2+u0(x-at)/2, u(0,x)=u0, u0=u0(x)/2+u0(x)/2.
Задача Коши на полупрямой utt=a2uxx (1) u(0,x)=u0(x) (2) ut(0,x)=u1(x) (3) u(t,0)=0 (4),t>0,x>0,ux(t,0)=0 (5),(1)-(4)-задача где на левом конце задана функция (1)-(3),(5)- задача где на левом конце задано значение производной.
Лемма1.Функции u0 и u1 нечетные. Тогда u(t,0)=0 t0 где u-решение задачи Коши (1)-(3). Доказательство U(t,0)=u0(at)/2+u0(-at)/2+1/2a-atatu1() d=0
Лемма 2 Функции u0 и u1 четные. Тогда ux(t,0)=0 t0. Доказательство продифференцируем по х формулу Даламбера и воспользуемся четностью функции u1 и нечетностью производной от u0.
U0(x)=(система) u0(x) x0, -u0(-x) x<0,U1(x)= (система)u1(x) x0, -u1(-x) x<0, U(t,x)=U0(x+aat)/2+U0(x-at)/2+1/2a-atatU1()d на основании леммы 1 u(t,x)=U(t,x) t0,x>0 решение задачи (1)-(4),Utt=a2Uxx ,U(0,x)=U0(x) ,Ut(0,x)=U1(x)
U(t,x)=(система) u0(x+at)/2+u0(x-at)+ 1/2ax-atx+atu1()d, x/a>t,u0(x+at)/2-u0(x-at)+ 1/2a0x+atu1()d+1/2ax-at0 -u1(-)d= u0(x+at)/2-u0(x-at)+ 1/2ax-atx+atu1()d, x/a<t.
Задача Коши для уравнения теплопроводности. Формула Пуассона.
ut=uxx (1), u(0,x)=u(x) x(-,+), t>0
Условие 1. |(x)|Mea|x|, M,a- const>0.Решение задачи (1),(2) дается формулой Пуассона U(t,x)=1/2t-+exp[-(-x)2/4t]*()d (3)
Утверждение 1. При выполнении условия (1) интеграл (3) сходится t,xП(0,+) и |u(t,x)|2Mexp[a2t]*ea|x| (4). Доказательство: U(t,x)=1/2t-+exp[-(-x)2/4t]*()d ( условие 1) M/2t-+exp[-(-x)2/4t]* ea|| d(|||-x+x|) M/2-+exp[-(-x)2/2t]+a|x|+a|-x|2t/2t d(-x)/2t=[(-x)/2t=]= M*ea|x|/-+exp[-2/2t] d=(подынтегральная функция четная)= 2M*ea|x|/-+exp[-2+2at-a2t]*exp[a2t] d= 2M*ea|x|*exp[a2t]/ -+exp[-(-at)2] d<[-at=z]< 2M*ea|x|*exp[a2t]/-+exp[-z2] dz=2M*exp[a2ta|x|]
Утверждение 2. При выполнении условия 1 функция u(t,x) бесконечно дифференцируема по t,x при t>0. При этом m+nu(t,x)/tmxn=1/2-+()m+n( exp[-(-x)2/4t ]/td /tmxn (5).Интеграл в правой части (5) сходится в прямоугольнике R[t0,T,r]={t0tT,|x|<r}. Рассмотрим подинтегральное выражение (5) ()- известна оценка, оценим [ ]. m+n( exp[-(-x)2/4t ]/t/tmxn= конечная C*(-x)степень*exp[-(-x)2/4t]/tстепень, выражение можно записать как полином по переменной с коэффициентами зависящими от t,x=p(t,x,)* exp[-(-x)2/4t ], где p(t,x,)- полином степени к ( к зависит от n,m) |p(t,x,)|C(t0,T,r)(1+||k)=( в силу |ax2+bx+c|d*2(1+x2)= d*2(1+x2), |p(t,x,)|C(t0,T,r)(1+||k)( в силу ez=1+z+z2/2+……zk/k!+…, 1+||k<e*k!)N(t0,T,r)e, N(t0,T,r)=k!C(t0,T,r), |()|P(t,T,r)~Me(a+1)||,[]exp[-(-x)2/4t], ~Me(a+1)|| Заметим: aa+1;;MM*N(t0,T,r) и по утверждению 1, этот интеграл сходится, причем равномерно. По теореме о дифференцирования несобственных интегралов функция u имеет производную по m+nu/tmxn при этом выполняется соотношение (5).
Утверждение 3. Функция u(t,x) заданная соотношением (3) является решением уравнения (1) при t>0, x(-,+). Доказательство: из (5) следует ut-uxx=-+,{[exp(--x)2/4t]t-2[exp(--x)2/4t]/2tx2}()d. Легко проверить что {} 0 t>0,xE1 ut-uxx=0.
Утверждение 4. Пусть х0-точкинепрерывности функции (х). Тогда limt+0,xx0u(t,x)=(x0) (6).Доказательсво: P (t,x) P |u(t,x)-(x0)|< , (x0)=1/2t*-+exp[(--x)2/4t]*(x0)d, |u(t,x)-(x0)|=| 1/2t*-+exp[(--x)2/4t]*[()-(x0]d| 1/2t*-+exp[(--x)2/4t]*[()-(x0]d=[(--x)/2t=]=1/*-+exp[-2]*[(x+2t *)-(x0]d (*)=1/(--N ++N++-NN exp[-2]*[(x+2t *)-(x0]d. Интегралы обозначим I1,I2,I3. Зафиксируем 0 и пусть Р прямоугольник {(t,x),0<tt0, |x-x0|<} I1: подынтегральная функция удовлетворяет условию такому же как и функция . По утверждению 1 (*) сходится равномерно по (t,x) из P .Следовательно при достаточно большом N I1</3 (t,x)P при большом выборе N аналогично I2</3 Рассмотрим I3=-NN exp[-2]*[(x+2t *)-(x0]d, если t,x близки к 0 и х0 соответственно .Пусть |x+2t-x0|<(x), [-N,N], |(x+2t)-(x0)|</6N . Тогда I3</6N-NNd=/3 таким образом |u(t,x)-(x0)|< при выполнении (*). Доказали следующую теорему : При выполнении условий утверждения 1 интеграл Пуассона (3) есть решение задачи Коши (1)-(2). Начальный условия имеют смысл в точке непрерывности.
Свойство 1. Если |(x)|<M то и |u(t,x)|M T>0,xE1, |u(t,x)| 1/2t*-+exp[(--x)2/4t]d*M=1*M
Свойство 2. Пусть (х)-нечетная функция функция, тогда u(t,0)= 1/2t*-+exp[-2]()d=0 подынтегральная функция нечетная. Следовательно можем решать задачу система ut=uxx (7), u(0,x)=u0(x) (8) x>0, u(t,0)=0 (9) Решение ищется в первой четверти. Продолжим U0=(система) u0(x), x0, -u0(-x),x<0 . Построим функцию пуассона для u(t,x)= 1/2t*-+exp[(--x)2/4t]u0()d.
Принцип максимума для параболических функций.
L(u)=f (1) где L(u)=i,j=1naij(t,x)2u/xixj+i=1nbi(t,x)u/xi+c(t,x)-ut, f=f(t,x), aij,bi,c- коэффициенты оператора f. Заданы в QT{(t,x)|0<t<T,x},En. Все коэффициенты и правая часть ограничены в QT |aij(t,x)|+|bi(t,x)|+|c(t,x)|+(f(t,x)|K, 0<i,jnaji(f,x)ij (2) (t,x)QT,(1…n)En.
Уравнение параболического типа.
*Функция u(t,x)C(QT)C1,2(QT) называется классическим решением уравнения (1) если она удовлетворяет уравнению (1) в QT.
C1,2(QT)={f1,ft,fxi,fxi,xjC(QT)}
Замечание 1. Заменой можно добиться того, что бы коэффициент c(t,x) будет строго отрицательным в Q( мы предпологали что все коэффициенты в QТ ограничены) v(t,x)=u(t,x)*e-t, где -const, u=vet, et i,j=1naij(t,x)vxixj+ et i=1nbi(t,x)vxi+с(t,x) et v-vt et -v et =f,
~L(v) i,j=1naij(t,x)vxixj+ i=1nbi(t,x)vxi+(с-)v-vt=f/et .Если удовлетворяет условию с< (t,x) QT. Обозначим c-=~c(с волной). Получим уравнение ~L(v)=~F где ~F=fe-t. Если функция f имеет постоянный знак в Q то и ~F=fe-t имеет тот же знак.
Теорема 1.
Пусть функция u непрерывна в QT все ее производные входящие в оператор L непрерывны в QT\ГТ и выполняется неравенство L(u(t,x))0 в QT\ГТ (3) u(t,x)0 на ГТ (4).Пусть коэффициент c(t,x) ограничена сверху некоторой постоянной M[ c(t,x)<M (t,x) QT] Тогда u(t,x)0 в QT.
Лемма 1. Пусть выполнено условие (2) , пусть функция u(t,x)C(QT) C1,2(QT\ГТ) достигает своего минимума в точке (t0,x0) QT\ГТ . Тогда i,j=1naij(t0,x0)uxixj(t0,x0)0 (5) .Без доказательства.
Доказательство теоремы. Пусть утверждение теоремы не выполняется следовательно u(t,x) принимает отрицательное значение в QT на ГТ u(t,x)0 она принимает отрицательные значение в QT\ГТ так как функция u(t,x) непрерывна в QT то она достигает своего максимума и минимума в QT следовательно min<0 и достигается он в QT\ГТ Пусть это точки (t0,x0). Рассмотрим L(u(t0,x0))=i,j=1naij(t0,x0)uxixj (t0,x0)+i=1nbi(t0,x0)uxi(t0,x0)+с(t0,x0) u(t0,x0)-vt (t0,x0) =A1+A2+A3+A4 .По лемме А10, А2=0 uxi=0 в точке минимума. A4: Если (t0,x0) QTА4 =ut(t0,x0)=0 (точка локального минимума). В точке (t0,x0): ut(T,x0)=limtT-0 u(T,x0)- u(t,x0)/T-t0, A3=c(t0,x0)u(t0,x0) 1.Предположим что c(t,x)<0 в QT тогда c(t0,x0)<0, u(t0,x0)<0, A3>0,L(u(t0,x0))>0, (t0,x0) QT\ГТ. противоречие. Значит теорема доказана в случае когда C<0. 2.Пусть c(t,x)<M тогда замена u=veMt приводит нас к рассмотрению оператора ~L(v)0 причем v0 на ГТ, c(t,x)<0 в QT. По доказанному v(t,x)0 в QT. А Значит и u=veMt0 в QT.
Введем N,q неотрицательные константы, С0- строго положительная константа.
Теорема 2. Пусть функция u(t,x) непрерывна в QT удовлетворяет в QT\ГТ уравнению (1) и |u(t,x)|Гтq. Пусть f ограниченная функция а коэффициент c(t,x) не положителен |f(t,x)|N, c(t,x)<0 (t,x) QT. Тогда всюду в QT |u(t,x)|q+Nt (6)
Доказательство: Рассмотрим 2 формулы w+(t,x)=q+Nt+u(t,x), w-(t,x)=q+Nt-u(t,x) Нетрудно проверить выполнение условий теоремы (Пусть функция u непрерывна в QT все ее производные входящие в оператор L непрерывны в QT\ГТ и выполняется неравенство L(u(t,x))0 в QT\ГТ (3) u(t,x)0 на ГТ (4).Пусть константа c(t,x) ограничена сверху некоторой постоянной M[ c(t,x)<M (t,x) QT] Тогда u(t,x)0 в QT) для w+-, действительно w+-|Гт=q+Nt+-u(t,x)|ГтNt0 1 условие теоремы выполнено, 2 условие: L(w+-)=L(q+Nt)+-L(u)=cq+cNt-N-+f 0 по условию Th1 имеем неравенство w+-(t,x)=q+Nt+-u(t,x)0, q+Ntu(t,x)-(q+Nt) (6) t,xQT