数学物理方程复习.docVIP

数学物理方程复习 一．三类方程及定解问题 方程 波动方程（双曲型） Utt = a2Uxx +f; 0xl,t0 U(0,t)= Φ1（t）; U(l,t)= Φ2（t）; U(x,0)= Ψ1（x）; Ut(x,0)=Ψ2（x）。 热传导方程（抛物型） Ut = a2Uxx +f; 0xl,t0 U(0,t)= Φ1（t）; U(l,t)= Φ2（t）; U(x,0)= Ψ1（x）. 稳态方程（椭圆型） Uxx +Uyy =f; 0xa;0yb;t0. U(0,x)= Φ1（x）; U(b,x)= Φ2（x）; U(y,0)= Ψ1（y）; Ut(y,a)=Ψ2（y）。 解题的步骤 建立数学模型，写出方程及定解条件 解方程 解的实定性问题（检验） 写方程的定解条件 微元法：物理定理 定解条件：初始条件及边界条件 解方程的方法 分离变量法（有界区域内） 行波法（针对波动方程，无界区域内） 积分变换法（Fourier变换Laplace变换） Fourier变换:针对整个空间 奇：正弦变换 偶：余弦变换 Laplace变换：针对半空间 Green函数及基本解法 Bessel函数及Legendre函数法 例一：在弦的横震动问题中，若弦受到一与速度成正比的阻尼，试导出弦阻尼振动方程。 解：建立如图所示的直角坐标系，设位移函数为U（x,t）,取任意一小段△x进行受力分析，由题设，单位弦所受阻力为b Ut（b为常数），在振动过程中有△x所受纵向力为：（T2COSa2-T1COSa1）横向力为：（T2SINa2-T1SINa1-b Ut(x+n△x））(0n1). T2,T1为△x弦两端所受的张力，又因为弦做横振动而无纵振动，由牛顿定律有T2COSa2-T1COSa1=0，T2SINa2-T1SINa1-b(x+n△x)Ut=p Utt(x+n△x)△x 在小的振动下SINa1≈TANa1=Ux(x,t), SINa2≈TANa2=Ux(x+△x,t), COSa2≈COSa1≈1,T=T1=T2.(ρ是密度) 即(T/ρ)[ Ux(x+△x,t)- Ux(x,t)]/ △x-(b/ρ) Ut(x+n△x，t) 即令△x?0时有：Utt+ aUt=a2Uxx 例二：设扩散物质的源强（即单位时间内单位体积所产生的扩散物质）为F（x,y,z,t），试导出扩散方程。 解：设U（x,y,z,t）为粒子的浓度（单位体积内的粒子数），在空间内画出一个立方体，体积△V=△X△Y△Z,考虑在△t内△V内的粒子流动情况。 由扩散定律知： 流入X方向的流粒子数为：[qx(x,t)- qx(x+△x,t)] △t△y△z, 流入Y方向的流粒子数为: [qY(y,t)- qY(y+△y,t)] △t△x△z, 流入Z方向的流粒子数为: [qz(z,t)- qz(z+△z,t)] △t△x△y. 而源强产生的粒子数为：F（x,y,z,t）△t△x△y△z. 由质量守恒定律为： [qx(x,t)- qx(x+△x,t)] △t△y△z+[qY(y,t)- qY(y+△y,t)] △t△x△z+[qz(z,t)- qz(z+△z,t)] △t△x△y+ F(x,y,z,t）△t△x△y△z= [U（x,y,z,t+△t）- U（x,y,z,t）] △t△x△y△z. 令△t△x△y△z?0时有：（＠是求偏导） -＠qx/＠x-＠qy/＠y-＠qz/＠z+ F（x,y,z,t）= Ut 由自由扩展定律得： ＠(D＠u/＠x）/＠x+＠（D＠u/＠y）/＠y+＠（D＠u/＠z）/＠z+F= Ut 若扩散粒子是均匀的： Ut= a2△U. 二．线性偏微分方程 （一）二阶线性偏微分方程 LU=a11Uxx+2a12Uxy+a22Uyy+b1Ux+b2Uy+c+f 1.主要部分：a11Uxx+2a12Uxy+a22Uyy 2.判别式△= a212- a11a △0 双曲线方程 △=0 抛物型方程 △0 椭圆方程 3.特征方程 a11（-dy/dx）2-2a12（-dy/dx）+a22=0 特征根：dy/dx=（a12±△1/2）/ a11 特征曲线：y=[（a12+△1/2）/ a11]x+C1 y=[（a12-△1/2）/ a11]x+C2 新旧变量关系：ζ=y+λ1x，η= y+λ2x 令Q=省略 例一：把方程x2Uxx+2xyUxy-3y2Uyy-2xUx+4yUy+16x4U=0改成标准形式，并判断类型。 例二：x2Uxx+2xyUxy+y2Uyy=0 例三：化简2aUxx+2aUxy+aUyy+2bUx+2cUy+U=0,并判断类型。a≠0 （二