《从非线性动力学到复杂系统》
- 格式:docx
- 大小:92.40 KB
- 文档页数:8
非线性动力学中的生态系统变化随着环境问题的严重性日益凸显,生态系统的变化也成为了人们关注的焦点。
生态系统本身就具有高度的复杂性和非线性特征,因此非线性动力学的理论和方法在生态系统研究中得到了广泛应用。
本文将介绍非线性动力学在生态系统变化研究中的应用和意义。
一、生态系统变化的非线性特征生态系统是以生物群体、环境、营养物质等为要素的生态整体,具有高度的复杂性和非线性特征。
实际上,生态系统中的各个要素之间相互作用,彼此牵引,共同构成一个巨大的生态网。
在这个生态网中,各种生物之间相互制约、相互利用、相互竞争,对环境进行各种适应和调节。
这些复杂的相互作用一旦发生变化,就会引起不可逆的生态系统变化。
由于生态系统中各种生物之间的相互作用是动态和非线性的,因此生态系统表现出了很强的非线性特征。
非线性动力学正是一种用于描述这种复杂系统动力学行为的数学方法。
而当我们尝试建立一些简单的数学模型来描述生态系统时,这种非线性特征就会显现出来。
比如典型的食物链模型,这其中每个物种的增长速率都是受到其食物来源(进食物种)的限制。
但是,由于食物来源的数量也是受到其天敌的控制的,当这个网络中出现某个环节崩塌的时候,或是某个环节增加过大,都会引发一系列连锁反应。
这些反应的无效时间也很长,同时也都受到每个种群所受到的生态环境的影响与制约,还原整个系统的稳定状态变得异常复杂,以至于在这种非线性动力学的控制下,很难对系统进行精细地调节和控制。
二、基于非线性动力学方法的生态系统模拟与传统生态学建模方法相比,基于非线性动力学方法的生态系统模拟更加接近于自然生态系统的真实情况。
在该模拟过程中,我们可以用数学模型描述生态系统中不同生物之间相互作用的动态关系,并通过数值模拟来模拟生态系统的未来演化。
例如,在模拟森林群落动态变化时,可以建立一个基于生境影响的种群增长模型,然后使用不同的参数来描述环境因素对每个种群的影响。
再加上对其生殖行为、迁徙和竞争等众多因素进行系统建模,就可以在不同情境下,对系统的稳定性和动态性进行建模和预测,从而为生态管理和保护提供理论支持。
非线性动力学方法在系统科学中的应用
随着科学研究的不断深入,人们对于简单线性系统进行研究的
时代已经过去了。
现如今,具有非线性复杂性特征的系统科学越
来越得到人们的关注。
因此,非线性动力学方法应运而生,它可
以帮助我们更好地理解和建模这些复杂系统。
何谓非线性动力学?
首先,我们需要明确非线性动力学的定义。
非线性动力学是一
种研究非线性现象这种现象在某种程度上可能表现出如下结构特征:当系统的初始状态发生微小变化时,系统的响应也会发生显
著的变化,表现为不可逆、不确定、异构、分形、混沌、奇点等。
非线性动力学和系统科学
非线性动力学方法在系统科学中的应用越来越普遍,因为它可
以帮助我们更好地了解和解释复杂系统的行为。
作为一种独立的
学科,非线性动力学涵盖了许多领域。
它可以被应用于系统科学
的许多方面,包括生物学、化学、物理学、经济学、社会科学等。
应用于实践
非线性动力学方法的应用远不止于理论研究,它也为实践活动
提供了很多实用的方法。
例如,非线性动力学方法可以用于智能
控制、自适应控制、优化设计。
还可以被用于金融预测、管道损
伤检测、信号处理、图像分析等领域。
总结和展望
在系统科学中,非线性动力学方法正逐渐成为一种重要的工具。
它可以帮助我们更好地理解和掌握复杂系统的行为,为科学家和
工程师们提供全新的视角。
未来,随着复杂系统研究的深入,非
线性动力学方法的应用将会变得越来越广泛,对于解决复杂问题
和实现科学进步将具有重要的作用。
数学中的复杂系统和非线性动力学随着科学技术的不断发展,人们对于复杂系统的研究越来越感兴趣。
复杂系统是由许多相互联系和作用的元素组成,这些元素之间的关系可能会随时间变化而产生变化。
但是,这些系统往往比较难以理解,因为它们的行为具有不可预测性和不可确定性。
在数学中,复杂系统的研究和分析是一门重要的学科,而这门学科的一个核心内容就是非线性动力学。
在数学中,动力学是一门研究物体运动规律和光、电、热能的传递规律的学科。
线性动力学研究的是线性系统,这些系统中的元素之间的关系是可以通过线性方程来描述的。
然而,非线性动力学研究的是非线性系统,这些系统中的元素之间的关系无法通过线性方程来描述。
非线性系统往往具有复杂的行为和不可预测性。
例如,天气系统就是一个典型的非线性系统。
虽然我们可以通过观测气温、湿度、气压等指标来预测未来的天气趋势,但是由于天气系统是一个非线性系统,即使是微小的扰动也可能导致非常不同寻常的天气现象的出现。
这也是为什么我们在气象预报中经常听到关于概率和范围的预测,而不是绝对准确的预测。
其他的非线性系统包括心脏病理系统、人口系统和经济系统等。
非线性动力学的研究很大程度上集中在混沌理论的探究上。
混沌理论研究的是一种决定论的系统,这种系统虽然是确定性的,但是由于初始条件的微小差异,也可能会产生完全不同的演化结果。
换句话说,混沌系统可能会在看似随机的行为模式之间不断跳跃,这就是所谓的混沌现象。
混沌现象可以用一种叫做“迭代函数系统”的方式来描述。
迭代函数系统使用一系列非线性方程来描述非线性系统,这些方程通常是相互联系的,并且会随着时间变化而产生变化。
在迭代函数系统中,我们可以通过输入初始条件并反复应用这些方程来计算系统的演化。
因此,迭代函数系统和计算机模拟合作是研究非线性系统和混沌现象的重要工具。
虽然非线性动力学和混沌理论可能听起来很复杂,但是它们也具有实际的应用。
例如,在金融市场中,非线性动力学可以用来预测市场趋势和风险。
物理学中的复杂系统与动力学在我们所生活的这个世界中,物理学无疑是一门揭示自然规律和奥秘的重要学科。
而在物理学的广袤领域里,复杂系统与动力学则像是一座神秘而迷人的迷宫,吸引着无数科学家去探索和求解。
复杂系统,顾名思义,是由大量相互作用的成分组成的系统。
这些成分之间的相互作用并非简单的线性叠加,而是呈现出复杂的非线性关系。
这种非线性使得系统的行为难以预测,充满了不确定性和多样性。
比如生态系统,其中包含了无数的生物物种,它们之间相互依存、竞争、共生,形成了一个复杂的网络。
再比如社会系统,由众多的个体和组织构成,个体之间的交流、合作与冲突,共同塑造了社会的发展和变迁。
动力学则关注的是系统随时间的变化和演化。
在复杂系统中,动力学的表现形式更是丰富多样。
有些系统可能会在一段时间内保持相对稳定的状态,但在某些条件的触发下,会突然发生剧烈的变化,进入一个全新的状态。
这种从一种稳定态到另一种稳定态的转变,被称为相变。
想象一下,水在不同的温度和压力条件下,会从固态的冰转变为液态的水,再转变为气态的水蒸气。
这就是一个典型的相变过程。
在这个过程中,温度和压力等条件的微小变化,都可能导致系统状态的巨大改变。
而在复杂系统中,类似的相变现象更是屡见不鲜。
复杂系统的动力学行为往往具有自组织的特点。
也就是说,系统在没有外部指令的情况下,能够通过内部成分之间的相互作用,自发地形成某种有序的结构或模式。
比如鸟群在空中飞行时,没有一个统一的指挥者,但它们却能自发地形成整齐的队列;再比如大脑中的神经元网络,通过相互之间的电信号传递和连接,能够产生复杂的思维和认知活动。
那么,研究复杂系统与动力学对于我们有什么意义呢?首先,它有助于我们更好地理解自然界中的各种现象。
从气候变化到生物进化,从地震的发生到星系的形成,几乎所有的自然现象都可以看作是复杂系统的动力学过程。
通过研究这些过程,我们能够更深入地认识自然规律,为预测和应对自然灾害提供理论依据。
复杂系统的非线性动力学分析复杂系统是指由许多相互作用的要素构成的系统,在很多领域都有广泛的应用,例如生物、天气、社会和经济等。
复杂系统具有许多特点,其中一个重要的特点就是非线性。
在非线性系统中,当系统受到外部扰动或内部变化时,系统的响应是不可预测的,因此非线性动力学的研究受到越来越多的关注。
非线性动力学研究的重点是在非线性系统中探索混沌、周期和奇异吸引子等动力学现象。
为了更好地理解这些现象,我们需要掌握一些基本的概念和数学工具。
1. 相空间和相轨道相空间是指由系统的所有可能状态组成的空间,在非线性系统中,相空间通常是高维的。
在相空间中,一个状态可以用一个点来表示。
当系统的初始状态确定后,系统将按照一定的动力学规律在相空间中运动,这条轨道就是相轨道。
2. 相平面和自由度相平面是相空间的一个子空间,它只包含与某个变量有关的状态。
在非线性系统中,自由度指的是系统中独立的有意义的变量的个数。
3. 动力学方程动力学方程是描述系统运动规律的数学公式。
在非线性系统中,由于存在相互作用和非线性效应,动力学方程通常是高阶的、非线性的微分方程组。
4. 混沌吸引子混沌吸引子是非线性系统中的一种特殊的吸引子,指的是在相空间中具有分维的奇异吸引子。
混沌吸引子的出现是混沌现象的重要特征,它是一种长期的、不可预测的振荡行为。
5. 可控性和可观测性可控性和可观测性是非线性动力学研究中的重要概念。
可控性指的是系统是否能被外部控制,可观测性指的是系统是否能被观测。
这两个概念对于控制和预测非线性系统的行为至关重要。
6. 分岔理论分岔理论是非线性动力学的重要分支之一,它研究的是系统在参数变化时的运动规律。
在分岔理论中,分岔点是一个重要的概念,指的是当系统的参数变化到一定程度时,系统的运动规律发生突然变化。
总之,非线性动力学是一门重要而复杂的学科,它涉及到数学、物理、化学、生物和工程等多个领域。
随着科技的不断发展,非线性动力学的研究将在许多领域发挥重要的作用。
复杂性科学理论综述在当代科学领域中,复杂性科学逐渐成为一门独立的学科。
复杂性科学研究的是那些由大量相互作用的个体组成的系统,这些系统表现出非线性、自组织和难以预测的特征。
它的研究对象包括自然界中的生态系统、大脑、气候系统,以及社会经济系统等。
本文将综述复杂性科学的理论发展、应用领域以及未来的研究方向。
复杂性科学的理论发展可以追溯到20世纪40年代的系统论研究。
系统论强调整体性思维,将系统看作一个整体,由各个子系统相互作用而成。
然而,系统论的主要局限是对复杂系统的刻画过于简单,缺乏对系统内部的动力学和复杂性的深入理解。
为了克服这一限制,复杂性科学成为了一个新兴的领域。
复杂性科学的核心理论之一是复杂网络理论。
复杂网络理论从网络的结构、性质和动态演化等方面研究网络系统的特点。
网络由节点和连接边组成,可以用来描述身份关系网络、社交网络以及脑神经网络等。
复杂网络理论通过度分布、聚集系数和小世界结构等指标来研究网络的特性,揭示了网络系统的规模自相似性和无标度特性。
另一个重要的理论是非线性动力学,它研究的是复杂系统中的非线性行为。
复杂系统常常表现出非线性响应,这意味着系统的行为是非线性的,并且可能出现周期性、混沌和自组织等特征。
通过非线性动力学的方法,可以揭示系统内在的关联和相互作用,预测系统的行为,并解释系统中的复杂现象。
另外,复杂性科学还借鉴了信息论和统计物理学的方法。
信息论提供了熵、互信息和复杂度等指标,用于度量和量化系统的复杂性。
统计物理学则将统计学的方法引入到复杂系统的研究中,通过模拟和建模来解析系统的行为。
这些方法使得研究人员可以通过收集和分析大量数据来揭示系统的内在规律和特征。
复杂性科学的应用领域广泛。
在生态学领域,复杂性科学被应用于生态系统的保护和管理中。
研究人员通过对生物群落结构、物种相互作用和食物网等复杂网络的研究,揭示了物种灭绝的模式和传染病的传播机制。
在社会科学领域,复杂性科学可以帮助我们理解城市的增长与发展、社交网络的形成和演化。
非线性动力学系统深度研究深度研究非线性动力学系统引言:非线性动力学系统是一类常见的复杂系统,广泛应用于物理、化学、生物学等领域。
与线性系统相比,非线性系统具有更为复杂的行为和动力学特性。
本文将对非线性动力学系统进行深度研究,探讨其定义、模型、稳定性和混沌等关键概念。
一、非线性动力学系统的定义和基本概念非线性动力学系统是指系统中的状态变量和控制参数之间的关系是非线性的系统。
其基本概念主要包括状态变量、动力学方程和相空间等。
1. 状态变量:状态变量是系统的内部变量,它们描述了系统在不同时间的状态。
通常采用向量形式表示,例如(x1, x2, ..., xn)。
2. 动力学方程:动力学方程是描述系统演化规律的数学方程。
对于非线性动力学系统,动力学方程通常是一组非线性微分方程或差分方程。
3. 相空间:相空间描述了非线性动力学系统的所有可能状态的集合。
在相空间中,每个状态被表示为一个点,而系统的演化则对应于在相空间中的运动轨迹。
二、非线性动力学系统的模型与常见例子非线性动力学系统的模型通常采用一组微分方程或差分方程来描述。
下面给出两个常见的非线性动力学系统模型。
1. Lorenz系统:Lorenz系统是一个三维非线性动力学系统,由爱德华·洛伦兹发展而来,主要用于描述大气环流的运动。
Lorenz系统的动力学方程如下:dx/dt = σ(y - x)dy/dt = ρx - y - xzdz/dt = xy - βz其中,x、y、z分别表示系统的三个状态变量,σ、ρ、β分别为控制参数。
2. Van der Pol振荡器:Van der Pol振荡器是一个二阶非线性动力学系统,广泛应用于电子工程和生物学中。
其动力学方程如下:d²x/dt² - μ(1 - x²)dx/dt + x = 0其中,x表示系统的状态变量,μ为控制参数。
三、非线性动力学系统的稳定性分析在研究非线性动力学系统时,稳定性是一个关键问题。
复杂系统中的非线性动力学我们身处的世界是一个充满复杂性的系统,无论是自然界还是人类社会,都存在着众多相互作用的要素。
在这个复杂系统中,非线性动力学起着至关重要的作用。
本文将从理论和实践的角度探讨复杂系统中的非线性动力学,并分析其对我们理解和应对现实世界的意义。
非线性动力学是研究系统行为的一门学科,它关注的是复杂系统中各个要素之间的相互作用及其对整体系统行为的影响。
与线性系统不同,非线性系统的响应不能简单地通过线性关系来描述。
相反,非线性系统的行为是由其内部相互作用及外部环境的影响所决定的。
在自然界中,非线性动力学能够帮助我们更好地理解各种自然现象。
例如,气象学中的天气系统往往是非线性的,气象变量之间互相影响,因此我们无法准确预测天气的未来发展。
地球上的生态系统也是一个复杂而非线性的系统,各种物种之间的相互关系以及环境因素的变化都会对整个生态系统产生深远的影响。
非线性动力学的应用可以帮助我们更好地理解这些自然现象,并且为环境保护和资源管理提供指导。
在人类社会中,非线性动力学同样发挥着重要的作用。
社会系统是一个巨大而复杂的网络,各种社会要素之间的相互作用导致了社会的演化和变化。
研究社会系统中的非线性动力学可以帮助我们更好地理解经济、人口、文化等方面的问题。
例如,经济学中的经济波动往往具有非线性的特点,各种经济变量之间的相互作用导致了经济的不稳定性和突发性的变化。
非线性动力学的方法可以帮助我们对经济现象进行建模和预测,从而更好地应对经济风险。
实际上,非线性动力学的应用还延伸到了许多其他领域。
生物学、化学、物理学等学科都在研究非线性系统中的动力学行为。
例如,神经科学领域中的脑电波和神经网络的研究都需要考虑非线性动力学的影响。
此外,还有许多工程领域使用非线性动力学的理论和方法来研究各种工程系统的行为,包括航空航天、能源、交通运输等。
然而,非线性动力学也带来了挑战和困难。
非线性系统的行为往往复杂且难以预测,需要借助数学模型和计算方法来进行研究。
复杂系统的非线性动力学模型引言复杂系统是由大量相互作用的组件构成的系统,其行为可能表现出非线性特征。
非线性动力学模型是描述复杂系统中非线性行为的数学工具。
本文将介绍复杂系统的特点、非线性动力学模型的基本原理以及其在实际应用中的重要性。
复杂系统的特点复杂系统具有以下几个特点:1.由多个相互作用的组件组成:复杂系统由许多相互作用的组件构成,例如生态系统中的物种、社交网络中的个体等。
2.非线性关系:复杂系统中的组件之间存在非线性关系,即系统的整体行为不能简单由各个组件的线性叠加得到。
3.自组织性:复杂系统具有自组织的能力,即系统中的组件可以通过相互作用形成新的结构和行为。
4.非平衡状态:复杂系统处于非平衡状态,其内部和外部的能量和物质交换使得系统的状态不断变化。
非线性动力学模型的基本原理非线性动力学模型是描述复杂系统中非线性行为的数学模型。
它基于以下几个基本原理:1.状态变量和演化规律:非线性动力学模型使用状态变量来描述系统的状态,并通过演化规律描述状态变量随时间的演化过程。
2.非线性关系:非线性动力学模型中的演化规律包含非线性关系,即状态变量之间的相互作用具有非线性特征。
3.启动条件和边界条件:非线性动力学模型需要指定适当的启动条件和边界条件来确定系统的起始状态和外部约束。
非线性动力学模型可以使用不同的数学方法来描述,包括微分方程、差分方程、离散映射等。
其中,常用的非线性动力学模型包括洛伦兹系统、范德波尔系统、Logistic映射等。
非线性动力学模型的应用非线性动力学模型在各个领域的应用非常广泛,包括自然科学、社会科学、工程技术等。
以下是一些应用实例:1. 生态系统建模非线性动力学模型在生态学中被广泛应用于描述物种的动态演化过程。
例如Lotka-Volterra模型描述了捕食者和猎物之间的相互作用,通过非线性关系描述了食物链中的生态系统行为。
2. 经济系统建模非线性动力学模型在经济学中被用于描述市场的波动和不确定性。
非线性动力学理论及其预测复杂系统演化规律非线性动力学理论是研究非线性系统行为的数学分析方法,并且被广泛应用于各个科学领域,如物理学、化学、生物学、经济学等。
非线性系统具有复杂的演化规律,这使得其行为难以预测和理解。
然而,非线性动力学理论在揭示这些复杂系统演化规律方面发挥着重要的作用。
非线性动力学理论的核心概念之一是混沌现象。
混沌是一种看似无序却具有确定性的动力学行为。
混沌系统对初始条件极其敏感,微小的变化可能导致系统的完全不同演化。
这使得预测混沌系统的未来行为非常困难,但非线性动力学理论可以帮助我们了解和描述混沌系统中的规律。
在非线性动力学中,还存在着很多其他重要的现象,如周期运动、分岔现象和吸引子等。
周期运动是指系统在某个轨道上周期性地运动,而分岔现象则指系统参数或初始条件微小改变时系统行为突变的现象。
而吸引子是指系统演化过程中的某些稳定状态,吸引子可以是点、线、面甚至是复杂的分形结构。
非线性动力学理论的研究方法包括数学模型的建立、动力学方程的求解和稳定性分析等。
数学模型是描述系统行为的基础,可以是连续的或离散的。
动力学方程是描述系统演化的数学表达式,可以是常微分方程、偏微分方程或差分方程等。
通过求解动力学方程,我们可以获得系统的轨道和稳定性信息,从而了解系统的演化规律。
利用非线性动力学理论,我们可以预测和理解复杂系统的演化行为。
复杂系统是由多个相互作用的元素组成,其整体行为难以通过简单的线性关系推导。
非线性动力学理论为我们提供了一种分析和预测复杂系统行为的框架。
例如,在生物学中,非线性动力学可以帮助我们理解生物系统中的自组织行为和生物钟的运作机制。
在物理学中,非线性动力学理论被应用于研究混沌系统和相变现象,如液滴的形成和磁性材料的相变。
在经济学中,非线性动力学理论可以帮助我们理解市场波动和金融危机的爆发机制。
尽管非线性动力学理论在理论和实践中发挥了重要作用,但仍存在一些挑战和限制。
首先,非线性动力学的数学工具较为复杂,需要较高的数学背景和计算能力。
《从非线性动力学到复杂系统》
段法兵
系统理论博士生课程
第一讲动态系统的发展
系统是一些相互关联的客体组成的集合,动态(动力dynamica)系统是系统状态变量,比如温度、位移、价格、信号幅值等,随着时间变化的。
它的描述可以用微分方程或者离散方程。
微分方程历史悠久,可追溯到牛顿、伽利略、欧拉、雅克比等人,用以描述行星的运动轨迹。
研究中发现即使满足牛顿引力定律的三体运动也非常复杂,其微分方程是非线性的,非线性是指不满足叠加定律的方程, 解无法利用已知函数进行描述,如果能够描述的我们称为显式解。
因此,庞加莱在1880年-1910年期间,试图利用解的拓扑几何性质来解释动态系统的运动规律,发现即使确定性系统,其运动规律也会出现随机性态,非常复杂(确定性系统是指其外力是确定的不随机,只要知道初始条件和演化方程,其运动是可预先确定的)。
非线性系统运动的复杂性:李雅普诺夫研究了系统平衡点?的稳定性?问题,随后本迪尔松等发现系统的解包含(1)平衡态(静止不动);(2)周期运动(比如行星)(3)拟周期,就是几个频率不可公约周期之和。
接着1975年Li和Yorke提出了混沌的概念,即系统的解是非周期的一种类似随机运动的现象,这其中就包含了洛伦兹提出的“蝴蝶效应”,根源在于这类
非线性动力系统对于初始条件的极其敏感性,初始条件的微小变化导致了系统状态的巨大改变,从此有关非线性科学的发展异常迅速,形成了现代动力学理论,其最重要的贡献是揭示了一个简单的模型可能蕴含了无比复杂的动力学性态。
例子:Van der Pol (范德波尔)方程
1920年Van der Pol利用电子震荡管研究心脏的跳动问题,比如人工心脏起
电子震荡管电路
如上图所示,利用电路分析,变量变换之后得到方程
这里■■就是电路的固有震荡频率。
Matlab code 这是内嵌的,无需自己定义。
fun cti on outl = vdp1(t,y)
outl = [y(2); (1-y(1)A2)*y(2) - y(1)];
解释:首先定义微分方程函数,因为方程(1-1)可以用变量代换
X = u(1 _ X )x -,X
设向量y = 丫1 = x ,
-1,那么
丿2 一X '
因此上面的函数就定义了向量微分方程的导数部分。
求解
[t,y] = ode45('vdp1',[0 20],[3; 0]);
plot(t,y(:,1),'-',t,y(:,2),'-.')
ode45表示了求解常微分方程的一种精度达到 4-5阶的算法,微分方程vdp1定义
的,[0 20]表示起止时间,向量[3;0]表示了微分方程的边界条件,即
搏器
P 厂、
2 2 x -u(1「X )x • X = 0 (1-1)
对—X ] -31
Van der Pol (范德波尔)方程解轨迹
可以看出x (实际是电压)呈现周期性震荡现象,如果将
7=;之间关系画出来,如下图
相图(x,dx/dt)
这种没有时间轴,只是几何表示系统解的几何图形就是相图1
f
1 -I J J 1 I I 1 1 ■' 9 1 1 [ r i * /
2 : > 1 \ 1. 1
f II 1
* 1 1 1 1 \ 10 20 30 40 50 y i
-1
-2
60
-30 x 与其导数,就是 二维图叫相图,
维叫相空间
从图中可以看到一个,系统的解渐渐地趋向一个封闭的环
-极限环。
上面是初 >> hold on >> [t,y] = ode45('vdp1',[0 60],[0.2; 0]);
>> plot(y(:,1),y(:,2),'-g')
相图(x,dx/dt)
见绿色部分,同样系统的解渐渐地趋向一个封闭的环 -极限环。
因此,不论给定 初始电压是小还是大,系统都慢慢地归结为一个震荡环,这个震动环的频率与固 有频率一致,称为自激震荡。
补充:对于
dt 初始条件■
我们希望在固定的步长tn 上求解
Un P ), r )
— 0,1,****
匸J 3,如果取
设
h =切+1 —切
严
y(t + h)= y(t) + J f(s,y(s))ds.
Euler认为方程左侧的导数近似为
心严=f (t, y)
h
那么自然有迭代形式
1/n+l = Un"H hf (t n t y n),
S+l = + h.
这种公式也称为矩形公式,即差分形式。
为了获得好的精度,必须采用很小的步
长,这种Euler方法没有误差估计方法,无法自动确定步长来达到期望的精度。
类比于积分中的中点公式和梯形公式,在中点h/2处估计函数的值,然后在计算
导数,
S f n ? Vn ),
.( h h \
s2 = f(如+ 厅,弘+ 7T51 I ,
y«+l = l/n十h g、
®+l = 十h*
两位德国数学家卡尔•龙格和马丁•威尔海姆•库塔于1900年左右发明提出了种方法Runge-Kutta法
Si = f ( tn l/n)?
V h h \ = J 1
1 2/n + ~S 1 J I £(
h h \ 53 = J I + ^2 * i/n + 齐$2 ) i
«4 = f^n + h, y n 十 ^53),
h
咖+1 = J/n + ~{s l + 2舅 2 + 2^3 + 百4人 b
b+1 =t n + h.
它每一步要用到4个函数估值。
Felhberg 又进行了改进,每个步长内进行了 6次 求值,形成了 5阶公式,因此精度又进行了提高,并且定义了一个误差向量来调 节步长的大小,平滑的阶段步长大些,急剧变化的导数处步长小些,这样形成了
变步长的Runge-Kutta-Felhberg 算法。
此即Matlab 中ode45算法,ode 就是常微分 方程,45就是比较4阶和5阶公式的误差进行变步长计算的意思。
刚性方程:所谓刚性方程,就是说存在两
(多)重尺度,一个尺度比另外一个尺度大很多。
所导致的麻烦就是在计算中很难兼顾两者。
例如下面的方程:
dx/dt=-100 x-100.1 y
dy/dt=100.1 x-100 y 两个特征值 lambda_1 = -200.1 ,lambda_2 = -0.1,所以解为 a1*exp(-200.1*x) + a2*exp(-0.1*x), 无论你用什么样的尺度(单一尺度)都不能很好刻画解的行为。
一个是快变行为,一个是慢变 行为。
所有这样的方程计算时候,稳定性条件比较苛刻。
实际情况要比这还复杂得多。
fun ction out1 = vdp1000(t,y) out1 = [y(2); 1000*(1-y(1)A 2)*y(2) - y(1)];
求解还是用ode45就比较慢
tic
[t,y] = ode45('vdp1000',[0 2000],[0.2; 0]); toe I
— * x_ - 〔1000*(1 —y :)y2—y
非常慢,如果采用ode15s专门求解刚性方程的算法,那么会非常快。
tie;[t,y] = ode15s('vdp1000',[0 2000],[0.2; 0]);toe
例子:Lorenz方程,起源于气象学,里面的参数都具有物理意义,和大气热对流系数,粘性系数,导温系数,Rayleigh数等
x = TOx 10y
y = 28x _ y _ xz
,z= -舟z+ xy
20 -40
xyz相空间
如何理解上述方程的行为,周期是固定的吗?是唯一的吗?这些都属于非线
性动力学的范畴,进入第二讲动力系统的形态分析。