fun.m  文件里面定义微分方程组

function dx = fun(t, x)

dx = zeros(2, 1);

dx(1) = -x(1).^3 - x(2);

dx(2) = x(1) - x(2).^3;

end

运行命令:

[t,x] = ode45(@fun, [0,30], [1;0.5]);

[0, 30] 是 t 的区间,[1; 0.5] 是初值。

得到:一列 t, 两列 x。

精度参数:

options = odeset('RelTol', 1e-4, 'AbsTol', [1e-4, 1e-5]);

x(1), x(2) 的总的相对误差(relative error)是 1e-4,绝对误差(absolute error)分别是1e-4, 1e-5。

[t,x] = ode45(@fun, [0,30], [1;0.5], options);

matlab ode45相关推荐

  1. matlab ode45例子,matlab的ode45

    matlab ode45 解微分方程在用 odesolver(ode45, od... 在控制系统仿真中,常用的求微分方程数值解的函 数是ode23和ode45. 2 1. ode23 在MATLAB ...

  2. matlab生成有向网络,matlab ode45和矩阵生成有向网络图

    Matlab中解常微分方程的ode45 ode是专门用于解微分方程的功能函数,他有ode23,ode45,ode23s等等,采用的是Runge-Kutta算法.ode45表示采用四阶,五阶runge- ...

  3. matlab ode45 二阶微分,matlab关于ode45解二阶微分方程的困惑

    matlab关于ode45解二阶微分方程的困惑 matlab关于ode45解二阶微分方程的困惑 一个二阶微分方程: y''+y'+y=sin(t) 初始条件为y(0)=5,y'(0)=6. 过程: 先 ...

  4. matlab ode45求解齿轮动力学,ode45求解多自由度动力学方程实例.doc

    ode45求解多自由度动力学方程实例 Ode45函数调用形式如下:[T,Y]=ode45(odefun,tspan,y0)相关参数介绍如下:参数名称参数说明odefun用于存放待求解的方程的m文件名, ...

  5. matlab ode45使用,ODE45函数的使用——翻译

    <ODE45函数的使用--翻译>由会员分享,可在线阅读,更多相关<ODE45函数的使用--翻译(7页珍藏版)>请在人人文库网上搜索. 1.在Matlab中使用ode45简介Ma ...

  6. matlab ode45设置步长,MATLAB中用ode45求解微分方程,如何设置最大步长?

    如果你用过simulink里的ode45配置,我觉得你就会发现高赞就是在扯淡,你给的时间序列只是采样点,根本不是设置步长用的.我最近偷懒不想用simulink就研究了下ode45的函数配置项,在mat ...

  7. matlab ode45求解齿轮动力学,Matlab拟合动力学参数遇到问题(ode45)

    function [rho2,F] = rho2_F(k,r,ss,Ne,Np) sy = sum(r.^2); rho2 = 1 - ss./sy F = (sy - ss).*(Ne-Np)./( ...

  8. matlab ode45 初值,为什么设定的初值在ode45里无效?

    function dx= exercisesub(t,x,u,d) J=2.0; d=sin(t); x1dot=x(2); x2dot=(1/J)*u+(1/J)*d; dx=[x1dot; x2d ...

  9. matlab定步长ode,[转载]matlab ode45 函数传自定义参数用法及定步长ode

    要用的时候总是忘记,这回给把它写在这里! %%程序1 arg1 = 2; arg2 = 1; [T,Y] = ode45('vdp1000',[0 10],[2 0], [], arg1, arg2) ...

  10. matlab ode45的使用

    一元微分方程 在使用ode45函数时,要先把方程变形一下. 比如方程x¨(t)+x˙(t)+x(t)=0\ddot{x} \left(t\right)+\dot{x} \left(t\right)+x ...

最新文章

  1. 爬虫之requests模块中cookies参数的使用
  2. es 调整gc_实际中进行GC调整
  3. android 上下偏差怎么写_详解 Android 热更新升级如何突破底层结构差异?
  4. Luogu P3731 [HAOI2017]新型城市化
  5. linux mount挂载大小,Linux中mount挂载问题小结
  6. sql2005通用分页存储过程
  7. java equals和==以及hashmap原理和区别
  8. Flink技术到底是什么?Flink原理及深度解析
  9. eclipse maven 打war包的几种方式
  10. vue的自定义指令的坑
  11. yuicompressor java_YUI Compressor使用配置方法 JS/CSS压缩工具
  12. 2012年3月份第2周51Aspx源码发布详情
  13. Python安装pyqt4
  14. OMNeT++学习---TicToc(1)
  15. 全自动反冲洗叠片过滤器
  16. 震撼人心的战争类背景音乐
  17. Python画樱花树代码和小猪佩奇
  18. Linux基线检查( 一)
  19. 不用的MacBook这样处理比较好
  20. BIUT Blockchain Explorer功能介绍

热门文章

  1. 用代码拆掉你思维的墙
  2. EGE基础:高级绘图篇
  3. python飞机大战项目描述_飞机大战项目梳理(一)
  4. 我的十年编程路 2013年篇
  5. php+js+背景特效,基于canvas+html5炫酷星空背景动画特效
  6. 树莓派学习笔记003-超声波距离传感器
  7. Win2D 入门教程 VB 中文版
  8. stm32低功耗模式下mcu的实际功耗(ST芯片手册)
  9. 黑作坊灯光大师完美破解版
  10. 光圈叶片数量与衍射导致的品质和成像问题