一、模型的常微分方程及初值

二、ode45求解函数

function r=hudie3(t,x)

global a;

global b;

global c;

a=8/3;b=10;c=28;

r=zeros(3,1);

r(1)=-a*x(1)+x(2)*x(3);

r(2)=-b*x(2)+b*x(3);

r(3)=-x(1)*x(2)+c*x(2)-x(3);

end

三、绘制混沌系统三维曲线图

a=8/3;b=10;c=28;

t0=[0,100];f0=[0,0,1e-10];

[t,x]=ode45('hudie3',t0,f0);

plot3(x(:,1),x(:,2),x(:,3));

title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');

grid on

四、绘制混沌系统三个平面上的投影图

a=8/3;b=10;c=28;

t0=[0,100];

f0=[0,0,1e-10];

[t,x]=ode45('hudie3',t0,f0);

subplot(2,2,1)

plot3(x(:,1),x(:,2),x(:,3));

title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');

grid on

subplot(2,2,2)

plot(x(:,1),x(:,2));

title('Lorenz模型X-Y平面图');xlabel('x');ylabel('y');

grid on

subplot(2,2,3)

plot(x(:,1),x(:,3));

title('Lorenz模型X-Z平面图');xlabel('x');ylabel('z');

grid on

subplot(2,2,4)

plot(x(:,2),x(:,3));

title('Lorenz模型Y-Z平面图');xlabel('y');ylabel('z');

grid on

五、修改初值观察混沌系统图像的变化

以下是给出的几组任意的初值:

其图像如下:

六、探究各初值对图像的影响的参考代码

a=8/3;b=10;c=28;

for k=0:50:400

t0=[0,100];f0=[k,0,1e-10];

[t,x]=ode45('hudie3',t0,f0);

subplot(3,3,(k+50)/50)

plot3(x(:,1),x(:,2),x(:,3));

title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');

grid on

end

上面程序的图像如下:

探究第一个初值的变化对图像的影响的时候将变量K放在f0向量的第一个参数上,如上程序所示,探究第二个初值的变化对图像的影响的时候将变量K放在f0向量的第二个参数上,以此类推。

参数的取值范围通过修改 for k=0:50:400语句实现,但是注意在修改该语句时对应的 subplot(3,3,(k+50)/50)语句中(k+50)/50也要修改,比如将for k=0:50:400改为for k=0:500:4000,则subplot(3,3,(k+50)/50)应该相应的改为subplot(3,3,(k+500)/500)

java绘制蝴蝶效应曲线代码_利用MATLAB的常微分方程数值解绘制蝴蝶效应(混沌系统)三维曲线图、平面上的投影图...相关推荐

  1. python绘制繁花曲线代码_使用python和pygame绘制繁花曲线的方法

    前段时间看了一期<最强大脑>,里面各种繁花曲线组合成了非常美丽的图形,一时心血来潮,想尝试自己用代码绘制繁花曲线,想怎么组合就怎么组合. 真实的繁花曲线使用一种称为繁花曲线规的小玩意绘制, ...

  2. python绘制繁花曲线代码_使用python和pygame绘制繁花曲线

    前段时间看了一期<最强大脑>,里面各种繁花曲线组合成了非常美丽的图形,一时心血来潮,想尝试自己用代码绘制繁花曲线,想怎么组合就怎么组合. 真实的繁花曲线使用一种称为繁花曲线规的小玩意绘制, ...

  3. matlab低通滤波器库函数代码_利用Matlab filterDesigner 工具生成FIR滤波器函数,并调用实现低通滤波...

    本文使用的开发环境为:Win10 Matlab2018a 版本. 在matlab命令窗口输入:filterDesigner命令,即可打开filterDesigner设计工具. 按照下图调整FIR低通滤 ...

  4. python画pr曲线代码_利用Python中的numpy包实现PR曲线和ROC曲线的计算

    闲来无事,边理解PR曲线和ROC曲线,边写了一下计算两个指标的代码.在 python 环境下,sklearn里有现成的函数计算ROC曲线坐标点,这里为了深入理解这两个指标,写代码的时候只用到numpy ...

  5. python绘制繁花曲线代码_做个绘图仪,画出美丽的繁花曲线 | 上周实验回顾

    原标题:做个绘图仪,画出美丽的繁花曲线 | 上周实验回顾 在刚刚过去的周末,北京.上海.广州.深圳.杭州.南京的孩子们又在科学队长实验室玩嗨啦!一起来回顾一下吧~ 开眼时间 认识繁花曲线 纸币上的背景 ...

  6. python绘制繁花曲线代码_繁花曲线

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 #include #include #include // 引用图形库头文件 #include #include #include #define pi ...

  7. 【通信原理课程设计】利用MATLAB实现PCM编码的语音基带传输系统

    目录 一.摘要 二.设计任务 项目作业分析: 三.设计内容 1.整体设计方案 2.设计内容 四.结果与分析 这个项目在我的B站上有专门的视频演示:[通信原理课程设计]利用MATLAB实现PCM编码的语 ...

  8. 在matlab中实现累乘,如何利用matlab设计一个线性相位FIR带通滤波器,并在FPGA上实现...

    设计要求 利用matlab设计一个线性相位FIR带通滤波器,并在FPGA上实现. 1.滤波器指标:过渡带带宽分别为100~300HZ,500~700HZ,阻带允许误差为0.02,通带允许误差为0.01 ...

  9. matlab相关性分析频谱_利用Matlab绘制正弦信号的频谱图并做相关分析[共6页]

    <利用Matlab绘制正弦信号的频谱图并做相关分析[共6页]>由会员分享,可在线阅读,更多相关<利用Matlab绘制正弦信号的频谱图并做相关分析[共6页](6页珍藏版)>请在人 ...

最新文章

  1. SubVersion和Subclipse的简单使用方法
  2. VTK:Qt之BorderWidgetQt
  3. ​冲刺最后一公里——音视频场景下的边缘计算实践
  4. 如何控制product search attribute支持的操作类型
  5. 關於@Override
  6. NodeJS解决跨域问题:Access-Control-Allow-Origin
  7. 曝光原理_泰国精戈咖啡效果反馈 作用原理曝光
  8. iOS下数据存储的方式
  9. CIRD主站与DOPAS构建笔记#1
  10. [开源]quakeIII(雷神之锤3)源码
  11. 流媒体应用程序Mobdro或存在安全隐患
  12. StringBuilder常用函数
  13. ANSYS SPEOS VRXPERIENCE-基于物理特性的智能驾驶传感器高精度仿真
  14. python多行写入文件_写入文件
  15. it,GitHub与GitLab的区别
  16. [转]尺度不变特征变换(SIFT算法)Matlab程序代码测试例子的说明(Lowe的代码)
  17. PHP自学---黑马程序员笔记【持续更新】
  18. 为什么手机信号满格,但网速却还是那么慢?这4点原因是关键
  19. 流程图中参与者_探索流程参与者的移动流程门户
  20. 用多备份将网站数据备份到百度云,七牛云存储,阿里云OSS,亚马逊S3,金山云等云存储服务上

热门文章

  1. 升降摄像头电机静音驱动器TMC2300
  2. Object的wait、notify和notifyAll方法
  3. 【真实点评】迪士尼幼儿科学大揭秘全套12册 儿童成长故事百科绘本故事书 3-6岁
  4. 石鼓篆书第一人:吴昌硕数字藏品限量发行!
  5. Android11取消强制加密,悲催!Android 6.0设备强制要求开启全盘加密
  6. 华为eNSP模拟器中路由器设置密码+空闲时间
  7. 1688接口 -电商代购
  8. Hive UDF 开发手册
  9. 快速掌握EDIUS的几个常用技巧
  10. 大学计算机基础实践教程百度云,易学宝Office 2010微视频系列教程(300集 )