| 查看: 371 | 回复: 1 | |||
| 当前主题已经存档。 | |||
nickler荣誉版主 (知名作家)
|
[交流]
泛涵微分简介(转载）
|
||
|
考查下列泛函: 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. 上述思想可用于非常复杂的情况.比如说弯曲能决定的膜的平衡形状.下次 介绍膜的问题. |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
liu_dong88
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3021.7
- 红花: 1
- 帖子: 511
- 在线: 100.9小时
- 虫号: 90311
- 注册: 2005-09-01
- 专业: 结构陶瓷
2楼2006-04-30 08:52:17
25