图片中是转子碰磨故障的动力学模型,应用龙格库塔法求解,在画分岔图和poincare截面图,画出来不对,不知道为什么,下面是我自己编的程序,请大侠们帮忙看看我错在哪里了?

global omega;              %以转速为变化画分岔图

range=[1.1:0.001:1.35];

k=0;

YY1=[];

step=2*pi/100;  %步长。

for omega=range

period=2*pi/omega;

x0=[0;0;0.001;0];

omega

k=k+1;

% discard the first 60 periodic data;

%除去前面60个周期的数据,并将最后的结果作为下一次积分的初值

tspan=[0:step:60*period];

[t,x]=ode45(@fangcheng,tspan,x0);

y0=x(end,:);

j=1;

for i=60:200

tspan=[i*period:step:(i+1)*period];

[t,x]=ode45(@fangcheng,tspan,x0);

YY1(k,j)=x(end,1);   % get the omega data from every period end

j=j+1;               %取出每一个周期内的第一个解的最后一个值。

y0=x(end,:);

end

<

matlab画虫口模型分岔图,Matlab画图片碰磨转子系统中的分岔图,poincare映射图???...相关推荐

  1. matlab画虫口模型分岔图,[转载]Matlab——混沌分岔图

    经过近期的研究发现,目前对于系统单参数分岔图的计算共有以下的几种方法: 1)最大值法 即对系统微分方程(组)进行求解,对求解的结果用getmax函数进行取点,并绘图. 2)Poincare截面法 对系 ...

  2. 【数学建模】基于matlab三维海浪模型仿真【含Matlab源码 1159期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab三维海浪模型仿真[含Matlab源码 1159期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付 ...

  3. 现在有一个未分库分表的系统,未来要分库分表,如何设计才可以让系统从未分库分表动态切换到分库分表上?

    问题 现在有一个未分库分表的系统,未来要分库分表,如何设计才可以让系统从未分库分表动态切换到分库分表上? 分析 你看看,你现在已经明白为啥要分库分表了,你也知道常用的分库分表中间件了,你也设计好你们如 ...

  4. 异构系统中的UUV、异构系统模型统一化及单应矩阵在异构系统中的应用

    异构系统中的UUV.异构系统模型统一化及单应矩阵在异构系统中的应用 UUV发展趋势及关键技术 异构编队系统编队 异构系统模型统一化 单应矩阵在异构系统中的应用 UUV发展趋势及关键技术 上一篇多智能体 ...

  5. matlab画平面心型线,如何用matlab画出心形线

    心形线,是一个圆上的固定一点在它绕着与其相切且半径相同的另外一个圆周滚动时所形成的轨迹线.下面就简单讲解一下如何用matlab画出心形线. 1.心形线的数学定义 2.编制的绘制心形线的matlab程序 ...

  6. matlab画y x 2,怎么用MATLAB画x^2+y^2=4*x图像

    第一问:不能用通常的三维命令来绘制椭圆球体,应用专用的绘制椭圆球体命令ellipsoid().更改后 由热心网友提供的答案1: 给你图片吧,纯代码常常会被百度认为不符合要求,你自己打这个方程应该是单叶 ...

  7. 基于AR模型的功率谱matlab,基于AR模型法功率谱估计的Matlab实现

    2005iF 12月 笫 21卷第 6期 武警工程学 院学报 jOURNAL OF ENGG COLLEGE OF ARMED POLICE FORCE DeC.2005 Vo1.21 No.6 [计 ...

  8. 用matlab画大数据曲线_基于MATLAB的大数据分析

    王媚 摘要:传统计算机模式与MATLAB软件技术相比较,传统软件运行起来较为复杂.以此基于MATLAB软件下的网络数据技术,它以高速化.关联化的优势成为人们眼中的焦点.本文针对传统网络软件模式中出现的 ...

  9. matlab 画狗头彩蛋,我用matlab画制动力分配曲线,求程序代码

    matlab画函数曲线 c = 3.00e8;h = 6.63e-34;k = 1.38e-23;lambda = (1:10:5000).*1e-9; figure(1)subplot(2,1,1) ...

最新文章

  1. Control Channel Element (CCE)
  2. 告别低分辨率网络,微软提出高分辨率深度神经网络HRNet | CVPR 2019
  3. JS开发利器-IxEdit傻瓜式JavaScript开发工具(附下载、汉化版、视频教程)
  4. 射频MOS管和三极管优缺点对比
  5. [Oracle]UNIX与Windows 2000上Oracle的差异(I)
  6. 修改el-popover和el-select样式
  7. centos下nginx安装和配置
  8. 成都阿里技术一面后感,给广大码友的一些小建议
  9. 疯狂代码 写给WEB2.0的站长
  10. WIN10电脑自动关机命令
  11. android imageview 获取bitmap缩放大小,android – Imageview缩放方法“centercrop”作为代码...
  12. HTTP 接口测试的流程
  13. PCB的制作工艺流程
  14. dx11学习笔记-3.三维空间变换(自己整理,基础详尽)
  15. HTML如何把输入框变成必填值,required输入框为必填项
  16. 在软件开发的早期阶段为什么要进行可行性研究?应该从哪些方面研究目标系统的可行性?...
  17. android 拨打带有分机号的电话
  18. 如何用计算机蓝牙发送文件,怎么用蓝牙传文件
  19. java软件测试经典案例,java语言编程案例 - Mrsjjl的个人空间 - 51Testing软件测试网 51Testing软件测试网-软件测试人的精神家园...
  20. 计算机毕业设计,vue+springboot的农产品溯源系统,内附源码

热门文章

  1. f12弹出debug_调试Javascript代码(浏览器F12及VS中debugger关键字)
  2. JS基础面试题一(原型)
  3. Mac - 查看/关闭 端口占用的方法
  4. office无法显示联机服务器,Microsoft Office Outlook出现无法打开服务器的解决方法
  5. imgui绘图框架使用
  6. 图像平滑(边缘检测)
  7. Java Lambda 连接字符串
  8. Python3网络爬虫:使用Beautiful Soup爬取小说
  9. linux服务器基线配置
  10. java与plc进行通迅