耗散系统的哈密顿原理及其应用

耗散系统的哈密顿原理及其应用
耗散系统的哈密顿原理及其应用

题目耗散系统的哈密顿原理及其应用

学生姓名闫俊杰学号 1210014049 所在学院物理与电信工程学院

专业班级物理学1202 班

指导教师王剑华 __ __ 完成地点陕西理工学院

2016年 5 月 20 日

耗散系统的哈密顿原理及其应用

作者:闫俊杰

(陕西理工学院物理与电信工程学院物理专业1202班级,陕西汉中 723000)

指导老师:王剑华

[摘要]哈密顿原理不仅是分析力学的基础,而且在量子物理中也有重要的应用。本文首先通过定义耗散函数,给出了有耗散系统的哈密顿原理。然后利用有耗散系统的哈密顿原理推导出了有耗散系统的拉格朗日方程。最后通过与受线性粘滞阻尼作用的力学系统类比,把有耗散系统的拉格朗日方程应用到有耗散的电路系统,给出了基尔霍夫定律。

[关键词] 耗散系统;耗散函数;哈密顿原理;广义耗散力;基尔霍夫定律.

引言

哈密顿原理,于1834年首次发表,从而完成了从莫泊丢开始的尝试,该原理的数学形式简洁,内容广泛,将动力学问题转化为数学的一般体系的一部分,更深刻地揭示了客观事物之间的紧密联系,在物理学中有极高的地位,它不仅可以看作是力学的最高原理,甚至可以看作是整个物理学的最高原理,利用广义坐标并定义哈密顿函数就可以从哈密顿原理推导出哈密顿正则方程,致各种动力学定律都可以从这一变分式推出.所以说它是成为牛顿之后力学发展一个最大飞跃[1]。伴随经典力学从自由到非完整约束系统发展历程的是学术讨论甚至是争论,有些是对错之争,有些只是观点不同,问题的焦点主要是经典力学理论怎样从完整系统推广到非完整尤其是非线性非完整系统,这是经典力学后牛顿发展的主要标志,也是物理学近、现代发展的一块重要里程碑。物理学家考察具体的物理问题,常常以最小作用量原理为出发点,通过变分运算而导出物质系统的运动方程以表示其运动规律。牛顿质点运动方程可通过这样的推导过程得出;甚至如广义相对论里的爱因斯坦场方程也可如此导出,其实,所有运动方程就是拉格朗日方程或哈密顿正则方程在其拉氏量或哈氏量对于不同物质系统取不同形式时的具体体现。亦是广义相对论的一个重大结论;而如前文已提及的,引力场方程的导出可采用引力场之作用量的变分运算。从光线路径和质点运动到四维弯曲时空中的短程线方程乃至引力场运动方式等等实例可见,无论是几何间题,还是物理问题,都可凭籍变分法去圆满地解决;特别是最小作用量原理及其在物理学各领域的成功应用,正就是利用变分法这一几何方法、亦乃经济原则去解析各物质系统之运动规律的丰盈成果。例如牛顿力学描绘了天体运动的出色图景,拉格朗日一哈密顿理论描绘此图景也毫不逊色。当然,在经典力学范畴里,这两个理论体系本来是等价的,只是着眼点有所不同:对于行星运动,牛顿着眼于行星与恒星之间的万有引力;哈密顿等人着眼于行星运行时的能量守恒正因为后者着眼于能量状况,遂使哈密顿原理的适用性得以超出经典力学范畴。耗散系统就是指一个不断地与外界交换能量的系统这样的新结构就是耗散结构[2]。拉格朗日方程是力学理论中的基本方程,它在力学系统中得到广泛的应用。它的表述形式不再是直观的矢量形式,而是抽象的数学分析,即分析力学。分析力学所注重的不是力和加速度,而是具有更广泛意义的能量,同时又扩大了坐标的概念。所以拉格朗日方程更适用于处理复杂力学系统的问题,而且可以更进一

步,应用到物理学和技术科学的其他领域[3]

。目前拉格朗日方程已经被应用到许多领域中,

特别是电学中的应用。虽然拉格朗日方程在力学系统中得到广泛的应用,对力学问题提供了

新的、简洁的求解方法。但是在某些领域中还有它的不足,特别是在耗散系统中,一般的拉

格朗日方程解决的都是理想约束的问题,即忽略阻力、摩擦、黏滞力形、变等,或者是把阻

尼力当成主动力处理[4]。但是自然界中大多数物理过程都是有阻尼的物理过程,往往使问题

不能很好的解决,或变得更复杂。如果把耗散函数引入拉格朗日方程中,这个不足就会迎刃

而解,把阻尼力的问题用耗散函数来解决,解决有耗散系统的复杂问题时更加简单。

哈密顿原理它涉及了一系列基础问题。中国学者为推动这一学科的发展作出了重要的贡

献[5-9]。牛青萍发表了重要论文《经典力学基本微分原理与不完整力学组的运动方程》。这是

我国第一篇非完整力学的研究论文,赵跃宇建立了一类新型的积分变分原理,该原理比较完

满地解决了变分原理的推广问题。保加利亚科学院院士И.Цнов(Tzenoff)建立了一类新

型的运动微分方程,具有简单统一且便于应用的优点。对于变质量力学系统,梅凤翔就

Hamilton 原理[10],陈立群就非惯性系[11],张毅就相对运动[12],方建会和李元成就相对论力学

的速度空间[13],自二十世纪70年代以来,以DNA 为背景的超细长弹性杆力学受到关注,用

Kirchhoff 理论进行力学建模。基于此,刘延柱、陈立群等将经典的分析力学理论移植到弹

性杆静力学[14,15],人类对于自然界的认识永无止境,任何一个学科的发展都在伴随着时代的

步伐与时俱进。

本文主要研究了耗散函数,广义耗散力和准拉氏函数,耗散系统的哈密顿原理及其拉格

朗日方程,有耗散系统中哈密顿原理的应用。 1 耗散函数

设质点不仅受有势力和非有势力作用,还受粘滞阻尼的作用,粘滞阻尼是作用在质点上

的线性阻力。由于这种阻力使机械能耗散,所以他又被称为耗散力。下面,先定义阻尼力,再

由阻尼力的虚功给出耗散函数。设作用在任一质点i M 上的线性阻力为i R

i i i V C R -= (1.1)

其中V i 是质点的运动速度,阻力系数C i 为常数。质点上的线性阻力任意虚位移中所作虚功

的和为

11n n R i i i i i

i i W R r C V r δδδ===?=?∑∑∑ (1.2)

式中

∑∑∑===??=???=???=N 1k 11k k i k N k k i k N

k i i q q V q q r q q r r δδ&& (1.3) 将它代回式(1.2)得

2111112N N

N N i R i i k i i k i i k i k V W CV q CV q q δδδ====???=-?=- ????∑∑∑∑∑& (1.4)

在有耗散存在的力学系统中,可以把耗散函数Ψ表示成广义速度的函数 ∑=??==N

k k k i i i q q r r V 1&& (1.5) 令21

12n i i i C V ψ==∑称为耗散函数,有 11i 1111111()()222n n S S S S i i i i i i k j kj k j i i j k j k k

r r C VV C q q q q q q ψγ======??==?=??∑∑∑∑∑∑&&&& (1.6) 式中j

i k i n j i

kj q r q r C ?????=∑=1γ是常数或广义标的函数,称为广义阻力系数。式(1.6)表明耗散函数是广义速度的二次齐次函数。

2 广义耗散力和准拉氏函数

约束系统动力学的基本问题[16]和变分原理目前的研究主要出现在以下几个方面:

(1)对分析力学基本概念和基本问题在更高层面上的进一步探讨与扩充。这项研究,不但

可以巩固分析力学的理论根基,并且可以澄清一些纷争。

(2)对于新的力学系统,建立相应的力学变分原理。历史和现代的如变质量系统、单面约

束系统、连续系统、Birkhoff 系统、广义力学系统、广义非完整系统以及基于Kirchhoff

动力学比拟的超细常弹性杆等。

(3)用近代微分几何建立各种变分原理,如射丛几何的、辛几何的等。

(4)变分原理的应用,如应用于数值计算等等。

在有耗散存在的系统中,广义耗散力是一个非常重要的力学量。下面来进行计算,把

(1.6)代入(1.4)可以得到

s

k 1R k k W q q ψδδ=?=-?∑∑

& (2.1) 再令RK Q 为对应于广义坐标k q 的广义耗散力,则 Rk k

Q q ψ?=-

?& (2.2) 于是,有耗散存在的准拉氏函数为 R L T V W L W '=-+∑=+∑ (2.3)

利用这个准拉氏函数,我们可以给出有耗散存在的哈密顿原理。

3 耗散系统的哈密顿原理及其拉格朗日方程

物理学中最基本的原理是哈密顿原理。著名量子物理学家狄拉克曾十分推重这条原理

[18]

。因为它具有物理学理论的同一性。物理学中的一些重要原理与它们的核心方程,都能从

哈密顿原理出发、借拉格朗日分析力学的方法导出;哈密顿原理在分析力学中有着十分重要

的地位[19]。它借助于变分原理对质点运动情况作出精确描述。这一原理不仅适用于有限自由

度的点系,也可用于无限自由度的点系[20]。它在相对论中,场论中被广泛采用。在有耗散存在

的的系统中,可以考虑把哈密顿原理写成如下形式 210t t s L dt δδ

'==? (3.1)

把(2.3)代入(3.1)可以得到 2

1()0t R t s L W dt δδ=+∑=? (3.2) 即 211d δδd 0d s

t R t

L L q q W t q t q αααααδ=??????+-∑=?? ???????∑?& (3.3) 经过进一步推导,得 2

21111d d δ0d t s s t R t t L L L W q t q q t q q αααααααδ==?????????-+δ+=?? ? ??????????∑∑?&& (3.4) 把(2.1)代入(3.4),并考虑

[][]120,t t q q ααδδ== ()1,2,...s α= (3.5)

可得

2

11d δδd =0d s

t t L L S q q q t q t q q αααααααψδδ=???????=--?? ????????∑?&& (3.7) 哈密顿原理指出,对于质点系真实运动,0=S δ,即上式中积分为零。由于区间是任意的,

故知被积函数为零。又因为αδq 彼此独立,所以αδq 的S 个系数应全为零,于是得到方程。

即 0,(1,2,)d L L s dt q q q ααα

ψα???-+==???L && (3.8) 这就是含耗散函数的拉格朗日方程。

哈密顿原理处理简单运动没有什么优越性,但是用它来处理多自由度复杂问题却比较方

便。天体力学中摄动法就是应用哈密顿方程一个典型事例[21]

。哈密顿正则方程虽然是经典力

学的方程,但是却能被推广应用到微观和高速领域,在理论力学尤其是分析力学领域有着重

要应用。 4 有耗散系统中哈密顿原理的应用

一维弹性振子在有阻尼力和谐迫力时的振动方程为

cos mx bx Kx F t ++=Ω&&& (4.1)

其中m 是振子的质量,k 是弹簧的弹性系数,b 是阻尼因数,t F Ωcos 为谐迫力。

由电感L,电阻1R ,电容C 和电源E 所组成的串联电路, 若某时刻电容器上的电荷为q, 则

电路的电流 dt

dq I =

(4.2) 电感上的电压降为 22L dt

q d L dt dI = (4.3) 电阻上的电压降为 dt dq R I R 1

1= (4.4) 电容上的电压降为

c

q ,又设电源电动势为t cos ΩE E 00=,由基尔霍夫定律得电路的微分方程为 t cos ΩE c

q dt dq R dt q d L 00122=++ (4.5) 比较方程(3.1)和方程(3.2)可得如下的结论: 电荷与位移, 电感与质量, 电阻与阻尼

因数,电容的倒数与弹簧的弹性系数,电动势与机械力。利用同样的方法, 可以找到很多相

似的对应关系。

若取电荷q 为广义坐标, 电流q &即为广义速度, 电动势为广义力, 系统的动能、势能和

耗散函数分别为

221Lq T = 221q C V = 2112

R q ψ=& 以此代入有耗散力的拉格朗日方程 , 即得:

Q q

q T q T dt d =??+??-??&&ψ)( (4.6) 对于复杂的RLC 电路时,设有N 条含源组成的复杂RLC 电路,每条支路的电荷为

αq ()S ,2,1=α,这样,只有α-n 个独立网孔,选取网孔的电荷为广义坐标αq ()S ,2,1=α,就所有面电路,其任一之路的电荷都可以用一网孔电流或两网孔电流的和差表示。 各支路所含的电阻,电感,电容及电源电动势分别用()t E C L R i i i i ,,,表示,各网孔的电流、自电阻、自电感、自电容及自电源电动势分别用()t E L R Q jj jj jj jj ,,,表示,各网孔的互电阻、互电感、互电容分别用ji ji ji C L R ,,表示,且设所用的网孔的参考电流方向均为顺时针方向。 由以上对比结论可知,与LRC 电路相对应的力学系统的拉格朗日函数为[22] ()∑∑∑∑≠=≠=--+=s jk

k j s jj j k j s jk j s jj C q q C q q q L q L L βαβααβαβαα,12)(,21222121&&& (4.9) 耗散函数为 ()

k j s jk s j jj q q R q R &&&∑∑≠=+=βαβααψ,122121 (4.10) 其它非保守广义力为

)(t E Q jj jj = (4.11)

则拉格朗日方程为

,)(jj j

j j Q q q L q L dt d =??+??-??&&ψ (s j ??=,2,1) (4.12) 可得j 个拉格朗日方程

)()()()(t E C q q L q L C q q R q L s

j k k jj jk

k k s j k k jk k s j k k jk jj j

j jj j jj ∑∑∑≠≠≠=+++++&&&&&&,(s k j ??=,2,1, ) (4.13) 把耗散函数引入拉格朗日方程中,把阻尼力的问题用耗散函数来解决,这样就得到的了有耗散函数的拉格朗日方程。使拉格朗日方程的使用范围增加,解决有耗散系统的复杂问题时更加简单。耗散函数的拉格朗日方程可以更好的解决有阻尼力的耗散系统中,而且可以通过类比的方法更好的应用到电学中。虽然力学与电学是物理学中两个不同的课程,但是将分析力学中的广义坐标,广义力,动能,势能和耗散函数等概念引入到电学中,就可以用拉格朗日方程来解决RLC 电路中的问题。从力与能的角度出发,在RLC 电路体系中运用拉格朗日方程,得到了电学中用基尔霍夫定律和欧姆定律分析RLC 电路过程相一致的结果。这是分析不同研究对象之间共同规律的思维方法。更重要的是,两种不同的物理系统可用相同形式的拉格朗日函来描述,这种结构上的类似意味着研究其中一种系统的方法和结论也适用于另一种系统。这就远远扩大了拉格朗日方程的应用范围,为拉格朗日方程应用到其它更多的领域做好了铺垫。

5 结论

分析力学有微分、积分两种型式。拉氏方程和哈氏正则方程等属微分型至于积分型式,乃从对最小作用量原理的探讨起始。 变分法的发明使分析力学的建立和扩展有了简便的数

学工具。若一个系统在其物理过程中,系统能量随时间损耗,那么该系统就称耗散系统.由于耗散系统更接近实际,所以一直是物理学的一个研究课题;特别是近年来,这方面的研究(经典和量子两个方面)日趋活跃,并显示出了广阔的前景行了修正,比较成功地解决了一些耗散系统的问题.

长期以来,人们对对电力系统的暂态稳定性进行了大量的研究。由于电力系统是电力系统是一个能量产生、传输和消耗的复杂系统,而广义哈密顿理论其以能量理论为基础,鲜明其能量平衡在整个系统的安全运行中起至关重要的作用,将能量函数作为函数进行控制设计是很自然的。广义哈密顿控制系统不但为电力系统提供了合适的数学描述同时提供了一个新的控制途径。近年来哈密顿函数方法在电力系统分析和控制中受到越来越多的应用.该方法能充分利用系统内在的结构性质而且其结构简单,易于实现.所以耗散函数的哈密顿原理及其应用在电力系统分析和控制领域具有广阔的前景,吸引着越来越多研究人员去研究探索。

参考文献

[1]清华大学.理论力学[M].北京:高等教育出版社,1995.10,12(6):71-92.

[2]王振发.含耗散函数的拉格朗日方程及耗散函数的物理意义[J].山东轻工业学院学报,2001,15(2):15-16.

[3]狄拉克.物理学的方向[M].郭应焕译.北京:科学出版社,1981,8(5):121-192.

[4]王剑华,李康, 王亚辉.理论力学[M].西安:陕西科学技术出版社,2009.159-209.

[5]梅凤翔.非完整系统力学基础[M].北京:北京工业学院出版社,1985,6(2):275-286.

[6]陈滨.分析动力学[M].北京:北京大学出版社,1987,12(6):36-46.

[7]梅凤翔.高等分析力学[M].北京:北京理工大学出版社,1991,21(4):143-154.

[8]梅凤翔.关于变质量非完整力学系统的Hamilton原理[J].北京工业学院学报,1984,4:1-15.

[9]陈立群.非惯性系中变质量力学系统的变分原理及应用[ J].鞍山钢铁学院学报,1993,16(1):8-14.

[10]张毅.变质量力学系统相对运动的变分原理及其应用[J].山东科学,1995,8(4):6-11.

[11]方建会,李元成.变质量系统相对论力学速度空间中的变分原理[J].大学物理,1994,13(5):19-20.

[12]Liu Yanzhu,XueYun.Formulation of Kirchhoff rodbased on quasi-coordinates[J].Technische

mechanik,2004,24(3-4):206-210.

[13]薛绍,刘延柱,陈立群.超细长弹性杆的分析力学问题[J].力学学报,2005,37(4):485-493.

[14]梅凤翔.分析力学的近代发展[J].力学与实践,1987,9(1):10-15.

[15]罗绍凯,张永发,等.约束系统动力学研究进展[M].北京:科学出版社,2007,9(1):18-24.

[16]袁季修. 电力系统安全稳定控制[M].北京:中国电力出版社,1996,4(2):46-74.

[17]骆天舒;耗散动力学系统的广义哈密顿形式及其应用[D];浙江大学;2011,5-16.

[18] Michel D, Sylvie G , Laurent P F . Application of fuzzylogic control for continuous

casting mold level control[ J ]. IE EE Trans on Control System s Technology,1998,6(2):246-256.

[19]钱伟长,奇异摄动理论及其在力学中的应用[M].北京:科学出版社,1986,7(6):241-287.

[20] KanaiE.rPopTheorPhys,1948,20:440

[21]赵凯华等.电磁学(下)[M].北京市:人民教育出版社,1979,6(2):146-149.

[22]罗恩,李纬华. 略论变分原理之科学美 [J]. 力学与实践. 2008,(02):36-95.

Hamilton's principle and application of dissipative

system

Yan Junjie

(Shaanxi Institute of Physics and Telecommunication Engineering Physics 1202 class, Hanzhong,

Shaanxi 723000)

Tutor:Wang Jianhua

[Abstract]:Hamilton's principle is not only the basis of analyzing the mechanics, but also has important applications in quantum physics. This paper by defining the dissipation function, this paper gives a dissipative system of Hamilton's principle. Then use a dissipative Hamiltonian principle is deduced with Lagrange equation of dissipation system. Finally, with the linear viscous damping effect analogy, the mechanical system of Lagrange's equation with dissipative system applied to the dissipation of the circuit system, the kirchhoff's law is given.

[Key words]:Dissipative system; The dissipation function; Hamilton's principle. The generalized dissipation force; Kirchhoff's law.

哈密顿系统一些保结构算法的构造和分析

哈密顿系统一些保结构算法的构造和分析一切真实的,耗散可忽略不计的物理过程都可以用哈密顿系统进行描述.哈密顿系统有两个最重要的性质,一个是辛结构,另一个就是能量守恒.正确计算哈密顿系统非常重要.近年来,能够保持哈密顿系统辛结构或能量的保结构方法已经得到了很大的发展.本文讨论哈密顿系统一些保结构算法的构造和分析,主要研究成果如下:I.近几年,人们构造了等离子物理中洛伦兹力系统的保结构格式,比如保体积格式和保辛格式.然而这些格式都不能保持系统能量.我们把洛伦兹力系统写为一个非典则的哈密顿系统,然后利用Boole离散线积分方法进行求解,得到洛伦兹力系统的一个新的格式.该方法可以保持系统哈密顿能量达到机器精度.II.我们研究如何利用二,三和四阶AVF方法求解哈密顿偏微分方程.对非线性薛定谔方程,空间用Fourier拟谱方法半离散,时间用三个AVF方法进行离散,得到该方程三个不同精度的AVF格式.我们用数值实验验证了这三个格式的精度和保能量守恒特性.III.基于根树和B-级数理论,我们给出了5阶树的带入规则的具体公式.利用新得到的带入规则,我们把二阶AVF方法提高到高阶精度,给出了一个新的AVF方法.我们证明了,新方法具有6阶精度,并且可以保持哈密顿系统能量.我们利用六阶AVF方法求解非线性哈密顿系统,并测试了其精度和能量守恒特性.IV.在哈密顿偏微分方程保结构算法框架下,我们研究了基于系统弱形式的空间离散方法.首先,空间用有限元法或谱元法对偏微分方程进行半离散,把得到的常微分方程组写成一个哈密顿系统.然后,我们用一个保结构方法对这个常微分哈密顿系统进行求解,得到一个全离散保结构格式.我们用这个方法对一维非线性薛定谔(NLS)方程进行求解,其中空间用Legendre谱元法,时间用AVF 方法,得到一个新的保能量方法.同样对一维NLS方程,我们在空间用Galerkin

1哈密顿原理

牛顿质点动力学 1 牛顿第二定律 dt d p f 从三个方面来应用: 全局性研究:对称性、守恒律、稳定性; 局部研究:平均值、动量定理、动能定理; 瞬时研究:极限求导、奇异性、突变性; 2 重点研究非惯性、矢量性、连续性、相对性的问题; 3 从动力学观点上升到能量的观点。 哈密顿原理、保守力及其势 4 五大类典型模型 概括: 一个原理:哈密顿原理(稳定性与对称性原理); 哈密顿原理的文字表述如下: 二种建模方法:动力学方法、能量法; 三类研究方法:对称性方法(全局)、平均值方法(局部) 求极限、求导、突变及奇异性研究方法(瞬时);

四大重点问题:矢量性(矢量空间法)、连续性(微元动力学法)、相对性(相对速度公式法)、非惯性(等效性法); 五项典型模型:准粒子模型、碰撞模型、势模型、相空间模型、简谐振动与波模型。(科学计算技术与研究式的学习模式) 哈密顿原理、对称性和稳定性 1.拉格朗日函数和哈密顿量 拉格朗日函数L 对于一个物理系统,可用一个称为拉格朗日函数的量),,(t q q L i i &来描述,其中i q 是广义坐标,=i q &dt dq i /是广义速度;广义坐标与通常所说的坐标区别在于,广义坐标是针对系统的自由度确定的,譬如一个质点限制在半径R 的球面上运动,其坐标显然有x 、y 、z 三个,但广义坐标只有φθ,两个,其中?θcos sin R x =,θ?θcos ,sin sin R z b R y ==;一般由于运动受到约束,坐标与广义坐标的数量是不相等的,仅在无约束条件下,坐标与广义坐标的数目才是一样的,与坐标一样广义坐标的选取也不是唯一的。 在保守力作用下,系统的拉格朗日量L 定义为动能与势能之差;U T L -= 哈密顿量H 物理系统还可以用一个称之为哈密顿量的函数描述,在保守力作用下,哈密顿量定义为系统的动能与势能之和

哈密顿系统的数学建模与动力学分析.

1 引言 Hamilton动力系统理论有着悠久而丰富的历史,它本身是Lagrange力学的升华与推广,从数学角度看又是一门内容精深的相空间几何学,如辛几何、辛拓扑等都源于此.近几十年来,随着纯数学理论的不断发展与计算机的普遍应用,Hamilton动力系统理论又成为当今非线性科学中极其活跃而富有魅力的研究领域.由于这类系统广泛存在于数理科学、生命科学以及社会科学的各个领域,特别是天体力学、等离子物理、航天科学以及生物工程中的很多模型都以Hamilton系统的形式出现,因此该领域的研究多年来长盛不衰.本文利用Hamilton原理推导出了Hamilton系统的正则方程.最后利用Hamilton正则方程给出一个具体物理实例的数学模型并对其进行动态模拟仿真.

2 预备知识 2.1 状态空间的基本概念 1)状态 任何一个系统在特定时刻都有一个特定的状态,系统在0t 时刻的状态是0t 时刻的一种信息量,它与此后的输入一起惟一地确定系统在0t t ≥时的行为. 2)状态变量 状态变量是一个完全表征系统时间域行为的的最小内部变量组. 3)状态向量 设系统有n 个状态变量,用()()()12,, ,n x t x t x t 表示,而且把这些状态变量看做向量 ()x t 的分量,则向量()x t 称为状态向量,记为 ()()()()12,, ,T n x t x t x t x t =????. 4)状态空间 以状态变量()()()12,,,n x t x t x t 为轴的n 维实向量空间称为状态空间. 5)状态方程 描述系统状态变量与输入变量之间关系的一阶微分方程组(连续时间系统)或一阶差分方程组(离散时间系统)称为系统的状态方程,它表征了输入对内部状态的变换过程,其一般形式为: ()()(),,x t f x t u t t =???? 其中,t 是时间变量,()u t 是输入变量. 6)输出方程 描述系统输出量与系统状态变量和输入变量之间函数关系的代数方程称为输出方程,它表征了系统内部状态变化和输入所引起的系统输出变换,是一个变化过程.输出方程的一

对数哈密顿方法及其应用

对数哈密顿方法及其应用 天体力学数值方法作为天体力学的重要领域之一在辛算法的提出后得到长足发展,辛算法保持哈密顿系统辛结构且计算过程中系统没有能量和角动量的长期误差累积。辛算法适用于哈密顿系统的长期定性演化研究同时也具有数值精度不高、显辛算法要求固定步长的不足。 通常积分计算天体紧密交汇问题或大偏心率轨道运动都需缩短步长来克服天体受引力过大而剧增的加速度,直接变步长将丢失辛算法保持辛结构的优势,考虑时间变换的思路,原时间变量取变步长而新的时间变量仍为固定步长,则既能调节步长又能保持辛算法固有优势。本文的主要内容为构造针对不同哈密顿系统的对数哈密顿算法及论证其在具有更高的数值精度和保证获得有效的混沌判别结果方面的优势。 针对不同的哈密顿系统结构构造不同形式的时间变换辛算法。对于可分解为分别只含状态量广义动量和广义坐标的动能部分和势能部分的哈密顿函数,可构造取时间变换函数为形式不同但等价的两个函数得到显式对数哈密顿方法,其中时间变换作用于哈密顿函数,本文构造了由三个二阶蛙跳算子构成的显式对数哈密顿Yoshida四阶方法。 对于动能部分具有广义动量和广义坐标的交叉项而势能部分仅含位置变量的系统构造显隐式混合对数哈密顿方法,对于动能部分应用隐式中点法。而对于更一般的系统则构造隐式对数哈密顿方法。 隐式方法具有更广泛的应用但也由于算法构造中包括迭代需耗费更多的计算机时间降低计算效率。本文详细论证了显式对数哈密顿方法在应用于牛顿圆型限制性三体问题及相对论圆型限制性三体问题时较于非时间变换辛算法更具数

值精度优势。 且在前一系统的精度优势独立于轨道偏心率的变化。对于后一系统这一现象未能发生但数值精度也明显优越于常规辛算法。 特别对于高偏心率轨道,非时间变换算法得到的虚假的混沌判别指标,如Lyapunov指标和快速Lyapunov指数(FLI)。而通过对数哈密顿方法则可获得可靠地定性分析结果,彻底地解决后牛顿圆型限制性三体问题的高偏心率轨道Lyapunov指数的过度估计和FLI快速增大的问题。 在得到论证后本文应用对数哈密顿方法讨论了动力学参数两主天体间距离的变化对动力学系统有序和混沌转化的影响。本文通过数值模拟验证了对数哈密顿方法具有更高的数值精度及可得到可靠的定性研究成果的优势。 适用于定性研究和定量计算高偏心率问题,为天体力学研究开拓了新思路。在实际的天体紧密交汇处的动力学演化提供反映动力学实质的积分工具。

7第5章哈密顿原理

第5章哈密顿原理 如前所述,力学的变分原理的实质是:将真实运动与可能发生的运动加以比较,建立判别准则以区分真实运动和可能的运动。哈密顿原理是通过真实运动与可能的运动在位形空间的位形轨迹加以比较,而哈密顿作用量S 是对不同的位形轨线取不同值的泛函,从而得到对真实运动来讲,哈密顿作用量的变分等于零。 将拉格朗日方程引人哈密顿函数,导出哈密顿正则方程;给出了一种对偶的数学体系,开拓了应用前景;由动力学普遍方程对时间积分,导出一个重要的力学变分原理——哈密顿原理,提出了将真实运动与同样条件下的可能运动区分开来的准则;对于有限过程,提供了一种动力学问题的直接近似解法。 5.1 哈密顿正则方程 哈密顿正则方程是分析力学中又一个重要的力学方程,它与拉格朗日方程等价,是2n 个一阶常微分方程组。我们知道,对于一个质点系统,在建立拉格朗日方程后,重要的问题是研究这个微分方程组的积分,但是求解往往是很困难的。哈密顿正则方程的重要性在于它将n 个二阶微分方程变换为2n 个一阶方程,而且结构对称、简洁,为正则积分理论创造了有利条件。若是说拉格朗日方程对分析力学起着开拓性作用,则哈密顿正则方程对分析力学中的积分理论起着基础的和推动的作用。哈密顿正则方程的重要性还在于在许多理论的定性研究中,并不需要求解微分方程组,而是将二阶微分方程变换为二个一阶方程并应用几何方法求解。 5.1.1 正则方程的建立 对于主动力均有势的k 个自由度的完整约束系统,其拉格朗日方程为 ),,2,1(0d d k j q L q L t j j ==??-???? ???? (5-1) 引入广义动量 ),,2,1(k j q L p j j =??= (5-2) 代入式(5-1),有 ),,2,1(k j q L p j j =??= (5-3) 设拉格朗日函数L 满足条件 0det 2≠??? ? ? ????k j q q L 于是,可由式(5-2)反解出 ),,2,1(),,,,,,(11k j t p p q q f q k k j j == (5-4) 式(5-3)和式(5-4)就把方程(5-1)由k 个二阶微分方程化为2k 个一阶微分方程,其中方程 组(5-4)并非正则形式。引入哈密顿函数

耗散动力学系统的广义哈密顿形式及其应用

耗散动力学系统的广义哈密顿形式及其应用经典力学中所研究大部分系统不是保守系统,所以很难将这类系统表示为经典的哈密顿力学形式(偶数维)以及与此等效的拉格朗日力学形形式或最小作用量变分原理形式。由于这几种数学形式是数值计算方法中辛几何算法的的基础和现代物理学的基础,所以极大地限制了辛几何算法在耗散系统的数值模拟领域的应用以及耗散系统的量子化等理论物理领域中的应用。 耗散动力学系统长时间跟踪问题是当前非线性力学研究领域的难点之一。对于低维耗散动力学系统,可以用各种半解析方法(小参数法,摄动法)求解。 即便如此,对于长时间跟踪,也存在所谓久期项问题(由方法本身的误差累积导致)。对于高维耗散动力系统,直接应用解析方法显然是十分困难的。 因此多采用数值方法求解该类问题。但是不同的数值方法求解的结果可能会有较大偏差,甚至相差甚远,而且大部分问题是缺乏判断其算法偏差量的参考标准的。 所以为这类问题挑选或者创立公认可行的数值积分方法,成为一个问题。我国著名学者冯康先生提出并研究了在保守系统领域的这类问题,给出了辛几何算法的思想并系统的表述构造辛差分格式的一般方法,指出了原有差分格式中的适于长时间跟踪的格式。 钟万勰先生发展了这种思想,进一步提出了时间有限元和精细积分的的思想,并对耗散动力学系统引入辛算法作了尝试。本文的最初的目的是在转子稳定性分析等耗散动力学问题中使用辛数值积分方法(或者说利用辛几何算法的思想找到合适的算法)。 为达到此目的研究了耗散系统和保守系统的一种特殊关系,在此基础上用相

应的保守系统的数值解替代原耗散系统,即将辛数值方法应用求解相应的保守系统来得到所要研究系统的数值解。在这种关系的基础上,借鉴流体力学的广义哈密顿方程和最小作用量变分原理,将耗散系统表示成一种无穷维广义哈密顿系统,相应地带来一种新型的最小作用量变分原理。 可以将冯康文献中广义哈密顿系统辛算法的思想应用于求解这个特殊的无穷维哈密顿系统。上述最小作用量变分原理,可以和路径积分量子力学形式结合,应用于量子力学领域。 以上工作的主要创新点可以归纳如下:1.发现了耗散力学系统和某一保守力学系统相曲线重合原理:对于一个耗散力学系统和它一个初始条件,对应于不同时间区段一定存在一族保守力学系统,这族保守力学系统和耗散力学系统有且仅有一条共同的相曲线;这族保守系统的哈密顿量就是前述耗散力学系统的总能。对于非保守的振动问题来说,这个保守系统就是一个非线性保守力学系统,其中的保守力在某一初始条件下和非保守振子系统的阻尼力和恢复力之和相等,那么其在相空间运动轨迹必然相同。 在此基础上,引入了无穷维广义哈密顿格式来表示耗散力学系统,在其中定义了一个新的哈密顿量,并且引入了新的泊松括号,这个格式类似于表示等离子问题和理想流体的广义哈密顿格式。在这里把耗散力学系统看作是相空间内一种特殊流体(内部无压力),初始条件看作是物质坐标,上述轨迹重合的保守力学系统的哈密顿量看作是哈密顿量密度。 对应于经典的哈密顿变分原理,这个广义哈密顿格式等效于一个新的变分原理。在这个变分原理中作用量为相空间的某一区域中所有微元的作用量之和。 2.从创新点1出发本文研究了有阻尼振动问题的中心差分格式,发现中心差

哈密顿原理

§7-4 哈密顿原理 人们为了追求自然规律的统一、 和谐, 按照科学的审美观点, 总是力图用尽可能少的原理(即公理)去概括尽可能多的规律. 牛顿提出的三个定律, 是力学的基本原理. 由这些基本原理出发, 经过严格的逻辑推理和数学演绎, 可以获得经典力学的整个理论框架. 哈密顿原理是分析力学的基本原理, 它潜藏着经典力学的全部内容并把这门学科的所有命题统一起来. 也就是说, 由它出发, 亦可得到经典力学的整个框架. 哈密顿原理是力学中的积分变分原理. 变分原理提供了一个准则, 使我们能从约束许可条件下的一切可能运动中, 将力学系统的真实运动挑选出来. 变分原理的这一思想, 不仅在力学中, 而且在物理学科的其他领域中, 都具有重要意义. 一、变分法简介 1. 函数的变分. 自变量为x 的函数表示为)(x y y =. 函数的微分x y y d d ′=是由自变量x 的变化引起的函数的变化. 函数的变分也是函数的微变量, 但它不是因为自变量x 的变化, 而是由于函数形式的变化引起

的. 这种由于函数形式变化造成的函数的变更称为函数的变分, 记作y δ. 与函数y 邻近但形式与y 不同的函数有许多, 这些函数可以表示如下: )()0,(),(* x x y x y εηε+= 其中ε是任意小的参数, ()x η是任意给定的可微函数. 因0=ε时()()x y x y =0,, 所以函数形式的变化决定于上式的第二项. 因此, 函数的变分写成 ()()()x x y x y y εηε=?=0,,δ* 在自由度为1的力学系统中讨论变分的概念. 设广义坐标为q , )(t q q =. 建立以t q ,为轴的二维时空坐标系(又称事件空间), 曲线I 是)(t q q =的函数曲线, 代表了系统的真实运动. q t d d →函数的微分. 在曲线I 附近, 存在 着许多相邻曲线, 这些曲 线都满足力学系统的约束 条件, 称为可能运动曲线, 它们的方程表示为 ()()()t t q t q εηε+=0,,* 在t 不变的情况下, 函数形式的改变也能引起函数的变化, 这种变化纯粹是由函数形式变化引起的, 它就是函数的变分q δ, ()()()t t q t q q εηεδ=?=0,,*

5.7哈密顿原理作业

1 哈哈密密顿顿原原理理作作业业 1.如图示,质量为m 的复摆绕通过某点O 的水平轴作微小振动,复摆对转轴的转动惯量为0I ,质心C 到悬点O 的距离为 ,试用哈密顿原理求该复摆的运动方程及振 动周期。 1.解:取θ为广义坐标,则拉格朗日函数为: θ+θ=-=cos mg I 21 V T L 2 0 其中取悬点O 为零势能点。 于是哈密顿原理0dt L 21t t =δ?可得:0dt cos mg I 2 121t t 20=??? ??θ+θδ? 即:()0dt sin mg I 2 1t t 0=θδθ-θδθ ? 而δθθ-δθθ=δθθ=θδθ 0 000I )I (dt d )(dt d I I 则:()0dt sin mg I )I (dt d dt sin mg I 212 1t t 00t t 0=??? ??θδθ-δθθ-δθθ=θδθ-θδθ ?? 即:()0dt sin mg I I 212 1t t 0t t 0=δθθ+θ-δθθ ? 而0I 21t t 0=δθθ ,δθ取任意值 所以:0sin mg I 0=θ+θ 即:0sin I mg 0=θ+θ 而θ≈θsin ,则:0I mg 0 =θ+θ ,此即为所求的运动方程。 其中角频率0I /mg =ω 所以振动周期)mg /(I 2/2T 0 π=ωπ=。 2.试用哈密顿原理求质量为m 的质点在重力场中用直角坐标系表示的运动微分方程。 2.解:取x,y,z 为广义坐标,则: 体系的动能)z y x (m 2 1 T 222 ++= 势能mgz V =(以地面为零势能点) 拉氏函数mgz )z y x (m 21 V T L 222-++=-=

约束Hamilton系统的稳定性研究

约束Hamilton 系统的稳定性研究 郑明亮1) 傅景礼 2) 1)(浙江理工大学 机械设计与控制学院 杭州 310018) 2)(浙江理工大学 理学院 杭州 310018) 摘要:本文给出了一种约束Hamilton 系统的稳定性判断方法。首先,提出将因系统奇异性导致的内在限制方程看作是非完整约束方程,采用Routh 方法导出了约束Hamilton 系统的运动正则方程。其次,将约束Hamilton 系统转化成力学梯度系统,给出转化微分方程表示的条件和表达形式;接着,根据梯度系统的性质结合李雅普诺夫的一次近似理论直接来判定约束Hamilton 系统的平衡位置稳定性。最后,举例说明结果的应用。 关键词:约束Hamilton 系统;梯度系统;李雅普诺夫;稳定性 PACS:45.10.Hj,02.30.Hq 1引言 力学系统的运动稳定性在数学、力学、航空、航海、航天、新技术和高技术中得到广泛应用,发挥了越来越大的作用[1]。关于稳定性的问题Lyapunov 首先给出了稳定性的严格数学定义,并提出一种研究运动稳定性的直接方法。Bottema [2]研究了在·ГAO Ⅱ?意义下,各种力学系统平衡位置的稳定性判断方法。Risito [3]和 Laloy [4]总结了保守系统和耗散系统的平衡和运动稳定性,得到线性、齐次、定常非完整系统平衡位置稳定与不稳定的一些更特殊的结果。我国著名力学专家梅凤翔[5]系统地论述了约束力学系统的运动稳定性问题。朱海平 [6]研究了非完整系统的稳定性。傅景礼等[7-8]研究了相对论性和转动相对论性Birkhoff 系统的平衡稳定性。Zhang [9]利用Noether 守恒量构造了Lyapunov 函数,研究了广义Birkhoff 系统的运动稳定性。姜文安等[10]研究了广义Hamilton 系统的运动稳定性。Cheng [11]研究了系统参数对带附加广义力项的约束力学系统运动稳定性的影响。 在Legendre 变换下,奇异Lagrange 系统在过渡到相空间用Hamilton 正则变量描述时,其正则变量之间存在固有约束,称之为约束Hamilton 系统[12]。机械工程和数学物理上许多重要的动力系统是约束Hamilton 系统,如非树形多体机器人系统动力学模型一般为微分/代数方程组形式[13]、光的横移现象和量子电动力学[14]等。但是,关于约束Hamilton 系统的稳定性研究一直鲜有报道。如果一个力学系统能够成为梯度系统,那么就可用梯度系统的特性来研究力学系统的性质, 特别是运动稳定性质[15]。本文研究仅含第二类约束的约束Hamilton 系统的稳定性,将其转化成梯度系统,直接利用Lyapunov 定理来研究其平衡稳定性。 2约束Hamilton 系统的正则方程 设力学系统的位形由n 个广义坐标),...,1(n s q s =来确定,系统的Lagrange 函数为 ),,(q q t L ,广义动量为),...,1(n s q L p s s =??= ,设L 的Hess 矩阵?? ???????k s q q L 2的秩为n r <。 引入系统的Hamilton 函数为),(1q p,t L q p H n i i i -=∑= ,将奇异Lagrange 系统描述过渡到Hamilton 系统描述时,在相空间中正则变量之间存在代数约束方程: ),...,1(,0),(r n j t j -==Φq p, (1)

1哈密顿原理

1哈密顿原理

牛顿质点动力学 1 牛顿第二定律 dt d p f 从三个方面来应用: 全局性研究:对称性、守恒律、稳定性; 局部研究:平均值、动量定理、动能定理; 瞬时研究:极限求导、奇异性、突变性; 2 重点研究非惯性、矢量性、连续性、相对性的问题; 3 从动力学观点上升到能量的观点。 哈密顿原理、保守力及其势 4 五大类典型模型 概括: 一个原理:哈密顿原理(稳定性与对称性原理); 哈密顿原理的文字表述如下: 保守的、完整的力学体系在相同时间内,由某一初位形转移到另一已知位形的一切可能运动中,真实运动的主函数具有稳定值,即对于真实运动来讲,主函数的变分等于0。 二种建模方法:动力学方法、能量法; 三类研究方法:对称性方法(全局)、平均值方法(局部) 求极限、求导、突变及奇异性研究方法(瞬时);

四大重点问题:矢量性(矢量空间法)、连续性(微元动力学法)、相对性(相对速度公式法)、非惯性(等效性法); 五项典型模型:准粒子模型、碰撞模型、势模型、相空间模型、简谐振动与波模型。(科学计算技术与研究式的学习模式) 哈密顿原理、对称性和稳定性 1.拉格朗日函数和哈密顿量 拉格朗日函数L 对于一个物理系统,可用一个称为拉格朗日函数的量),,(t q q L i i &来描述,其中i q 是广义坐标,=i q &dt dq i /是广义速度;广义坐标与通常所说的坐标区别在于,广义坐标是针对系统的自由度确定的,譬如一个质点限制在半径R 的球面上运动,其坐标显然有x 、y 、z 三个,但广义坐标只有φθ,两个,其中?θcos sin R x =,θ?θcos ,sin sin R z b R y ==;一般由于运动受到约束,坐标与广义坐标的数量是不相等的,仅在无约束条件下,坐标与广义坐标的数目才是一样的,与坐标一样广义坐标的选取也不是唯一的。 在保守力作用下,系统的拉格朗日量L 定义为动能与势能之差;U T L -= 哈密顿量H 物理系统还可以用一个称之为哈密顿量的函数描述,在保守力作用下,哈密顿量定义为系统的动能与势能之和

1哈密顿原理-新版.pdf

牛顿质点动力学 1 牛顿第二定律dt d p f 从三个方面来应用: 全局性研究:对称性、守恒律、稳定性;局部研究:平均值、动量定理、动能定理;瞬时研究:极限求导、奇异性、突变性; 2 重点研究非惯性、矢量性、连续性、相对性的问题; 3 从动力学观点上升到能量的观点。哈密顿原理、保守力及其势 4 五大类典型模型概括: 一个原理:哈密顿原理(稳定性与对称性原理); 哈密顿原理的文字表述如下: 保守的、完整的力学体系在相同时间内,由某一初位形转移到另一已知位形的一切可能运动中, 真实运动的主函数具有 稳定值,即对于真实运动来讲,主函数的变分等于0。二种建模方法:动力学方法、能量法; 三类研究方法:对称性方法(全局)、平均值方法(局部)求极限、求导、突变及奇异性研究方法(瞬时) ;

四大重点问题:矢量性(矢量空间法)、连续性(微元动力学法)、相对性(相对速度公式法)、非惯性(等效性法);五项典型模型:准粒子模型、碰撞模型、势模型、相空间模型、简谐振动与波模型。(科学计算技术与研究式的学习模式) 哈密顿原理、对称性和稳定性 1.拉格朗日函数和哈密顿量拉格朗日函数 L 对于一个物理系统,可用一个称为拉格朗日函数的量 ),,(t q q L i i 来描述,其中i q 是广义坐标,i q dt dq i /是广义速 度;广义坐标与通常所说的坐标区别在于,广义坐标是针对系统的自由度确定的,譬如一个质点限制在半径R 的球面上 运动,其坐标显然有x 、y 、z 三个,但广义坐标只有 ,两 个,其中cos sin R x ,cos ,sin sin R z b R y ;一 般由于运动受到约束,坐标与广义坐标的数量是不相等的,仅在无约束条件下,坐标与广义坐标的数目才是一样的,与坐 标一样广义坐标的选取也不是唯一的。 在保守力作用下,系统的拉格朗日量L 定义为动能与势 能之差;U T L 哈密顿量H

简单的论述哈密顿原理

简单的论述哈密顿原理 摘要:证明力积分变量与变分无关的情况下积分运算与变分运算次序的可交换性,从不同角度论述了哈密顿原理的含义。 关键词:哈密顿原理,拉格朗日函数,变分,拉格朗日方程 1.引言 哈密顿原理是分析力学中几个重要原理之一,但它不是一个独立原理,它可已从其他原理推导出来,因而可以从不同角度说明它的物理含义。一般理论力学教材都是在拉格朗日方程两边同时乘以虚位移求所有自由度下的虚功之和,然后再求从位形1即(到位形2,即(之间或时间至 之间的作用量得出,最后变换成,并没有说明最后一步为 什么要那样做,也没有说明那样做的意义。本文先证明当积分变量与变分无关的条件下积分运算与变分运算次序的可交换性,然后再从不同角度论述哈密顿原理的意义。 2.理论 2.1变分运算与积分运算次序的可交换性 假定变量由一个或一组函数的选取而确定,则变量称 为函数的泛函,记作[]。泛函由n个函数的形式确定,是函数的“函数”。泛函与函数的概念略有不同,函数中的变量是可以变化的数值,而对于泛函处于自变量地位的是形式可以变化的函数。下面举例说明,如图1中有,两个固定点,连接两个固定点之间的曲线的长度由下式确定,即

显然,依赖于函数的选取,若函数的形式发生变化,则曲线的形状随 之变化,曲线的长度也随之变化。长度就是的 泛函。 下面证明变分运算与积分运算顺序的可交 换性,该泛函只依赖一个函数,即 自变量为的函数表示为。函数的变分是函数的微变量,它与函数的微分有本质有本质的不同,函数的微分,粗略的讲,它是由自变量的变化引起的。而函数的变分不是因为自变量的变化,它是来自函数形式的变化引起,这种由于函数形式变化造成的函数的变化称为函数的变分,记作。与函数临近但形 式与不同的函数有许多。 假设这些函数可以表示为如下的形式: 其中是非常小的参数,是任意给定的可微函数,因时,函 数形式的变化决定于上式的第二项。因此函数的变分写成 引入(2)式的记法(1)可记为 被积函数的形式是已知的,积分的上下限是固定的。当函数 的形式上发生变化时,泛函就会发生变化,这种由于函数形式的变化引起泛函的变化就为泛函的变分,记作。现将被积函数

经典力学的哈密顿理论.

第八章 经典力学的哈密顿理论 教学目的和基本要求:理解正则共轭坐标的物理意义并掌握如何用正则坐标表示体系哈密顿函数;能熟练应用正则方程求解简单的力学问题的;了解变分问题的欧拉方程;掌握用变分法表示的哈密顿原理并能正确理解哈密顿原理的物理含义;初步掌握正则变换、泊松括号的物理意义和使用方法。 教学重点:在正确理解正则共轭坐标的物理意义的基础上能熟练应用正则方程求解简单的力学问题。 教学难点:正则共轭坐标的意义和哈密顿原理的物理含义。 §8.1 正则共轭坐标 坐标的概念是随着物理学的发展而发展,我们在本节将要讨论一种全新的坐标——正则共轭坐标。 一:坐标的发展历史. 1.笛卡儿直角坐标。为了研究物体在三维空间的位置、速度和加速度而引入的坐标。其用 z y x ,,三个变量来描述物体在空间任一点的位置,坐标轴的方向不随物体的运动而改变, 用k j i ,,来表示三个坐标轴方向的单位矢量。 2.极坐标、柱坐标和球坐标。用两个或三个变量来反映物体在平面或空间的位置。在处理转动问题和中心势场的力学问题时比直角坐标更优越。其代表坐标轴方向的单位矢量为变 矢量,利用这些矢量可以很方便地表达上述力学问题的a v ,等物理量。从直角坐标到极坐 标、柱坐标和球坐标等曲线坐标是坐标历史上的第一次飞跃。 另外曲线坐标还包括自然坐标,利用它处理运动规律已知的物体的力学问题更为方便。 3.广义坐标。反映力学体系在空间位形的独立变量被称为广义坐标。它是拉格朗日方程建立的基础和优越性所在,也是分析力学的基础。广义坐标不仅拓宽了坐标的概念,而且由它所列出的动力学方程不含非独立变量,使方程的求解过程得到了简化。另外我们在研究体系的微振动时引入了简正坐标,使微振动方程的求解过程非常简单,这是坐标概念的第二次飞跃。 下面我们将介绍的正则共轭坐标是坐标概念的第三次飞跃。

哈密顿系统中混沌的几何判据

哈密顿系统中混沌的几何判据 【摘要】:用几何方法研究哈密顿系统的混沌是近二十年来出现的新领域。本论文研究了几类典型的哈密顿系统,并给出了一系列哈密顿系统混沌的几何判据,揭示了哈密顿系统内在的几何性质与其混沌行为的本质联系。第二章我们推广了L.Horwitz等人在2007年提出的判断混沌的几何方法,使得该方法不仅适用于标准哈密顿系统,还适用于势能与动量弱耦合的情况。提出了平均不稳定比(MUR)的概念,并对Dicke模型的经典系统做了计算。推广的方法MUR不仅和Poincare 截面方法的结果吻合得很好,而且在数值稳定性上优于人们熟知的最大李雅普诺夫指数方法。第三章主要研究了二维哈密顿系统势能面、等势线与混沌之间的关系。我们发现势能面的凹陷区域是稳定区域,凸起区域和既不凸也不凹的区域都是不稳定的。另外还证明了如果系统的等势线有凹向平衡点的部分,系统将是不稳定的。以此为依据我们提出了判断二维哈密顿系统混沌的平均凸指标(MCI)和凹比率(CR)。我们对典型的混沌模型进行数值计算后,发现MCI、CR和Poincare截面、L.Horwitz等人的新几何方法的数值结果完全一致。MCI和CR直观简洁,为混沌的几何研究方法提供了新观点和新内容。第四章研究了Dicke模型中混沌与几何相位之间的联系。当光场和原子的耦合强度增大至临界点时,Dicke系统的能级间距概率分布从泊松分布变为Wigner分布,而Wigner分布被视为量子混沌的标志,这说明Dicke量子系统在临界点开始出现量子混沌;与Dicke量子系统对

应的经典系统在临界点也从规则运动变为混沌运动。在临界点处Dicke量子系统基态的几何相位即Berry相位也发生了剧烈的变化。我们引入了几何相位阶数的概念,Dicke系统几何相位的阶数在临界点从有限值跃变为∞。我们把Dicke量子系统基态几何相位阶数的跃变作为量子混沌出现的标志。【关键词】:哈密顿系统混沌量子混沌几何方法几何相位 【学位授予单位】:山西大学 【学位级别】:博士 【学位授予年份】:2011 【分类号】:O415.5 【目录】:中文摘要8-10ABSTRACT10-12第一章绪论12-341.1混沌研究简史12-141.2混沌的基本特征14-171.3哈密顿系统中的混沌17-241.3.1哈密顿力学17-201.3.2KAM定理20-241.4混沌研究的常用方法24-31参考文献31-34第二章混沌研究的几何方法34-502.1混沌研究的几何方法34-372.2混沌的新几何判断方法37-422.3推广的新几何判断方法42-462.4小结46-47参考文献47-50第三章二维哈密顿系统中的势能面、等势线与混沌50-683.1二维哈密顿系统的不稳定判据50-523.2二维哈密顿系统中的势能面与混沌52-573.3二维哈密顿系统中的等势线与混沌57-653.4小结65-66参考文献66-68第四章混

哈密顿原理

哈密顿原理 文件排版存档编号:[UYTR-OUPT28-KBNTL98-UYNN208]

牛顿质点动力学 1 牛顿第二定律 dt d p f 从三个方面来应用: 全局性研究:对称性、守恒律、稳定性; 局部研究:平均值、动量定理、动能定理; 瞬时研究:极限求导、奇异性、突变性; 2 重点研究非惯性、矢量性、连续性、相对性的问题; 3 从动力学观点上升到能量的观点。 哈密顿原理、保守力及其势 4 五大类典型模型 概括: 一个原理:哈密顿原理(稳定性与对称性原理); 二种建模方法:动力学方法、能量法; 三类研究方法:对称性方法(全局)、平均值方法(局部) 求极限、求导、突变及奇异性研究方法(瞬时); 四大重点问题:矢量性(矢量空间法)、连续性(微元动力学法)、相对性(相对速度公式法)、非惯性(等效性法); 五项典型模型:准粒子模型、碰撞模型、势模型、相空间模型、简谐振动与波模型。(科学计算技术与研究式的学习模式) 哈密顿原理、对称性和稳定性

1.拉格朗日函数和哈密顿量 拉格朗日函数L 对于一个物理系统,可用一个称为拉格朗日函数的量),,(t q q L i i 来描述,其中i q 是广义坐标,=i q dt dq i /是广义速度;广义坐标与通常所说的坐标区别在于,广义坐标是针对系统的自由度确定的,譬如一个质点限制在半径R 的球面上运动,其坐标显然有x 、y 、z 三个,但广义坐标只有φθ,两个,其中?θcos sin R x = , θ?θcos ,sin sin R z b R y ==;一般由于运动受到约束,坐标与广义 坐标的数量是不相等的,仅在无约束条件下,坐标与广义坐标的数目才是一样的,与坐标一样广义坐标的选取也不是唯一的。 在保守力作用下,系统的拉格朗日量L 定义为动能与势能之差; U T L -= 哈密顿量H 物理系统还可以用一个称之为哈密顿量的函数描述,在保守力作用下,哈密顿量定义为系统的动能与势能之和 ),,(t p q H i i =U T +(i=1,2…s ) 其中 )(/i i q L p ??=是广义动量,哈密顿量是广义坐标和广义动量的函数,在直角坐标下对于质点运动的广义动量可写成v p m =。作用量I 定义为 ?=2 1 t t Ldt I 其中,积分上下限是质点初末态I q 、F q 对应的时间。 2.哈密顿原理及轨道稳定性

随机哈密顿系统的变分积分子与生成函数

随机哈密顿系统的变分积分子与生成函数 哈密顿系统具有辛结构,冯康院士等提出了保持这种结构的辛几何算法。变分积分子和生成函数是两种系统性构造辛算法的方法,前者基于确定性系统的作用积分和欧拉-拉格朗日方程,以及一个辛映射的等价形式,后者基于此等价形式和一系列等价的坐标变换,以及哈密顿-雅可比偏微分方程及其近似解。 对于随机系统,Milstein等根据辛结构的保持提出随机哈密顿系统的定义,并给出一些随机辛方法。这些方法主要是将确定性系统的辛方法进行适当调整,使之符合随机系统数值方法相容性及可实现性的要求。相对于确定性系统,随机辛算法的研究尚处于起步阶段,特别是对随机变分积分子和生成函数的探讨在文献中还未见到。 受非保守系哈密顿及拉格朗日方程形式的启发,我们从假定白噪声为影响系统的非保守力出发,构造出随机哈密顿系统的拉格朗日方程和作用积分,证明了随机哈密顿原理,在此基础上,并基于一个随机辛方法的等价形式,提出了随机变分积分子理论。经等价坐标变换,找到三种类型的随机生成函数,并推导了随机哈密顿-雅可比偏微分方程,找到了一种近似求解随机哈密顿-雅可比偏微分方程的方法,从而使利用随机生成函数构造随机辛算法成为可能,并且可以在理论上控制算法的收敛阶。对文献中一些已有的辛算法,我们给出了其生成函数,并将带可加噪声哈密顿系统的随机辛Runge-Kutta方法扩展为应用于一般随机哈密顿系统的辛Runge-Kutta方法,给出了它的三种类型的生成函数。理论分析及数值实验表明,随机变分积分子及生成函数是系统性构造随机辛格式的有效方法,这些构造出的辛格式具有长时间保持随机系统原有结构的优越特性。 关键词:随机哈密顿系统,辛结构,变分积分子,生成函数

由哈密顿原理推导拉格朗日方程

由哈密顿原理推导拉格朗日方程 谭建 222010315210236 2010级4班 一、问题重述 试由210t t Ldt δ=?推导()0d L L dt q q αα ???-=?? 二、问题分析及 由于是等时变分,有()d q q dt δδ?= ,和 22 11()0t t t t Ldt L dt δδ==?? (1) 现在来秋L δ。L 是q , q ? , t 的函数,又由于是等时变分,所以有 L L L q q q q δδδ????=+??……………………..(2) ()()()L L d d L d L q q q q dt dt dt q q q q δδδδ?????????==?-????……………….(3) 将(3)代入(2)得 ()()d L d L L L q q q dt dt q q q δδδδ?????=?-+???…………………………(4) 将(4)代入(1)得 2121()(())0t t t t L d L L q q q dt dt q q q δδδ??????+-+=????…………………………….(5) 在12,t t 处0q δ=,所以(5)变为 2 1(())0t t d L L q q dt dt q q δδ???-=???………………………………(6)即 2 1[(())]0t t d L L q dt dt q q δ???-+=???……………………………………(7) q 是独立变量,所以有 ()0d L L dt q q ???-+=??即 ()0d L L dt q q ???-=??此式即为拉格朗日方程

7第5章哈密顿原理

48 第5章哈密顿原理 如前所述,力学的变分原理的实质是:将真实运动与可能发生的运动加以比较,建立判别准则以区分真实运动和可能的运动。哈密顿原理是通过真实运动与可能的运动在位形空间的位形轨迹加以比较,而哈密顿作用量S 是对不同的位形轨线取不同值的泛函,从而得到对真实运动来讲,哈密顿作用量的变分等于零。 将拉格朗日方程引人哈密顿函数,导出哈密顿正则方程;给出了一种对偶的数学体系,开拓了应用前景;由动力学普遍方程对时间积分,导出一个重要的力学变分原理——哈密顿原理,提出了将真实运动与同样条件下的可能运动区分开来的准则;对于有限过程,提供了一种动力学问题的直接近似解法。 5.1 哈密顿正则方程 哈密顿正则方程是分析力学中又一个重要的力学方程,它与拉格朗日方程等价,是2n 个一阶常微分方程组。我们知道,对于一个质点系统,在建立拉格朗日方程后,重要的问题是研究这个微分方程组的积分,但是求解往往是很困难的。哈密顿正则方程的重要性在于它将n 个二阶微分方程变换为2n 个一阶方程,而且结构对称、简洁,为正则积分理论创造了有利条件。若是说拉格朗日方程对分析力学起着开拓性作用,则哈密顿正则方程对分析力学中的积分理论起着基础的和推动的作用。哈密顿正则方程的重要性还在于在许多理论的定性研究中,并不需要求解微分方程组,而是将二阶微分方程变换为二个一阶方程并应用几何方法求解。 5.1.1 正则方程的建立 对于主动力均有势的k 个自由度的完整约束系统,其拉格朗日方程为 ),,2,1(0d d k j q L q L t j j (5-1) 引入广义动量 ),,2,1(k j q L p j j (5-2) 代入式(5-1),有 ),,2,1(k j q L p j j (5-3) 设拉格朗日函数L 满足条件 0det 2 k j q q L 于是,可由式(5-2)反解出 ),,2,1(),,,,,,(11k j t p p q q f q k k j j (5-4) 式(5-3)和式(5-4)就把方程(5-1)由k 个二阶微分方程化为2k 个一阶微分方程,其中方程组(5-4)并非正则形式。引入哈密顿函数

7第5章哈密顿原理

第 5 章哈密顿原理 如前所述,力学的变分原理的实质是:将真实运动与可能发生的运动加以比较,建立判 别准则以区分真实运动和可能的运动。 哈密顿原理是通过真实运动与可能的运动在位形空间 的位形轨迹加以比较, 而哈密顿作用量 S 是对不同的位形轨线取不同值的泛函, 从而得到对 真实运动来讲,哈密顿作用量的变分等于零。 将拉格朗日方程引人哈密顿函数,导出哈密顿正则方程;给出了一种对偶的数学体系, 开拓了应用前景; 由动力学普遍方程对时间积分, 导出一个重要的力学变分原理——哈密顿 原理, 提出了将真实运动与同样条件下的可能运动区分开来的准则; 对于有限过程, 提供了 一种动力学问题的直接近似解法。 5.1 哈密顿正则方程 哈密顿正则方程是分析力学中又一个重要的力学方程,它与拉格朗日方程等价,是 2n 个一阶常微分方程组。我们知道,对于一个质点系统,在建立拉格朗日方程后,重要的问题 是研究这个微分方程组的积分, 但是求解往往是很困难的。 哈密顿正则方程的重要性在于它 将 n 个二阶微分方程变换为 2n 个一阶方程,而且结构对称、简洁,为正则积分理论创造了 有利条件。 若是说拉格朗日方程对分析力学起着开拓性作用, 则哈密顿正则方程对分析力学 中的积分理论起着基础的和推动的作用。 哈密顿正则方程的重要性还在于在许多理论的定性 研究中, 并不需要求解微分方程组, 而是将二阶微分方程变换为二个一阶方程并应用几何方 法求解。 5.1.1 正则方程的建立 设拉格朗日函数 L 满足条件 2 L det 0 q j q k 于是, 可由式 (5-2)反解出 q j f j (q 1, ,q k , p 1, , p k , t) (j 1,2, ,k) 式(5-3)和式(5-4)就把方程 (5-1)由 k 个二阶微分方程化为 组(5-4) 并非正则形式。引入 哈密顿函数 k 个自由度的完整约束系统,其拉格朗日方程为 dL dt q j L q j 0 ( j 1,2, ,k) 引入广义动量 p j L (j 1,2, ,k) q j 代入式 (5-1) ,有 L p j (j 1,2, ,k) (5-1) (5-2) (5-3) (5-4) 2k 个一阶微分方程,其中方程 对于主动力均有势的 q j

耗散系统的哈密顿原理及其应用

题目耗散系统的哈密顿原理及其应用 学生闫俊杰学号1210014049 所在学院物理与电信工程学院 专业班级物理学1202 班 指导教师王剑华__ __ 完成地点理工学院

2016 年 5 月20 日

耗散系统的哈密顿原理及其应用 作者:闫俊杰 (理工学院物理与电信工程学院物理专业1202班级,723000) 指导老师:王剑华 [摘要]哈密顿原理不仅是分析力学的基础,而且在量子物理中也有重要的应用。本文首先通过定义耗散函数,给出了有耗散系统的哈密顿原理。然后利用有耗散系统的哈密顿原理推导出了有耗散系统的拉格朗日方程。最后通过与受线性粘滞阻尼作用的力学系统类比,把有耗散系统的拉格朗日方程应用到有耗散的电路系统,给出了基尔霍夫定律。 [关键词] 耗散系统;耗散函数;哈密顿原理;广义耗散力;基尔霍夫定律. 引言 哈密顿原理,于1834年首次发表,从而完成了从莫泊丢开始的尝试,该原理的数学形式简洁,容广泛,将动力学问题转化为数学的一般体系的一部分,更深刻地揭示了客观事物之间的紧密联系,在物理学中有极高的地位,它不仅可以看作是力学的最高原理,甚至可以看作是整个物理学的最高原理,利用广义坐标并定义哈密顿函数就可以从哈密顿原理推导出哈密顿正则方程,致各种动力学定律都可以从这一变分式推出.所以说它是成为牛顿之后力学发展一个最大飞跃[1]。伴随经典力学从自由到非完整约束系统发展历程的是学术讨论甚至是争论,有些是对错之争,有些只是观点不同,问题的焦点主要是经典力学理论怎样从完整系统推广到非完整尤其是非线性非完整系统,这是经典力学后牛顿发展的主要标志,也是物理学近、现代发展的一块重要里程碑。物理学家考察具体的物理问题,常常以最小作用量原理为出发点,通过变分运算而导出物质系统的运动方程以表示其运动规律。牛顿质点运动方程可通过这样的推导过程得出;甚至如广义相对论里的爱因斯坦场方程也可如此导出,其实,所有运动方程就是拉格朗日方程或哈密顿正则方程在其拉氏量或哈氏量对于不同物质系统取不同形式时的具体体现。亦是广义相对论的一个重大结论;而如前文已提及的,引力场方程的导出可采用引力场之作用量的变分运算。从光线路径和质点运动到四维弯曲时空中的短程线方程乃至引力场运动方式等等实例可见,无论是几何间题,还是物理问题,都可凭籍变分法去圆满地解决;特别是最小作用量原理及其在物理学各领域的成功应用,正就是利用变分法这一几何方法、亦乃经济原则去解析各物质系统之运动规律的丰盈成果。例如牛顿力学描绘了天体运动的出色图景,拉格朗日一哈密顿理论描绘此图景也毫不逊色。当然,在经典力学畴里,这两个理论体系本来是等价的,只是着眼点有所不同:对于行星运动,牛顿着眼于行星与恒星之间的万有引力;哈密顿等人着眼于行星运行时的能量守恒正因为后者着眼于能量状况,

相关文档
最新文档