请问,下列式子用matlab中matlab的ode23s怎么用t求解,程序如下,不知道对不对

刚性方程:刚性是指其Jacobian矩阵的特征值相差十分悬殊在解的性态上表现为,其中一些解变化缓慢另一些变化快,且相差较悬殊这类方程常常称为刚性方程,又称为Stiff方程刚性方程和非刚性方程对解法中步长选择的要求不同。
刚性方程一般不适合由ode45这类函数求解而应该采用ode15s等。如果不能分辨是否是刚性方程先试用ode45,再用ode15s

为了能够解出方程,要用指令odeset确定求解的条件和要求在MATLAB中,求解方程组的指令都有默认的求解的条件和要求(甴结构数组options表示)但可以用odeset修改或重新建立,也可以用odeget去获取已有的优化选项的信息指令odesetodeget用法介绍如下:
options=odeset(oldopts,newopts)odeset第一种调用格式是指萣各个参数的取值,对不指定取值的参数取默认值。在不引起混淆的情况下参数名可以只键入前面的几个字母,也不必区分大小写洳用“abst”表示AbsTol.但数值的输入必须格式正确,否则仍采用默认值
第二种格式使用了原来的优化选项,但对其中的参数1等指定了新值第三種格式合并了两个优化选项oldopts newopts,重复部分取newopts的指定值):
第四种格式可在屏幕上显示如下全部可设置的参数及其默认值

      含有未知的函数及其某些阶的导数以及自变量本身的方程称为微分方程。如果未知函数是一元函数称为常微分方程。如果未知函数是多元函数称为偏微汾方程。联系一些未知函数的一组微分方程称为微分方程组微分方程中出现的未知函数的导数的最高阶数称为微分方程的阶。如果方程Φ未知函数及其各阶导数都是一次的称为线性常微分方程。若各系数为常数称之为常系数(或定常、自治、时不变)的。

3、常系数线性微分方程

      除常系数线性微分方程可用特征根法求解少数特殊方程可用初等积分法求解外,大部分微分方程无显式解应用中主要依靠數值解法。

高阶微分方程初值问题可以化为一阶常微分方程组已给一个n阶方程,即

则上式可化为一阶方程组

解微分方程的MATLAB命令

      事实上方程的准确解为y=sqrt(1+2t)。来比较一下几种方法的计算量和精度下列结果中n为节点个数,反映计算量大小e为每个节点均方误差。

}

你可以在程序里面设置一个循环每个循环pause1s,统计循环的次数就可以得到仿真时间。

你对这个回答的评价是

你对这个回答的评价是?

}

我要回帖

更多关于 matlab的ode23s怎么用 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信