This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
keynote:2011-lesson06 [2011/06/02 20:48] 11021038 [6.3 可动边界的变分问题] |
keynote:2011-lesson06 [2023/08/19 21:02] (current) |
||
---|---|---|---|
Line 1: | Line 1: | ||
+ | . | ||
+ | ====== 第六课 Calculus of Variations(变分法)====== | ||
+ | ====6.1 变分概念与变分法基本引理==== | ||
+ | ==6.1.1变分法引例== | ||
+ | 泛函分析是研究拓扑线性空间到拓扑线性空间之间满足各种拓扑和代数条件的映射的分支学科。 | ||
+ | |||
+ | 泛函的定义域是一个函数集,而值域是实数集或者实数集的一个子集,泛函就是从任意的向量空间到标量的映射。 | ||
+ | |||
+ | 泛函是从函数空间到数域的映射。 | ||
+ | |||
+ | 变分学研究的主要内容是泛函的极值。因此,凡是求泛函极值的问题,都可以称作变分问题。 | ||
+ | |||
+ | |||
+ | **引例1,降线问题 | ||
+ | ** | ||
+ | 如图所示,小球自由从O点到B点滑落(仅受重力影响),问:当O点与B点之间的曲线是何种轨迹时,小球的滑落时间最{短? | ||
+ | |||
+ | {{:keynote:例1_初始示例_.jpg|}} | ||
+ | |||
+ | **解: | ||
+ | ** | ||
+ | ①建立坐标系 | ||
+ | |||
+ | {{:keynote:例1_坐标系_.jpg|}} | ||
+ | |||
+ | ②设曲线函数 | ||
+ | |||
+ | y=y(x) | ||
+ | |||
+ | ③假设在图中N点的速度是v,分析从O到N的过程,根据能量守恒定律,可得: | ||
+ | |||
+ | mgh=1/2mV^2. | ||
+ | v=√2gy | ||
+ | |||
+ | ④设A为计算曲线弧长的起点,则有 | ||
+ | |||
+ | ds/dt=v=√2gy | ||
+ | |||
+ | ds=√(dx^2+dy^2)=√[dx^2(1+dy^2/dx^2)]=dx√(1+y'^2) | ||
+ | |||
+ | 故dt=ds/v=[√(1+y'^2)/√(2gy)]dx | ||
+ | |||
+ | ⑤从O到B进行积分,总时间 | ||
+ | |||
+ | T=∫0→t dt=∫0→a √(1+y^2)/√(2gy) dx | ||
+ | |||
+ | 观察上式可以看出,泛函T的值与函数y=y(x)相关,在这种情况下,函数T是无法直接求其极值的。 | ||
+ | |||
+ | 因此,降线问题最后可以被描述为:求光滑曲线γ : y = y(x), (0 ≤ x ≤ 1),在其满足条件y(0) = 0, y(a) = b下使泛函T[y(x)]取最小值。 | ||
+ | |||
+ | |||
+ | **引例2,短程线问题 | ||
+ | ** | ||
+ | 设A(x0,y0,z0)和B(x1,y1,z1)是曲面∑:φ(x,y,z)=0上的两点,求在曲面∑上过A和B的长度为最短的曲线。所求曲线称为短程线或侧地线。 | ||
+ | |||
+ | **解: | ||
+ | ** | ||
+ | 设经过两点的曲线Γ的方程为 | ||
+ | |||
+ | y=y(x) | ||
+ | |||
+ | z=z(x) x0≤x≤x1 | ||
+ | |||
+ | 则曲线Γ的长度为L=∫x0→x1 √(1+y'^2+z'^2)dx | ||
+ | |||
+ | 泛函L的值与Γ相关。即L与y = y(x), z = z(x)有关。记L = L[y(x), z(x)],短程线问题可以归结为求曲线 | ||
+ | |||
+ | Γ: | ||
+ | y=y(x) | ||
+ | |||
+ | z=z(x) x0≤x≤x1 | ||
+ | |||
+ | 在其满足边界条件y0 = y(x0),z0 = z(x0),y1 = y(x1),z1 = z(x1),及约束条件ϕ(x, y(x), z(x)) = 0下,使泛函L取得最小值。 | ||
+ | |||
+ | **引例3,等周问题 | ||
+ | ** | ||
+ | 求长为定值l的平面封闭曲线Γ,使其所围成的平面区域D的面积最大。 | ||
+ | |||
+ | 解: | ||
+ | |||
+ | 设曲线Γ的方程为: | ||
+ | |||
+ | x = x(t), | ||
+ | |||
+ | y = y(t),t0 ≤ t ≤ t1, | ||
+ | |||
+ | 由弧长公式知:l=∫t0→t1√(x'^2(t) + y'^2(t))dt | ||
+ | |||
+ | 由Green公式,Γ所围成的面积A为 | ||
+ | |||
+ | {{:keynote:例3公式.jpg|}} | ||
+ | |||
+ | 泛函A的值与x(t), y(t)有关。本例的问题归结为在上述条件下使得泛函A取最大值。 | ||
+ | |||
+ | **引例4,最小旋转面问题 | ||
+ | **在XOY平面内求一条边界固定的曲线,使其绕着横轴旋转所产生的空间曲面面积最小。 | ||
+ | |||
+ | **解: | ||
+ | ** | ||
+ | |||
+ | 设过点A(x0, y0), B(x1, y1)的曲线方程为y = y(x), x0 ≤ x ≤ x1, | ||
+ | |||
+ | 则旋转曲面的面积S为 | ||
+ | |||
+ | S = 2π∫x0→x1 y√(1 + y'^2)dx. | ||
+ | |||
+ | 泛函S的值与y=y(x)有关,记S=S[y(x)] 于是最小旋转曲面问题归结为:求曲线Y=y(s)(x0≤x≤x1)在其满足边界条件 | ||
+ | |||
+ | y(x0)=y0,y(x1)=y1 | ||
+ | |||
+ | 下面使泛函S取极小值 | ||
+ | |||
+ | |||
+ | **引例5,悬链形状问题 | ||
+ | ** | ||
+ | 求长度为l两段系于A和B两点,绝对柔软而不伸长的匀质链的形状。 | ||
+ | |||
+ | **解: | ||
+ | ** | ||
+ | 设链的方程y = y(x), (0 ≤ x ≤ x1),链的线密度ρ,按照位能原理,链的形状应该使链具有最小位能。 | ||
+ | |||
+ | 小弧段ds的位能为:dU = ρds ˙ gy = ρgy√(1 + y'^2)dx | ||
+ | |||
+ | 链的总位能为:U = ρg∫y√(1+y'^2)dx | ||
+ | |||
+ | 泛函U的值与曲线y=y(x)有关。记PU = U[y(x)],于是悬链问题归结为在上述条件下求泛函U的极小值。 | ||
+ | |||
+ | **综上可以看出共通点: | ||
+ | |||
+ | ①这些例子都是求泛函的极值问题。 | ||
+ | ②边界都是固定的,端点有不变的值。 | ||
+ | ③除边界条件外,有些再无其他条件,有些则另外有附加条件。有的条件是一般函数,有的条件也是泛函。其中,固定边界且无条件的泛函是最基本最简单的泛函。 | ||
+ | **\\ | ||
+ | |||
+ | ==6.1.2 变分概念与引理== | ||
+ | |||
+ | 变分问题即泛函的极值问题。 | ||
+ | |||
+ | 1.变分概念 | ||
+ | |||
+ | **定义1: | ||
+ | ** | ||
+ | 设 J[y(x)] 是定义在函数集合<jsm>Y = { y(x) } </jsm>上的泛函,称 y(x) 为 J[y(x)] 的宗量,Y 是 J[y(x)] 的定义域(或容许曲线类(簇))。\\ | ||
+ | 引例1中,<jsm>J[y(x)] = \int_0^1 y(x) dx</jsm>是定义在<jsm>Y = \left\{y(x), y(x) \in C[0,1]\right\}</jsm>上的**泛函**,y(x) 是该**泛函的宗量**,Y 为**容许曲线簇**。\\ | ||
+ | 通俗地理解,宗量即自变量,为了与函数定义中的自变量区分开来引入的概念,在讨论泛函的连续以及其他性质时需要此概念。宗量和容许曲线簇之间没有区别,都是指泛函的自变量。 | ||
+ | |||
+ | **定义2: | ||
+ | **泛函 J[y(x)] 的宗量与另一宗量<jsm>y_0(x)</jsm>之差,称为宗量 y(x) 的变分,记为**<jsm> δy = δy(x) = y(x) - y_0(x)</jsm>.** | ||
+ | |||
+ | **定义3: | ||
+ | **若<jsm>max|y(x)-y_0(x)|</jsm>很小,则称 y(x) 与<jsm>y_0(x)</jsm>具有零阶接近,称函数集合$\{ {y(x)} | {|y(x)-y_0(x)| < δ} \}$ 为 <jsm>y_0(x)</jsm> 的**零阶 δ- 邻域**。\\ | ||
+ | 设 k 为正整数,若<jsm>max\{ |y(x)-y_0(x)|, |y'(x)-y_0'(x)|,...,|y^{k}(x)-y_0^{k}(x)| \}</jsm>很小,则称 y(x) 与<jsm>y_0(x)</jsm>有k阶接近,称函数集合<jsmath>\{ y(x) | max[|y(x)-y_0(x)|,|y'(x)-y_0'(x)|,...,|y^{k}(x)-y^{k}_0(x)|] < δ \}</jsmath> | ||
+ | 为<jsm>y_0(x)</jsm>的** k 阶 δ- 邻域**。 | ||
+ | |||
+ | **定义4: | ||
+ | **如果 ∀ε > 0,使对<jsm>y_0(x)</jsm>的 k 阶 δ- 邻域中的任何y(x),都有<jsm>| J[y(x)]-J[y_0(x)]| < ε</jsm>,则称 J[y(x)] 是在<jsm>y_0(x)</jsm>处具有 k 阶接近度的连续泛函。 | ||
+ | |||
+ | **定义5: | ||
+ | **设 F(x, y, y') 是关于三个变元 x, y, y' 的二阶连续可微函数,x 任意固定,η(x) 是任意可微函数,ε 是无穷小参数,则 F(x, y, y') 的增量为<jsmath>ΔF = F(x, y+εη, y'+εη') - F(x, y, y')</jsmath>按Taylor公式,有<jsm>{ΔF = \frac{∂F}{∂y}εη + \frac{∂F}{∂y'}εη' + R_1}</jsm>,其中<jsm>R_1</jsm>是较 ε→0 时高阶的无穷小,称<jsm>{δF = \frac{∂F}{∂y}εη + \frac{∂F}{∂y'}εη'}</jsm>为函数 F(x, y, y') 的变分。\\ | ||
+ | <note>变分符号与微分符号可以互换,且<jsm>δF = \frac{∂F}{∂y}δy + \frac{∂F}{∂y'}δy'</jsm>.\\ | ||
+ | |||
+ | 并且满足<jsm>δ(F_1 + F_2) = δF_1 + δF_2</jsm>和<jsm>δF_1F_2 = F_1δF_2 + F_2δF_1</jsm> | ||
+ | </note> | ||
+ | 证明过程: | ||
+ | 令 F(x, y, y') = y,则 δy = εη;令 F(x, y, y')= y',则 δy = εη';因此有(δy)' = (εη)' = εη' = δy';令 F(x, y, y') 依次等于<jsm>{y^{(2)},y^{(3)},...}</jsm>,由此可得<jsm>{(δy)^{(k)} = δy^{(k)} }</jsm> | ||
+ | |||
+ | **定义6(泛函的变分第一种定义): | ||
+ | **如果 J[y(x)] 的增量 ΔJ = J[y + δy] - J[y] 可以表示为<jsmath>ΔJ=J[y,δy] + β(y,δy)max|δy|,</jsmath>其中,J[y,δy] 对于δy而言是线性的,且当 δy→0 时,β(y,δy)→0,则称 J[y,δy] 为泛函 J[y] 的变分,记为δJ,即 **δJ = J[y,δy].**\\ | ||
+ | |||
+ | **定义7(泛函的变分较弱的定义): | ||
+ | **考虑含参数 α 的一系列容许曲线 y(x)+ αδy,将 y(x),δy 任意固定,泛函 J[y+ αδy] 看成是参数 α 的函数,记为Φ(α) = J[y+ αδy]。如果<jsm>Φ'(0) = \frac{∂}{∂α}J[y+ αδy]|_{α=0}</jsm>存在,则称 Φ'(0) 为泛函 J[y] 的变分,仍记为 δJ,即<jsmath>δJ = Φ'(0) = \frac{∂}{∂α}J[y+ αδy]|_{α=0}.</jsmath> | ||
+ | <note>考察两种定义下的变分关系\\ | ||
+ | **定理1. | ||
+ | **对于泛函 J[y] ,若按定义 6 的变分存在,则在定义 7 意义下 J[y] 的变分也存在并且两者相等.\\ | ||
+ | 变分问题转换为泛函的一阶求导 | ||
+ | </note> | ||
+ | 证明. | ||
+ | 设 δJ = J[y,δy] + β(y,δy)max|δy|,则 | ||
+ | <jsmath> | ||
+ | Φ'(0) = \lim_{α\rightarrow 0}\frac{Φ(α)-Φ(0)}{α} | ||
+ | = \lim_{α\rightarrow 0}\frac{J[y + αδy]-J[y]}{α} | ||
+ | = \lim_{α\rightarrow 0}\frac{J[y, αδy] + β(y,αδy)max|αδy|}{α} | ||
+ | = \lim_{α\rightarrow 0}\{J[y, δy] \pm β(y,αδy)max|δy|\}. | ||
+ | </jsmath> | ||
+ | 当 α → 0 时,αδy → 0,从而,β(y,αδy) → 0,故 | ||
+ | <jsmath> | ||
+ | Φ'(0) = J[y, δy] = ΔJ. | ||
+ | </jsmath>证明完毕. | ||
+ | |||
+ | **定义8: | ||
+ | **设<jsm>y_0(x)</jsm>是泛函 J[y] 的容许曲线簇 Y 中的某一函数,若对 ∀y ∈ Y,都有<jsm>J[y(x)] ≤ J[y_0(x)](或J[y(x)] ≥ J[y_0(x)])</jsm>,则称泛函 J[y] 在<jsm>y_0(x)</jsm>处达到**极大(小)值(或绝对极大(小)值)**,并称<jsm>y_0(x)</jsm>为** J[y] 的极大(小)值曲线.**\\ | ||
+ | 若对于<jsm>y_0(x)</jsm>的零阶 δ- 邻域内的所有函数 y(x), 都有<jsm>J[y(x)] ≤ J[y_0(x)](或J[y(x)] ≥ J[y_0(x)]),</jsm>则称泛函 J[y] 在<jsm>y_0(x)</jsm>处达到**强极大(小)值.**\\ | ||
+ | 若对于<jsm>y_0(x)</jsm>的一阶 δ-邻域内的所有函数 y(x),都有<jsm>J[y(x)]≤J[y_0(x)](或J[y(x)]≥J[y_0(x)]),</jsm>则称泛函 J[y] 在<jsm>y_0(x)</jsm>处达到**弱极大(小)值**,此时<jsm>y_0(x)</jsm>是 J[y]的**极大(小)值曲线**。 | ||
+ | <note> | ||
+ | **推论1 | ||
+ | **强极值必是弱极值,但反之不成立.\\ | ||
+ | **推论2 | ||
+ | **绝对极值必是强极值.\\ | ||
+ | **推论3 | ||
+ | **泛函 J[y] 达到绝对极值的必要条件也是达到弱极值的必要条件,更是达到强极值的必要条件.\\ | ||
+ | </note> | ||
+ | <note> | ||
+ | **定理2. | ||
+ | **若具有变分的泛函 J[y(x)] 在<jsm>y=y_0(x)</jsm>达到极值,则<jsm>δJ|_{y=y_0(x)} = δJ[y_0(x)] = 0</jsm>\\ | ||
+ | 证明.<jsm>Φ(α) = J[y_0(x) + αδy]</jsm>在 α = 0 时达到极值,故 Φ'(0) = 0. 即<jsm>δJ[y_0(x)] = 0.</jsm> | ||
+ | </note> | ||
+ | |||
+ | **6.1.3 变分法基本引理** | ||
+ | |||
+ | **1.基本引理1. | ||
+ | ** | ||
+ | 设<jsm> Φ(x) ∈ C[x_0, x_1]</jsm>,如果对<jsm>∀η(x) ∈ C'[x_0, x_1], η(x_0) = η(x_1) = 0,</jsm>都有 | ||
+ | <jsmath> | ||
+ | \int^{x_1}_{x_0}Φ(x)η(x)dx = 0, | ||
+ | </jsmath> | ||
+ | 则在区间<jsm>[x_0, x_1]</jsm>上,Φ(x)≡ 0. | ||
+ | \\ | ||
+ | 证明:恒等于风格的问题使用反证法8-)\\ | ||
+ | 假设 Φ(x) 在<jsm>[x_0, x_1]</jsm>中的某点 ξ 处不等于0,不妨设Φ(ξ) > 0. 由 Φ(x) 的连续性知,必存在区间<jsm>[ξ_0, ξ_1]</jsm>,使 ξ ∈ <jsm>[ξ_0, ξ_1]</jsm> ⊂ <jsm>[x_0, x_1]</jsm>,并且当 x ∈ <jsm>[ξ_0, ξ_1]</jsm> 时,Φ(x) > 0. 此时η(x)的取值情况如下:\\ | ||
+ | *当<jsm>x_0 ≤ x ≤ ξ_0</jsm>时,η(x) = 0.\\ | ||
+ | *当<jsm>ξ_0 ≤ x ≤ ξ_1</jsm>时,η(x) = <jsm>(x - ξ_0)^2(x - ξ_1)^2.</jsm>\\ | ||
+ | *当<jsm>ξ_1 ≤ x ≤ x_0</jsm>时,η(x) = 0.\\ | ||
+ | 则<jsm>η(x_0) = η(x_1) = 0</jsm>,η'(x) 在<jsm>[x_0, x_1]</jsm>上连续,即 η(x) 适合引理1的条件,但是 | ||
+ | <jsmath> | ||
+ | \int^{x_1}_{x_0}Φ(x)η(x)dx = \int^{ξ_1}_{ξ_0}Φ(x)(x - ξ_0)^2(x - ξ_1)^2dx > 0, | ||
+ | </jsmath> | ||
+ | 与引理结论矛盾,所以假设不成立,故 Φ(x)≡ 0. | ||
+ | |||
+ | **2 基本引理2. | ||
+ | **设 D 为平面区域,D 的边界为 Γ,Φ(x,y) ∈ C(D), 若对任一在 D + Γ 上连续可微且在 Γ 上取零值的函数 η(x,y), 都有 | ||
+ | <jsmath> | ||
+ | \int\int_{D}Φ(x,y)η(x,y)dxdy = 0 , | ||
+ | </jsmath> | ||
+ | 则在 D 上 Φ(x,y) ≡ 0.\\ | ||
+ | 证明: 用反证法。假设 ∃(ξ,ξ) ∈ D,使 Φ(ξ,ξ) ≠ 0,不妨设 Φ(ξ,ξ) > 0. 由 Φ(x,y) 的连续性知,∃ρ > 0, 使 Φ(x,y) 在圆<jsm>S_ρ: (x - ξ)^2 + (y - ξ)^2 < ρ</jsm>内仍为正,且<jsm>S_ρ ⊂ D</jsm>。构造函数:\\ | ||
+ | *当<jsm>(x - ξ)^2(y - ξ)^2 ≥ ρ^2</jsm>时,η(x,y) = 0.\\ | ||
+ | *当<jsm>(x - ξ)^2(y - ξ)^2 < ρ^2</jsm>时,η(x,y) = <jsm>[(x - ξ)^2 + (y - ξ)^2 - ρ^2]^2</jsm>.\\ | ||
+ | 函数 η(x,y) ∈ <jsm>C^1(D)</jsm>, 且 η(x,y)<jsm>|_Γ</jsm> = 0. 但是:\\ | ||
+ | <jsmath> | ||
+ | \int\int_{D}Φ(x,y)η(x,y)dxdy = \int\int_{D}Φ(x,y)[(x - ξ)^2 + (y - ξ)^2 - ρ^2]^2dxdy | ||
+ | </jsmath> | ||
+ | 与引理结论矛盾,故假设不成立,有Φ(x,y) = 0.\\ | ||
+ | <note> | ||
+ | Section**6.1.2** edited by 冉清(11021038)2011/06/01 | ||
+ | </note> | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | ====6.2 固定边界的变分问题==== | ||
+ | ==6.2.1 泛函极值的必要条件与Euler方程== | ||
+ | * **泛函极值的必要条件** | ||
+ | |||
+ | 最简单的泛函的极值问题:\\ | ||
+ | $ J[y(x)]= \int_{x_0}^{x_1} {F(x,y,y')dx,} $ (5) \\ | ||
+ | 其中,$F∈C^2$,容许曲线$y(x)∈C^2[x_0,x_1]$,且满足边界条件: \\ | ||
+ | $y(x_0)=y_0, y(x_1)=y_1.$ (6) \\ | ||
+ | 极值的必要条件:使泛函$J[y]=\int_{x_0}^{x_1} {F(x,y,y')dx}$达到极值的必要条件,是y(x)满足Euler方程: \\ | ||
+ | $F_y- \frac{d}{dx}F_y=0.$ (7) \\ | ||
+ | 证明:构造以α为参数的容许曲线族$\overline {y}(\alpha )=y(x)+\alpha \delta y$,\\ | ||
+ | 其中,$\delta y$为宗量y(x)的变分,满足$\delta y|_{x=x_0}=\delta y|_{x=x_1}=0$.\\ | ||
+ | 所以,$\overline {y}|_{x=x_0}=y_0,\overline {y}|_{x=x_1}=y_1,\overline {y}\in C^2[x_0,x_1]$.\\ | ||
+ | 假设,当α=0时,$\overline {y}=y(x)$是使泛函(18)达到极值的曲线,将$\overline {y}$代入方程(18),得到\\ | ||
+ | $\Phi (\alpha )=J[y+\alpha \delta y]=\int_{x_0}^{x_1} {F(x,y+\alpha \delta y,y'+\alpha \delta y')dx.}$\\ | ||
+ | 由于$\Phi (\alpha )$在α=0时取得极值,由极值的必要条件,有\\ | ||
+ | $\Phi '(0)=\frac{\partial }{\partial \alpha }J[y+\alpha \delta y]|_{\alpha =0}=\int_{x_0}^{x_1} {(F_y\delta y+F_y\delta y')dx}=0.$ (8)\\ | ||
+ | 由分部积分法,并注意到$\delta y|_{x=x_0}=\delta y|_{x=x_1}=0$,有\\ | ||
+ | $\int_{x_0}^{x_1} {F_{y'}}\delta y'dx=\int_{x_0}^{x_1} {F_{y'}}d(\delta y)=F_{y'}\delta y|_{x_0}^{x_1}-\int_{x_0}^{x_1} {\delta y\frac{d}{dx}F_{y'}dx}=-\int_{x_0}^{x_1} {\delta y\frac{d}{dx}F_{y'}dx}$.\\ | ||
+ | 将上式代入方程(8),得到\\ | ||
+ | $\Phi '(0)=\delta J=\int_{x_0}^{x_1} {[F_y-\frac{d}{dx}F_{y'}]\delta ydx}=0.$\\ | ||
+ | 由基本引理1知,使泛函J[y]达到极值的函数y(x)必满足Euler方程7.\\ | ||
+ | 说明:注意到满足Euler方程7仅是泛函取到极值的必要条件,Euler方程7的解曲线仅仅是可能的极值曲线。但是,在实际问题中,往往可以事先确定泛函的极值是否存在。在确定了极值存在的情况下,Euler方程的解曲线就是极值曲线.\\ | ||
+ | |||
+ | |||
+ | |||
+ | ==6.2.2 含有多个未知函数的变分问题== | ||
+ | |||
+ | * **含两个未知函数的泛函极值的必要条件** | ||
+ | |||
+ | ☆问题: \\ | ||
+ | $J[y(x),z(x)]=\int_{x_0}^{x_1}{F(x,y,z,y',z')}dx,$ (11) \\ | ||
+ | 其中,F关于所含变量具有二阶连续领导,容许曲线$y(x)$,$z(x)$∈$C^2[x_0,x_1]$,且满足边界条件: \\ | ||
+ | $\left\{ {\begin{array}\ {y(x_0)=y_0, y(x_1)=y_1,} \\ {z(x_0)=z_0, z(x_1)=z_1}. \\ \end{array} } \right.$ \\ | ||
+ | ☆极值的必要条件:使泛函$J[y(x),z(x)]=\int_{x_0}^{x_1}{F(x,y,z,y',z')}dx$达到极值的必要条件是,y(x),z(x)满足Eular方程组,即: \\ | ||
+ | $\left\{ {\begin{array}\ {F_y-\frac{d}{dx}F_{y'}=0,} \\ {F_z-\frac{d}{dx}F_{z'}=0}. \\ \end{array} } \right.$ (12)\\ | ||
+ | |||
+ | ☆例1.求泛函$J[y,z]=\int_{0}^{\frac{π}{2}{(y'^2+z'^2+2yz)}}dx$满足边界条件: \\ | ||
+ | $\left\{ {\begin{array}\ {y(x_0)=0, y(\frac{π}{2})=1,} \\ {z(x_0)=0, z(\frac{π}{2})=-1}. \\ \end{array} } \right.$ \\ | ||
+ | 的极值曲线。 \\ | ||
+ | 解:$F=y'^2+z'^2+2yz$,故Eular方程组为:\\ | ||
+ | $\left\{ {\begin{array}\ {2z-\frac{d}{dx}(2y')=0,} \\ {2y-\frac{d}{dx}(2z')=0}. \\ \end{array} } \right.$\\ | ||
+ | 得: \\ | ||
+ | $\left\{ {\begin{array}\ {z-y^{(2)}=0,} \\ {y-z^{(2)}=0}. \\ \end{array} } \right.$\\ | ||
+ | 由此方程组消去z,得$y^{(4)}-y=0$,其通解为:\\ | ||
+ | $\left\{ {\begin{array}\ {y=c_1e^x+c_2e^{-x}+c_3cosx+c_4sinx,} \\ {z=y''=c_1e^x+c_2e^{-x}-c_3cosx-c_4sinx}. \\ \end{array} } \right.$\\ | ||
+ | 由边界条件可确定$c_1=c_2=c_3=0$,$c_4=1$,故,所求极值曲线为: \\ | ||
+ | $\left\{ {\begin{array}\ {y=sinx,} \\ {z=-sinx}. \\ \end{array} } \right.$\\ | ||
+ | |||
+ | | ||
+ | ★例2 求泛函$J[y(x),z(x)]=\int_{x_0}^{x_1}{F(y',z')dx,}$ 的极值曲线,其中假设 F_{y'y'}-F^2_{y'z'} ≠0 \\ | ||
+ | 解:因为F_y=F_z=0,故Euler方程组为 \\ | ||
+ | $\left\{ {\begin{array}\ {\frac{d}{dx}(F_{y'})=0,} \\ {\frac{d}{dx}(F_{z'})=0}. \\ \end{array} } \right.$\\ | ||
+ | 得:\\ | ||
+ | $\left\{ {\begin{array}\ {F_{y'y'}y"+F_{y'z'}z"=0,} \\ {F_{y'z'}y"+F_{z'z'}z"=0}. \\ \end{array} } \right.$\\ | ||
+ | 根据假设条件,该方程只有零解y"=0,z"=0,即 \\ | ||
+ | $\left\{ {\begin{array}\ {y=c_1x+c_2,} \\ {z=c_3x+c_4}. \\ \end{array} } \right.$\\ | ||
+ | 是所求的极值曲线。\\ | ||
+ | * **含n(n>2)个未知函数的泛函极值的必要条件**\\ | ||
+ | 泛函$J=\int_{x_0}^{x_1}{F(x,y_1,y_2,…,y_n,y'_1,y'_2,…y'_n)dx,}$在满足边界条件 \\ | ||
+ | $\left\{ {\begin{array}\ {y_i(x_0)=y_{i0},} \\ {y_i(x_1)=y_{i1}},i=1,2,…n, \\ \end{array} } \right.$\\ | ||
+ | 下,取得极值的必要条件是y_1(x),y_2(x),…,y_n(x),满足Euler方程组\\ | ||
+ | {{:keynote:xu3.png|}}\\ | ||
+ | 上述方程组通解中的常数,可以由所给的边界条件确定。\\ | ||
+ | 8-)<note> Section**6.2** edited by 徐磊(11021041)2011/05/31</note> | ||
+ | |||
+ | ====6.3 可动边界的变分问题==== | ||
+ | ==6.3.1 最简单的可动边界问题== | ||
+ | 固定边界的变分问题:在未知函数 y(x) 的边界点<jsm>A(x_0, y_0), B(x_1, y_1)</jsm>固定的情况下,求泛函<jsm>J[y(x)] = \int_{x_0}^{x_1}F(x,y,y')dx</jsm>的极值曲线。 | ||
+ | |||
+ | 当两边界点中的一个或两个是可以变动的 → 可动边界的问题。可动边界情况下的容许曲线簇是**包含**固定边界情况下的容许曲线簇。因此,当函数 y(x) 使可动边界的泛函达到极值,也能够使固定边界的泛函达到极值。此时的求解条件为:\\ | ||
+ | 1. Euler方程<jsm>F_y - \frac{d}{dx}F_{y'} = 0.</jsm>\\ | ||
+ | 2. 斜截条件.\\ | ||
+ | |||
+ | **斜截条件: | ||
+ | **y(x)有一个固定的左端点<jsm>A(x_0, y_0)</jsm>,右端点<jsm>B(x_1, y_1)</jsm>在某曲线 ω(x, y) = 0上移动。**考察<jsm>x_1</jsm>满足的特征.**\\ | ||
+ | 考察当右边界点<jsm>B(x_1, y_1)</jsm>移动到<jsm>\bar{B}(x_1+δx_1, y_1+δy_1)</jsm>的时,泛函 J 的增量为 | ||
+ | <jsmath> | ||
+ | ΔJ = \int^{x_1+δx_1}_{x_0}F(x, y+δy, y'+δy')dx - \int^{x_1}_{x_0}F(x, y, y')dx = \int^{x_1+δx_1}_{x_1}F(x, y+δy, y'+δy')dx + \int^{x_1}_{x_0}[F(x, y+δy, y'+δy')-F(x, y, y')]dx ---(19) | ||
+ | </jsmath> | ||
+ | (19)式中第一项,由积分中值定理及F的连续性,有: | ||
+ | <jsmath> | ||
+ | \int^{x_1+δx_1}_{x_1}F(x, y+δy, y'+δy')dx = F|_{x=x_1 + θδx_1}· δx_1 = F(x,y,y')|_{x=x_1}· δx_1 + ε_1δx_1, 0 < θ < 1; δx_1 → 0, ε_1 → 0. ---(20) | ||
+ | </jsmath> | ||
+ | (19)式中第二项,由Taylor展开有: | ||
+ | <jsmath> | ||
+ | \int^{x_1}_{x_0}[F(x, y+δy, y'+δy')-F(x, y, y')]dx = \int^{x_1}_{x_0}[F_y(x, y, y') + F_{y'}(x, y, y')δy']dx + R_1. ---(21) | ||
+ | </jsmath> | ||
+ | <jsm>R_1</jsm>是较 δy 和 δy' 的无穷小量。利用分部积分法以及 Euler方程,且 <jsm>δy|_{x=x_0} = 0,</jsm> 可得: | ||
+ | <jsmath> | ||
+ | \int^{x_1}_{x_0}[F_yδy + F_{y'}δy']dx = [F_{y'}δy]|^{x_1}_{x_0} + \int^{x_1}_{x_0}[F_y - \frac{d}{dx}F_{y'}]δydx = [F_{y'}δy]|_{x=x_1}. ---(22) | ||
+ | </jsmath> | ||
+ | {{:keynote:math01.png|}} | ||
+ | |||
+ | 如上图所示,<jsm>δy_1</jsm>是当边界点移到<jsm>(x_1+δx_1, y_1+δy_1)</jsm>位置时<jsm>y_1</jsm>的增量,<jsm>δy|_{x=x_1}</jsm>是当通过<jsm>(x_0, y_0)</jsm>和<jsm>(x_1, y_1)</jsm>两点的极值曲线移到通过<jsm>(x_0, y_0)</jsm>和<jsm>(x_1+δx_1, y_1+δy_1)</jsm>两点的极值曲线时,在点<jsm>x_1</jsm>处纵坐标的增量。\\ | ||
+ | 在图中, <jsm>BD = δy|_{x=x_1}, FB' = δy_1, EB'\approx{y'(x_1)δx_1}</jsm>,所以 BD = FB' - EB',即<jsm>δy|_{x=x_1} \approx{δy_1 - y'(x_1)δx_1}</jsm>。此时,近似等式与精确等式相差一个较<jsm>δx_1</jsm>高阶的无穷小量,将(22)式代入(21)式得到: | ||
+ | <jsmath> | ||
+ | \int^{x_1}_{x_0}[F(x, y+δy, y'+δy')-F(x, y, y')]dx = [F_{y'}δy]|_{x=x_1} \approx F_{y'}|_{x=x_1}(δy_1-y'(x_1)δx_1).---(23) | ||
+ | </jsmath> | ||
+ | 综合以上结论: | ||
+ | <note> | ||
+ | <jsmath> | ||
+ | \int^{x_1+δx_1}_{x_1}Fdx \approx F|_{x=x_1}δx_1. | ||
+ | </jsmath> | ||
+ | <jsmath> | ||
+ | \int^{x_1}_{x_0}[F(x, y+δy, y'+δy')-F(x, y, y')]dx = [F_{y'}δy]|_{x=x_1} \approx F_{y'}|_{x=x_1}(δy_1-y'(x_1)δx_1). | ||
+ | </jsmath> | ||
+ | 上述两个近似等式与精确等式相差较<jsm>δx_1</jsm>和<jsm>δy_1</jsm>的高阶无穷小。 | ||
+ | </note> | ||
+ | 代回到(19)式,得到: | ||
+ | <jsmath> | ||
+ | δJ = F|_{x=x_1}δx_1 + F_{y'}|_{x=x_1}(δy_1-y'(x_1)δx_1) = (F - y'F_{y'})|_{x=x_1}δx_1 + F_{y'}|_{x=x_1}δy_1 | ||
+ | </jsmath> | ||
+ | 所以**极值的基本必要条件** <jsm>δJ = 0</jsm>就成为: | ||
+ | <jsmath> | ||
+ | (F - y'F_{y'})|_{x=x_1}δx_1 + F_{y'}|_{x=x_1}δy_1 = 0. ---(24) | ||
+ | </jsmath>\\ | ||
+ | 如果<jsm>δx_1</jsm>和<jsm>δy_1</jsm>无关 → <jsm>(F - y'F_{y'})|_{x=x_1} = 0, </jsm><jsm>F_{y'}|_{x=x_1}=0.</jsm> ---(25)\\ | ||
+ | 如果<jsm>δx_1</jsm>和<jsm>δy_1</jsm>相关,设右端点沿某一曲线<jsm>C_1: w(x, y)=0</jsm>移动,此时有 <jsm>ω_{x_1}δx_1 + ω_{y_1}δy_1 = 0. ---(26)</jsm>\\ | ||
+ | 假设<jsm>ω_{y_1} ≠ 0,</jsm>, <jsm>δy_1 = -\frac{ω_{x_1}δx_1}{ω_{y_1}}.</jsm>\\ | ||
+ | 代入(24)式有,<jsm>[(F - y'F_{y'})ω_{y} - F_{y'}ω_{x}]|_{x=x_1} = 0. ---(27)</jsm>\\ | ||
+ | (27)时为**斜截条件(贯截条件)**。表明当<jsm> y = y(x), (x_0 ≤ x ≤ x_1)</jsm>为变动右端点的泛函极值曲线,则右端点满足此式子。\\ | ||
+ | |||
+ | <note> | ||
+ | 三种特殊情形: \\ | ||
+ | * 右端点 B 可以沿曲线<jsm>C_1: y = φ(x)</jsm>移动,此时有<jsm>ω = y - φ(x),ω_x = -φ'(x), ω_y = 1;</jsm> 斜截条件为:<jsm>[(φ'(x)-y')F_{y'} + F]|_{x=x_1}.</jsm>\\ | ||
+ | * 右端点 B 可以沿直线 <jsm>x = x_1</jsm>上移动,此时称 B 为自由端点,而<jsm>δx_1 = 0</jsm>. 斜截条件成为: <jsm>F_{y'} | {x=x_1} = 0.</jsm> 表示<jsm>F_{y'}</jsm>沿着 y = y(x)在 <jsm>x = x_1</jsm>处的值为零。此时的斜截条件也称为 **J[y] 的自然边界条件**.\\ | ||
+ | * 右边界点B在直线 <jsm>y = y_1</jsm>上移动,于是<jsm>δy_1 = 0.</jsm> 斜截条件为: <jsm>[F - y'F_{y'}] = 0.</jsm>\\ | ||
+ | 当左右边界均可动时可以推导出相应的斜截条件. | ||
+ | </note> | ||
+ | |||
+ | |||
+ | **例2. | ||
+ | **求联结平面上两条曲线 x + y + 1 = 0 与 xy = 1 的最短曲线.\\ | ||
+ | 解: \\ | ||
+ | 设所求平面曲线为 y = y(x),现求泛函 <jsm>J[y] = \int_{x_0}^{x_1}\sqrt{1 + {y'}^2}dx</jsm>的极值曲线,且,两个端点A和B都是变动的: \\ | ||
+ | <jsmath>C_1: y = ψ(x) = -1 - x, A(x_0,ψ(x_0))∈ C_1; </jsmath> | ||
+ | <jsmath>C_2: y = φ(x) = \frac{1}{x}, B(x_1, φ(x_1)) ∈ C_2</jsmath> | ||
+ | 易求得 J[y] 的 Euler 方程的解为 <jsm>y = c_1x + c_2.</jsm>,所以:\\ | ||
+ | <jsmath>y' = C_1, F = \sqrt{1 + {y'}^2}, F_{y'} = \sqrt{y'}{1+y'} = \frac{c_1}{1+c_1^2}</jsmath> | ||
+ | 联立斜截条件和所给边界条件 | ||
+ | <jsmath>[(φ'(x)-y')F_{y'} + F]|_{x=x_1} = 0</jsmath> | ||
+ | <jsmath>[(ψ'(x)-y')F_{y'} + F]|_{x=x_0} = 0</jsmath> | ||
+ | <jsmath>y(x_0) = ψ(x_0), y(x_1) = φ'(x_1)</jsmath> | ||
+ | 即 | ||
+ | {{:keynote:math02.png|}} | ||
+ | |||
+ | 解上述方程得到: <jsm>x_0 = -\frac{1}{2}, x_1 = 1, c_1 = 1, c_2 = 0.</jsm>\\ | ||
+ | 若x > 0,那么所求曲线为 y = x, (-1/2 ≤ x ≤1)。 | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | <note> | ||
+ | Section**6.3.1** edited by 冉清(11021038)2011/06/01 | ||
+ | </note> | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | <note important> 本节编撰作者(请大家在这里报到): | ||
+ | * [[zjuwlw@zju.edu.cn|王礼文]] (ID: 11021035), 编写了... | ||
+ | * [[xxxx@xxx.xxx|AuthorName2]] (ID: xxxxxxxxx), 编写了... | ||
+ | * [[ranqing@zjucadcg.cn|冉清]] (ID: 11021038), 编写了6.1.2变分概念和基本问题.6.3.1最简单的可动边界问题 | ||
+ | * [[cnxxxx@xxx.xxx|AuthorName3]] (ID: xxxxxxxxx), 编写了... | ||
+ | * [[kiriyame@163.com|张娜]] (ID: 11021040), 编写了6.1变分法引例 | ||
+ | * [[leixu@zju.edu.cn|徐磊]] (ID: 11021041), 编写了6.2固定边界的变分问题 | ||
+ | |||
+ | 浙江大学2008-2010版权所有,如需转载或引用,请与[[zhx@cad.zju.edu.cn | 作者联系]]。 | ||
+ | </note> |