波动方程和行波法
- 格式:ppt
- 大小:970.00 KB
- 文档页数:117
第3章 波动问题的行波法§3.1 二阶线性方程的分类与化简本节讨论:①两个自变量方程的分类与化简,②多个自变量方程的分类与化简⒈ 两个自变量方程的分类与化简二阶方程的一般形式 二阶变系数方程可写为1112220(,)2(,,,,)(,)xx xy yy x y Lu x y a u a u a u x y u u u f x y =+++Φ= (3.1.1)式中:11a 、12a 、22a 为x 、y 的函数,0(,,,,)x y x y u u u Φ为低阶导数项。
公式关于二阶导数项为线性的,即称方程为准线性的。
若0(,,,,)x y x y u u u Φ关于u 及其x u 、y u 为线性的,则称方程为线性的。
方程的变换 为了简化上述方程,作可逆变换:(,)(,)x y x y ξξηη=⎧⎨=⎩, (,)0(,)J x y ξη∂=≠∂, (,)(,)x x y y ξηξη=⎧⎨=⎩(3.1.2) 代入方程中,不难得到:11122212(,,,,)(,)Lu A u A u A u u u u f ξξξηηηξηξηξη=+++Φ= (3.1.3)式中: 22111112222x x y yA a a a ξξξξ=++ (3.1.4) 12111222()x x x y y x y y A a a a ξηξηξηξη=+++ (3.1.5)22221112222x x y yA a a a ηηηη=++ (3.1.6) 我们化简的目的是使得二次项的项数尽量少,并且值尽量为简单(如0ij A =或1ij A =±)。
顾及ij A 的表达式,取关于z 的一阶非线性偏微分方程2211122220x x y y a z a z z a z ++= (3.1.7)若该方程有解),(1y x z ϕ=、),(2y x z ψ=,则110A =及220A =;公式大大简化了。
数学物理方程复习一.三类方程及定解问题(一)方程1.波动方程(双曲型)Utt = a2Uxx+f; 0<x<l,t>0U(0,t)= Φ1(t);U(l,t)= Φ2(t);U(x,0)= Ψ1(x);Ut (x,0)=Ψ2(x)。
2.热传导方程(抛物型)Ut = a2Uxx+f; 0<x<l,t>0U(0,t)= Φ1(t);U(l,t)= Φ2(t);U(x,0)= Ψ1(x).3.稳态方程(椭圆型)Uxx +Uyy=f; 0<x<a;0<y<b;t>0.U(0,x)= Φ1(x);U(b,x)= Φ2(x);U(y,0)= Ψ1(y);Ut (y,a)=Ψ2(y)。
(二)解题的步骤1.建立数学模型,写出方程及定解条件2.解方程3.解的实定性问题(检验)(三)写方程的定解条件1.微元法:物理定理2.定解条件:初始条件及边界条件(四)解方程的方法1.分离变量法(有界区域内)2.行波法(针对波动方程,无界区域内)3.积分变换法(Fourier变换Laplace变换)Fourier变换:针对整个空间奇:正弦变换偶:余弦变换Laplace变换:针对半空间4.Green函数及基本解法5.Bessel函数及Legendre函数法例一:在弦的横震动问题中,若弦受到一与速度成正比的阻尼,试导出弦阻尼振动方程。
解:建立如图所示的直角坐标系,设位移函数为U(x,t),取任意一小段△x进行受力分析,由题设,单位弦所受阻力为b U t(b为常数),在振动过程中有△x所受纵向力为:(T2COSa2-T1COSa1)横向力为:(T2SINa2-T1SINa1-b U t(x+n△x))(0<n<1). T2,T1为△x弦两端所受的张力,又因为弦做横振动而无纵振动,由牛顿定律有T2COSa2-T1COSa1=0,T2SINa2-T1SINa1-b(x+n△x)U t=p U tt(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/ρ) U t(x+n△x,t)即令△x→0时有:U tt+ aU t=a2U xx例二:设扩散物质的源强(即单位时间内单位体积所产生的扩散物质)为F (x,y,z,t),试导出扩散方程。