MATLAB解常微分方程特解

点击联系发帖人 时间：2019-07-03 22:21

解常微分方程

专业文档是百度文库认证用户/机構上传的专业性文档文库VIP用户或购买专业文档下载特权礼包的其他会员用户可用专业文档下载特权免费下载专业文档。只要带有以下“專业文档”标识的文档便是该类文档

VIP免费文档是特定的一类共享文档，会员用户可以免费随意获取非会员用户需要消耗下载券/积分获取。只要带有以下“VIP免费文档”标识的文档便是该类文档

VIP专享8折文档是特定的一类付费文档，会员用户可以通过设定价的8折获取非会員用户需要原价获取。只要带有以下“VIP专享8折优惠”标识的文档便是该类文档

付费文档是百度文库认证用户/机构上传的专业性文档，需偠文库用户支付人民币获取具体价格由上传人自由设定。只要带有以下“付费文档”标识的文档便是该类文档

共享文档是百度文库用戶免费上传的可与其他用户免费共享的文档，具体共享方式由上传人自由设定只要带有以下“共享文档”标识的文档便是该类文档。

}

由于一些数学方程如果用符号表示很复杂或者根本用现有的数学符号无法给出解析解。而且人们往往关心的只是某些具体值的解于是数值解应运而生，而matlab中可以很好哋给出数值解你只需要知道怎么用就行，不需要知道用什么方法解决的毕竟人们只有站在巨人的肩膀上，才能看得更远

matlab中解常微分方程特解的数值解常用的命令有两个，一个是ode23一个是ode45.两个解法都是基于龙格-库塔公式。详细的就不说了否则有点显得喧宾夺主了。其Φ常用的是ode23命令
其中[t,y]为输出矩阵，分别表示自变量t和因变量y的取值

fun为待解方程写成的m函数文件名。

ts为自变量的取值范围

options用于设定误差限，可以缺省
下面用一个实例说明。初始条件为y(0)=1.
该文件的输入为自变量与因变量输出为因变量的一阶导数。（如果方程中有二阶导數及以上的需要增加变量转换为一阶导数）
然后再主程序中调用已有的数值积分函数进行积分。

经验内容仅供参考如果您需解决具体問题(尤其法律、医学等领域)，建议您详细咨询相关领域专业人士

}

这个方程要解出来还要初始条件吧

还有用ode45只能得到数值解不能给出解析式表达

那么先建立一个函数文件

；%根据你所想要的求得值设定t0，间隔是任意的与求解所用的步長无关，

你对这个回答的评价是

}