庞加莱映射

先对庞加莱映射作一简介,为了更清楚地了解运动的形态,庞加莱对连续运动的轨迹用一个截面(叫庞加莱截面)将其横截,那么根据 轨迹在截面上穿过的情况,就可以简洁地判断运动的形态,由此所得图像叫庞加莱映像。在截面图上,轨迹下一次穿过截面的点 可 以看成前一次穿过的点 的 一种映射

?? (n=0,1,2,…)

这个映射就叫庞加莱映射。它把一个连续的运动化为简洁的离散映射来研究

在庞加莱映射中的不动点反映了相空 间的周期运动,如果运动是二倍周期的,则庞加莱映射是两个不动点,四倍周期则有四个不动点等

绘制庞加莱映射是在普通的相平面上进行,它不是像画相轨道那样随时间变化连续地画出相点,而是每隔一个外激励周期( ) 取一个点,例如取样的时刻可以是t=0,T,2T…相应的相点记为 , , … 这些离散相点就构成了庞加莱映射

自激振动

1.实验题目

研究范·德·波耳(Van der pol) 方程

??????????????????????? (2.19.1)

所描述的非线性有阻尼的自激振动系统,其中 是 一个小的正的参量, 是 常数。下面简称范·德·波耳方程为VDP方程

在VDP方程中,增加外驱动力 项 所得到的方程

????????????????? (2.19.2)

称强迫VDP方程,其中外驱动力的振幅,角频率分别是V和 ,试研究强迫VDP方程的行为

2.实验目的和要求

⑴演示VDP方程所描述的系统在非线性能源供给下,从任意初始条件出发都能产生稳 定的周期性运动

⑵采用庞加莱映像,演示强迫VDP方程在不同参数下所存在四种吸引子,即周期1吸引子,周期2吸引子,不变环面吸引子和奇怪吸引子

⑶对于强迫VDP方程,在V和 为 定值条件下,逐渐增大 值, 将出现周期倍分岔和混浊现象

3.解题分析

自激系统是一个非线性有阻尼的振动系统,在运动过程中伴随有能量损耗,但系统存在一种机制,使能量能够由非振动的能源通过系统本身的反馈 调节,及时适量地得到补充,从而产生一个稳定的不衰减的周期运动,这样的振动称为自激振动

对VDP方程,可从机械振动角度理解, 是 阻尼系数,它是变化的,如果 , 则阻尼系数为正,系统将受阻尼,能量将逐渐减少,但如果 , 则发生负阻尼,意味着不仅不消耗系统的能量,反而给系统提供能量。此系统能通过自动的反馈调节,使得在一个振动过程中,补充的能量正好等于消耗的能量,从 而系统作稳定的周期振动

取方程中的 , , (这 些值可适当调整)。给出任一初始条件,通过计算机数值求解可以证明它的相轨道都将趋向于一条闭合 曲线,这一条闭合曲线,成为极限环,极限环以外的相轨道向里盘旋,而极限环以内的相轨道则向外盘旋,都趋向极限环(如图2.36所示),说明不论初始情况如何,系统最终都到达以极限环描述的周期性运动。由于这段程序较简单,我们没 有专门编写,事实上,只要将下面编写的关于强迫VDP方程的程序中令V=0, 再 取不同的初始条件,就能看到这个现象

?

?

下面研究强迫VDP方程的行为,我们同时采用时间历程图,相图,庞加莱映像图来研究系统在不同参数条件下的动力学行为,可以看到存在不同的吸引子,即周期1吸引子,周期2吸引子,不变环面吸引子和奇怪吸引子

先对庞加莱映射作一简介,为了更清楚地了解运动的形态,庞加莱对连续运动的轨迹 用一个截面(叫庞加莱截面)将其横截,那么根据轨迹在截面上穿过的情况,就可以简洁地判断运动的形态,由此所得图像叫庞加莱映像。在截面图上,轨迹下一次 穿过截面的点 可 以看成前一次穿过的点 的 一种映射

?? (n=0,1,2,…)

这个映射就叫庞加莱映射。它把一个连续的运动化为简洁的离散映射来研究

在庞加莱映射中的不动点反映了相空间的周期运动,如果运动是二倍周期的,则庞加莱映射是两个不动点,四倍周期则有四个不动点等

绘制庞加莱映射是在普通的相平面上进行,它不是像画相轨道那样随时间变化连续地画出相点,而是每隔一个外激励周期( ) 取一个点,例如取样的时刻可以是t=0,T,2T…相应的相点记为 , , … 这些离散相点就构成了庞加莱映射

设 , , 则(2.19.2)式可 化为

??????????????????????????????????????

?????????????? ? (2.19.3)

取 , , 进行以下数值计算研究

⑴在 ,V=1, 条 件下,存在周期1吸引子,它的周期等于外激励的周期,代表主谐波运动,如图2.37所示

⑵在 ,V=1, 条 件下,存在周期2吸引子,它的周期等于外激励的整数倍, 代表次谐波运动,如图2.38所示

⑶在 ,V=1, 条 件下,存在不变环面吸引子,它代表准周期(拟周期)运动,如图2.39所示

⑷ ,V=1, 条 件下,存在奇怪吸引子,它代表混浊运动,如图2.40所示

⑸保持V和 为

庞加莱映射导数matlab,庞加莱映射.doc相关推荐

  1. 微分matlab程序,数值微分及其matlab程序.doc

    数值微分及其matlab程序.doc 第八章数值微分1282一阶导数的数值计算及其MATLAB程序821差商求导及其MATLAB程序例821设215SINXF(1)分别利用前差公式和后差公式计算的 ...

  2. matlab数值分析拟合实例,数值分析函数拟合matlab代码.doc

    数值分析函数拟合matlab代码.doc 第一题MATLAB代码用SPLINE作图XI0204060810YI098092081064038X10012Y1NEWTON3XI,YI,X源代码见M文件Y ...

  3. matlab 二阶非线性微分方程组,二阶非线性常微分方程的打靶法matlab实现.doc

    二阶非线性常微分方程的打靶法matlab实现.doc 二阶非线性常微分方程的打靶法1.问题:试用打靶法求二阶非线性常微分方程亮点边值的数值解:要求用Matlab编程计算,请给出一些例子,验证你的算法与 ...

  4. 升降压斩波电路matlab,升降压斩波电路matlab仿真.doc

    <升降压斩波电路matlab仿真.doc>由会员分享,可在线阅读,更多相关<升降压斩波电路matlab仿真.doc(4页珍藏版)>请在人人文库网上搜索. 1.升降压斩波电路ma ...

  5. matlab分类器设计,简单分类器的MATLAB实现.doc

    简单分类器的MATLAB实现.doc 简单分类器的MATLAB实现 摘要:本实验运用最小距离法.Fisher线形判别法.朴素贝叶斯法.K近邻法四种模式识别中最简单的方法处理两维两类别的识别问题,最后对 ...

  6. matlab微分的语句格式,偏微分差分四种格式的matlab程序.doc

    偏微分差分四种格式的matlab程序.doc 偏微分差分四种格式的matlab程序 (1)Lax-Fridrichs格式 function u = LaxFridrichs(a,dt,n,x1,x2, ...

  7. matlab转向梯形优化设计,转向梯形优化设计matlab程序.doc

    转向梯形优化设计matlab程序.doc -WORD格式-可编辑-优化计算 MATLAB 程序首先,将目标函数写成 M 文件,其程序语句如下: function f = fun (x)global K ...

  8. 自动控制matlab实验,自动控制matlab实验.doc

    自动控制matlab实验.doc 实验一MATLAB的基本应用1.试验目的熟悉MATLAB的基本知识,会利用其与控制系统有关的函数分析和设计控制系统.2.实验内容(1)系统模型如下所示,判断系统的稳定 ...

  9. mlse均衡器matlab,基于LMS自适应均衡器matlab仿真.doc

    基于LMS自适应均衡器matlab仿真.doc 毕 业 设 计 (2014届) 题 目 一种基于OpenCV的摄像机标定方法 学 院 物理电气信息学院 专 业 电子信息工程 年 级 2010 学生学号 ...

最新文章

  1. 再见,工资!程序员工资统计平均14404元,网友:又跌了!
  2. android从放弃到精通 第一天 重拾项目开发
  3. shell中字符串基本用法
  4. python流行的爬虫框架_Python爬虫相关框架
  5. JavaScript中getBoundingClientRect()方法详解
  6. html匹配属性正则表达式,正则表达式匹配html标签的属性值
  7. 关于计算机的英语作文八百字,小学英语作文800字(精选10篇)
  8. java大数据组件HBase
  9. 扩展CheckBoxList实现选中绑定
  10. python 文本向量化_TFIDF向量化文档时出现问题?
  11. Python 修改图片的时候抗锯齿
  12. 文后参考文献著录规则 GB/T 7714-2005
  13. DL for Scratch 读书笔记
  14. 【LeetCode-面试算法经典-Java实现】【012-Integer to Roman(数字转罗马字符)】
  15. 卡巴斯基实验室解析勒索软件的发展与攻防
  16. 域用户登陆显示计算机图标,PC用户或域用户登陆后任务栏无显示桌面的小图标...
  17. NLP中Tokenizers总结(BPE、WordPiece、Unigram和SentencePiece)
  18. 如何复制权限受限PDF文件的内容(亲测有效,Microsoft Edge打开pdf文件)
  19. 女孩子适合软件测试这个行业吗?【工作内容、薪资、加班、怎么转行、职业规划】全面解析女生适不适合软件测试。
  20. python k线斜率计算公式_python – 计算Numpy(或Scipy)中的斜率

热门文章

  1. IPIC2020物界科技田民:Physical Internet物流和供应链领域的颠覆性创新
  2. IOS端K线系列之K线、OHLC线、分时线基础知识
  3. 使用Maya2018快速制作气球的简单漂浮动力学
  4. 安信可PB-01/02蓝牙模组实现远程OTA无线升级功能介绍,剖析整个实现原理和代码介绍。
  5. ESP8266 + MQTT固件(四博智联、安信可)
  6. 口腔ct重建服务器原理,基于口腔CT重建数据的牙齿分割与可视化研究
  7. 吃鸡专用计算机,死得不明不白 推荐你一台吃鸡专用显示器
  8. PHP 开发 APP 接口 学习笔记与总结 - 静态缓存
  9. SSM在线购物商城——电子产品商城(SSM毕业设计)
  10. Linux总结(2)