时滞反馈作用下ENSO充电振子模型的分岔分析-应用数学和力学.pdf

　 ， ３９ １０ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１８年１０月１ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３９，Ｎｏ．１０，Ｏｃｔ．１，２０１８ 文章编号：１０００⁃０８８７（２０１８）１０⁃１１２８⁃０９ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 时滞反馈作用下 ＥＮＳＯ充电振子模型的 分 岔 分 析∗ 刘宇丹，　 杜智远，　 赵　 强 （北京大学 物理学院大气与海洋科学系，北京 １００８７１） 摘要：　 通过数学变换将一类含有时滞反馈机制的ＥＮＳＯ充电振子模型转换成时滞 Ｖａｎ ｄｅｒ Ｐｏｌ⁃ Ｄｕｆｆｉｎｇ方程，并以此为基础来研究该ＥＮＳＯ系统的零解稳定性、Ｈｏｐｆ分岔和极限环等动力学特征． 用平均法分析了其零解的稳定性与时滞效应的强度，其和时间都有关系，讨论了时滞负反馈对ＥＮ⁃ ＳＯ振荡的影响并通过简单数值模拟验证理论分析的结果． 关　 键　 词：　 ＥＮＳＯ；　 时滞Ｖａｎ ｄｅｒ Ｐｏｌ⁃Ｄｕｆｆｉｎｇ方程；　 Ｈｏｐｆ分岔 中图分类号：　 Ｐ４３３　 　 　 文献标志码：　 Ａ ＤＯＩ：１０．２１６５６／ １０００⁃０８８７．３８０３３２ 引　 　 言 ～ Ｅｌ Ｎｉｎｏ（厄尔尼诺）／ 南方涛动（ＥＮＳＯ）是大尺度海气相互作用（特别是热带海气耦合相 互作用）的集中表现，是年际气候变化中的最强信号，因此其成为认识年际气候变化规律和预 测年际气候变化的重要突破口．为了揭示ＥＮＳＯ现象的物理机制及其变化规律，科学家们在进 ［１⁃５］ 行观测资料统计和诊断分析、数值模拟的同时，也开展了ＥＮＳＯ机理的动力学理论研究 ．相 对于复杂的全球海⁃气耦合数值模式，通过简化海气非线性相互作用物理过程所得到的振子 概念模型能够更容易地刻画海气耦合运动本质和物理机理，从而更成功地模拟ＥＮＳＯ 的某些 重要物理现象．因此许多学者提出研究ＥＮＳＯ 的各种振子形式的动力系统模型，例如“时滞振 子” “平流⁃反射振子” “充电⁃放电振子”“西太平洋振子”，以及“时滞微分方程”理论来探索 ［６⁃１２］ ＥＮＳＯ动力学机理 ．另一方面，笔者也注意到ＶａｎｄｅｒＰｏｌ⁃Ｄｕｆｆｉｎｇ振子作为概念模型已应用 于耦合海冰模式、模拟高分辨率的冰芯资料、古气候动力学理论等研究［１３⁃１５］．由于时滞效应对 非线性系统稳定性的影响不可忽视［１６⁃１７］ ，本文中，用一个受到时滞反馈作用的Ｖａｎ ｄｅｒ Ｐｏｌ⁃ ［６］ Ｄｕｆｆｉｎｇ振子来研究时滞反馈效应对ＥＮＳＯ的影响．基于Ｊｉｎ 所提出的ＥＮＳＯ充电振子理论模 型，增加考虑时滞反馈作用机制，该ＥＮＳＯ振子模型的动力学控制方程组可以写成如下形式： ｄＴ ì ３ ï ＝ ＋ － ＋ ＣＴ Ｄｈ εＴ Δ， ïｄｔ 　 　 í （１） ｄｈ ï ＝－ － ＥＴ Ｒ ｈ， ïｄｔ ｈ î ∗ 收稿日期：　 ２０１７⁃１２⁃２２；修订日期：　 ２０１８⁃０２⁃２６ 基金项目：　 国家自然科学基金 作者简介：　 刘宇丹（１９８３—），女，硕士（通讯作者． Ｅ⁃ｍａｉｌ：ｖｅｔａ＿ｅｌｌｅ＠ｐｋｕ．ｅｄｕ．ｃｎ）．