www.rtmj.net > MAtlAB使用oDE45解微分方程组的基本问题,下面有一...

MAtlAB使用oDE45解微分方程组的基本问题,下面有一...

建立myfun.m文件,把以下代码复制进去 function dx=myfun(t,x)dx=zeros(2,1);dx(1)=4*x(1)+x(2).^2;dx(2)=8*x(1)+2*x(2).^2+5*x(2);运行以下代码 t=[0 0.1]; x0=[1 1]; [tt,xx]=ode45('myfun',t,x0); plot(tt,xx(:,1),'-',tt,xx(:,2),'-.')%画出...

在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dxdt=zhidao_rk4_5(t,x) global a b c d r x1=x(1); x2=x(2); x3=x(3); x4=x(4); dxdt=[ a*(x2-x1)+x4; d*x1-x...

这个方程要解出来还要初始条件吧 不妨设(x(0),y(0))=(A,B); 还有用ode45只能得到数值解,不能给出解析式表达 那么先建立一个函数文件 test_fun.m function dx=tsst_fun(x,t) dx(1)= c-b*x(1)-w*x(1)*x(2); dx(2)=w*x(1)*x(2)-(b+m)*x(2); 那么调...

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

首先,不要求你多么的精通matlab,但能够理解你的代码中的那几条注释语句就可以了,否则就不要继续往下看了。。。 测试结果:代码没问题,没有语法错误;算法没问题,可以用45阶龙格库塔求解;精度没问题,该方程可以使用ode45默认精度;运行没...

ode45函数 求一阶常微分方程的初值问题 [t,y] = ode45(@(t,y)y-2*t/y,[0,4],1); plot(t,y); 求解 y’ – y + 2*t / y且初值y(0) = 1的常微分方程初值问题,返回自变量和函数的若干个值. 若不写返回值,则会自动作出函数随自变量的变化图像. ode45(...

使用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) ...

相关常数都没给,这里随便设置一些数值: 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))

分段函数的一个简单方法就是分了几段就写几个函数: 定义一个主调函数、两个子函数 最后一个主调程序调用 function dx=fun(t,x) dx=zeros(2,1); f=; dc(t)=diff('f','t');%f就是你的c(t) dx(1)=x(2)*dc(t); dx(2)=x(1)+6*t^2;%0

function myode2(wx,wy,wz) clc; x0=[0,1,1];%x0表示三个初始值 tspan=[0:0.1:12]; [t,x]=ode45(@(t,x,w) weifen(t,x),tspan,x0,[]); disp(' 时间 x y z ') disp([t,x]) function dx=weifen(t,x) dx=zeros(3,1); dx(1)=x(2)*x(3); dx(2)=-x(1)*x...

网站地图

All rights reserved Powered by www.rtmj.net

copyright ©right 2010-2021。
www.rtmj.net内容来自网络,如有侵犯请联系客服。zhit325@qq.com