matlab中的ode45怎样设置定步长解微分方程?
matlab中的ode45是变步长解法,定步长在simulink里有ode4和ode5,但是在matlab里面好像没有现成的函数,听说ode45可以实现定步长的输出,如何在matlab中(不是simulink)实现定步长解微分方程,
function [x,y]=runge_kutta1(ufunc,y0,h,a,b) %参数表顺序依次是微分方程组的函数名称,初始值向量,步长,时间起点,时间终点,
n=floor((b-a)/h); %求步数
x(1)=a; %时间起点
y(:,1)=y0; %赋初值,可以是向量,但是要注意维数
for ii=1:n
x(ii+1)=x(ii)+h;
k1=ufunc(x(ii),y(:,ii));
k2=ufunc(x(ii)+h/2,y(:,ii)+h*k1/2);
k3=ufunc(x(ii)+h/2,y(:,ii)+h*k2/2);
k4=ufunc(x(ii)+h,y(:,ii)+h*k3);
y(:,ii+1)=y(:,ii)+h*(k1+2*k2+2*k3+k4)/6; %按照龙格库塔方法进行数值求解
end
这是给你编的定步长龙哥库塔,即ode4
-
点赞数:0回答数:1
-
点赞数:0回答数:2
-
点赞数:0回答数:1
-
点赞数:0回答数:1
-
点赞数:0回答数:2
-
点赞数:0回答数:1
-
点赞数:0回答数:1
-
点赞数:0回答数:3
-
点赞数:0回答数:2
-
点赞数:0回答数:1

最新问答: 男生人数占女生人数的八分之七,男生人数占全年级人数的几分之几?女生人数占全年级人数的几分之几? 小林买了3本练习本,付了5元.每本练习本多少元 你是我的小生命用英语最短的怎么说 Geographical Location是什么意思 已知椭圆E的中心在原点,焦点在x轴上,椭圆的焦距为2,离心率e=1/2,直线l:y=k(x-1)(k≠0) 太阳是固体吗无 设F(X)是可导的偶函数,且f'(0#存在.证明f'#0#=0求大神帮助 证明德摩根法则RT 是数学集合中的概念 又叫摩根定律.....我问的是 怎么证明... 为什么飞翔鸟类胸肌细胞线粒体的数目比不上不飞翔鸟类的多? 某研究性学习小组在练习使用打点计时器的实验中,所打出的纸带如图所示,纸带上相邻两点对应的时间间隔为0.02s.在纸带上选 touch是不是感官动词哦 将6个长2厘米,宽1厘米的小长方开摆放三角形,拼成的图形的周长是多少厘米 已知正三角形ABC边长为a,当一点P在三角形ABC的外接圆上的劣弧AB(AB上面有一弧)上移动时,求S三角形PAC +S 一腰上的中线也是一腰上的高的等腰三角形是等边三角形.怎么证?急 果糖和葡萄糖性质上有哪些区别? Would you like a trip ( ) the zoo with me? There is a bus under the tree.复数形式 对于自变量在取值范围内一个确定的值,函数都有唯一确定的值与之对应,这个对应值叫做______ 什么动物经常滑到?iq题 某同学身高1.8m,体重为70㎏,在迎奥运田径运动会上他参加背越式跳高比赛,起跳后身体横着越过了1.8m高的杆,请问: