偏微分方程







二維熱傳導方程式的解
































偏微分方程英语:partial differential equation,缩写作PDE)指含有未知函数及其偏导数的方程。描述自变量、未知函數及其偏导數之間的關係。符合這個關係的函数是方程的解。


偏微分方程分為線性偏微分方程式與非線性偏微分方程式,常常有幾個解而且涉及額外的邊界條件。




目录






  • 1 記號及例子


    • 1.1 拉普拉斯方程


    • 1.2 泊松方程


    • 1.3 波動方程式


    • 1.4 熱傳導方程式




  • 2 分类


    • 2.1 一阶偏微分方程


    • 2.2 二阶偏微分方程


    • 2.3 混合形式方程




  • 3 偏微分方程有关问题


    • 3.1 适定问题




  • 4 解析法解偏微分方程


    • 4.1 分离变量法


    • 4.2 特征线法


    • 4.3 积分变换


    • 4.4 变量变换


    • 4.5 基本解


    • 4.6 叠加原理




  • 5 数值法解偏微分方程


  • 6 参考文献





記號及例子


方程式中常以u為未知數及偏微分,如下:


ux=∂u∂x{displaystyle u_{x}={partial u over partial x}}{displaystyle u_{x}={partial u over partial x}}

uxy=∂2u∂x∂y{displaystyle u_{xy}={partial ^{2}u over partial x,partial y}}{displaystyle u_{xy}={partial ^{2}u over partial x,partial y}}

用於空間偏微分的梯度運算子=(∂x,∂y,∂z){displaystyle nabla =({partial over partial _{x}},{partial over partial _{y}},{partial over partial _{z}})}{displaystyle nabla =({partial  over partial _{x}},{partial  over partial _{y}},{partial  over partial _{z}})}


時間偏微分=∂u∂t{displaystyle {dot {u}}={partial u over partial t}}{displaystyle {dot {u}}={partial u over partial t}},線性偏微分方程式的例子如下:



拉普拉斯方程


uxx+uyy+uzz=0{displaystyle u_{xx}+u_{yy}+u_{zz}=0}u_{xx} + u_{yy} + u_{zz} = 0

適用于重力场问题的求解



泊松方程


uxx+uyy+uzz=f(x,y,z){displaystyle u_{xx}+u_{yy}+u_{zz}=f(x,y,z)}u_{xx} + u_{yy} + u_{zz} = f(x,y,z)

適用於所有物質或電荷的重力場或靜電場。



波動方程式


未知函數u(x,y,z,t):



utt=c2(uxx+uyy+uzz){displaystyle u_{tt}=c^{2}(u_{xx}+u_{yy}+u_{zz})}u_{tt} = c^2( u_{xx} + u_{yy} + u_{zz} )

=c2∇2u{displaystyle {ddot {u}}=c^{2}nabla ^{2}u}ddot u=c^2nabla^2u



熱傳導方程式


ut=k(uxx+uyy+uzz){displaystyle u_{t}=k(u_{xx}+u_{yy}+u_{zz})}u_t = k ( u_{xx} + u_{yy} + u_{zz} )

其中k代表該材料.



分类


一些线性二阶偏微分方程可以分为:抛物线方程,双曲线方程和椭圆方程。其他的像Euler–Tricomi方程在不同应用领域中也有不同的形式。这种分类便于在解偏微分方程时寻找初始条件提供依据。



一阶偏微分方程



一阶偏微分方程是只和未知數的一階導數有關的偏微分方程,表示式为:


Aux+Buy+⋯=0,{displaystyle Au_{x}+Bu_{y}+cdots =0,}{displaystyle Au_{x}+Bu_{y}+cdots =0,}

其中参数A,B是x,y的變數。



二阶偏微分方程


表示式为:


Auxx+2Buxy+Cuyy+⋯=0,{displaystyle Au_{xx}+2Bu_{xy}+Cu_{yy}+cdots =0,}Au_{xx} + 2Bu_{xy} + Cu_{yy} + cdots = 0,

其中参数A,B,C是x,y的變數。如果在xy平面上有A2+B2+C2>0{displaystyle A^{2}+B^{2}+C^{2}>0}A^2 +B^2 + C^2 > 0,该偏微分方程在该平面上为二阶偏微分方程。二阶偏微分方程類似以下的圓錐方程:


Ax2+2Bxy+Cy2+⋯=0.{displaystyle Ax^{2}+2Bxy+Cy^{2}+cdots =0.}Ax^2 + 2Bxy + Cy^2 + cdots = 0.

该二阶偏微分方程可分类为:抛物线方程,双曲线方程和椭圆方程,其分类方式为:




  1. B2−AC<0{displaystyle B^{2}-AC,<0}B^2 - AC , < 0 : 椭圆方程;


  2. B2−AC=0{displaystyle B^{2}-AC=0,}B^2 - AC = 0, : 抛物线方程;


  3. B2−AC>0{displaystyle B^{2}-AC,>0}B^2 - AC , > 0 :双曲线方程。



混合形式方程


如果偏微分方程的系数不是一个常数,该偏微分方程可能不属于以上几种类别之一,而可能是混合形式方程。一个简单的例子为Euler–Tricomi方程:


uxx=xuyy{displaystyle u_{xx},=xu_{yy}}<br />
u_{xx} , = xu_{yy}

该方程称为椭圆双曲线方程。因为当x < 0时是椭圆形式,当x > 0时是双曲线形式。



偏微分方程有关问题



适定问题


偏微分方程解中任意函数的出现必然产生解的各种差异,考虑到几乎不知道这些解的详情,在大多数问题中惯常的目标是找满足合适的和确定的条件(例如在空间的边界处和某固定时刻)的那些解,要求这些条件可以确定唯的解是自然的要求。


而且还有更进一步的考虑,即这些条件的大小或量的微小改变在解本身也带来相应地小的改变。


法国数学家阿达马强调后一方面,当解不连续地依赖于原始数据变化时称此问题是不适定的或提得不正确的


  • 不适定的例子

对于双变量的Laplace方程:


2z∂x2+∂2z∂y2=0(y>0){displaystyle {frac {partial ^{2}z}{partial x^{2}}}+{frac {partial ^{2}z}{partial y^{2}}}=0(y>0)}frac{partial^2 z}{partial x^2}+frac{partial^2 z}{partial y^2}=0 (y>0)


在边界条件


z(x,0)=0{displaystyle z(x,0)=0}z(x,0)=0z(x,0)∂y=1ncos⁡nx{displaystyle {frac {partial z(x,0)}{partial y}}={frac {1}{n}}cos nx}frac{partial z(x,0)}{partial y}=frac{1}{n}cos nx


之下,符合条件的解为


z(x,y)=1n2sinh⁡(ny)cos⁡(nx){displaystyle z(x,y)={frac {1}{n^{2}}}sinh(ny)cos(nx)}z(x,y)=frac{1}{n^2}sinh (ny) cos (nx)


n→+∞{displaystyle {begin{smallmatrix}nrightarrow +infty end{smallmatrix}}}begin{smallmatrix} nrightarrow + infty end{smallmatrix}<br />
其数据在y=0{displaystyle {begin{smallmatrix}y=0end{smallmatrix}}}begin{smallmatrix} y=0 end{smallmatrix}<br />
z{displaystyle {begin{smallmatrix}zend{smallmatrix}}}begin{smallmatrix} z end{smallmatrix}<br />
z∂y{displaystyle {begin{smallmatrix}{frac {partial z}{partial y}}end{smallmatrix}}}begin{smallmatrix} frac{partial z}{partial y} end{smallmatrix}<br />
的指定值趋于0,而z(x,y){displaystyle {begin{smallmatrix}z(x,y)end{smallmatrix}}}begin{smallmatrix} z(x,y) end{smallmatrix}<br />
的值在无穷大的范围内震荡,所以这个解不适定。



解析法解偏微分方程


一些有效的解析法解偏微分方程方法:



分离变量法



通过分离变量法减少偏微分方程中的变量,将一个偏微分方程分解成若干个常微分方程。



特征线法



沿着一阶偏微分方程的特征线,偏微分方程简化为一个常微分方程。沿着特征线求出对应常微分方程的解就可以得到偏微分方程的解。



积分变换


利用积分法,将偏微分方程变换为可分离的偏微分方程,方便求解。一般为傅里叶变换分析。



变量变换


通过适当的变量变换,可以简化偏微分方程的求解。一个典型的例子为Black–Scholes方程:


V∂t+12σ2S2∂2V∂S2+rS∂V∂S−rV=0{displaystyle {frac {partial V}{partial t}}+{frac {1}{2}}sigma ^{2}S^{2}{frac {partial ^{2}V}{partial S^{2}}}+rS{frac {partial V}{partial S}}-rV=0} frac{partial V}{partial t} + frac{1}{2}sigma^2 S^2frac{partial^2 V}{partial S^2} + rSfrac{partial V}{partial S} - rV = 0

可以简化为热力方程:


u∂τ=∂2u∂x2{displaystyle {frac {partial u}{partial tau }}={frac {partial ^{2}u}{partial x^{2}}}} frac{partial u}{partial tau} = frac{partial^2 u}{partial x^2}

通过如下变换:



V(S,t)=Kv(x,τ){displaystyle V(S,t)=Kv(x,tau ),} V(S,t) = K v(x,tau),

x=ln⁡(S/K){displaystyle x=ln(S/K),} x = ln(S/K),

τ=12σ2(T−t){displaystyle tau ={frac {1}{2}}sigma ^{2}(T-t)} tau = frac{1}{2} sigma^2 (T - t)

v(x,τ)=exp⁡(−αx−βτ)u(x,τ).{displaystyle v(x,tau )=exp(-alpha x-beta tau )u(x,tau ).,} v(x,tau)=exp(-alpha x-betatau) u(x,tau).,



基本解


非齐次偏微分方程可通过寻找基本算子,然后通过带有初始条件的卷积来解答。
该法常用于信号处理中通过冲激响应来求解滤波器。



叠加原理


因为一个线性齐次偏微分方程解的重叠也可看做一个解,所以可以通过交叉重叠这些解得到偏微分方程的一个解。



数值法解偏微分方程


在众多求解偏微分方程的数值方法中,三种应用最广的方法为有限元法(Finite Element Method, FEM)、有限体积法(Finite Volume Method, FVM)和有限差分法(Finite Difference Method, FDM)。其中,有限元法占主要地位,尤其是它的高效高阶版本—hp-FEM英语hp-FEM。其它版本的有限元法还有:广义有限元法(Generalized Finite Element Method, FFEM)、扩展有限元法(eXtended Finite Element Method, XFEM)、无网格有限元法(Meshfree Finite Element Method)、离散迦辽金有限元法(Discontinuous Galerkin Finite Element Method, DGFEM)等。



参考文献







Comments

Popular posts from this blog

Information security

Lambak Kiri

章鱼与海女图