Рассмотрим одномерную задачу теплопроводности на отрезке от α до β. Уравнение теплопроводности (1.8∂t∂T(M,t)=a2⋅ΔT(M,t)+f(M,t),) в этом случае принимает вид
∂t∂T(x,t)=a2⋅∂x2∂2T(x,t)+f(x,t),
(2.23)
Краевые условия Дирихле (1.13T(M,t)=Φ(M,t),M∈S.) в одномерном случае имеют вид
⎩⎨⎧T(x,0)=T0(x),T(α,t)=Φα(t),T(β,t)=Φβ(t)
(2.24)
Для избавления от неоднородности введем T(x,t)=T(x,t)+U(x,t), выберем функцию U(x,t)=a1(t)+a2(t)⋅x, и вычислим коэффициенты, в соответствии с (2.2T(M,t)+U(M,t)=Φ(M,t),T(M,t)=0⇒U(M,t)=Φ(M,t),M∈S.)
Для избавления от неоднородности введем T(x,t)=T(x,t)+U(x,t), выберем функцию U(x,t)=a1(t)⋅x+a2(t)⋅x2, и вычислим коэффициенты, в соответствии с (2.3∂n∂T(M,t)+∂n∂U(M,t)=Φ(M,t),∂n∂T(M,t)=0⇒∂n∂U(M,t)=Φ(M,t),M∈S.)
Тогда уравнение (2.23∂t∂T(x,t)=a2⋅∂x2∂2T(x,t)+f(x,t),) и краевые условия (2.28⎩⎨⎧T(x,0)=T0(x),∂x∂T(x,t)x=α=Φα(t),∂x∂T(x,t)x=β=Φβ(t)) примут вид
Понятно, что из первого уравнения c1=0 и решение имеет смысл только тогда, когда c2=0, что возможно только тогда, когда sin[γ⋅(β−α)]=0
γ⋅(β−α)=π⋅n,n∈(1..∞).
(2.35)
Нулевое и отрицательные значения n не годятся: при n=0 имеем γ=0, и Ψ(x)=c2⋅sin[0⋅(x−α)]≡0 — тождественный нуль, не являющийся собственной функцией; отрицательным же n отвечает то же собственное значение γ2 (ведь γ−n=−γn), а sin[−γn⋅(x−α)]=−sin[γn⋅(x−α)] линейно зависимо с решением для положительного n и нового вклада не даёт. Поэтому n∈(1..∞).
Таким образом, собственные значения и собственные функции имеют вид
где c2n — произвольный постоянный множитель: собственная функция определена с точностью до него. Положим c2n=1 — этот множитель всё равно вбирается в коэффициенты разложения T0n и fn, как и в общем решении.
Для разложения функций в ряд Фурье по Ψn(x) необходимо вычислить норму ∥Ψn(x)∥2, вес для декартовых координат ρ(x)=1.
Теперь подставим решение в граничные условия (2.31⎩⎨⎧∂t∂T(x,t)=a2⋅∂x2∂2T(x,t)+f(x,t),T(x,0)=T0(x)−a1(0)⋅x−a2(0)⋅x2,∂x∂T(x,t)x=α=0,∂x∂T(x,t)x=β=0,) и получим
Понятно, что из первого уравнения c2=0 и решение имеет смысл только тогда, когда sin[γn⋅(β−α)]=0
γn⋅(β−α)=π⋅n,n∈(0..∞).
(2.39)
Здесь, в отличие от задачи Дирихле, мода n=0 сохраняется: при γ0=0 собственная функция Ψ0(x)=c10⋅cos[0⋅(x−α)]=c10 — ненулевая постоянная, то есть полноценная собственная функция (постоянная мода). Поэтому индексация начинается с нуля.
Таким образом, собственные значения и собственные функции имеют вид
где c1n — произвольный постоянный множитель: собственная функция определена с точностью до него. Положим c1n=1 — этот множитель всё равно вбирается в коэффициенты разложения T0n и fn, как и в общем решении.
Норму собственных функций вычислим так же, как для случая с синусом, понизив степень косинуса по формуле половинного угла. Для n≥1
поскольку интеграл от косинуса по целому числу полупериодов равен нулю: sin[2⋅γn⋅(β−α)]=sin(2⋅π⋅n)=0. Особый случай n=0 (при γ0=0 собственная функция Ψ0(x)=1) вычислим отдельно
∥Ψ0(x)∥2=∫αβΨ0(x)2dx=∫αβ1dx=β−α.
(2.42)
Финальное решение одномерной краевой задачи Дирихле с неоднородными граничными условиями примет вид
В задаче Неймана, в отличие от Дирихле, спектр содержит нулевую моду (γ0=0, Ψ0(x)=1) — случай γn=0, который в общем решении (2.18T(M)=n=1∑∞∥Ψn(M)∥2Ψn(M)⋅a2⋅γn21⋅∭Gρ(M)⋅f(M)⋅Ψn(M)dG,γn=0.) был оставлен для отдельного рассмотрения; именно она даёт первое (незатухающее) слагаемое. Финальное решение одномерной краевой задачи Неймана с неоднородными граничными условиями примет вид
Чтобы получить стационарные решения, устремим в найденных решениях время к бесконечности (t→∞), как это сделано для общего случая в (2.18T(M)=n=1∑∞∥Ψn(M)∥2Ψn(M)⋅a2⋅γn21⋅∭Gρ(M)⋅f(M)⋅Ψn(M)dG,γn=0.). Слагаемое с начальным условием занулится, поскольку e−a2⋅γn2⋅t→0; источники и граничные условия перестают зависеть от времени, а интеграл по времени для γn=0 даёт множитель a2⋅γn21:
Стационарное решение одномерной краевой задачи Неймана не такое простое. В отличие от задачи Дирихле, чисто неймановская стационарная задача разрешима не всегда. Приведённая функция T удовлетворяет однородным условиям Неймана T′(α)=T′(β)=0 и стационарному уравнению a2⋅T′′(x)+f(x)=0. Проинтегрируем его по отрезку [α,β]:
Граничный член обратился в нуль в силу однородных условий Неймана, и осталось условие разрешимости (2.22∭Gf(M)dM=−a2⋅∬SΦ(M)dS,): ∫αβf(ω)dω=0 — суммарный приведённый источник должен обращаться в нуль (иначе тепло накапливается, среднее неограниченно растёт и стационара нет). При выполнении этого условия решение определяется лишь с точностью до произвольной аддитивной постоянной C: ей отвечает нулевая мода (γ0=0, Ψ0(x)=1), амплитуда которой стационарным уравнением не фиксируется. С этой оговоркой решение примет вид