Math 版 (精华区)

发信人: micheal (平凡的世界), 信区: Math
标  题: 泛涵微分简介(转载）
发信站: 哈工大紫丁香 (Fri Apr  7 13:21:05 2000), 转信

发信人: space (排骨教主), 信区: mathematics
标  题: 泛涵微分简介
发信站: BBS 曙光站 (Tue May 19 14:24:58 1998)


考查下列泛函:

S(q)=Integrate[L(q(t),q'(t),t), (t,t1,t2)]

S就是关于函数q(t)的泛函,因为积分的结果完全依赖于函数q(t)的形式. L是一个任意
的三元函数.估计你的情行比这应简单.

所谓泛函微分是指当q(t)的型式有小变化时S取值的改变. 假设h(t)是一个小函数,即
在定义域内值接近零.另外假设导数h'(t)也是小函数. 考查:

S(q+h)=Integrate[L(q(t)+h(t),q'(t)+h'(t),t), {t,t1,t2}]

由于h(t)(h(t)称为一个变分)和h'(t)很小,可对L作Talor展开:

L(q(t)+h(t),q'(t)+h'(t),t)=L(q(t),q'(t),t)+D[L,q]*h(t)+D[L,q']*h'(t)

所以:

\Delta S=S(q+h)-S(q)=Integrate[L(q(t)+h(t),q'(t)+h'(t),t)
         -L(q(t),q'(t),t), {t,t1,t2}]
        =Integrate[D[L,q]*h(t)+D[L,q'(t)]*h'(t), {t,t1,t2}]

估计你们是要求极值问题,这就要求\Delta S=0. =>

Integrate[D[L,q]*h(t)+D[L,q'(t)]*h'(t), {t,t1,t2}]=0

考虑到:
Integrate[D[L,q'(t)]*h'(t) dt,{t1,t2}]=Integrate[D[L,q'(t)] dh(t),{h(t1),h(t2)
]
=h(t2)*D[L,q'(t2)]-h(t1)*D[L,q'(t1)]-Integrate[h(t)*(d/dt)D[L,q'(t)] 
dt,{t1,t2].

(以上用的是分部积分法)

由于问题中积分边界上一般有限制: 边界固定. 所以: h(t1)=h(t2)=0.于是:

Integrate[D[L,q'(t)]*h'(t) dt,{t1,t2}]
=-Integrate[h(t)*(d/dt)D[L,q'(t)] dt,{t1,t2]

Hence:

Integrate[D[L,q]*h(t)+D[L,q'(t)]*h'(t), {t,t1,t2}]=0=>

Integrate[D[L,q]*h(t)-h(t)*(d/dt)D[L,q'(t)] dt, {t,t1,t2}]=0=>
由于问题中积分边界上一般有限制: 边界固定. 所以: h(t1)=h(t2)=0.于是:

Integrate[D[L,q'(t)]*h'(t) dt,{t1,t2}]
=-Integrate[h(t)*(d/dt)D[L,q'(t)] dt,{t1,t2]

Hence:

Integrate[D[L,q]*h(t)+D[L,q'(t)]*h'(t), {t,t1,t2}]=0=>

Integrate[D[L,q]*h(t)-h(t)*(d/dt)D[L,q'(t)] dt, {t,t1,t2}]=0=>

Integrate[h(t)*(D[L,q]-(d/dt)D[L,q'(t)]) dt, {t,t1,t2}]=0

由于h(t)是任意小函数,上式成立必要求

D[L,q]-(d/dt)D[L,q'(t)]=0

这是关于t的常微分方程. 由于L是给定的函数,该方程可求出函数q(t).即:

使泛函S(q)取极值的函数形式q(t).

上述方程在物理中称为: Lagrange equation.

上述思想可用于非常复杂的情况.比如说弯曲能决定的膜的平衡形状.下次
介绍膜的问题.



--
You own me a bottle of beer!                 
                      
                    

※ 来源:．BBS 曙光站 bbs.ncic.ac.cn．[IP: 159.226.41.99]※ 修改:·a24 於 Sep 16 14:10:23 修改本文·[FROM: bbs.ndc.neu.edu]

--
☆ 来源:．哈工大紫丁香 bbs.hit.edu.cn．[FROM: sillystone.bbs@smth.]