数理方程习题解答

竖立方程习题解答 习题一 1 习题二 3 习题三 5 习题四 错误！未定义书签。 习题五 错误！未定义书签。 1 习题 1.1.4 1．导出弦受阻力的波动方程 其中阻力与速度成正比， 为常数． 解 我们考虑弦的一个微元。令 为端点处的张力，如 教材 民兵爆破地雷教材pdf初中剪纸校本课程教材衍纸校本课程教材排球校本教材中国舞蹈家协会第四版四级教材 图 1.1 所示，沿锤直方向作 用在这个微元上的力是 ,阻力为 ,由牛顿(Newton)第二定律1， 此合力等于质量乘以加速度．因此 (1) 其中 是密度， 是微元弦的弧长．因为运动弦的斜率是很小的，故有 s xΔ ≈ Δ . 因角 和 也很小，所以我们有， ，于是(1)式变成 （2） 但由微积分学我们知道，在时刻 , ,于 是，方程(2)便可写成 令 取极限，我们求得 （3） 其中 2. 设长度为的均匀弹性杆的线密度为 ，杨氏模量为 ，试列出杆的微小纵振动方程。 解 考虑杆在无外力作用下的振动。取杆的一端为原点，干的方向为 轴建立坐标系： 则杆上各点 在时刻 的位移是 。 在杆上任取一段，其两端点静止时的坐标为 ,此小杆段在时刻 的相对伸长 为： ，令 得 点在时刻 的相对伸长为 ( , )xu x t ，由Hooke 定 律知张力为 ，再此小杆段上用Newton 第二定律得 两边同除 并令 得： 若杨氏模量为 为常数则得： 。 1 牛顿(Newton)第二定律与动量守恒定律等价，也可以用动量守恒定律来见方程，见《数学物理方程 讲义》 （姜礼尚、陈亚浙）P1 习题 1.2.4 3 习题 1.2.4 1 设悬浮粒子由重力引起的沉淀速度 是不变的，又假定在同一水平面上粒子的浓度 是相同的，试给出悬浮粒子的浓度 所满足的方程． 解 取竖直向下的方向为 轴，考虑介于平面 之间，截面积为常数 的柱体。质 量守恒关系为 其中 为这段时间内柱体内粒子质量的增加，而 为 这段时间内由于扩散作用经由柱体的上、下底面进入柱体内的粒子质量， 是由于 沉淀经柱体的上、下底面进入柱体的粒子质量。 其中 是扩散系数。 从而 注意到 与 的任意性，由上式立即得 2． 没有一厚为 l的无限平面板，在其表面与温度为 的周围介质发生热交换，如果 板的温度不随其厚度而变化(即在垂直于板面的直线上的点的温度均相同)，试给出板 冷却的初边值问题。 解 取平面 上任意一区域 ，在从 到 ( )这段时间内考察柱体 中热量的平衡关系： 其中 为从 到 这段时间内， 中温度的变化所吸收的热量，而 与 则分别为这 段时间内，通过 由板内其他部分流入 的热量以及通过上、下板面与周围 介质的热交换所获得的热量。不难计算得 数学物理方程习题解答 4 从而 所以 其中 。此外还有初始条件 于是得二维Cauchy问题2 2 取垂直于板平面的方向为 轴，则在温度 依赖于 的情况下，所讨论的是由两平面 所界的无界区城内的问题，边界条件给在两平面上，此时边值问题为 在温度 u 不依赖于 z 的情况下，不能简单地由三维方程得出 u 满足二维热传导方程， 因为若将所给问题当作二维问题，则此时上、下板面与周围介质的热交换不能再当作 边界条件来处理，而应考虑到方程中去。也就是说，边界条件和方程不是截然分开的 两个不同的东些，而是同一事物的不同表现方式。在齐次方程和非齐次方程、齐次边 界条件和非齐次边界条件之间的转换时，就会从数学上遇到此问题。 习题 1.3.3 5 习题 1.3.3 1．(《数学物理方程讲义》 姜礼尚、陈亚浙P34,17) 设 ( )2 2 21 1( ) | | ( )2 2J v v v dx a x v ds fvdx gvdsΩ ∂Ω Ω ∂Ω= ∇ + + − −∫ ∫ ∫ ∫ 其中 ( ) 0a x ≥ 。考虑以下三个问题： 问题 I(变分问题) ：求 1 ( )u M C∈ = Ω 使得 ( ) min ( ) u M J u J v∈= 问题 II：求 1 ( )u M C∈ = Ω 使得它对于任意v M∈ 都满足 ( ) ( )( ) 0u v u v fv dx a x uv gv ds Ω ∂Ω ∇ ⋅∇ + ⋅ − + − =∫ ∫ 问题 III (第三边值问题)：求 2 1( ) ( )u C C∈ Ω ∩ Ω 满足以下边值问题 ( ) u u f x u a x u g x n −Δ + = ∈Ω⎧⎪∂⎨ + = ∈∂Ω⎪∂⎩ (1) 证明问题 I 与问题 II等价． (2) 当 2 1( ) ( )u C C∈ Ω ∩ Ω 时，证明问题 I、II 、III 等价． 证明 设问题 I(变分问题)成立，即 1 ( )u M C∈ = Ω 使得 ( ) min ( ) u M J u J v∈= ，则对 v M u v Mε ε∀ ∈ ∀ ∈ + ∈R � ，从而 ( )2 2 2 1( ) ( ) | ( ) | ( ) 2 1 ( )( ) ( ) ( ) 2 j J u v u v u v dx a x u v ds f u v dx g u v ds ε ε ε ε ε ε ε Ω ∂Ω Ω ∂Ω = + = ∇ + + + + + − + − + ∫ ∫ ∫ ∫ 在 0ε = 取的最小值，所以 (0) 0j′ = ，即 ( )(0) ( ) 0j u v uv dx a x uvds fvdx gvds Ω ∂Ω Ω ∂Ω ′ = ∇ ⋅∇ + + − − =∫ ∫ ∫ ∫ （1） 所以问题 I⇒问题 II。由于 ( )j ε 是一个二次函数，故 (0)j 取最小值等价于 (0) 0j′ = ，也就是问题 I⇔问题 II。 当 2 1( ) ( )u C C∈ Ω ∩ Ω 时由高斯公式得 6 ( ) uu vdx v u v u dx v dS v udx nΩ Ω ∂Ω Ω ∂∇ ⋅∇ = ∇ ∇ − Δ = − Δ⎡ ⎤⎣ ⎦ ∂∫ ∫ ∫ ∫ 所以 ( ) ( ) 0uu u f vdx a x u g vds v M nΩ ∂Ω ∂⎛ ⎞−Δ + − + + − = ∀ ∈⎜ ⎟∂⎝ ⎠∫ ∫ 由变分引理得 ( ) u u f x u a x u g x n −Δ + = ∈Ω⎧⎪∂⎨ + = ∈∂Ω⎪∂⎩ 所以问题 II⇒问题 III。而问题 III⇒问题 II是显然的，所以，问题 II⇔问题 III。 习题 2.1.3 1 将下列方程化成标准型 ( ) ( ) 2 2 2 2 2 2 2 (1) 0, 0; (2) 0, 0, 0; (3) 4 4 2 0, (4) 0, 0; (5) sin 2 sin 0; (6) 1 1 0. xx yy xx yy xx xy yy y xx yy xx xy yy xx yy x y u x u x y u x u x y u u u u u yu y x u y x u y u x u y u xu yu + = > + = > > + + − = + = ≠ ⋅ − ⋅ + = + + + + + = 解（1）由特征方程 解得两簇共轭特征线 ，做变换 ，则由链式法则原方程化为 （2）由特征方程 解得两簇共轭特征线 ，做变换 ，则由链式法则原方程化为 （3）由特征方程 解得一簇特征线 ，做变换 ，则由链式法则原方程化为 （4）由特征方程 2 2 0dy ydx+ = 当 0y > 时解得一簇共轭特征线2 y xi c± = ，做变换 2 ,y xξ η= = ，则由链式 法则原方程化为 1 0u u uξξ ηη ξξ+ − = 7 当 0y < 时解得两簇特征线2 y x c− ± = ，做变换 2 , 2x y x yξ η= + − = − − ，则由链式法则原方程化为 ( )1 0 2( ) u u uξη ξ ηξ η+ − =− （5）由特征方程 2 2 2 2sin 2 sin 0xdy y xdxdy y dx− + = 解得一簇特征线 2 xytg c= ，做变换 , 2 xytg yξ η= = ，则由链式法则原方程化为 2 2 2 0u uηη ξξ ξ ξ η− =+ （6）由特征方程( ) ( )2 2 2 21 1 0x dy y dx+ + + = 解得两簇共轭特征线 ( ) ( )2 2ln 1 ln 1x x i y y c+ + ± + + = ，做变换 ( )2ln 1x xξ = + + ， ( )2ln 1y yη = + + ，则由链式法则原方程化为 0u uξξ ηη+ = 。 2 证明常系数方程 0xy x yu au bu cu+ + + = 必可通过未知函数变换化为 1 0xyv c v+ = 。 证明 令 ( ) 2ax byu ve− += ，则 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( )( ) ,bx ay bx ay bx ay bx ayx x y y bx ay bx ay bx ay bx ay xy xy x y xy x y bx ay bx ay bx ay bx ay xy x y bx ay bx ay bx ay bx ay bx x y u v e bve u v e ave u v e av e bv e abve u au bu cu v e av e bv e abve a v e bve b v e ave cve − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − = − = − = − − + ∴ + + + = − − + + − + − + ( ) ( ) ( ) ( ) ( ) 0 0 ay bx ay bx ay xy xy v e c ab ve v c ab v + − + − += + − = ∴ + − = 习题 2.2.3 1 证明常系数椭圆型方程必能通过自变量与未知函数的变换化为 0u cuΔ + = 的形 式。 证明： 的系数矩阵 是正定的。取 8 正交矩阵 使得 。做变换 则原方程化为 。在令 可得 。最后做变换 ，则得 。 2 证明常系数双曲型方程必能通过自变量与未知函数的变换化为 的形式。 证明： 的系数矩阵 的特征值全部 为零，且有一个与其 n 个异号。取正交矩阵 使得 。做变换 则原方程化为 。在令 可得 。最 后做变换 ，则得 。 习题 2.3.1 1 求下列方程的特征方向 1 1 2 2 3 3 4 4 (1) (2) (3) x x x x x x x x tt xx yy zz t xx yy u u u u u u u u u u u + = + = + + = − 解 （1）设 ( )1 2 3 4, , ,α α α α α= ，则 2 2 2 21 2 3 4α α α α+ = + ，取单位特征方向， 2 2 2 2 1 2 3 4 1α α α α+ + + = 。所以， 2 2 2 21 2 3 4 12α α α α+ = + = 。记 1 2 3 4 1 1 1 1cos , sin , cos , sin 2 2 2 2 α θ α θ α ϑ α ϑ= = = = ，则 1 1 1 1cos , sin , cos , sin 2 2 2 2 α θ θ ϑ ϑ⎛ ⎞= ⎜ ⎟⎝ ⎠。 （2）设 ( )1 2 3 4, , ,α α α α α= ，则 2 2 2 21 2 3 4α α α α= + + ，取单位特征方向， 2 2 2 2 1 2 3 4 1α α α α+ + + = 。所以， 2 2 2 21 2 3 4 12α α α α= + + = 。记 1 2 3 4 1 1 1 1, sin sin , sin cos , cos 2 2 2 2 α α θ β α θ β α θ= ± = = = ，则 1 1 1 1, sin sin , sin cos , cos 2 2 2 2 α θ β θ β θ⎛ ⎞= ±⎜ ⎟⎝ ⎠。 9 （3）设 ( )1 2 3, ,α α α α= ，则 2 22 30 α α= − ，取单位特征方向， 2 2 21 2 3 1α α α+ + = 。 所以， 2 21 22 1α α+ = 。记 1 2 31 1cos , sin , sin2 2α θ α θ α θ= = = ± ，则 1 1cos , sin , sin 2 2 α θ θ θ⎛ ⎞= ±⎜ ⎟⎝ ⎠。 2 对波动方程 ( )2 0tt xx yyu a u u− + = 过直线 : 0, 2l t y x= = 的特征平面。 解 设单位特征方向 ( ), ,n α β γ= ，则 2 2 2 2 2 2 2 ( ) 0 (1), 1 (2) aα β γ α β γ − + = + + = 易见直线 l 过了原点的。设所求特征平面为 0t x yα β γ+ + = ，直线 l 过点(0,1,2) 所以 2 0 (3)β γ+ = 。由（1）、（2）、（3）得特征平面为 5 2 0at x y± + − = 。 习题 3.1.7 1 求 的一般解解。（提示：令 ） 解3 令 ，则 化为 ， 其通解为 ，从而原方程的一般解为 。 2 求 2 2 2 2 u ux x x x y ∂ ∂ ∂⎛ ⎞ =⎜ ⎟∂ ∂ ∂⎝ ⎠ 一般解。 （提示： ( , ) ( , )v x y xu x y= ） 解 令 ( , ) ( , )v x y xu x y= ，则 2 2 2 2 u ux x x x y ∂ ∂ ∂⎛ ⎞ =⎜ ⎟∂ ∂ ∂⎝ ⎠ 化为 2 2 2 2 0 v v x y ∂ ∂− =∂ ∂ ，通解 为 ( , ) ( ) ( )v x y x y x yϕ ψ= + + − ，从而原方程的一般解为 3 通过未知函数的变换求解，习题2 思想一样。问题是如何寻找这种变换，这是非常困难的。 10 ( )1( , ) ( ) ( )u x y x y x y x ϕ ψ= + + − 。 3 求 一般解。 （提示： ） 解4 特征方程为 ，解得两簇特征线 , 。做变换 ，则原 方程化为 所以 ,故 。 4 求 2 2 2 2 2 2 22 0 u u u u ux xy y x y x y x yx y ∂ ∂ ∂ ∂ ∂− + + + =∂ ∂ ∂ ∂∂ ∂ 的一般解。 （提示： xy yξ η= = ） 解5 特征方程为 dy y dx x = − ，解得特征线为 xy c= ，令 xy yξ η= = ， 则原 方程化为 2 2 1 1 ,u u v uv vη η η η η η ⎛ ⎞∂ ∂ ∂ ∂= − ⇒ = − =⎜ ⎟∂ ∂ ∂ ∂⎝ ⎠ 所以 ln ln ln ( )v η ϕ ξ= − + 。从而， ( )u ϕ ξη η ∂ =∂ 。于是 ( , ) ( ) ln ( )u x y xy y xyϕ ψ= + 5 求初值问题的解 （提示： ） 解 做变换 ，则原方程化为 。其一般解为 。从而 ，由初始条解得 解得 4 化方程为标准型的思想的应用。 5 与习题3思想一样，这里化简后的方程是一个含参数的特殊的二阶常微分方程. 11 于是 . 6 求解问题 2 2 2 2 2 2 sin sin 2cos sin sin 0 ( ) ( ) y x y y x u u u ux x x x y yx y u g x u h x = = ⎧∂ ∂ ∂ ∂+ − − =⎪ ∂ ∂ ∂∂ ∂⎪⎪ =⎨⎪ ′ =⎪⎪⎩ 解 由特征方程解得特征两簇线 sin sinx y x c x y x c+ − = − + = ，做变换 sin sinx y x x y xξ η= + − = − + ，则原方程化为 2 0uξ η ∂ =∂ ∂ 。其一般解为 ( , ) ( ) ( )u ξ η ϕ ξ ψ η= + 。从而 ( , ) ( sin ) ( sin )u x y x y x x y xϕ ψ= + − + − + ， 由初始条件得 ( ) ( ) ( ) ( ) ( ) ( ) x x g x x x h x ϕ ψ ϕ ψ + =⎧⎨ ′ ′− =⎩ 解得 0 0 1 1( ) ( ) ( ) ( ) ( ) ( ) 2 2 x x x g x h x dx x g x h x dxϕ ψ⎡ ⎤ ⎡ ⎤= + = −⎢ ⎥ ⎢ ⎥⎣ ⎦ ⎣ ⎦∫ ∫ 从而 sin sin ( sin ) ( sin ) 1( , ) ( ) 2 2 x x y x x y g x y x g x y xu x y h x dx − + + − + − + − += + ∫ 7 求解定解问题 （提示：令 ） 解 由特征方程解得特征两簇线 ，令 ，则原方程化为 。其一般解为 。从而 ，由初始 条件得 12 从而， 。从而解得 所以 . 8. 求解Goursat 问题 其中 为已知常数．并画图指出定解条件的决定区域。 解6 方程的一般解为 带入定界条件解得 9. 求解Cauchy 问题 解 由 Dalembert 公式解得 i ) . ii ) 10. 求解非齐次方程的Cauchy问题 解 由推广的Dalembert公式解得 6 Goursat问题也称特征编制问题，其定解条件是给在特征线上的。 13 或者令 则原问题化为 再用Dalembert 公式解之。 11. 求解半无界问题 解 (i) 代入第一类边界条件的半无界问题的公式可得 (ii) 代入第一类边界条件的半无界问题的公式可得 另解7：原方程的通解为 ，由初始条件得 解得 ,所以 12. 求解半无界问题 7 为何没有用到边界条件？这是因为函数x，xcosx本身是奇函数，便界条件是自然满足的。若将xcosx 改为xsinx，则就不能这样解。 14 解8 原方程的通解为 ，由初始条件得 解得 ，又由边界条件得9 由于 ，所以 ，从而我们有 故原问题的解为 13. 求无限长理想传插线上传播着的电压和电流．设初始电压分布为 ；初始电流分布为 。 解 (1) 根据传输线方程和电学上的知识以及所给条件，我们得关于传轴线上电压 所满足的定解问题 （2）同理，电流 应满足的定解问题为 由Dalembert 公式解得] (1) ; (2) 8 这是一个特殊解法，一般地说，应该先将边界条件齐次化。 9利用这个条件来对g做延拓 15 14. 求解Cauchy问题 解 由推广的Dalembert公式得 . 另解：令 ,则原问题化为 ， 由Dalenbert 公式解得 . 15． 求解以下Cauchy问题 （1） ( ) 2 3 0 2 2 2 0 0 0 ( , , , ) (0, ) 0 0 tt tt t u a u x y z t u r R u u r x y z r R= = ⎧ − Δ = ∈ × ∞⎪ <⎧⎨ = = = + +⎨⎪ ≥⎩⎩ R (2) ( ) 0 , , 0 ( , ,0) 3 2 , ( , ,0) 0 , tt xx yy t u u u x y t u x y x y u x y x y − + = −∞ < < ∞ >⎧⎪⎨ = + = −∞ < < ∞⎪⎩ (3) ( ) ( ) ( ) 0 , , 0 ( , , ,0) , ( , , ,0) 3 , , tt xx yy zz a x y z a x y z t u u u u x y t u x y z e u x y z ae x y z+ + + + − + + = −∞ < < ∞ >⎧⎪⎨ = = −∞ < < ∞⎪⎩ 解 （1）由Poisson 公式得 00 0 0 ( ) ( ) ( ) ( )( , ) , ( ) 2 0 u r Rr at u r at r at u r atu r t u r r r R <⎧− − + − −= = ⎨ ≥⎩ （2）取 3 2v x y= + ，则 3 2v x y= + 满足初始条件，令u w v= + ，则w满足Cauchy 问题， ( ) 0 , , 0 ( , ,0) 0, ( , ,0) 0 , tt xx yy t w w w x y t w x y w x y x y − + = −∞ < < ∞ >⎧⎪⎨ = = −∞ < < ∞⎪⎩ （*） 易见（*）的解为 0w ≡ ，所以原问题的解为 ( , , ) 3 2u x y t x y= + 。 （3）取 ( )a x y z atv e + + += ，则 ( )a x y z atv e + + += 满足初始条件，令u w v= + ，则w 满足Cauchy问题， ( ) 0 , , 0 ( , ,0) 0, ( , ,0) 0 , tt xx yy zz t w w w w x y t w x y w x y x y − + + = −∞ < < ∞ >⎧⎪⎨ = = −∞ < < ∞⎪⎩ （*） 易见（*）的解为 0w ≡ ，所以原问题的解为 ( )( , , , ) a x y z atu x y z t e + + += 。 16. 验证形式为 ( , , , ) ( )u x y z t f x y z atα β γ= + + + 的函数(其中 , ,α β γ 为任意 16 满足 2 2 2 1α β γ+ + = 的实常数， 2f C∈ 满足方程 2 ( ) 0tt xx yy zzu a u u u− + + = 。 这种形式的波动方程的解称为平面波解。 解 由 ( , , , ) ( )u x y z t f x y z atα β γ= + + + 得 ( ) ( ) 2 2 2 2 2 2 2 2 2 2 ( ) , ( ) ( ) , ( ) 0 xx yy zz tt tt xx yy zz u f x y z at f u f x y z at f u f x y z at f u f x y z at a f u a u u u a f a f α β γ α α β γ β α β γ γ α β γ α β γ = + + + = = + + + = = + + + = = + + + = ∴ − + + = − + + = 习题 3.2.6 1 求下列函数的Fourier变换 （1） 解 记 ，则 ， 从而 (2) (3) 解 ， 所以 用幂级数法解此一阶线性常微分方程Cauchy问题得 另解： 17 同理可得 2 解定解问题 （1） 解 由热传到方程的Cauchy 问题的解的Piossion 公式计算得 （2） 解 先作奇延拓再由热传到方程的Cauchy 问题的解的Piossion 公式计算得 18 （3） Hint: 偶延拓 解 先作偶奇延拓再由热传到方程的Cauchy问题的解的Piossion 公式计算得 3 求解在上半平面 的静电场的电势，已知， 解 因为电势函数 满足方程Laplace 方程，从而有 对 进行富里哀变换10，有 , 解得 由 得 1 0c = ，由 得 ，所以 从而由 得 4 求解Cauchy 问题 其中， 为常数。 解 先作变换把方程化成标准型。令 ，其中 为待定常数。代入原方程 计算得 取 ，则原问题化为 10 这里先只考虑了边界条件 ( ,0) ( )u x f x= ，在用Fourier 变换法解出像函数时，再考虑无穷远边 界条件。 19 由热传导方程Cauchy 问题Poisson 公式解得 从而 5 两种不同 材料 关于××同志的政审材料调查表环保先进个人材料国家普通话测试材料农民专业合作社注销四查四问剖析材料 的均匀细杆的热传导现象可归结为定解问题 2 0 2 0 0, 0, 0 (1) ( ), 0 (2) 0, 0 , 0 (3) ( ), 0 (4) ( 0, ) ( 0, ), ( 0, ) ( 0, ) (5) t l xx t t r xx t l x r x u a u x t u x x v a v x t v x x u t v t k u t k v t ϕ ψ = = ⎧ − = −∞ < < >⎪ = −∞ < <⎪⎪ − = < < ∞ >⎨⎪ = < < ∞⎪⎪ − = + − = +⎩ 其中 ,l rk k 分别是两种材料的导热系数。另外设 (0) (0)ϕ ψ= ，是求解上述定解问题。 解 设法将问题化为（1）或（3）的Cauchy问题来讲解决。设 ( , ), ( , )u x t v x t 为问题 （1）和(3)的解，则 0x > 时， ( , ), ,r l au x t v x t a ⎛ ⎞− ⎜ ⎟⎝ ⎠ 均为（1）的解，令 ( , ) ( , ) , 0r l au x t u x t v x t x a α β ⎛ ⎞= − + >⎜ ⎟⎝ ⎠ % 其中 ,α β 为待定常数。 显然， 0x > 时 ( , )u x t% 是（1）的解。令 20 ( , ) 0, 0 ( , ) ( , ) 0, 0 u x t x t U x t u x t x t > >⎧= ⎨ < >⎩ % 取 ,α β 使得 ( , ), ( , )xU x t U x t 在 0x = 处连续。即 ( 0, ) ( 0, ), ( 0, ) ( 0, )x xu t u t u t u t− = + − = +% % （6） 由 ( , ) ( , ) , 0r l au x t u x t v x t x a α β ⎛ ⎞= − + >⎜ ⎟⎝ ⎠ % 得 ( 0, ) ( 0, ) ( 0, ), ( 0, ) ( 0, ) ( 0, )rx x x l au t u t v t u t u t v t a α β α β+ = − + + + = − + +% % 由（5）、（6）得 1 1l r r l k a k a α β α β + =⎧⎪⎨− + =⎪⎩ （7） 于是， ( , )U x t 满足 2 0 0, 0, 0 ( ) 0 ( ) 0 t l xx r lt U a U x t ax x x U a x x αϕ βψ ϕ = ⎧ − = ≠ >⎪ ⎧ ⎛ ⎞⎪ − + >⎨ ⎪ ⎜ ⎟= ⎨ ⎝ ⎠⎪ ⎪⎪ <⎩⎩ 由Poisson 公式解得 ( ) ( ) ( ) 2 20 2 2 2 2 0 1( , ) ( ) exp exp 4 42 1 exp 42 l ll r l ll x x U x t d a t a ta t xa d a a ta t ξ ξϕ ξ α ξπ ξψ ξ ξπ −∞ ∞ ⎡ ⎤⎛ ⎞ ⎛ ⎞− +⎢ ⎥⎜ ⎟ ⎜ ⎟= − + −⎜ ⎟ ⎜ ⎟⎢ ⎥⎝ ⎠ ⎝ ⎠⎣ ⎦ ⎛ ⎞−⎛ ⎞ ⎜ ⎟+ −⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠ ∫ ∫ 于是，得 ( , ) ( , ) 0 1( , ) , , , 0l l r r u x t U x t x a av x t U x t u x t x a a α β β = < ⎛ ⎞ ⎛ ⎞= − − >⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠ 21 习题 3.3.7 1 一根均匀弦固定于 x=0和 x=l两端，在开始的瞬间，它的形状是一条以过 2x l= 点 的铅直线为对称轴的抛物线，假定没有初速度，确定弦上诸点距离平衡状态的位移。 解 原问题归结为初边值问题 直接用用变量分离法解之得 2 在一端固定而另—端自出的条件下，园柱状杆的纵向微小振动的方程式问题归结为: 其中 是坐标为 的截面的位移，是杆长， 为单位长度的质量， 是杨氏系数。 解 直接用变量分离法解之得 其中 3 求解初边值问题 解 直接用变量分离法解之得 22 其中 是 的解 4 求单位圆内的调和函数，他在边界上的值为 解 在极坐标下研究，原边值问题化为 以分离变数形式的试探解 代入泛定方程得 上式左边是 的函数，与 无关；右边是 的函数，与 无关．两边不可能相等，除非 两边实际上是同一个常数．把这常数记作 λ，这就分解为两个常微分方程 (1) (2) 常微分方程（1）隐含着一个附加条件。事实上，一个确定地点的极角可以加减 的 整倍数，而电势 在确定的地点应具确定数值，所以 ，即 亦即 (3) 这叫做自然的周期条件。常微分方程（1）与周期条件（3）构成本征值问题。求解得 本征值和本征函数为 以本征值代入常微分方程（2）得 这是欧拉型常微分方程，作代换 ，即 方程化为 其解为 这样，由于 u在D内有界， 所以D＝0，于是分离变数形式的解是 23 拉普拉斯方程是线性的，它的一般解应是所有本征解的叠加，即 (4) 为确定（4）中的系数，把（4）代入原边值问题的边界条件．先代入齐次边界条件得 从而， 代入（4）计算得 ,Poisson 公式 5 弦在一端 是固定的，而在另一端 受到谐和的扰动力，因而引起该端的位 移 。研讨这个弦的横向振动． 解 问题归结为： 1) 分解问题：令 ，则得 2) 用分离变量法求(I)的一个特解， 3) 解(II)应用分离变量法可得 6 试证拉普拉斯方程(用球坐标) 24 有特解 ，其中 为任意常数， 是 阶 Legendre11多项式。 证明 考虑球对称的情形，即与φ无关的解12。由球对称原方程化为 令 代入方程得 令 得 令 代入（1）得 ，解得 ，从而（1） 的两个线性无关解为 。 令 ，则（2）化为， （Legendre 方程） 他有解 （ 为 次勒让德多项式）。从而原方程由特解为 7 设有一弹簧一端固定，另一端在外力作用下作周期运动，当振动进行很多次后，初 始条件产生的影响因阻力而消失，它与初始条件无关，问题归结为： 2 0 (0, ) 0 ( , ) sin tt xxu a u u t u l t A tω ⎧ − =⎪ =⎨⎪ =⎩ 试求它的周期解。 解 设 ( , ) ( ) i tU x t X x e ω= ，满足 (0) 0, ( )X X l A= = ，则 ( , )U x t 的虚部就是本问 题的解13，将 ( , )U x t 代入原定解问题得 勒让德（1752～1833）Legendre，Adrien-Marie 法国数学家。1752年9月18日生于巴黎 ，1833 年1 月10日卒于同地。1770 年毕业于马萨林学院 。1782 年以外弹道方面的 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 获柏林科学院奖。1783 年被选为巴黎科学院助理院士，两年后升为院士。1795年当选为法兰西研究院常任院士。1813年继 任 J.-L.拉格朗日在天文事务所的职位。 勒让德的主要研究领域是 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 学(尤其是椭圆积分理论)、数论、初等几何与天体力学，取得了许多 成果，导致了一系列重要理论的诞生。勒让德是椭圆积分理论奠基人之一。在L.欧拉提出椭圆积分加 法定理后的40年中，他是仅有的在这一领域提供重大新结果的数学家。但他未能像N.H.阿贝尔和C.G.J. 雅可比那样洞察到关键在于考察椭圆积分的反函数 ，即椭圆函数。在关于天文学的研究中，勒让德 引进了著名的“勒让德多项式” ，发现了它的许多性质 。他还研究了B函数和Γ函数，得到了Γ函 数的倍量公式。他陈述了最小二乘法，提出了关于二次变分的“勒让德条件”。 勒让德对数论的主要贡献是二次互反律，这是同余式论中的一条基本定理。他还是解析数论的先 驱者之一，归纳出了素数分布律，促使许多数学家研究这个问题。 12 实际上，完全球对称的解应该是只与 r有关，而与无φ、θ关的解，这里选的是与φ无关的特解。 考虑只与θ无关的解会得到什么结果？ 13 这不是一个标准的定解问题，因为从常规来看，它差一个定解条件-----初始条件。但是，题目明确 说初始条件因长时间运动而消除了。由于振动纯属边界振型造成，所以弦振动的角频率也应该为ω， 25 2 2 0 (0) 0, ( ) X X a X X l A ω⎧ ′′ + =⎪⎨⎪ = =⎩ 解得 1 2( ) cos sinX x c x c xa a ω ω= + 由边界条件得 1 20, sinc c l Aa ω= = 所以 ( ) sin sin AX x x al a ω ω= 从而 ( , ) sin sin i tAU x t xe al a ωω ω= 于是 ( , ) Im ( , ) sin sin sin Au x t U x t x t al a ω ωω= = 。 8 求解没有初始条件的半无界杆的热传导问题 的周期解。 解 设14 ，则 ，于是有 特征方程为： 解得特征根为 所以 由 得 ，由 得 ，于是 所以可用此特殊的变量分离法。 14这不是一个标准的定解问题，因为从常规来看，它差一个定解条件-----边界条件。但试题上明确说明 是求周期解，方程是齐次的，所以周期解得周期应该与初始函数的周期相同。特殊的变量分离法，这 里用到一个特殊的边界条件lim ( ) x X x→∞ =有界 . 26 。从而原问题的解为 9 求解振动向题 解 用叠加原理、齐次化原理和变量分离法解得 其中 10 求解热传导方程定解问题 解 先作变换 ，然后再用变量分离法解之得 其中 ，这里 是特征函数， 而 则是权函数。 27 2 2 2 2 2 2 2 4 2 1 2 0 2 2 , [0, ] ( , ) sin 2 ( ) sin ( ) sin na b bt xl a a n n bt x a n b x a n b x athey are orthogonal on l with weight function nu x t e e x l nwhere x e xdx l l nNote the eigen functions are X x e x l e π πϕ πϕ ϕ π ⎡ ⎤⎛ ⎞− +⎢ ⎥∞ ⎜ ⎟⎝ ⎠⎢ ⎥⎣ ⎦ = − − = ⋅ = = ∑ ∫ l1. 求解细杆导热问题秆长 l，初始温度均匀为 u0。两端分别保持温度u1和 u2。 解 数学模型为 2 0 1 2 0 ( ,0) (0, ( , ) t xxu a u u x u u t u u l t u ⎧ − =⎪ =⎪⎨ =⎪⎪ =⎩ 作变换 2 11 u uu v u x l −= + + 将边界条件其次化再用变量分离法解之得 ( )2 2 1 1 1 0 2 1 1 ( , ) ( )2 ( 1) 1 ( 1) sin na l tn n n u uu x t u x l u u u u n x e n n l ππ π π ∞ − = −= + − −⎧ ⎫⎡ ⎤+ − − + − ×⎨ ⎬⎣ ⎦⎩ ⎭∑ 12 解一个具有热绝缘的侧表面均匀杆冷却问题，假定这个汗的初始温度分布情况是 ( ,0) ( )u x xϕ= ，它的一端热绝缘，而另一端永远保持为 0u ，问题归结为， 2 0 0 0 , 0 ( ,0) ( ) 0 (0, ) 0, ( , ) t xx x u a u x l t u x x x l u t u l t u ϕ ⎧ − = < < >⎪ = ≤ ≤⎨⎪ = =⎩ 解 作变换 0u u v= + ，然后用分离变量法解得 22 1 2 0 0 0 2 2 1 ( 1) 4 2 1( , ) ( )cos cos 2 (2 1) 2 n l na t l n n u nv x t e d x l l n l π ϕ ξ πξ ξ ππ +⎛ ⎞∞ −⎜ ⎟⎝ ⎠ = ⎡ ⎤⎛ ⎞+ − +⎢ ⎥= −⎜ ⎟+⎢ ⎥⎝ ⎠⎣ ⎦ ∑ ∫ 13 有一个半径为R 的均匀球体，它的中心就是坐标的原点，已知球内任意一点的初 始温度只与该点到球心的距离有关，球的表面温度永远保持零度，确定球内任一点在 t>0 时的温度。解 由球对称性可知方程为 2 2 2 2 0 u u ur a r t r r ⎡ ⎤∂ ∂ ∂− + =⎢ ⎥∂ ∂ ∂⎣ ⎦ ，从而有 定解问题 28 2 2 2 2 0 ( , ) 0, (0, ) ( ,0) ( ) u u ur a r t r r u R t u t u r f r ⎧ ⎡ ⎤∂ ∂ ∂− + =⎪ ⎢ ⎥∂ ∂ ∂⎣ ⎦⎪⎪ =⎨⎪ =⎪⎪⎩ 有界 作变换 ( , ) ( , )v r t ru r t= 得 2 2 2 0 ( , ) 0, (0, ) ( ,0) ( ) v va t r v R t v t v r rf r ⎧∂ ∂− =⎪∂ ∂⎪⎪ =⎨⎪ =⎪⎪⎩ =0 应用分离变量法解得 2 1 0 ( , ) ( )sin sin n la t R n n nv r t e f d r R R π π πρ ρ ρ ρ ⎛ ⎞∞ −⎜ ⎟⎝ ⎠ = ⎡ ⎤⎛ ⎞⎢ ⎥= ⎜ ⎟⎢ ⎥⎝ ⎠⎣ ⎦ ∑ ∫ 14．在圆域 aρ < 上求解 4uΔ = − ，边界条件是 0 a u ρ = = 。 解 作变换 2 2( )u v x y= − + 则 v 满足Lapulace 方程，在极坐标系下用变量分离法解 之得 2 2u a ρ= − 15. 在圆域 aρ < 上求解 u xyΔ = − ，边界条件是 0 a u ρ = = 。 解 作变换 3 31 ( ) 12 u v x y xy= − + 则 v 满足Lapulace 方程，在极坐标系下用变量分离 法解之得 ( )2 2 21 sin 2 24 u aρ ρ ϕ= − 16. 在矩形城0 , 2 2 b bx a y≤ ≤ − ≤ ≤ 上求解 2uΔ = − ，且 u 在边界上的值为零。 解 作变换 ( )u v x x a= + − 后再用变量分离法解之得 [ ] [ ] ( ) 2 3 3 0 (2 1) sin (2 1)8( , ) ( ) (2 1) 1 2n ch n y a n y aau x y x x a n ch n b a π π π π ∞ = + −= − − + +⎡ ⎤⎣ ⎦∑ 17 求解 2 2 1 1 ( , ) 4 xx yy x y u u u x y + = + =⎧⎪⎨ =⎪⎩ 解 u v w= + ，其中 29 2 2 1 1 ( ) ( , ) 0 xx yy x y v v i v x y + = + =⎧⎪⎨ =⎪⎩ 2 2 1 0 ( ) ( , ) 4 xx yy x y w w ii w x y + = + =⎧⎪⎨ =⎪⎩ 意见(i)有解 2 2 1 4 4 x yv += − ，(ii)有解 4w = 。所以 2 2 15( , ) 4 4 x yu x y += + 18 求解 2 2 3 1 1 ( , ) 9 xx yy zz x y z u u u u x y + + = + + =⎧⎪⎨ =⎪⎩ 解 u v w= + ，其中 2 2 2 1 1 ( ) ( , ) 0 xx yy zz x y z v v v i v x y + + = + + =⎧⎪⎨ =⎪⎩ 2 2 2 1 0 ( ) ( , ) 9 xx yy zz x y z w w w ii w x y + + = + + =⎧⎪⎨ =⎪⎩ 意见(i)有解 2 2 2 1 6 6 x y zv + += − ，(ii)有解 9w = 。所以 2 2 2 53( , , ) 6 6 x y zu x y z + += + 19 求解 2 2 2 2 1 ( , ) 2 xx yy x y u u x y u x y + = ⎧ + = +⎪⎨ =⎪⎩ 解 u v w= + ，其中 2 2 2 2 1 ( ) ( , ) 0 xx yy x y v v x y i v x y + = ⎧ + = +⎪⎨ =⎪⎩ 2 2 1 0 ( ) ( , ) 2 xx yy x y w w ii w x y + = + =⎧⎪⎨ =⎪⎩ 意见(i)有解 ( )3 22 2 1 9 9 x y v += − ，(ii)有解 2w = 。所以 ( )3 22 2 17( , )