当前位置:
文档之家› 3[1].变分法与Hamilton原理
3[1].变分法与Hamilton原理
;
泛函Φ 的自变量为
定义域 ,自变量的变分
,
数的无穷小改变。
,即函
小参量法的定义为
,其中 为任意无穷小量, 为任意连续有界函数。类似于数
学分析中的 ‐ 语言,这是严格的数学定义。下面的叙述我们不追求数学上的严格性。
3. 泛函的变分
函数的微分
,,,
,
,,
,,, .
泛函的变分 Φ
Φ
Φ.
1 / 40
例Φ
解出 , ,得
⁄
0.562551 0.422487 π θ θ 0.211530 π θ θ 0.0604328 π θ θ
⁄
0.0173699 π θ θ 0.00202953 π θ θ
0.000561428 π θ θ
⁄
0.00164894 π θ θ
1 绳长不变,得约束条件
0
1
由于有约束, 不独立。引进拉氏乘子,得 0
10 / 40
1 1 1 1
1 1
1 1
1 1
1
1
1
1
01
0 1
可以直接求解方程,也可以利用“广义能量积分”,
1 1
1
,
1 1
1
ln cosh
3 个待定常数(2 个积分常数,1 个来自拉氏乘子)由代数方程
,
,
1
1 sinh
2 cosh
点为原点,设 0, 0。 , 两点之间用曲线
连接,一个质点被束缚在曲线上运
动,在重力作用下自由下降,初速为 0。什么样的曲线形状可以使质点从 到 所花的时间最少?
解 机械能守恒, 通过这段弧所需的时间为
2 ;弧长
1,
从 到 所花的时间为
1 2
1 2
取极值的必要条件为
,,
0.
5 / 40
注意上式无需对端点成立,因为端点处
13 / 40
在微分方程中令
0
0
0
0
/2,得隐藏的另外一个边界条件,
2
0
2
2
0 2
其几何意义是,上表面正中间的点切面和水平面夹角为 0。
求解微分方程后,代入约束条件
cos
,可确定参数 λ 。
图 1 用 RUNGE‐KUTA 法做数值计算,求出的水银滴在水平玻璃上的形状。计算中取 水银与玻璃的接触角为 ,水银的质量为 0.1 克。左图中横轴 的单位是度,竖轴 的形状,单位是毫米。
1
2
1 1
21
1
|
|
另一种计算方法:
12
Φ
12
1
1
2
1
2
1
3 / 40
2 例S
1 √1
,,
S
1
2
121Βιβλιοθήκη 216. 变分的运算规则
与微分法则类似,可以证明,
ΦΦ
Φ
ΦΦ
Φ Φ Φ Φ,
Φ (线性)
Φ
Φ Φ ΦΦ
Φ
Φ
Φ
Φ Φ,
Φ , ,Φ
7. 变分可以与微分、积分交换次序
按定义,
Φ Φ
lim
lim
lim
, , Δ, , lim
Δ
, ,, ,
,, ,
泛函的导数定义为
Φ Φ
2 / 40
按习惯, 一般写成 。
例Φ
,Φ
Φ
Φ
,
Φ Φ
0
这里δ~ ,由 Lagrange 引理,
例Φ
,
Φ
例 Φ, Φ 0
Φ
其中 是阶跃函数(Heaviside function)。
例Φ
12
Φ
1
2
2 1
1
1 1 12
上面用了 函数的性质
则水银的表面可以用广义坐标
描述,其中 是表面上的点 P 到原点的距离, 是 P 点与原
点的连线对桌面的夹角, 0, /2 。
按虚功原理,水银的形状应该使得势能最低,同时还必须满足体积为 的约束条件,
sin
2
势能分为重力势能和表面能两部分,
2
cos
3
cos
。重力势能为
sin
sin
1
2
sin cos
4
总表面能为液气表面积 和液固表面积 的表面能之和
lim
这一步交换求和顺序不严格,对病态函数可能不成立
因此对虚位移,
,在分析力学中称为等时变分或简单变分。
分析力学中还有其它方式定义的变分,一般说来等价于 与平移的混合运算,此时Δ
Δ。
对于积分,
lim
Δ
lim
Δ lim
Δ
lim
Δ 这一步交换求和顺序不严格,对某些病态函数可能不成立
8. 泛函的高阶变分
4 / 40
1
。把第一式代入第二式,
1
sinh
sinh
.
,
cosh
;
cosh
,
,
11 / 40
(2) 水银滴的表面形状 一滴水银静止于水平桌面上,求它的表面形状。
设水银的密度为 ,体积为 ;水银与空气之间的表面张力系数为 ,水银与桌面间的表面 张力系数为 ,桌面与空气之间的表面张力系数为 。
由对称性,水银的表面旋转对称。取坐标原点为桌面上水银滴的中心点,桌面为 ‐ 平面,
0 Euler 方程
,
0,应视为独立的变分,得自然边界条
0,
0
(3) EULER方程的首次积分
①
, ,“广义能量”
constant
②
, ,“广义动量” constant
③
, ,得代数方程 0,可直接解出 。
7 / 40
例 最速降线的另一种解法(沈 p219)设
,
有“广义能量积分”,……
(4) 多变量泛函的固定边界驻值
2 sinh 2
确定。 注:将坐标原点平移,总可以把悬链线方程写成
cosh 。
沈、金教材上由于没有严格地按条件极值来处理,得到的是特解。如果以绳子中心点即最
低点为坐标原点,则绳子的形状为
cosh 1 ,显然没有包含于课本给出的解中。
此问题也可以用牛顿力学来解。设张力为 ,水平方向力的平衡:
0
1
;垂直方向力的平衡:
11. 固定边界的泛函极值
在历史上,对最速降线问题(the brachistochrone problem)的分析导致了变分法(calculus of variation)的发明。(Fermat 原理比 Bernoulli 的最速降线更早了几十年)
(1) 最速降线
例 (J. Bernoulli, 1696 年)垂直平面上有两个固定的点 , 。以水平方向为 轴,向下方向为 轴,
sin ,
2 22
sin
arcsin
1
0,又 0,所以只能取正号。曲线方程为
1
arcsin
1
这个曲线方程可以改成参数形式,arcsin
sin ,取
,
是摆线方程。
sin 1 cos
这里我们只验证了驻值条件,还需进一步验证是否有二次变分 0,
6 / 40
1
21
1
1
21
确实是极小值。
另一种解法是设
,见沈 p219.
左边
12
1 23
1 12
右边 0
设 0时,
,
0,1
,
1
14 / 40
代入上式得
1
1
0
考虑到左右对称,令
1
2
,
0.
2
1 2!
代入微分方程,并将微分方程左边按
1
⁄1
3!
4!
的级数展开,得
2
4
,
,
3
保留到 项,再由
23 3 3
13
,
4 12
5
103 16 48
,
45
5
2
2 2
2 2
2
2,
3
cos
0.0001
例Φ
sin Φ
0.337979.
1
泛函可类比于多变量函数,
,,, ~ ,,, ,
Φ~
的定义域 .
这里的 相当于前面多变量函数的自变量 ,
,
;泛函Φ 是以不可数无穷多个变
量,
, 作为自变量的函数。
例 复合函数可以看成是一族泛函,
,,
,其中 是参数。
2. 变分
泛函可类比于多变量函数,
多变量函数 , , , ,自变量的微分为 , , , ,
Φ ,, ,
,, , ,, ,
,
,,
,
,
Euler 方程为
不指定边值时,自然边界条件为
0,
1,2, , .
(5) 多重积分型泛函的固定边界驻值
Φ
, , ,, ,
0
,
, 边界值
Euler 方程
已指定。
0
注意左边第二项 是对整个括号中式子求偏导数,
。
不指定边值,可得自然边界条件
0
,,
(6) 含高阶导数泛函的固定边界驻值 以 2 阶为例计算。……
07-12-06, 09:55
第 3 章 变分法与HAMILTON原理
一、 泛函与变分
1. 泛函
普通函数 :
,
泛函是普通函数概念的推广,Φ: 集合
自变量的集合通常取为函数,Φ
。
例
2 是函数;Φ
1 是泛函,并且有
Φ sin sin 1 0.84, Φ exp e 2.718, Φ 3.