文献综述题目：线性谐振子相图研究姓名：学号：系别：物理与电子信息工程系专业：物理学年级：指导教师：2009年2月7 日文献综述一、前言线性谐振子是量子力学中可以精确求解的有限几个事例之一[1]，其中最简单的线性谐振子是简谐振子。

自然界中任何一个力学系统，只要某一个物理量在其稳定平衡点附近作微小振动，便可以用简谐振子模型来描述，例如：复摆的振动、分子的振动、晶格的振动、原子核表面振动以及辐射场的振动等。

在选择适当的坐标系之后，复杂的运动往往可以分解成若干彼此独立的一维简谐振动（simple harmonic vibration ）。

简谐振动作为一种最简单最基本的振动，往往还是复杂运动的初步近似，是研究振动的基础。

因此研究它在理论上和应用上都有重大的意义。

其中从相空间的角度来研究振动系统的力学问题如今已经成为一个研究趋势。

因为相图里包含着完整的力学系统的全部信息，无须去解复杂的运动方程[2]。

计算机技术软硬件的飞速发展，为此研究趋势提供了现实条件。

本论文从简谐振子的基本定义出发，在Fortran 90条件下进行数值模拟并在Origin75 软件下获得简谐振子的相图。

二、主体2.1简谐振动的定义定义一： 物体只在弹性力或准弹性 （线性回复力）作用下发生的运动，即动力学方程为的运动为简谐振动[2]。

定义二： 在无外来强迫力作用下, 物体相对于平衡点的位移随时间按余弦（或正弦）规律变化即 则称物体作简谐振动式即简谐振动的表达式[3]。

—振幅；—角频率；—相位；—初相位。

位移随时间的变化曲线称为振动曲线。

广义定义：某个物理量随时间的变化是按正弦或余弦规律，则可称该物理量做简谐动，可用表示 。

自然界中任何一个力学系统中，只要某一个物理量在其稳定平衡点附近作微小振动，便可以用这种简谐振子模型来描述，例如：复摆的振动、分子的振动、晶格的振动，原子核表面振动、辐射场的振动以及电磁场振动等等。

2.2简谐振动的基本特征及动力学特征简谐振动位移随时间的变化 cos()x A t ωφ=+222d d xx o tω+=()cos()x t A t ωφ=+()cos()x t A t ωφ=+物体作简谐振动时,速度为：物体作简谐振动时，加速度为：可见物体做简谐振动时，其速度、加速度都以同样的角频率作简谐振动，相位依次超前π/2。

根据牛顿第二定律Fm a= ，得22cos()F m a m A t m xωωφω==-+=-力与位移大小成正比，符号相反，这样的力就是线性回复力。

这是简谐振动的动力学特征。

2.3相空间和相迹的概念 H q P αα∙∂=∂(1,2,3)s α=H p q αα∙∂=-∂(1,2,3)s α=以上是哈密顿正则方程[4]，其中q α为广义坐标，q α∙为广义速度,p α广义动量。

而正则方程共有 2S 个相互联立的一阶常微分方程组 ，对这 2S 个方程求解，即可完全确定力学系统的运动状态。

在哈密顿方法里，我们引入相空间概念[2]。

相空间：是以S 个广义坐标qα 和 S 个广义动量p α为变数而构成的2S 维抽象空间称为力学系统的相空间。

任一瞬时力学系统的广义坐标和广义动量确定了相空间的一个点，我们称为相点，每个相点对应于系统的一个确定状态，当时间变化时，由于系统运动，这个相点也在相空间中运动，它在相空间中描画出的一条曲线称为相迹。

对于二阶系统，比如要研究简谐振动，它的状态变量只有两个，所以简谐振动的相空间即简化为二维相平面。

相平面：对于以运动物体的位移（x ）和速度(y)作为坐标参量构建的空间，就是的相平面。

相点：相平面每个点对应着系统的一个运动状态，这个点就称为相点，“相”是指物体的运动状态。

d sin()cos()d 2x A t A t tπωωφωωφ==-+=++v 222d d cos()d d x a A t ttωωφ===-+v相轨迹：相点随时间t 的变化在x x ∙- （公式）平面上描绘出的轨迹线，这种轨迹称为相轨迹，它表征了系统运动状态（相）的演变过程。

简而言之，相图上的每一条曲线表示在不同初始条件下，物体在相空间内的运动轨迹，曲线上的任意一点代表物体的某一运动状态。

2.4相轨迹作图方法相轨迹的作图方法可分为：解析法和图解法。

其中解析法主要针对相对比较简单的系统，比如简谐振动系统，可直接由动力学方程求出位移（x）与速度（x ∙）之间的关系的。

用解析法求相轨迹是比较麻烦的，特别是对非线性系统，有时可能无法求出相轨迹的解析表达式。

而图解法则主要针对不能直接由方程求出,x x ∙关系的系统，原则上说，此法对任何非线性系统都适用。

2.4.1用解析法求解简谐振子相图方程[5]广义定义：某个物理量随时间的变化是按正弦或余弦规律，则可称该振动为简谐振动，可用 表示。

速度(y)y x∙=2y x x ω∙∙∙==-2dy y xdxω=-2ydy xdxω=-方程两边积分得：22222yx Cω=-+ 简谐振子相轨迹方程：222x y Cω+=①式中ω为系统初始值，C 是由初始状态决定的常量 法二：cos()x A t ωφ=+cos()x A t ωφ=+d sin()d x A t tωωφ==-+v cos()x A t ωφ=+法一：位移（x ）222221x vAAω+=式中ω为系统初始条件，而常数22020nx x A ω +=，它是由初始状态决定的常量。

由简谐振动系统的相图方程式①②中，可以很显然地看出简谐振动的相轨迹是为闭合的椭圆。

2.5简谐振动能量以弹性系数为k ，质量为m 的弹簧振子为例,圆频率ω 满足2k m ω=2.6计算机辅助获得相图Fortran( Formula Translation System)可谓目前计算机运算中的程序语言之父，它是第一个能将数学公式转换成计算机程序的语言。

Fortran 擅长于数学函数运算，主要应用于数值分析、系统仿真及自动控制等领域。

自1954年，从第一个Fortran 程序诞生到现在接近50年啦，很长一段222211sin ()22k E m m A t υωωφ==+212k p E E E kA=+=消去tcos()x A t ωφ=+d sin()d x A t tωωφ==-+v 22211cos ()22p E kxk A t ωφ==+对于简谐振子系统来说，系统的总能量与振幅A2成正比。

对于系统来说，振动过程中任意时刻机械能守恒②时间内是科学计算语言的唯一选择，期间积累的大量的正确、可靠的程序，尽管Fortran标准改了多次，但由于其向下兼容，很多程序是招值即来，来之能战，并且相比于C语言，Fortran语言的纠错功能更强大，因而可算是个易掌握的好工具。

Origin是有着强大的数据分析和绘图功能的作图软件。

它的数据分析功能包括数据的排序、调整、计算、统计、频谱变换、曲线拟合等各种完善的数学分析功能，而且此软件操作简单，易于掌握。

计算机辅助获得相图的方法本质上是图解法获得相图。

从简谐振子的定义出发，在Fortran 环境下编程，理论上获得无数个模拟数值，然后用这些模拟数值在Origin软件下计算机描点作图。

相比与传统的图解获相图法，计算机辅助画相图效率更高，不仅能任意更改初始条件，而且所作的相图更精确[6]。

2.7用相空间研究完整系统的力学问题[11] [12] [13]每一条相轨迹对应在一段时间内力学系统的状态变化过程。

相轨迹是等能曲线，相应于一定能量的质点运动。

不同相轨迹，对应于不同能量的质点运动，不同能量的相轨迹是不可能相交的。

相空间的每一点称为相点，对应力学系统的一个瞬时状态。

而相图中给定一个相点，实际上等于给定了一组完备的初始定值，即给出完整的动力学系统。

用相空间的概念研究力学问题，实际上是从几何的方法出发研究力学。

它能把力学系统的全部信息在相图上表示出来，而无须去解复杂的运动方程。

三、总结简谐振动系统是一个可以精确求解的模型。

它的动力学系统对我们来说虽然说是完全已知的，但从相图角度研究简谐振动系统，仍展现了它具大的优势——无需解复杂的运动方程，更直观更形象。

本论文在Fortran 90编程环境下对简谐振子进行数值模拟，在Origin75 软件下画出相图。

给定任意的初始条件，画出不同相图。

四、参考文献[1] 刘明.线性谐振子问题研究[J].培训与研究——湖北教育学院学报,2004年9月,21(5)：15-16.[2] 陈建仁.用相空间研究完整系统的力学问题.河南教育学院学报(自然科学版) ,2004年3月,l3(1):25～27[3] 漆安慎.杜婵英.力学[M].高等教育出版社[M].2005年6月第2版[4] 周衍柏.理论力学教程[M].高等教育出版社,1986年3月第2版.282～330[5] 严燕来.叶庆好.大学物理拓展与应用[M].北京：高等教育出版社,2002年12月.49～54.[6] 杨正波.一个非线性力学问题的理论分析和数值模拟研究[J].襄樊学院学,2008年5月，第29卷(5):19～21[7] 丁泽军.Fortran77和90/95编程入门[M].中国科技大学天文与应用物理系,2001年10月[8] 刘卫国.FORTRAN 90 程序设计上机指导与习题选解[M].北京邮电大学出版社,2003年2月.165-18[9] 叶卫平.方安平.于本方.Origin 7.0科技绘图及数据分析[M].机械工业出版社,2004年1月第一版[10] 周建平.精通Origin 7.0[M].北京航空航天大学出版社,2004年3月第一版[11] 许雪芬.巢毅敏倒摆运动稳定性的相空间动力学分析[J].江苏油化学院学报,1999年3月,第11卷(1):43～45[12] 谢利民.弹簧振子运动的实际动力学分析[J].上海师范大学学报(自然科学版),2002年6月,第31卷(2):91～95[13] 克劳斯.迈因策尔.复杂性中的思维[M].中央编译出版社.。