matlab动力学解析程序详解
- 格式:doc
- 大小:1.30 MB
- 文档页数:11
刚体动力学方程matlab simscape multibody 在MATLAB中,可以使用Simscape Multibody(前身为SimMechanics)工具箱来建模和仿真刚体动力学系统。
Simscape Multibody提供了一种基于物理原理的方法,可以对多体系统进行建模,并通过求解动力学方程来模拟其运动。
以下是使用MATLAB Simscape Multibody进行刚体动力学建模和仿真的一般步骤:
1. 定义刚体:使用Simscape Multibody中的刚体块创建刚体对象。
可以指定刚体的质量、惯性矩阵、形状和连接关系等属性。
2. 建立连接:使用连接块将刚体连接起来,形成多体系统。
连接可以是旋转关节、平移关节、万向节等。
3. 添加约束:根据需要,可以添加额外的约束条件,例如限制关节的运动范围或强制两个刚体之间的相对位移。
4. 定义输入:根据问题需求,可以添加外部力或扭矩作用于刚体上。
5. 设置初始条件:指定系统的初始位置、速度和其他状态。
6. 定义仿真参数:设置仿真时间、积分算法和其他仿真选项。
7. 运行仿真:运行仿真以计算系统的运动轨迹和其他输出。
8. 分析结果:使用MATLAB中的绘图和分析工具,对仿真结果进行可视化和后处理。
Simscape Multibody提供了丰富的示例模型和文档,可以帮助用户快速上手并理解刚体动力学建模的基本原理。
可以通过MATLAB的官方网站或在线文档获取更多关于Simscape Multibody的信息和资源。
1。
matlab拟合动力学方程
MATLAB可以用于拟合动力学方程。
在MATLAB中,我们可以使用curve fitting工具箱来实现这个目标。
首先,我们需要收集我们的数据,并确定我们要拟合的动力学方程的类型。
例如,我们可以选择一阶动力学方程:dy/dt = -k*y,其中y是我们的输出变量,t是时间,k是动力学常数。
然后,我们可以使用MATLAB的curve fitting工具箱来拟合这个方程。
以下是一些步骤:
1. 导入数据:将我们收集的数据导入MATLAB环境。
确保数据已经存储为一个列向量,例如y和t。
2. 建立起始参数:根据我们的动力学方程,我们需要为k提供一个初始猜测值。
这个值可以根据我们的应用和经验来确定。
3. 建立模型:使用fittype函数创建一个模型对象,该对象表示我们要拟合的动力学方程。
4. 进行拟合:使用fit函数拟合我们的数据。
该函数将数据和模型作为参数,并返回包含拟合结果的对象。
5. 分析结果:我们可以通过访问拟合对象的属性来分析拟合结果,例如拟合参数的值和置信区间。
6. 可视化结果:使用plot函数绘制原始数据和拟合结果的图像,以便我们可以直观地评估拟合的质量。
通过这些步骤,我们可以使用MATLAB拟合动力学方程,并从拟合结果中获得我们感兴趣的参数值。
sindy的matlab程序-概述说明以及解释1.引言1.1 概述Sindy是一种基于数据驱动的系统辨识方法,通过对系统的动态行为进行分析和建模,可以帮助我们更好地理解系统的运行机制和规律。
Matlab作为一种强大的科学计算工具,能够提供丰富的功能和工具,帮助我们进行数据处理、模型建立和结果分析。
本文将详细介绍Sindy在Matlab环境下的应用,探讨其在不同领域中的作用和价值。
通过对Sindy程序的优势和局限性进行分析,可以更全面地了解其在系统辨识方面的特点和适用范围。
最后,我们将总结Sindy 的Matlab程序的重要性,展望其未来在系统辨识领域的发展,并希望能为相关研究提供一定的参考和启发。
1.2 文章结构本篇文章主要分为三个部分:引言、正文和结论。
在引言部分,将会对文章的主题进行一定的概述,介绍Sindy的Matlab程序的背景和意义,以及对文章的结构进行简要的介绍。
在正文部分,将详细介绍Sindy的Matlab程序的相关内容,包括程序的介绍、应用领域、优势和局限性等方面。
最后,在结论部分,将总结Sindy的Matlab程序的重要性,展望其在未来的发展,并给出一些结束语,为全文画上一个完美的句号。
1.3 目的本文的目的是介绍Sindy的Matlab程序,探讨其在科学研究和工程领域中的应用情况。
通过对Sindy程序的介绍和分析,读者可以更深入地了解Sindy程序的原理和特点,以及其在系统辨识、动力系统建模等方面的重要性和价值。
同时,本文也将讨论Sindy程序的优势和局限性,对于读者在选择合适的程序工具时提供参考。
通过本文的阐述,旨在激发读者对于Sindy程序的兴趣,促进该程序在未来的发展和应用。
2.正文2.1 Sindy的Matlab程序介绍Sindy是一个用于系统辨识和模型推断的Matlab程序。
该程序的全称为Sparse Identification of Nonlinear Dynamics,意为稀疏非线性动力学识别。
齿轮故障动力学仿真matlab-概述说明以及解释1.引言1.1 概述齿轮是机械传动中常用的零部件,其在各种机械设备中起着至关重要的作用。
然而,由于工作环境的恶劣以及长期使用的磨损,齿轮可能出现故障,导致机械设备的性能下降甚至损坏。
为了更好地理解齿轮故障的动力学特性,可以通过仿真技术来模拟和分析齿轮系统的运行状态,并及时发现潜在的故障点。
本文将介绍齿轮故障动力学仿真在MATLAB中的应用,通过分析齿轮系统的动态特性,探讨不同故障模式对系统性能的影响,从而为齿轮故障诊断和预防提供有益的参考。
通过本文的研究,我们希望能够加深对齿轮故障动力学的理解,提高齿轮系统的可靠性和安全性。
1.2 文章结构文章结构部分的内容如下:文章结构包括以下几个部分:1. 引言:介绍文章的背景和研究意义,引出文章的主题和研究内容。
2. 正文:分为两个部分,分别是齿轮故障动力学简介和MATLAB在齿轮故障动力学仿真中的应用。
在齿轮故障动力学简介部分,将介绍齿轮故障动力学的基本概念和原理,为读者提供必要的背景知识。
在MATLAB 在齿轮故障动力学仿真中的应用部分,将详细介绍MATLAB在该领域的具体应用及其优势。
3. 结论:总结文章的主要内容和研究成果,对研究进行评价和展望未来的研究方向。
通过以上部分的内容安排,读者可以清晰地了解整篇文章的主要结构和内容安排,帮助他们更好地理解和阅读文章。
1.3 目的本文的主要目的在于探讨利用MATLAB进行齿轮故障动力学仿真的方法和技术。
通过对齿轮系统中可能出现的不同故障情况进行建模和仿真,我们可以更好地理解齿轮系统的运行机理,并且能够快速有效地诊断和解决齿轮故障问题。
同时,本文也旨在为工程师和研究人员提供一个基于MATLAB的齿轮故障动力学仿真平台,帮助他们更好地分析和优化齿轮系统的性能,推动齿轮传动技术的发展和应用。
通过本文的研究,我们希望能够为齿轮系统的设计、运行和维护提供更加有效的工程解决方案,提高齿轮系统的可靠性和稳定性。
MATLAB 动力学建模1. 引言动力学建模是一种描述物体运动和行为的数学建模方法。
在工程学和物理学中,动力学建模被广泛应用于设计、控制和优化系统。
MATLAB是一个强大的数值计算软件,可以用于动力学建模和仿真。
本文将介绍MATLAB在动力学建模中的应用。
2. 动力学建模基础动力学建模的基础是牛顿第二定律,即力等于质量乘以加速度。
根据这个定律,可以建立物体的运动方程。
在MATLAB中,可以使用符号计算工具箱来求解运动方程。
例如,考虑一个简单的弹簧振子系统,其中一个质量m通过一个弹簧与墙壁相连。
弹簧的劲度系数为k,质量m的加速度为a,弹簧的位移为x,墙壁的位置为0。
可以建立如下运动方程:m * a = -k * x在MATLAB中,可以使用符号计算工具箱来求解这个方程,并得到系统的运动方程。
3. 动力学建模方法在动力学建模中,有几种常用的方法可以用于建立系统的数学模型。
以下是一些常见的方法:3.1. 基于物理原理的建模基于物理原理的建模是一种常见的动力学建模方法。
这种方法基于系统的物理特性和力学原理,建立系统的数学模型。
例如,对于一个机械系统,可以根据质量、惯性、摩擦等物理特性,建立系统的动力学方程。
3.2. 系统辨识建模系统辨识建模是一种通过实验数据来建立系统模型的方法。
通过对系统进行实验观测,收集系统的输入和输出数据,然后使用系统辨识算法来估计系统的动力学模型。
MATLAB提供了多种系统辨识工具箱,可以用于建立系统的数学模型。
3.3. 仿真建模仿真建模是一种通过数值仿真来建立系统模型的方法。
通过使用数值计算方法和数学模型,可以模拟系统的运动和行为。
MATLAB提供了强大的仿真工具箱,可以用于建立系统的数学模型,并进行仿真研究。
4. MATLAB 动力学建模工具MATLAB提供了多种工具和函数,用于动力学建模和仿真。
以下是一些常用的工具和函数:4.1. 符号计算工具箱符号计算工具箱可以用于求解符号方程和符号运算。
matlab使用拉格朗日法计算动力学方程参数标题:用MATLAB使用拉格朗日法计算动力学方程参数导言:在工程和科学领域中,我们经常需要对系统进行动力学建模和分析。
动力学方程是描述系统运动的数学表达式,它们可以帮助我们理解和预测系统的行为。
而计算动力学方程的参数对于系统的设计、优化和控制具有重要意义。
本文将介绍如何使用MATLAB编程语言和拉格朗日法来计算动力学方程参数,通过数值求解和优化方法,为工程师和科学家提供一个有力的工具。
1. 动力学方程与参数计算在动力学中,系统可以通过一组微分方程来描述。
这些方程表示系统的行为和演变,通常包括质量、速度、加速度和其他相关因素。
对于复杂系统,计算动力学参数可能是一项繁琐且复杂的任务。
幸运的是,拉格朗日法可以简化这一过程,通过定义能量或运动方程来推导系统动力学方程。
2. 拉格朗日法介绍拉格朗日法是一种基于能量和运动方程的变分法。
通过将系统的动能和势能组合成拉格朗日函数,并使用欧拉-拉格朗日方程,可以得到系统的运动方程。
拉格朗日法不仅适用于具有广义坐标的连续系统,还适用于离散系统和有束缚条件的系统。
3. MATLAB编程实现在MATLAB中,我们可以使用符号计算工具包来处理复杂的数学表达式和符号计算。
我们需要定义系统的广义坐标、广义速度和拉格朗日函数。
通过欧拉-拉格朗日方程生成系统的运动方程。
我们可以使用数值求解和优化方法来计算动力学参数。
4. 示例:双摆系统为了更好地理解如何使用MATLAB实现拉格朗日法,我们将以一个双摆系统为例。
双摆系统由两个连杆组成,每个连杆上挂有一个质点。
我们需要计算系统的动力学参数,包括质点的位置、速度和加速度。
通过拉格朗日法,我们可以推导出系统的运动方程,并使用MATLAB 进行求解和优化。
5. 讨论与总结本文介绍了使用MATLAB编程语言和拉格朗日法计算动力学方程参数的方法。
通过拉格朗日法,我们可以简化动力学参数的计算,并为系统的设计和优化提供可行性的解决方案。
matlab动力学分析程序详解··1.微分方程的定义对于duffing 方程032=++x x xω ,先将方程写作--==3112221x x x x x ω function dy=duffing(t,x) omega=1;%定义参数f1=x(2);f2=-omega^2*x(1)-x(1)^3; dy=[f1;f2];2.微分方程的求解function solve (tstop) tstop=500;%定义时间长度 y0=[0.01;0];%定义初始条件[t,y]=ode45('duffing',tstop,y0,[]);function solve (tstop) step=0.01;%定义步长y0=rand(1,2);%随机初始条件tspan=[0:step:500];%定义时间范围[t,y]=ode45('duffing',tspan,y0);3.时间历程的绘制时间历程横轴为t ,纵轴为y ,绘制时只取稳态部分。
plot(t,y(:,1));%绘制y 的时间历程 xlabel('t')%横轴为t ylabel('y')%纵轴为y grid;%显示网格线axis([460 500 -Inf Inf])%图形显示范围设置4.相图的绘制相图的横轴为y ,纵轴为dy/dt ,绘制时也只取稳态部分。
红色部··分表示只取最后1000个点。
plot(y(end-1000:end,1),y(end-1000:end,2));%绘制y 的时间历程xlabel('y')%横轴为yylabel('dy/dt')%纵轴为dy/dt grid;%显示网格线5.Poincare 映射的绘制对于不同的系统,Poincare 截面的选取方法也不同对于自治系统一般每过其对应线性系统的固有周期,截取一次对于非自治系统,一般每过其激励的周期,截取一次例程:duffing 方程032=++x x x ω的poincare 映射function poincare(tstop) global omega; omega=1;T=2*pi/omega;%线性系统的周期或激励的周期step=T/100;%定义步长为T/100 y0=[0.01;0];%初始条件tspan=[0:step:100*T];%定义时间范围[t,y]=ode45('duffing',tspan,y0);for i=5000:100:10000%稳态过程每个周期取一个点plot(y(i,1),y(i,2),'b.'); hold on;% 保留上一次的图形 endxlabel('y');ylabel('dy/dt');Poincare 映射也可以通过取极值点得到 function poincare(tstop) y0=[0.01;0];tspan=[0:0.01:500];[t,y]=ode45('duffing',tspan,y0); count=find(t>100);%截取稳态过程 y=y(count,:);n=length(y(:,1));%计算点的总数··for i=2:n-1if y(i-1,1)+epsy(i+1,1)+eps % 简单的取出局部最大值plot(y(i,1),y(i,2),'.'); hold on end endxlabel('y');ylabel('dy/dt');6.频谱yy=fft(y(end-1000:end,1)); N=length(yy); power=abs(yy);freq=(1:N-1)*1/step/N;plot(freq(1:N/2),power(1:N/2)); xlabel('f(y)') ylabel('y')7.算例duffing 方程03=++x x x的时间历程,相图,频谱和poincare 映射。
matlab求解飞行器动力学方程
要在MATLAB中求解飞行器的动力学方程,你需要首先确定飞行器的动力学模型。
飞行器的动力学方程通常涉及到飞行器的质量、惯性、空气动力学力和控制力等因素。
一般来说,飞行器的动力学方程可以通过牛顿运动定律和欧拉运动方程来建立。
在MATLAB中,你可以使用数值求解方法来求解微分方程,常见的方法包括欧拉法、四阶Runge-Kutta法等。
你需要将飞行器的动力学方程转化为状态空间方程,然后使用MATLAB中的ode45函数或其他适合的数值求解函数来求解微分方程。
另外,你可能需要使用MATLAB中的符号计算工具箱来进行符号运算,以便得到飞行器的动力学方程的数值解。
你还可以使用MATLAB中的仿真工具箱来进行飞行器动力学仿真,以验证你得到的动力学方程的准确性。
总之,在MATLAB中求解飞行器的动力学方程需要对飞行器的动力学模型有深入的理解,以及熟练掌握MATLAB的数值计算、符号计算和仿真工具箱的使用方法。
希望这些信息能够帮助你开始在MATLAB中求解飞行器的动力学方程。
基于Newmark-β法的Matlab简单程序编写李坤明;黄菲【摘要】我们知道,对于结构动力学问题的直接求解往往是比较困难的,而Newmark-β法的出现,很巧妙地将这种困难化解.Newmark-β法是一种时程分析的方法,它将动力学微分方程的求解问题转化为若干个代数方程逐步求解,是一种隐式积分法,可以得到足够精确的解.本文基于Newmark-β法的理论知识,用Matlab 对一个简单动力学问题的求解进行编程运算.【期刊名称】《四川建材》【年(卷),期】2018(044)009【总页数】3页(P65-66,90)【关键词】Newmark-β法;编程;结构力学【作者】李坤明;黄菲【作者单位】西华大学土木建筑与环境学院,四川成都 610039;西华大学土木建筑与环境学院,四川成都 610039【正文语种】中文【中图分类】TP3131 逐步积分法在结构动力学问题的分析中,对于承受任意动荷载的线性结构,我们可以采用Duhamel、频域分析等方法,这些方法都可以很方便地求解出所需要的结果。
但是上面的两种方法都用到了叠加原理,所以它们只适用于求解线性结构体系,同时也必须要求这种体系的特性在反应过程中不能发生变化[1]。
然而,在另一方面,我们的实际生活中有很多重要的结构动力学问题,整个反应体系并不能单纯地认为是线性变化的,比如,在足以引起强烈破坏的地震作用下,很多建筑物的反应,我们都必须考虑非线性反应的影响。
所以,我们还必须要发展对于非线性结构的动力响应的分析方法。
对于非线性体系的动力响应问题,发展出了最有效的分析方法—逐步积分法[1]。
在这种方法中,我们采用一系列的短时间增量Δt来计算反应,而通常为了方便起见,我们将Δt取为等间距的时间步长。
在每个时间间隔的起点和终点建立动力平衡方程,并以一个假设的反应机理为根据,近似地计算在时间间隔内体系的运动(通常忽略在时间间隔内产生的不平衡),体系的非线性特性可以用每个时间间隔的起点所求得的当前变形状态的新特性来说明。
MATLAB 轴承动力学引言轴承是机械设备中广泛使用的部件,主要用于支撑和减少摩擦。
轴承动力学研究轴承在运转过程中的力学特性,包括轴承的负载、振动、磨损等。
MATLAB作为一种强大的工具,可以用于轴承动力学的建模、分析和仿真。
本文将介绍如何使用MATLAB进行轴承动力学的研究。
轴承建模输入参数轴承的建模需要输入一些基本参数,包括轴承的几何尺寸、材料性质、运转条件等。
这些参数将影响轴承的载荷分布、接触应力分布等。
轴承几何模型根据轴承的类型和尺寸,可以建立不同的几何模型。
常见的轴承类型包括球轴承、滚子轴承等。
根据轴承的尺寸和几何形状,可以确定轴承的内外圈半径、接触角度等参数。
轴承材料模型轴承的材料选择也对其力学性能有重要影响。
常见的轴承材料包括钢、陶瓷等。
根据轴承材料的力学性质,可以建立材料模型,包括弹性模量、泊松比等参数。
轴承载荷分析轴承在运转过程中承受着不同的载荷,包括径向载荷、轴向载荷、扭矩等。
这些载荷将导致轴承内部产生应力和变形。
轴承载荷计算根据轴承的运转条件和工作负荷,可以计算轴承承受的载荷大小。
对于复杂的工况,可以使用MATLAB进行多工况载荷计算,并对结果进行统计分析。
轴承应力分析根据轴承的载荷和材料模型,可以计算轴承内部的应力分布。
这对轴承的寿命和可靠性分析具有重要意义。
轴承振动分析轴承在运转过程中会产生振动,这对轴承的工作性能和寿命有一定影响。
轴承振动模态分析通过建立轴承的振动模型,可以计算轴承的固有频率和振动模态。
这对于减少轴承的共振现象和提高运转稳定性很重要。
轴承振动特性分析根据轴承的振动信号,可以对轴承的振动特性进行分析。
这包括频谱分析、时域分析等方法。
轴承磨损分析轴承在长时间运转过程中会发生磨损,这会影响轴承的工作性能和寿命。
轴承磨损模型根据轴承的工作条件和材料特性,可以建立轴承的磨损模型。
这包括磨损速率、磨损形式等参数。
轴承磨损预测通过分析轴承的磨损模型和运转条件,可以对轴承的磨损进行预测。
··
1.微分方程的定义
对于duffing 方程03
2=++x x x ω ,先将方程写作⎩⎨⎧--==31122
21x x x x x ω function dy=duffing(t,x)
omega=1;%定义参数
f1=x(2);
f2=-omega^2*x(1)-x(1)^3;
dy=[f1;f2];
2.微分方程的求解
function solve (tstop)
tstop=500;%定义时间长度
y0=[0.01;0];%定义初始条件
[t,y]=ode45('duffing',tstop,y0,[]);
function solve (tstop)
step=0.01;%定义步长
y0=rand(1,2);%随机初始条件
tspan=[0:step:500];%定义时间范围
[t,y]=ode45('duffing',tspan,y0);
3.时间历程的绘制
时间历程横轴为t ,纵轴为y ,绘制时只取稳态部分。
plot(t,y(:,1));%绘制y 的时间历程
xlabel('t')%横轴为t
ylabel('y')%纵轴为y
·· grid;%显示网格线
axis([460 500 -Inf Inf])%图形显示范围设置
4.相图的绘制
相图的横轴为y ,纵轴为dy/dt ,绘制时也只取稳态部分。
红色部分表示只取最后1000个点。
plot(y(end-1000:end ,1),y(end-1000:end ,2));%绘制y 的时间历程
xlabel('y')%横轴为y
ylabel('dy/dt')%纵轴为dy/dt
grid;%显示网格线
5.Poincare 映射的绘制
对于不同的系统,Poincare 截面的选取方法也不同
对于自治系统一般每过其对应线性系统的固有周期,截取一次 对于非自治系统,一般每过其激励的周期,截取一次
例程:duffing 方程03
2=++x x x ω
的poincare 映射 function poincare(tstop)
global omega;
omega=1;
T=2*pi/omega;%线性系统的周期或激励的周期
step=T/100;%定义步长为T/100
y0=[0.01;0];%初始条件
tspan=[0:step:100*T];%定义时间范围
[t,y]=ode45('duffing',tspan,y0);
for i=5000:100:10000%稳态过程每个周期取一个点
plot(y(i,1),y(i,2),'b.');
hold on;% 保留上一次的图形
end
·· xlabel('y');ylabel('dy/dt');
Poincare 映射也可以通过取极值点得到
function poincare(tstop)
y0=[0.01;0];
tspan=[0:0.01:500];
[t,y]=ode45('duffing',tspan,y0);
count=find(t>100);%截取稳态过程
y=y(count,:);
n=length(y(:,1));%计算点的总数
for i=2:n-1
if y(i-1,1)+eps<y(i,1) && y(i,1)>y(i+1,1)+eps % 简单的取出局部最大值
plot(y(i,1),y(i,2),'.');
hold on
end
end
xlabel('y');ylabel('dy/dt');
6.频谱
yy=fft(y(end-1000:end,1));
N=length(yy);
power=abs(yy);
freq=(1:N-1)*1/step/N;
plot(freq(1:N/2),power(1:N/2));
xlabel('f(y)')
ylabel('y')
7.算例
duffing 方程03
=++x x x
的时间历程,相图,频谱和poincare
映射。
function dy=duffing(t,x)
omega=1;%定义参数
f1=x(2);
f2=-omega^2*x(1)-x(1)^3;
dy=[f1;f2]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function duffsim(tstop)
step=0.01
y0=[0.1;0];
tspan=[0:step:500];
[t,y]=ode45('duffing',tspan,y0); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(2,2,1)
plot(t,y(:,1));%绘制y的时间历程
xlabel('t')%横轴为t
ylabel('y')%纵轴为y
grid;%显示网格线
axis([460 500 -Inf Inf])%显示范围设置%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(2,2,2)
plot(y(end-1000:end,1),y(end-1000:end,2));%绘制y的时间历程
xlabel('y')%横轴为y
ylabel('dy/dt')%纵轴为dy/dt
grid;%显示网格线%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(2,2,3)
··
yy=fft(y(end-1000:end,1));
N=length(yy);
power=abs(yy);
freq=(1:N-1)*1/step/N;
plot(freq(1:N/2),power(1:N/2));
xlabel('f(y)')
ylabel('y') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(2,2,4)
count=find(t>100);%截取稳态过程
y=y(count,:);
n=length(y(:,1));%计算点的总数
for i=2:n-1
if y(i-1,1)+eps<y(i,1) && y(i,1)>y(i+1,1)+eps % 简单的取出局部最大值
plot(y(i,1),y(i,2),'.');hold on;
end
end
xlabel('y');ylabel('dy/dt');
··
8.分岔图的绘制
随F变化的分岔图。
-
+
F
x
+
3.03=
x
x
x2.1
t
cos
function dy=duffing(t,x)
global c;
omega=1;%定义参数
f1=x(2);
f2=omega^2*x(1)-x(1)^3-0.3*x(2)+c*cos(1.2*t); dy=[f1;f2]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
global c; %定义全局变量
range=[0.1:0.002:0.9];%定义参数变化范围
··
k=0;
YY=[];%定义空数组
for c=range
y0=[0.1;0];%初始条件
k=k+1;
tspan=[0:0.01:400];
[t,Y]=ode45('duffing',tspan,y0);
count=find(t>200);
Y=Y(count,:);
j=1;
n=length(Y(:,1));
for i=2:n-1
if Y(i-1,1)+eps<Y(i,1) && Y(i,1)>Y(i+1,1)+eps % 简单的取出局部最大值。
YY(k,j)=Y(i,1);
j=j+1;
end
end
if j>1
plot(c,YY(k,[1:j-1]),'k.','markersize',3);
end
hold on;
index(k)=j-1;
end
xlabel('c');
ylabel('y');
··
随F变化的分岔图
F=0.20
··
··
F=0.27
F=0.275
F=0.2875
··
F=0.32
F=0.36
F=0.4
·
··
F=0.652
F=0.8。