不定期备考小tips[常微][2] #20210528

  • 保守系统(续)
    • 单摆
      • 画相图的典型示例
      • 周期大小
      • 圆柱上的相图
      • 带阻尼的情况
    • 考察保守性的例题
      • 进一步理解李雅普诺夫函数
  • 可逆系统
    • 例:随手例题
    • 例:力场
    • 可逆与保守

本专栏主要作个人笔记,有相关知识预备的同学也可作复习用。不保证无相应基础的人士能看明白。
万一考试考到了,或者对你的学习有较大帮助,一键三连不过分吧(斜眼笑)

保守系统(续)

单摆

画相图的典型示例

d2dt2θ=−glsinθ\frac{d^2}{dt^2} \theta = -\frac gl sin\thetadt2d2​θ=−lg​sinθ,取特征时间τ\tauτ满足ω=gl,τ=ωt\omega = \sqrt{\frac gl},\tau = \omega tω=lg​​,τ=ωt,则d2dτ2θ=θ¨=−sinθ\frac{d^2}{d\tau^2} \theta=\ddot \theta = -sin\thetadτ2d2​θ=θ¨=−sinθ
θ˙=v,v˙=−sinθ\dot \theta = v,\dot v = -sin\thetaθ˙=v,v˙=−sinθ,试画相图。

  1. 画nullcline,即所有θ˙=0\dot \theta=0θ˙=0的点和v˙=0\dot v = 0v˙=0的点。
  2. 标注每个区域处θ˙,v˙\dot \theta,\dot vθ˙,v˙的符号。
  3. 标注所有不动点,并对不动点作稳定性分析(可能手段:近似化,看图,看物理意义等)。
    举例:单摆例子中v=0,θ=πv=0,\theta=\piv=0,θ=π处显然不能稳定。(实际上是鞍点)
    对于线性近似后能求出特征方向的,可以求出特征方向,并辅助画图。
  4. 考察保守性、守恒量、闭轨。
    比如本例12θ˙2−cosθ\frac 12 \dot \theta^2-cos\theta21​θ˙2−cosθ守恒。(0,0)(0,0)(0,0)是12v2−cosθ\frac 12 v^2-cos\theta21​v2−cosθ的极小值点,根据保守场性质,闭轨存在。
    或者直接通过看等值线证明闭轨存在。
    或者通过对称性证明闭轨存在。
  5. 考察异宿轨(通过两个鞍点)、同宿轨(两个鞍点“合并”为一个时的异宿轨)等特殊轨道。对于本题,通过等值线或对称性都可说明通过(−π,0)(-\pi ,0)(−π,0)和(π,0)(\pi,0)(π,0)的异宿轨存在。
  6. 对称性在画相图中的作用
    确定相图形状,并且可以只画相图的一部分以画出其它部分。比如本题有中心对称、平移对称、θ\thetaθ和τ\tauτ同时反号的对称。这说明只需要考察第一象限0≤θ≤π0\le\theta\le\pi0≤θ≤π的情况即可画出整张图

周期大小

由12θ˙2−cosθ=C=−cosA\frac 12 \dot\theta^2-cos\theta=C=-cosA21​θ˙2−cosθ=C=−cosA得θ˙=2cosθ−2cosA\dot \theta =\sqrt{2cos\theta-2cosA}θ˙=2cosθ−2cosA​
利用对称性,容易计算得T(A)=C∫0Adθcosθ−cosAT(A)=C\int_0^A \frac{d\theta}{\sqrt{cos\theta-cosA}}T(A)=C∫0A​cosθ−cosA​dθ​,AAA是角度的振幅,CCC是一常数。
该积分的性质

  1. 当AAA很小,θ\thetaθ很小,线性近似下sinθ≈θsin\theta\approx \thetasinθ≈θ,原积分为C′∫0AdθA2−θ2=C′∫01dt1−t2=C′π/2C'\int_0^A \frac{d\theta}{\sqrt{A^2-\theta^2}}=C'\int_0^1 \frac{dt}{\sqrt{1-t^2}}=C'\pi/2C′∫0A​A2−θ2​dθ​=C′∫01​1−t2​dt​=C′π/2,则单摆有等时性,周期与振幅无关。
  2. 当AAA趋向于π\piπ时,相图中轨线趋向于异宿轨。回忆存在唯一性定理,推知“运动到平衡点”的所需时间发散(类比一维)。实际上容易验证A→πA\to \piA→π时原积分发散。

圆柱上的相图

由于θ\thetaθ变化2π2\pi2π在物理上表示一个点,可以将两侧粘合形成圆柱。
当能量E=12v2−cosθ>1E=\frac 12 v^2-cos\theta>1E=21​v2−cosθ>1,则物理意义是不断转圈(”进动“)不摆动,在平面相图上是不断向左或右一个方向运动,在圆柱上是不断绕柱作周期运动。
绕柱作周期运动的闭轨和前面平面上的闭轨具有不同的拓扑意义。
异宿轨在粘合后鞍点恰好合并了,成为同宿轨。

想象圆柱上的等能量线在设法弯曲圆柱后变为等高线,则直观看到:

E=0E=0E=0处有分岔现象。
E>0E>0E>0时看起来左右两情况旋转方向相同,但物理意义上旋转方向相反(把圆柱重新拉直可以看出)。

带阻尼的情况

θ¨+ϵθ˙+sinθ=0\ddot \theta + \epsilon \dot \theta +sin\theta=0θ¨+ϵθ˙+sinθ=0

在平面相图上出现螺旋。在柱面相图上的情况可以通过弯曲的柱直观想象。

考察保守性的例题

进一步理解李雅普诺夫函数

x˙=−y+ax3,y˙=x+ay3\dot x = -y+ax^3, \dot y=x+ay^3x˙=−y+ax3,y˙​=x+ay3
回顾:通过构造李雅普诺夫函数x2+y2x^2+y^2x2+y2,考虑其导数,证明a≠0a\neq 0a​=0时系统不保守。
并不是任何函数都能具有如此功能。比如重力场中高度也单调,同种电荷排斥的有心力场中角度也单调。
实际上,这里的x2+y2x^2+y^2x2+y2的导数真正说明的要点是:原点是排斥的或吸引的,与保守场性质矛盾。
具体地:例如a<0,ddt(x2+y2)=a(x4+y4)a<0, \frac d{dt} (x^2+y^2)=a(x^4+y^4)a<0,dtd​(x2+y2)=a(x4+y4),则x2+y2x^2+y^2x2+y2单减下有界,必有极限,且容易说明极限必须为0,则原点是吸引的。而若a>0a>0a>0就时间反演,证明原点是排斥的。

可逆系统

定义:x⃗˙=F⃗(x⃗)\dot {\vec x}=\vec F(\vec x)x˙=F(x),在变换t′=−t,x⃗′=R(x⃗)t'=-t,\vec x'=R(\vec x)t′=−t,x′=R(x)下不变,其中RRR是平面到平面的光滑双射且R(R(x⃗))=x⃗,∀x⃗R(R(\vec x))=\vec x,\forall \vec xR(R(x))=x,∀x(类比”对合矩阵“)

例:随手例题

x˙=−y+x3,y˙=x−y3\dot x = -y+x^3, \dot y=x-y^3x˙=−y+x3,y˙​=x−y3
容易求出零线y=x3y=x^3y=x3和x=y3x=y^3x=y3,不动点(0,0),(1,1),(−1,−1)(0,0),(1,1),(-1,-1)(0,0),(1,1),(−1,−1)
线性近似后:原点是中心。通过零线上的运动方向容易看出另外两个平衡点是鞍点。
线性近似前原点附近是否有闭轨?

  1. 可以直接交叉相乘相减x˙x−y3x˙+yy˙−x3y˙\dot x x-y^3\dot x+y\dot y-x^3\dot yx˙x−y3x˙+yy˙​−x3y˙​找守恒量,再利用保守场性质说明有中心、有闭轨。
  2. 可以换极坐标发现rr˙=r4cos2θ,(r−2)′=Ccos2θ,θ˙=1+小量r\dot r=r^4 cos2\theta,(r^{-2})'=Ccos2\theta,\dot\theta=1+小量rr˙=r4cos2θ,(r−2)′=Ccos2θ,θ˙=1+小量,根据rrr具有关于时间的周期性,容易说明附近不存在螺旋线。(实际上,线性近似后的中心在近似前的所有可能情况是能被完全分类的。若不存在螺旋,可以说明线性近似前是中心。具体阐述留待下一期)
  3. 或者可以利用体系关于x′=y,y′=x,t′=−tx'=y,y'=x,t'=-tx′=y,y′=x,t′=−t的对称性,直接考察y=xy=xy=x直线上点为起点的轨线,证明原点附近有一族闭轨。
  4. 和3.完全类似地:也可以利用关于y=−xy=-xy=−x直线的对称性。x′=−y,y′=−x,t′=−tx'=-y,y'=-x,t'=-tx′=−y,y′=−x,t′=−t后方程不变。现在考察x+y=0x+y=0x+y=0直线上点为起点的轨线,证明原点附近有一族闭轨。
    小复习:如何求关于直线y+x=0y+x=0y+x=0的反射变换:y+xy+xy+x变换后是−(y+x)-(y+x)−(y+x),x−yx-yx−y变换后不变,两式相加或相减发现xxx变为−y-y−y,yyy变为−x-x−x。
    注:方法3.在书写时的严格说法:若x=u(t),y=v(t)x=u(t),y=v(t)x=u(t),y=v(t)是解曲线,则x=v(−t),y=u(−t)x=v(-t),y=u(-t)x=v(−t),y=u(−t)也是解曲线,考察这两条解曲线就找到了闭轨。

例:力场

x¨=f(x)\ddot x=f(x)x¨=f(x)为奇函数,首先系统保守(可以写出x˙=y,y˙=f(x)\dot x=y,\dot y=f(x)x˙=y,y˙​=f(x)交叉相乘相减证明)。
t′=−tt'=-tt′=−t时,根据物理含义,y′=−yy'=-yy′=−y。t′=−t,y′=−yt'=-t,y'=-yt′=−t,y′=−y下方程不变,也就是物理规律不变

可逆与保守

  • 联系:都有中心与”一族闭轨“的结论。
    实际上可逆系统中,孤立的、线性近似后是中心的不动点周围都有一族闭轨。(重要补充说明:这里说到”周围“。其实这个命题是一个“局部性质”,要求不动点的邻域D映到自身,特别地,该孤立不动点在映射下不变。并不是所有满足原始定义的可逆系统都有此结论)
    具体证明略。根本思路是x⃗\vec xx运动到R(x⃗)R(\vec x)R(x)的解曲线”对称后“仍是解曲线。
  • 典型的可逆但不保守的例子:x˙=−2cosx−cosy,y˙=−2cosy−cosx\dot x = -2cosx-cosy, \dot y=-2cosy-cosxx˙=−2cosx−cosy,y˙​=−2cosy−cosx。构造该反例的思想:先构造出一类”足够大“的可逆系统,再从该类中试图找不保守的(比如找出其中有结点的系统)。实际上,容易知道x˙=f(x)+g(y),y˙=h(x)+i(y)\dot x=f(x)+g(y),\dot y = h(x)+i(y)x˙=f(x)+g(y),y˙​=h(x)+i(y),其中各f,g,h,if,g,h,if,g,h,i都是偶函数,则系统在t′=−t,x′=−x,y′=−yt'=-t,x'=-x,y'=-yt′=−t,x′=−x,y′=−y下不变。那么在各类偶函数中找出适当的,使得结点存在,就很简单。

不定期备考小tips[常微][2] #20210528相关推荐

  1. 不定期备考小tips[数模][0] #20210529

    不定期备考小tips[数模][0] #20210529 边值问题 S-L微分方程中两种内积的理解 一种内积:(u,v)=∫uvdx(u,v)=\int uvdx(u,v)=∫uvdx 另一种内积:(u ...

  2. 让你的网站更炫酷的一些小 tips

    本文讲的是让你的网站更炫酷的一些小 tips, 上周,我和一位老客户聊天,她说:"尼克,我觉得我的网站需要改进,但我不能确定我具体需要做什么." 然后我就去问了一圈,包括朋友.家人 ...

  3. tdm的应用计算机,2020计算机考研:TDM时分复用技术备考小知识点

    最近我们常听到一句话"得暑假者得考研",暑期将近,考研的你准备好了吗?暑假是考研的黄金备考期,很多同学已经忙着打基础或者赶进度.那么作为计算机考研的你,在备课计算机网络的过程中是否 ...

  4. LCD液晶屏的使用小tips

    LCD液晶屏和人类的生存发展越来越紧密,基本上家家户户都有一些液晶屏制品的存在,而那些还没有应用到人类生活中的"初始状态"液晶屏,你知道脆弱的它们,该如何小心使用吗?今天就对LCD ...

  5. 本人常用的一些编码小Tips(虽然不多,但很好用)

    本人常用的一些编码小Tips 我是可爱的目录 本人常用的一些编码小Tips 一.概述 二.Windows的一些常用的快捷键 (一).Win + others 类型 (二).Alt + others 类 ...

  6. php childnodes,小tips:HTML DOM中的children和childNodes属性

    childNodes 属性 标准的,childNodes 属性返回节点的子节点集合,以 NodeList 对象.包括HTML节点,所有属性,文本.可以通过nodeType来判断是哪种类型的节点,只有当 ...

  7. 小TIPS:合拼gridview的header

    http://www.cnblogs.com/jackyrong/archive/2006/08/15/477791.html 小TIPS:合拼gridview的header 在asp.net 2.0 ...

  8. 阿里高级技术专家邱小侠:微服务架构的理论基础 - 康威定律

    邱小侠 阿里高级技术专家 读完需要 10 分钟 速读仅需 4 分钟 邱小侠,阿里巴巴集团客户体验事业群高级技术专家,阿里花名肥侠.2014年加入阿里巴巴,现在负责客户体验驱动及创新中心有关商家业务的开 ...

  9. 小tips:JS之浅拷贝与深拷贝

    浅拷贝: function extendCopy(p) {var c = {};for (var i in p) {c[i] = p[i];}return c; } 深拷贝: function dee ...

  10. 商城前端模板_如何理解微信小程序和微商城,微信公众号以及APP之间的关系?一张图看懂了!...

    老张的一位粉丝,花了几天时间把知乎里面分享的一些关于微信小程序,微信商城,微信公众号,以及APP的相关介绍全看完了. 然后用他自己的话描述了微信小程序和微商城,微信公众号以及APP之间的关系,如下图所 ...

最新文章

  1. Java创建初始化List集合的几种方式
  2. linux kbhit扫描键盘,(转)检测按键(Linux中kbhit()函数的实现)
  3. jQuery 常用的效果函数(一)
  4. php 返回map,PHP Ds\Map get()用法及代码示例
  5. mvc php session,PHP Session入门教程
  6. Ask Me Anything #1 我是新晋CNCF TOC张磊,你有什么想问我的?
  7. c++图形中如何判断鼠标点击在一条直线上_中考数学常考题型精讲精练系列:函数图象上点的存在性问题中的距离与面积...
  8. angular input使用输入框filter格式化日期
  9. windows下使用cmake+mingw配置makefile(1)
  10. 《机器学习实战》KNN算法实现
  11. CentOS7.x编译安装nginx,实现HTTP2
  12. 充电速度公式_锂电充电时间计算公式
  13. android当电脑麦克风,电脑没有麦克风?让手机充当电脑麦克风!
  14. Linux下安装各种常用软件
  15. linux skb机制,skb 的分配细节
  16. 前端JS 实现将24位RGB颜色转换16位RGB颜色
  17. Django中引入bootstrap的方法
  18. 手机在信号好的地方一直无服务器,解决手机信号不好的几种方法
  19. tensorflow:自定义op
  20. 【图像识别】基于计算机视觉实现红绿灯识别含Matlab代码

热门文章

  1. cdr多页面排版_CDR排版须掌握三大功能 值得收藏
  2. cdr软件百度百科_cdr是什么软件?
  3. JAVA自动生成前缀加日期的编码
  4. 显卡、GPU和CUDA的概念整理
  5. HTML结构:自我简介网页
  6. 程序员3年工资从7千到2万2,晒出跳槽经历!网友:厉害
  7. font在html语言中是什么意思,font标签什么意思
  8. 结对编程-李晓彤+凌嘉文
  9. 蓝桥杯算法题--跑步锻炼(java实现)
  10. python zip 压缩文件夹