MATLAB解微分方程

用matlab时间也不短了,可是一直没有接触过微分方程。这次看看书,学习学习,记点儿笔记。

1.可以解析求解的微分方程。

dsolve()

调用格式为:

将yc的内容copy到latex中编译,得到结果。

关于Matlab的微分方程,直到今天才更新第2篇,实在是很惭愧的事——因为原因都在于太懒惰,而不是其他的什么。

在上一篇中,我们使用dsolve可以解决一部分能够解析求解的微分方程、微分方程组,但是对于大多数微分方程(组)而言不能得到解析解,这时数值求解也就是没有办法的办法了,好在数值解也有很多的用处。

数值分析方法中讲解了一些Eular法、Runge-Kutta 法等一些方法,在matlab中内置的ode 求解器可以实现不同求解方法的相同格式的调用,而不必太关心matlab究竟是用什么算法完成的。

这一回我们来说明ode45求解器的使用方法。

1.ode45求解的上手例子:

求解方程组

Dx=y+x(1-x^2-y^2);

Dy=-x+y*(1-x^2-y^2)

初值x=0.1;y=0.2;

先说明一下最常用的ode45调用方式,和相应的函数文件定义格式。

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

其中,Fun就是导函数,tspan为求解的时间区间(或时间序列,如果采用时间序列,则必须单调),x0为初值。

这时,函数文件可以采用如下方式定义

function dx=odefun(t,x)

对于上面的小例子,可以用如下的程序求解。

matlab循环求微分方程,MATLAB解微分方程相关推荐

  1. matlab求曲线极值程序,matlab函数求极值matlab函数求极值.ppt

    matlab函数求极值matlab函数求极值 * * 函数的极值 1.一元函数的极值 函数命令:fminbnd 调用格式:[x,feval,exitflag,output]=fminbnd(fun,x ...

  2. python求解析解,Python解微分方程

    Python解微分方程 微分方程回顾 微分方程:python 解析解(SymPy) 微分方程:python数值解(SciPY) 微分方程组:python数值解 微分方程回顾 微分方程是用来描述某一类函 ...

  3. matlab中求立方根,MATLAB基础入门

    MATLAB有许多使用方法,但最基本,也是入门时首先要掌握的是MATLAB命令窗口(Command Window)的使用方法. MATLAB命令窗口是用于输入数据,运行MATLAB函数和脚本,并显示结 ...

  4. matlab trapz求二重积分,matlab求积分(超详细,含int integral integral2/3 quad trapz

    matlab求积分 matlab求积分函数工具: int 用法1: 格式: int(fun,x,a,b) 功能: 计算定积分 用法2: 格式: int(f,x) 功能: 计算不定积分 注: 使用int ...

  5. 外推法程序matlab,急求用MATLAB用龙格库塔和外推法解一阶微分方程

    共回答了21个问题采纳率:81% f=inline('-y+x+1','x','y'); %微分方程的右边项 dx=0.05; %x方向步长 xleft=0; %区域的左边界 xright=10; % ...

  6. matlab分支与循环求函数,Matlab变量、分支语句和循环语句

    一.Matlab的变量 1.变量的命名 Matlab的变量名由数字.字母和下划线组成,但是只能由字母开头.大小写敏感,最大长度为63个字符,不能使用Matlab的关键字作为变量名,应当避免使用函数名作 ...

  7. matlab循环求 近似for,Matlab for 多个变量循环能不能这样啊 ,求教高手!!!!...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 for a=0.003:0.0005:1; b=0.002:0.0005:0.9; c=0.001:0.0005:0.8;d=0.0005:0.0005: ...

  8. matlab近似计算求积分,matlab 实验二 定积分的近似计算

    实验二定积分的近似计算 利用牛顿-莱布尼兹公式虽然可以精确地计算定积分的值,但它仅适用于被积函数的原函数能用初等函数表达出来的情形.如果这点办不到或者不容易办到,这就有必要考虑近似计算的方法.在定积分 ...

  9. matlab如何求极点,matlab求极点和零点

    让每个人平等地提升自我一.实验目的 1.利用 MATLAB 的 fdatool 观察传递函数 H(z)的零极点分布.幅度响应.脉冲响 应及阶跃响应,并观察零极点分布与系统稳定性...... 假设为一个 ...

  10. 在matlab中求协方差,matlab里面的求协方差函数

    一:matlab里面关于方差和协方差的几点需要注意的 求均值的函数:mean 求方差的函数:var 标准差: 方差: 求协方差的函数:cov 1 如果a,b是向量,则: cov(a,b)相当于是[D( ...

最新文章

  1. 物联网技术正颠覆传统医疗行业
  2. python中的set类型
  3. 小于n与n互质的所有数的和(欧拉函数+快速幂)
  4. SAP CRM和SAP Hybris的订单修改记录
  5. 不用ajax实现动态分页,详解ajax +jtemplate实现动态分页
  6. [react] React中在哪捕获错误?
  7. LNMP与CA认证的童话故事
  8. 使用.NET程序集作为Business Data Connectivity数据源(二)
  9. LVM+Xen虚拟化应用
  10. js中的逻辑与()和逻辑或(||)(转载)
  11. Eclipse 一直提示 loading descriptor for 的解决方法
  12. exchange无法收发邮件_MAPI session超过Exchange最大限制,邮件发送失败的终极解决方法...
  13. 现代信号处理——时频分析与时频分布(短时Fourier变换)
  14. java 接口 实验报告_java-接口练习实验报告
  15. 逆水寒捏脸服务器维护,《逆水寒》优化热门服务器排队体验 捏脸颜值评分机制优化...
  16. MTK_android11_以太网(内网)和WIFI(外网)共存
  17. 复利java_复利计算1.0,2.0,3.0(java)
  18. 项目实战:简易俄罗斯方块(附源码)
  19. 快速使用Vitamio框架播放网络视频
  20. 网上最全的系统服务想PF降到50以下(转)

热门文章

  1. excel+if函数+android,Excel中if函数多重条件的使用
  2. Windows系统如何mount盘符到目录
  3. 应届生计算机工程师岗位初试
  4. 已知两数的最大公因数和最小公倍数
  5. 【IoT】产品设计之市场概念:市场定位、产品定位、市场需求、产品需求
  6. WIN10删除微软拼音输入法,设置默认输入法为英文
  7. 阿里云虚拟主机内存占用大问题解决
  8. 如何向公众号添加的H5页面的链接(已解决)
  9. python爬取淘宝数据魔方_看我是怎么使用淘宝数据魔方--巧用数据魔方把握市场信息...
  10. python对seo有什么用_python对seo的帮助 – python对seo的帮助是什么?- 企业服务