Неоднородное уравнение с однородными краевыми и начальными условиями

Метод Фурье для уравнения теплопроводности

Содержание:

По этой ссылке вы найдёте полный курс лекций по математике:

Займемся решением первой смешанной задачи для уравнения теплопроводности: найти решение и(х, t) уравнения удовлетворяющее начальному условию и граничным условиям Начнем с простейшей задачи: найти решение u(x,t) однородного уравнения удовлетворяющее начальному условию и нулевым (однородным) граничным условиям Метод Фурье для уравнения теплопроводности.

Будем искать нетривиальные решения уравнения (4), удовлетворяющие граничным условиям (6), в виде Псдстаапя в форме (7) в уравнение (4), получим или откуда имеем два обыжювенных дифференциальных уравнения Чтобы получить нетривиальные решения и(х, *) вида (7), удовлетворяющие граничным условиям (6), необходимо найти нетривиальные решения уравнения (10), удовлетворяющие граничным условиям.

Таким образом, для определения фунмдои Х(х) мы приходим к задаче на собственные значения: найти те значения параметра А, при которых существуют нетривиальные решения задачи Эта задача была рассмотрена в предыдущей главе. Там было показано, что только при существуют нетривиальные решения При А = А„ общее решение уравнения (9) имеет вид удовлетворяют уравнению (4) и граничным условиям (6). Образуем формальный ряд.

Потребовав, чтобы функция и(х> t), определяемая формулой (12), удовлетворяла начальному условию , получим Ряд (13) представляет собой разложение заданной функции в ряд Фурье по синусам в интервале (О, I). Коэффициенты а„ разложения определяются по известным формулам Метод Фурье для уравнения теплопроводности Предположим, что Тогдаряд (13) с коэффициентами, определяемыми по формулам (14), будет сходиться к функции абсолютно и равномерно.

Так как при то ряд при также сходится абсолютно и равномерно.

Поэтому функция и(х, t) — сумма ряда (12) — непрерывна в области и удовлетворяет начальному и граничному условиям. Остается показать, что функция и(х, t) удовлетворяет уравнению (4) в области 0. Для этого достаточно показать, что ряды, полученные из (12) почленным дифференцированием по t один раз и почленным дифференцированием по х два раза, также абсолютно и равномерно сходятся при.

Но это следует из того, что при любом t > 0 если п достаточно велико. Единственность решения задачи (4)-(6) и непрерывная зависимость решения от начальной функции были уже установлены ранее. Таким образом, для t > 0 задача (4)-(6) поставлена корректно; напротив, для отрицательных t зада ча эта некорректна. Замечание.

В отличие отдомового уравнения уравнение неомметрично огноситн о времени t: если заменить t на -t, то получаем уравнение другого вида описывает необратимые процессы: Мы можем предсказать, каким станет данное и через промежуток времени данной t, но мы не можем с уверенностью сказать, какн м было это и за время t до рассматриваемого момента. Это раолич иемежду предсказание м и предысторией типично для параболического ура внения и не имеет места, например, для волнового уравн сния; в случае последнего заглянуть в прошлое так же легко, как и в будущее.

Возможно вам будут полезны данные страницы:

Пример:

Найти распределение температуры в однородном стерве длины ж, если начальная температура стержня и на концах стержня поддерживается нулевая температура. 4 Задача сводится к решению уравнения при начальном условии и граничных условиях Применяя метод Фурье, ищем нетривиальные решения уравнения (15), удовлетворяющие граничным условиям (17), в виде Подставляя u(x,t) в форме (18) в уравнение (15) и разделяя переменные, получим откуда Собственные значения задачи . собственные функции Хп(х) = мп пх.

При А = А„ общее решение уравнения (19) имеет вид Tn(t) = апе а п\ так что Решение задачи (15)—(17) ищем в виде ряда Потребовав выполнения начального условия (16), получим откуда . Поэтому решением исходной задачи будет фунхция 2. Рассмотрим теперь следующую задачу: найти решение гх(ж, t) неоднородного уравнения _ удовДстворя ющее начальному условию и однородным граничным услови м Предположим, что функци / непрерывна, имеет непрерывную производ-ную и при всех t > 0 выполняется условие .

Решение задач:

Решение задачи (1)-(3) будем искать в виде где определим как решение задачи а функци — как решение задачи Задача (8)—(10) рассмотрена в п. 1. Будем искать решение v(x, t) задачи (5)-(7) в виде ряда по собстве нным функциям < краевой задачи . Подсгааяяя t) в виде в уравнение (5), получим Разложим функцию /ОМ) в ряд Фурье по синусам, где Сравнивая два разложения (12) и (13) функции /(х, t) в ряд Фурье, получаем ! Пользуясь начальным условием для v(x, t).

Находим, что Решения уравнений (15) при начальных условиях (16) имеют вид: Подставляя найденные выражения для Tn(t) в ряд (11), получим решение Функция будет решением исходной задачи (1)-(3). 3. Рассмотрим задачу: найти в области решение уравнения при начальном условии и неоднородных граничных условиях Непосредственно метод Фурье неприменим из-за неоднородности условий (20).

Введем новую неизвестную функцию v(x, t), положив где Тогда решение задачи (18)—(20) сведется к решению задачи (1)-(3), рассмотренной в п. 2, для функции v(x, J). Упражнения 1. Задан бесконечный однородный стержень. Покажи те, что если начальная температура то влобой момент температура стержня 2. Ко|рцы стержня длиной ж поддерживаются при температуре, равной нулю. Начальная температура определяется формулой Определите температуру стержня для любого момента времени t > 0. 3.

Концы стержня длиной I поддерживаются при температуре, равной нулю. Начальная температура стержня определяется формулой Определите температуру стержня для любого момента времени t > 0. 4. Концы стержня длиной I поддерживаются при температуре, равной нулю. Начальное распределение температуры Определите температуру стержня для любого момента времени t > 0. Ответы

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

1 Метод Фурье для неоднородного параболического уравнения с однородными краевыми условиями второго рода.

Мария Огаркова 5 лет назад Просмотров:

1 1 Метод Фурье для неоднородного параболического уравнения с однородными краевыми условиями второго рода. 699 M. Решить неоднородную начально-краевую задачу для уравнения теплопроводности с однородными краевыми условиями второго рода. u t a u xx = fx, t, x,, t >, 1.1 u x, t = u x, t =, t >, 1. ux, = ϕx, x [, ]. 1.3 Шаг 1. Решение задачи Штурма Лиувилля. Рассмотрим задачу X x + λxx =, 1.4 X = X =. 1.5 Задача есть задача Штурма Лиувилля. Общее решение уравнения 1.4 имеет вид Xx = c 1 sin λ x + c cos λ x при λ > ; 1.6 Xx = c 1 e λ x + c e λ x при λ имеем из краевого условия X =, что c 1 =, Xx = c cos λ x X x = c λ sin λ x. Поэтому из второго краевого условия X = получаем, что λ = πk откуда имеем бесконечное множество собственных чисел задачи Штурма Лиувилля: πn λ n =, n N. 1.9 Им соответствует бесконечное множество собственных функций: X n x = cos, n N. 1.1 При λ 2 Шаг. Будем искать решение уравнения u t a u xx = fx, t с краевыми условиями u x, t = u x, t = в виде ux, t = X n T n t, где функции X n x имееют вид: n= X x 1, X n x = cos Заметим сразу, что каждое слагаемое приведённого ряда удовлетворяет краевым условиям 1., что достаточно если ряд допускает почленный переход к пределу при x +, x = для того, чтобы функция ux, t, определённая таким образом, также удовлетворяла краевым условиям 1.. Пусть функция fx, t разложена при каждом t [, T ] в ряд Фурье по косинусам fx, t = f t При этом, коэффициенты данного ряда Фурье ищутся по формулам: + cos f n t. 1.1 f n t = f, X n = fx, t cos dx Тогда уравнение 1.1 приобретает вид n= Xn xt nt a X n xt n t = f t Для его выполнения достаточно, чтобы + f n t cos. то есть T t = f t X n xt nt a X n xt n t = f n t cos T t = f t T nt + πna T n t cos = f n t cos для n = для n N, для n = для n N. Это заведомо выполнено, если T t = f t для n = 1.14 T nt + πna T n t = f n t для n N, 1.15 Итак, мы получили условия на функции T n t, достаточные для того, чтобы функция ux, t = T n t cos была если ряд «хороший» решением уравнения n= u t a u xx = fx, t с краевыми условиями u x, t = u x, t =. Шаг 3. Решаем задачу Из условий задачи мы ещё не использовали только начальные условия —

3 ux, = ϕx. Пусть функция ϕx, входящая в начальное условие, разлагается в ряд по косинусам Подставим функцию ux, t = T n t cos опять-таки в предположении, что ряд «хороший» в начальное условие: ϕx = ϕ + ϕ n cos, x [, ] где 1.16 n= n= ϕ n = ϕx cos dx T n cos = ϕ + ϕ n cos. Для выполнения этого равенства достаточно, чтобы T = ϕ для n = T n = ϕ n для n N. Таким образом, из 1.14, 1.15 и , для функций T n t имеем задачу Коши: T t = f t T = ϕ T nt + πna T n t = f n t для n = 1.18 T n = ϕ n для n N Эти задачи Коши имеют единственное решение при любых f n C[, T ] и любых значениях ϕ n R. При n = : T t = ϕ + 1 t f τdτ. 1. При n N: сначала решаем однородное уравнение: Его общее решение имеет вид: T nt + πna T n t =. T n t = c e πna t. Метод вариации постоянной позволяет нам искать решение уравнения 1.19 в виде T n t = cte πna t, = T nt = c t πna ct e πna t. Подставив эти равнества в 1.19, получим уравнение для нахождения ct: c t = f n te πna t, -3-

4 откуда, с учётом начального условия T n = ϕ n, ct = ϕ n + t f n τe πna τ dτ. 1.1 Таким образом, t T n t = ϕ n e πna t + f n τe πna t τ dτ. 1. Всё, что нам осталось сделать, это подставить 1., 1. в формулу Получаем ответ: ux, t = ϕ + 1 t f τdτ + ux, t = n= T n t cos. ϕ n e πna t + t f n τe πna t τ dτ cos 669 M. Решить неоднородную начально-краевую задачу для уравнения колебаний с однородными краевыми условиями второго рода.. u tt a u xx = fx, t, x,, t >,.1 u x, t = u x, t =, t >,. ux, = ϕx, x [, ]..3 u t x, = ψx, x [, ]..4 Шаг 1. Решение задачи Штурма Лиувилля. Этот шаг полностью повторяет Шаг 1. задачи 699 M. Шаг. Будем искать решение уравнения u tt a u xx = fx, t с краевыми условиями u x, t = u x, t = в виде ux, t = X n T n t, где функции X n x имееют вид: n= X x 1,.5 X n x = cos..6 Пусть функция fx, t разложена при каждом t [, T ] в ряд Фурье по косинусам fx, t = f t При этом, коэффициенты данного ряда Фурье ищутся по формулам: + cos f n t..7 f n t = f, X n = fx, t cos dx..8-4-

5 Тогда уравнение.1 приобретает вид n= Xn xt n t a X n xt n t = f t Для его выполнения достаточно, чтобы то есть T t = f t X n xt n t a X n xt n t = f n t cos + f n t cos. T t = f t T n t + πna T n t cos = f n t cos Это заведомо выполнено, если T t = f t для n = для n N, для n = для n N. для n =.9 T n t + πna T n t = f n t для n N,.1 Итак, мы получили условия на функции T n t, достаточные для того, чтобы функция ux, t = T n t cos была если ряд «хороший» решением уравнения n= u tt a u xx = fx, t с краевыми условиями u x, t = u x, t =. Шаг 3. Решаем задачу.1.4. Из условий задачи.1.4 мы ещё не использовали только начальные условия ux, = ϕx, u t x, = ψx. Пусть функции ϕx, ψx, входящие в начальные условия, разлагаются в ряд по косинусам Подставим функцию ux, t = T n t cos опять-таки в предположении, что ряд «хороший» в начальные условия: ϕx = ϕ + ϕ n cos, x [, ] где.11 ϕ n = ϕx cos dx..1 ψx = ψ + ψ n cos, x [, ] где.13 n= n= ψ n = ψx cos dx..14 T n cos = ϕ + ϕ n cos ; -5-

6 n= УМФ семинар Метод Фурье T n cos = ψ + ψ n cos. Для выполнения этих равенств достаточно, чтобы T = ϕ T = ψ для n = T n = ϕ n T n = ψ n для n N. Таким образом, из.9,.1 и.11.1, для функций T n t имеем задачу Коши: T t = f t T = ϕ для n =.15 T = ψ T n t + πna T n t = f n t T n = ϕ n для n N..16 T n = ψ n Эти задачи Коши имеют единственное решение при любых f n C[, T ] и любых значениях ϕ n R, ψ n R. При n = : T t = ϕ t + При n N: сначала решаем однородное уравнение: Его общее решение имеет вид: ψ + 1 τ T n t + πna T n t =. T n t = c 1 sin πnat f κdκ dτ c cos πnat. Метод вариации постоянной позволяет нам искать решение уравнения.16 в виде T n t = c 1 t sin πnat + c t cos πnat, где c 1, t есть решения системы c 1t sin πnat πna c 1 t cos πnat + c t cos πnat = ; c t sin πnat = fn t. откуда c 1t = πna f nt cos πnat, c t = πna f nt sin πnat. С учётом начальных условий T n = ϕ n, T n = ψ n окончательно получаем c 1 t = πna ψ n + t f n τ cos πnaτ dτ, c t = ϕ n t f n τ sin πnaτ dτ..18 πna πna -6-

7 Таким образом, T n t = ϕ n sin πnat + πna + ψ n πna sin πnat πnat cos + t f n τ cos πnaτ dτ cos πnat t f n τ sin πnaτ Всё, что нам осталось сделать, это подставить.17,.19 в формулу ux, t = T n t cos. n= dτ M. Классический способ. Решить неоднородную начально-краевую задачу для уравнения колебаний с однородными краевыми условиями второго рода. u tt a u xx = fx, x,, t >, 3.1 u, t = u, t =, t >, 3. ux, = β α x + α, x [, ]. 3.3 u t x, =, x [, ]. 3.4 Шаг 1. Решение задачи Штурма Лиувилля. Рассмотрим задачу 3.5 X x + λxx =, 3.6 X = X =. 3.7 Задача есть задача Штурма Лиувилля. Её решение нам уже известно: λ n =, X n x = sin, n N. Шаг. Будем искать решение уравнения u tt a u xx = fx, t с краевыми условиями u, t = u, t = в виде ux, t = X n T n t, где функции X n x имееют вид: X n x = sin, n N. 3.8 Пусть функция fx разложена в ряд Фурье по синусам так как в данном примере f не зависит от t, то f n тут просто константы, не зависящие от t fx = sin f n. 3.9 При этом, коэффициенты данного ряда Фурье ищутся по формулам: f n = f, X n = fx sin dx

8 Тогда уравнение 3.1 приобретает вид Xn xt n t a X n xt n t = Для его выполнения достаточно, чтобы X n xt n t a X n xt n t = f n sin то есть T n t + πna T n t sin = f n sin Это заведомо выполнено, если f n sin. для n N, для n N. T n t + πna T n t = f n для n N, 3.11 Итак, мы получили условия на функции T n t, достаточные для того, чтобы функция ux, t = T n t sin была если ряд «хороший» решением уравнения n= u tt a u xx = fx, t с краевыми условиями u, t = u, t =. Шаг 3. Решаем задачу Из условий задачи мы ещё не использовали только начальные условия ux, = ϕx = α β x α, u t x, = ψ =. Найдём разложение функций ϕx, ψx, входящих в начальные условия, в ряд по синусам ϕx = ψx = Подставим функцию ux, t = ϕ n sin, где ϕ n = ϕx sin ψ n sin, где ψ n = ψx sin T n t sin опять-таки в предположении, что ряд «хороший» в начальные условия: T n sin = dx. 3.1 dx T n sin = Для выполнения этих равенств достаточно, чтобы ϕ n sin ; ψ n sin. T n = ϕ n T n = ψ n для n N. Таким образом, из 3.11 и , для функций T n t имеем задачу Коши: T n t + πna T n t = f n T n = ϕ n для n N T n = ψ n -8-

9 Эти задачи Коши имеют единственное решение при любых f n R и любых значениях ϕ n R, ψ n R. Найдём ϕ n, ψ n из 3.1, 3.13 с учётом, что ϕx = β α x + α, ψ =. ϕ n = ϕx sin dx = + α β cos πn πn > << >β α x cos dx + α x= + x= sin dx = = 1n β α + α 1 1n = πn πn πn α 1n β, ψ n =. При n N: сначала решаем однородное уравнение: Его общее решение имеет вид: T n t + πna T n t =. T n t = c 1 sin πnat + c cos πnat. Простой вид правой части позволяет нам угадать частное решение неоднородного уравнения 3.14 в виде константы: f n. Поэтому общее решение 3.14 имеет вид: πna T n t = c 1 sin πnat + c cos πnat Из начального условия T n = ψ n = следует, что c 1 =. + f n πna. А второе начальное условие T n = ϕ n = πn α 1n β даёт нам c = πn α 1n β f n πna. Таким образом, T n t = πn α 1n β f n cos πnat + f n πna πna

10 Всё, что нам осталось сделать, это подставить 3.15 в формулу Получим ответ: ux, t = ux, t = T n t sin. πn α 1n β f n cos πnat sin + πna f n πna sin. 654 M. Короткий способ. Решить неоднородную начально-краевую задачу для уравнения колебаний с однородными краевыми условиями второго рода. u tt a u xx = fx, x,, t >, 4.1 u, t = u, t =, t >, 4. ux, = β α x + α, x [, ]. 4.3 u t x, =, x [, ]. 4.4 Шаг 1. Так как правые части всех равенств в этой задаче не зависят от времени, будем искать решение задачи в виде суммы ux, t = vx, t + wx. Найдём w = wx такую, чтобы 4.5 w tt a w xx = fx, x,, t >, 4.6 w, t = w, t =, t >. 4.7 Раз w = wx, то w tt =, и задача принимает более простой вид: Проинтегрируем уравнение 4.8 один раз: Проинтегрируем второй раз: w x = fx, a x,, 4.8 w = w =. 4.9 w x = 1 a wx = 1 a x y x fsds + c 1. fsdsdy + c 1 x + c. Из краевого условия w = получаем, что c =, а из w =, что = 1 a y fsdsdy + c 1, -1-

11 откуда c 1 = 1 a Итак, функция wx нам полностью известна: wx = x a y y fsdsdy. fsdsdy 1 a Тогда для vx, t = ux, t wx получается задача x y fsdsdy. 4.1 v tt a v xx =, x,, t >, 4.11 v, t = v, t =, t >, 4.1 vx, = β α x + α wx = ϕx, x [, ] v t x, =, x [, ] Такую задачу мы уже умеем решать см. номер 643. Её ответ: vx, t = sin A n cos где A n и B n задаются равенствами πna t + B n sin πna 4.15 t, 4.16 A n = ϕx sin dx; 4.17 B n = ψx sin aπn В нашем случае ψ =, а ϕ = β α x + α wx, откуда A n = β α x + α wx sin dx dx, B n = и функция v имеет вид vx, t = A n sin cos πnat Всё, что нам осталось сделать, это подставить в формулу найденные функции v и wиз 4.19 и 4.1. ux, t = vx, t + wx -11-

12 667. Решить неоднородную начально-краевую задачу для уравнения колебаний с однородными краевыми условиями второго рода. u tt a u xx = fx, t, x,, t >, 5.1 u, t = u x, t =, t >, 5. ux, =, x [, ]. 5.3 u t x, =, x [, ]. 5.4 Шаг 1. Решение задачи Штурма Лиувилля. Этот шаг мы проходили, когда решали задачу 649 M. Результат: бесконечное множество нетривиальных решений πn 1 πn 1x λ n =, X n x = sin, n N. Шаг. Будем искать решение уравнения u tt a u xx = fx, t с краевыми условиями u x, t = u x, t = в виде ux, t = X n T n t, где функции X n x имееют вид: n= 5.5 πn 1x X n x = sin. 5.6 Пусть функция fx, t разложена при каждом t [, T ] в ряд Фурье по функциям X n x: πn 1x fx, t = sin f n t. 5.7 При этом, коэффициенты данного ряда Фурье ищутся по формулам как мы убедились, решая 649 M : f n t = f, X n = πn 1x fx, t sin dx. 5.8 Тогда уравнение 5.1 приобретает вид Xn xt n t a X n xt n t = n= Для его выполнения достаточно, чтобы πn 1x X n xt n t a X n xt n t = f n t sin πn 1x f n t sin. то есть T n t + π n 1 a πn 1x πn 1x T 4 n t sin = f n t sin Это заведомо выполнено, если для n N, для n N. T n t + π n 1 a 4 T n t = f n t для n N,

13 Итак, мы получили условия на функции T n t, достаточные для того, чтобы функция ux, t = T n t cos была если ряд «хороший» решением уравнения n= u tt a u xx = fx, t с краевыми условиями u x, t = u x, t =. Шаг 3. Решаем задачу Из условий задачи мы ещё не использовали только начальные условия ux, =, u t x, =. Функции ϕx, ψx, входящие в начальные условия, разлагаются в ряд по функциям X n x πn 1x ϕx = ϕ n sin ψx = ψ + ψ n cos Подставим функцию ux, t = T n t sin «хороший» в начальные условия: πn 1x πn 1x T n sin = πn 1x T n sin =, x [, ] где ϕ n =, 5.1, x [, ] где ψ n = 5.11 опять-таки в предположении, что ряд πn 1x ϕ n sin = ; πn 1x ψ n sin =. Для выполнения этих равенств достаточно, чтобы T n = ϕ n = T n = ψ n = для n N. Таким образом, из 5.9 и , для функций T n t имеем задачу Коши: T n t + π n 1 a T 4 n t = f n t T n = n N. 5.1 T n = Эти задачи Коши имеют единственное решение при любых f n C[, T ] и любых значениях ϕ n R, ψ n R. сначала решаем однородное уравнение: Его общее решение имеет вид: T n t = c 1 sin T n t + π n 1 a 4 T n t =. πn 1at + c cos πn 1at. Метод вариации постоянной позволяет нам искать решение уравнения 5.1 в виде T n t = c 1 t sin πn 1at + c t cos πn 1at, где c 1, t есть решения системы -13-

14 откуда c 1t sin πn 1at πn 1a + c t cos πn 1at = ; c 1t cos πn 1at c t sin πn 1at = f n t. c 1t = πn 1a f nt cos πn 1at, c t = πn 1a f nt sin πn 1at. С учётом начальных условий T n = ϕ n =, T n = ψ n = окончательно получаем c 1 t = t πn 1aτ f n τ cos dτ, πn 1a Таким образом, T n t = c t = πn 1at t πn 1aτ sin f n τ cos dτ πn 1a cos πn 1at t πn 1aτ f n τ sin dτ. πn 1a 5.13 t f n τ sin Всё, что нам осталось сделать, это подставить 5.14 в формулу I. Найти решение ux, t задачи ux, t = T n t sin πn 1x. πn 1aτ dτ u tt a u xx = fx, t, x,, t >, 5.15 u, t = u, t =, t >, 5.16 ux, = ϕx, x [, ] u t x, =, x [, ] Решение: см. 654 M Классический способ, стр II. Найти решение ux, t задачи u t a u xx = fx, x,, t >, 5. u, t = u, t =, t >, 5.1 ux, = ϕx, x [, ]. 5. Решение: см. 654 M Короткий способ, стр