微分方程這樣設置邊界條件是否合理?

我想要解這樣一個偏微分方程:

[frac{partial y}{partial t}=afrac{partial ^2y}{partial x^2}-bfrac{partial y}{partial x}-c y,(a,b,cin mathbb{R})]

邊界條件是:

egin{equation}egin{cases}
y|_{x=0}=y_0+y_1 cos (	ext{$omega $t})\
y|_{x
ightarrow +infty }=y_2\
end{cases}
end{equation}

其中第二個邊界條件中:y_2
eq 0,且y_2是常數

用疊加原理把這個問題分解為兩個子問題:

egin{equation}egin{cases}
frac{partial y}{partial t}=afrac{partial ^2y}{partial x^2}-bfrac{partial y}{partial x}-c y\
y|_{x=0}=y_0\
y|_{x
ightarrow +infty }=y_2\
end{cases}
	ag{1}
end{equation}

egin{equation}egin{cases}
frac{partial y}{partial t}=afrac{partial ^2y}{partial x^2}-bfrac{partial y}{partial x}-c y\
y|_{x=0}=y_1 cos (	ext{$omega $t})\
y|_{x
ightarrow +infty }=0\
end{cases}
	ag{2}
end{equation}

那麼,現在我有四個比較困惑的地方:

1. 這樣用疊加原理分解問題的方法對嗎?如果不對的話,那麼我該怎麼做?

2. 現實的物理模型是,只有周期邊界條件y_1 cos (	ext{$omega $t})使得系統是非穩態的。如果沒有這個邊界條件,那麼系統就是穩態的。在這種條件下,我是否可以把子問題(1)改寫為:

egin{equation}egin{cases}
0=afrac{partial ^2y}{partial x^2}-bfrac{partial y}{partial x}-c y\
y|_{x=0}=y_0\
y|_{x
ightarrow +infty }=y_2\
end{cases}
	ag{3}
end{equation}

呢?如果不能的話,為什麼不能呢?

3. 假如我現在要解(3)的方程。這個方程實際上是一個常微分方程,得到通解:

egin{equation}y=C_1 exp left(frac{b+sqrt{4 a c+b^2}}{2 b}
ight)+C_2 exp left(frac{b-sqrt{4 a c+b^2}}{2 b}
ight)	ag{4}end{equation}

那麼我沒辦法處理方程的第二個邊界條件了,因為(4)在正無窮遠處,第一項趨於正無窮,第二項趨於0,無論我怎麼設置C_1C_2,邊界條件都沒辦法得到滿足。這是不是意味著,這樣的微分方程,我設置的第二個邊界條件是不合理的呢?是不是這種微分方程,我在正無窮遠點設置邊界條件,只能是y|_{x
ightarrow +infty }=0或者left.frac{partial y}{partial x}
ight|_{x
ightarrow +infty }=0呢?

4. 如果我的方程的邊界條件設置沒有錯誤,那麼糾結我的思路哪裡出了問題呢?請大家指教,謝謝!


這。。疊加原理還是有個大問題的,直接看到最後就好了。

我先給題主做個示範,怎麼通過疊加原理來特定邊界條件下的解偏微分方程。

首先,對於通解來說,邊界條件一定是齊次的。只有在齊次邊界條件下,通解中的任意係數的變化均不在邊界上產生影響,才可能進一步通過初始條件定下係數。

剩下的就是特解了。特解就需要使得通解的邊界條件齊次化。

我想不清楚應該如何措辭了,我上面說的通解是對應的齊次(邊界條件的)方程的通解,包含若干待定係數,包含0解。實際的包含待定係數的解(這個才是通常說的通解)由特解和之前說的通解疊加而成。

對於樓主的問題來說,偏微分方程(話說為啥b,c前面的係數是負的,形式一點都不好看,哼唧!)

[frac{partial y}{partial t}=afrac{partial ^2y}{partial x^2}-bfrac{partial y}{partial x}-c y,(a,b,cin mathbb{R})] (1)

的解y(x,t)可以通過疊加原理得到其通解,

y(x,t)=y_0(x,t)+z(x,t) (2)

我索性換一個符號表示,這樣就沒有歧義了。

其中

egin{equation}egin{cases}
frac{partial z}{partial t}=afrac{partial ^2z}{partial x^2}-bfrac{partial z}{partial x}-c z\
z|_{x=0}=0\
z|_{x
ightarrow +infty }=0\
end{cases}
end{equation}
(3)

egin{equation}egin{cases}
frac{partial}{partial t}y_0(x,t)=afrac{partial ^2}{partial x^2}y_0(x,t)-bfrac{partial}{partial x}y_0(x,t)-c y_0(x,t)\
y_0(x,t)|_{x=0}=y_0+y_1cos(omega t)\
y_0(x,t)|_{x
ightarrow +infty }=y_2\
end{cases}
end{equation}
(4)

這樣子,疊加原理就沒有問題了,方程(3)正常地分離變數就好了。分離變數的解是適用於非穩定情形的,因為分離變數的解是完備的,即便是非穩定解也能疊加出來。

方程(3)分離變數後的通解滿足z(x,t)=Z_t(t)Z_x(x)

frac{frac{partial}{partial x}Z_t(t)}{Z_t(t)}=frac{afrac{partial^2}{partial x^2}Z_x(x)-bfrac{partial}{partial x}Z_x(x)-cZ_x(x)}{Z_x(x)}=-kappa (5)

通解大概長成這個樣子

z(x,t)=e^{-kappa t}(C_1Z_{x1}(x)+C_2Z_{x2}(x)) (6)

其中有倆待定係數,需要通過初始條件來定,Z_{x1}(x)Z_{x2}(x)是適用(3)中的邊界條件得到的兩個線性獨立的解。

方程(4)中,需要找特解。這個特解也可以分成兩個部分,一個部分是穩定的,另一個部分是不穩定的,也就是題主想要找的原方程(1)(2)的兩個特解的和。

我在這裡就直接給出特解的形式好了,可假定特解的形式為

y_0(x,t)=alpha x+eta +e^{-gamma x}cos(omega t+kx) (7)

前兩項線性項是湊原方程(1)的邊界條件,第三項衰減的振動項是湊原方程(2)的邊界條件。

之後還有計算我就不算了。。

回到題主的問題。

問題(1),將方程拆分成原方程(1)(2)是沒有問題的,但是這兩個方程都具有非齊次邊界條件,仍然需要先做特解再給出帶有待定係數的通解,其中每個方程的特解可參考方程(7),其待定係數部分可參考方程(6)。此外,在疊加原方程(1)(2)的通解時,需要重新合併待定係數,因為這兩部分是線性相關的。

另外,通過疊加原理可以使我們只用計算一個好算的邊界條件下的偏微分方程,並不表示我們只用算一個好算的方程。題主的(2)(3)問題中,說好只改邊界條件的,怎麼把方程給改了呢。。

所以問題(4)實際上就是題主沒有意識到,疊加原理不改變方程的形式。反過來,恰恰是方程的形式才決定了我們可以使用疊加原理來解方程。


既然題主說清楚了,那我也說一點自己的想法吧。

如題主所言,這樣子拆應當是不對的,因為拆了之後的一個方程無解,理由也如題主所說。事實上,如果邊界條件不含時的話,這種拋物方程的解應當以橢圓常微分方程的解為極限,那麼這樣的話確實是無解的。

至於怎麼解,沒什麼想法。。無窮遠處非零導致固定t時的解不能做Fourier變換也不能算能量。。

摺疊我吧。


推薦閱讀:

是否存在一個常微分方程存在整體解的充要條件?
積分不是初等函數和微分方程無初等解法的證明是否關鍵在於群論?
導熱微分方程的推導中為什麼要用到泰勒展開,思路是什麼?
動力系統主要做的是哪幾個方向?
一階線性齊次微分方程中的齊次的含義?

TAG:數學 | 物理學 | 數學物理方程 | 微分方程 |