fkjj.net
当前位置:首页 >> mAtlAB程序oDE45 >>

mAtlAB程序oDE45

ode45表示采用四阶-五阶Runge-Kutta算法,它用4阶方法提供候选解,5阶方法控制误差,是一种自适应步长(变步长)的常微分方程数值解法,其整体截断误差为(Δx)^5。解决的是Nonstiff(非刚性)常微分方程。 ode45语法: [T,Y] = ode45(odefun,tspan,...

参考代码: % 参数定义m1 = 1; m2 = 2; m3 = 3;k1 = 4; k2 = 5; k3 = 6; k4 = 7;P0 = 8; w = 9; % 微分方程A = [k1+k2 -k1 0; -k2 k2+k3 -k3; 0 -k3 k3+k4];dx = @(t,x)[x(4:6); (-A*x(1:3)+[P0*sin(w*t);0;0])./[m1;m2;m3]]; % 初始条件(x4~x6...

这张图来自于MATLAB 技术论坛,这是一个非常有价值的网站,里面高手云集,希望你可以进去看看,里面有非常详细的关于ODE的介绍 http://www.matlabsky.com/thread-528-1-1.html

1、你贴出来的报错信息和代码对不上号:前面显示错误的那行代码和你贴出来的完全不一样;而后面的错误(Input argument 'u1' is undefined)也不可能是目前的代码所导致的——的确是有错,但错误应该是iL未定义才对。 2、使用ode*系列函数解常微分...

相关常数都没给,这里随便设置一些数值: m=1;c=2;k=3;y0=4;yp0=5;tf=10;dy=@(t,y)[y(2);-(c*y(2)+k*y(1))/m];[t,y]=ode45(dy,[0 tf],[y0 yp0]);plot(t,y(:,1))

使用ode45解微分方程组,可以用下列格式来求解。 [t,x]=ode45(odefun,tspan,x0) 式中:odefun——给定微分方程组的自定义函数;tspan——自变量t的范围,即【t0,tf】;x0——x的初值向量,即x1(0),x1'(0),x2(0),x2'(0),x3(0),x3'(0) ...

MATLAB提供了7个常微分方程求解器(solver),分别是ode45,ode23,ode113,ode15s,ode23s,ode23t,ode23tb,其中前3个适用于求解非刚性(Nonstiff)问题,后4个适用于刚性问题。所谓刚性问题,简单点说,就是系统包含多个相互作用但变化速度相差十...

function dx = rigid(t,x)dx = ones(6,1); dx(1) = 0.0515*x(2)-0.388*x(1)-0.136*x(1);dx(2) = 0.388*x(1)-0.0515*x(2);dx(3) = 0.136*x(1)+0.0254*x(6)-0.067*x(3)-0.0788*x(3);dx(4) = 0.0788*x(3)-0.0068*x(4)+0.001*x(5);dx(5) = 0.0068*x(...

您好,以下这张图来自于MATLAB 技术论坛,这是一个非常有价值的网站,里面高手云集,希望你可以进去看看,里面有非常详细的关于ODE的介绍: http://www.matlabsky.com/thread-528-1-1.html

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...

网站首页 | 网站地图
All rights reserved Powered by www.fkjj.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com