数理方程习题集综合

例 1.1.1 设v=v(线x,y),二阶性偏微分方程v xy =xy 的通解。

解 原方程可以写成 ð／ðx（ðv／ðy） =xy 两边对x 积分，得v y =¢（y ）+1/2 x 2Y,其中¢（y ）是任意一阶可微函数。

进一步地，两边对y 积分，得方程得通解为v （x,y ）=∫v y dy+f （x ）=∫¢（y ）dy+f （x ）+1/4 x 2y 2=f （x ）+g （y ）+1/4 x 2y 2其中f （x ）,g （y ）是任意两个二阶可微函数。

例1.1.2即 u(ξ,η) = F(ξ) + G(η),其中F(ξ),G(η)是任意两个可微函数。

例1.2.1设有一根长为L 的均匀柔软富有弹性的细弦，平衡时沿直线拉紧，在受到初始小扰动下，作微小横振动。

试确定该弦的运动方程。

取定弦的运动平面坐标系是O XU ,弦的平衡位置为x 轴，弦的长度为L ，两端固定在O,L 两点。

用u(x,t)表示弦上横坐标为x 点在时刻t 的位移。

由于弦做微小横振动，故u x ≈0.因此α≈0，cos α≈1,sin α≈tan α=u x ≈0,其中α表示在x 处切线方向同x 轴的夹角。

下面用微元法建立u 所满足的偏微分方程。

在弦上任取一段弧'MM ,考虑作用在这段弧上的力。

作用在这段弧上的力有力和外力。

可以证明，力T 是一个常数，即T 与位置x 和时间t 的变化无关。

事实上，因为弧振动微小，则弧段'MM 的弧长dx u xx xx ⎰∆++=∆21s ≈x ∆。

这说明该段弧在整个振动过程中始终未发生伸长变化。

于是由Hooke 定律，力T 与时间t 无关。

因为弦只作横振动，在x 轴方向没有位移，故合力在x 方向上的分量为零，即T(x+x ∆)cos α’-T(x)cos α=0.由于co's α’≈1，cos α≈1,所以T(X+∆x)=T(x),故力T 与x 无关。

于是，力是一个与位置x 和时间t 无关的常数，仍记为T.作用于小弧段'MM 的力沿u 轴方向的分量为 Tsin α’-Tsin α≈T(u x (x+x ∆,t)-u x (x,t)).设作用在该段弧上的外力密度函数为F （x,t ）那么弧段'MM 在时刻t 所受沿u 轴方向的外力近似的等于F(x,t)x ∆.由牛顿第二定律得T （u x (x+x ∆,t)-u x (x,t)+F(x,t)x ∆=ρx ∆tt u ,其中ρ是线密度，由于弦是均匀的，故ρ为常数。

这里tt u 是加速度tt u 在弧段'MM 上的平均值。

设u=u(x,t)二次连续可微。

由微分中值定理得Tu zz （x+θx ∆，t)x ∆+F(x,t)x ∆=ρtt u x ∆, 0<θ<1. 消去x ∆，并取极限x ∆→0得 Tu xx （x,t ）+F(x,t)=ρu tt , 即u tt =ɑ2u xx +ƒ(x,t), 0<x<L,t>0,其中常数ɑ2=T/ρ，函数ƒ（x,t ）=F(x,t)/ρ表示在x 处单位质量上所受的外力。

上式表示在外力作用下弦的振动规律，称为弦的强迫横振动方程，又称一维非齐次波动方程。

当外力作用为零时，即ƒ=0时，方程称为弦的自由横振动方程。

类似地，有二维波动方程u tt =ɑ2（u xx +u y y ）+ƒ（x.y.t ）, (x,y)Ω∈,t>0, 电场E 和磁场H 满足三维波动方程E c E 2222t ∇=∂∂和H c H 2222t∇=∂∂, 其中c 是光速和2222222x zy ∂∂+∂∂+∂∂=∆=∇⋅∇=∇。

例1.2.2设物体Ω在无热源。

在Ω中任取一闭曲面S （图1.2）。

以函数u(x,y,z,t)表示物体在t 时刻，M=M(x,y,z)处的温度。

根据Fourier 热传导定律，在无穷小时段dt 流过物体的一个无穷小面积dS 的热量dQ 与时间dt ，曲面面积dS 以及物体温度u 沿曲面的外法线n 的方向导数三者成正比，即,其中k=k(x,y,z)是在物体M(x,y,z)处的热传导系数，取正值。

我们规定外法线n 方向所指的那一侧为正侧。

上式中负号的出现是由于热量由温度高的地方流向温度低得地方。

故当0nu>∂∂时，热量实际上是向-n 方向流去。

对于Ω任一封闭曲面S ，设其所包围的空间区域为V ，那从时刻t 1到时刻t2经曲面流出的热量为1Q =dSdt nukS⎰⎰⎰∂∂21t t - 设物体的比热容为c(x,y,z)，密度为ρ(x,y,z)，则在区域V ，温度由u(x,y,z,1t )到u(x,y,z)所需的热量为[]dvdt tuc dv t z y x u t z y x u c t t VV∂∂=-=⎰⎰⎰⎰⎰⎰⎰21),,,(),,,(Q 122ρρ. 根据热量守恒定律，有12Q Q -=即[]dSst nukdv t z y x u t z y x t t S⎰⎰⎰⎰⎰⎰∂∂=-21),,,(),,,u c 12V（ρ 假设函数u(x,y,z,t)关于x,y,z 具有二阶连续偏导数，关于t 具有一阶连续偏导数，那么由高斯公式得0][21t =⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂∂∂-∂∂⎰⎰⎰⎰dvdt z u k z y u k y y u k x t u c t Vρ. 由于时间间隔[]21t ,t 及区域V 是任意的，且被积函数是连续的，因此在任何时刻t ，在Ω任意一点都有⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂=∂∂z u k z y u k y y u k y x u ρc(1.2.6)方程称为非均匀的各向同性体的热传导方程。

如果物体是均匀的，此时k,c 及ρ均为常数，令2a =ρc k，则方程(1.2.6)化为 u a z u y u xu ∆=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂=∂∂22222222a t u ,(1.2.7)它称为三维热传导方程若物体有热源，其热源密度函数为，则有热源的热传导方程为),,,(a u 2t t z y x f u +∆= (1.2.8)其中ρc F f =类似地，当考虑的物体是一根均匀细杆时如果它的侧面绝热且在同一截面上的温度分布相同，那么温度只与有关，方程变成一维热传导方程xx u 2t a u = （1.2.9）同样，如果考虑一块薄板的热传导，并且薄板的侧面绝热，则可得二维热传导方程 ）yy 2t u (u +=xx u a （1.2.10）（P16)例1.3.1一长为L 的弹性杆，一端固定，另一端被拉离平衡位置b 而静止，放手任其振动。

试写出杆振动的定解问题。

解 取如图1.3所示的坐标系。

O L L+b x泛定方程就是一维波动方程（杆的纵振动方程） u tt =a 2u xx , 0<x<L.在初始时刻（即放手之时），杆振动的速度为零，即u t (x,0)=0,0≤x ≤L.而在x=L 端拉离平衡位置，使整个弹性杆伸长了b 。

这个b 是来自整个杆各部分伸长后的贡献，而不是x=L 一端伸长的贡献，故整个弹性杆的初始位移为u|0=t =Lbx, 0≤x ≤L. 再看边界条件。

一端x=0固定，即该端位移为零，故有u(0,t)=0,0≤x ≤L.另一端由于放手任其振动时未受外力，故有u x (L,t)=0,t ≥0.所以，所求杆振动的定解问题为u tt =a 2u xx , 0<x<L,t>0, u(x,0)=Lbx, u t （x,0)=0, 0≤x ≤L, u(0,t)=0, u x (L,t)=0, t ≥0.（P17）例1.3.2 ：长为L 的均匀弦，两端x=0和x=L 固定，弦中力为T ，在x=x0处以横向力F 拉弦，达到稳定后放手任其振动。

试写出初始条件。

解：建立如图坐标系。

的两端绝热；（3）杆的一端为恒温零度，另一端绝热。

试写出该绝热传导问题在以上三种情况下的边界条件。

解：设杆的温度为u(x,t),则（1） u （x,t ）=0，u(L,t)=0. （2） 当沿杆长方向有热量流动时，由Fourier 实验定律得Lx x xu kq xukq ==∂∂-=∂∂=201,'其中q1,q2分别为x=0和x=L 处的热流强度。

而杆的两端绝热，这就意味着杆的两端与外界没有热交换，亦没有热量的流动，故有q1=q2=0和,0),0(=t x u 0),(=t L u x .(3)显然，此时有 0),(,0),0(==t L u t u x .例1.5.1求Poisson 方程Uxx +Uyy =X^2 +XY+Y^2的通解 解:先求出方程的一个特解V=V （x ，y),使其满足 Vxx +Vyy=X^2 +XY+Y^2 由于方程右端是一个二元二次齐次多项式，可设V （x ，y) 具有形式V(x,y)=aX^4 +bX^3 Y+cY^4,其中a,b,c 是待定常数 Vx=4aX^3+3bX^2 Y Vy=bX^3+4cY^3 Vxx=12aX^2+6bXY Vyy=12cY^2 得Vxx+Vyy=12aX^2 +6bXY+12cY^2=X^2 +XY+Y^2 比较两边系数，可得 a=1/12,b=1/6,c=1/12于是V （x,y)=1/12(X^4 +2X^3 Y+Y^4）下面求函数W=W(x,y),使其满足Wxx+Wyy=0.作变量代换e=x,n=iy(为d)Ue=du/de=du/dx=Ux Un=du/dn=du/dy *dy/dn=-iy Uee=dUe/de=Uxx Unn=-Uyy 可得Wee-Wnn=0 再作变量代换s=e+n,t=e-nUe=du/de(s,t)=Us+Ut Un=du/dn=Us-Ut Uee=dUe/de=d(Us+Ut)/de=Uss+Utt+2Ust Unn=dUn/dn=d(Us-Ut)/dn=Uss+Utt-2Ust 那么方程进一步化为Wst=0其通解为W=f(s)+g(t)=f(e+n)+g(e-n)=f(x+iy)+g(x-iy),其中f,g 是任意两个二阶可微函数。

那么根据叠加原理，方程的通解为u(x,y)=V+W=f(x+iy)+g(x-iy)+1/12(X^4+2X^3 Y+Y^4)(P32)例2.1.1 判断方程U xx +2U xy -3U yy +2U x +6U y =0（2.1.22）的类型，并化简。

解： 因为a 11= 1,a 12= 1,a 22= -3,所以 =a 212-a 11a 22=4>0,故方程为双曲型方程。