在matlab中用 ode45解方程 y'=x-y,y(0)=1,0

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/30 07:14:08

在matlab中用 ode45解方程 y'=x-y,y(0)=1,0
在matlab中用 ode45解方程 y'=x-y,y(0)=1,0

在matlab中用 ode45解方程 y'=x-y,y(0)=1,0
%%%%%%%%%%%%%%%%M-file 以下要单存一个文件在工作路径下%%%%%%%%%%%%%function dotx=myfun(t,y)
dotx=t-y;
%%%%%%%%%%%%%%%%M-file 以上要单存一个文件在工作路径下%%%%%%%%%%%%%在command 窗口
[t,y]=ode45('myfun',[0,3],1)
用plot(t,y) 绘图后,你要哪点值都可以,直接在图上看就可以了,即使用插值得到的也是近似值,因为你选择ode45已经是近似了,要想要精确值建议求解析解如下
y=dsolve('Dy=x-y','y(0)=1','x')
得到
y=x + 2/exp(x) - 1
把x=1,2,或3 带入即可

不知

在matlab中用 ode45解方程 y'=x-y,y(0)=1,0 matlab中用ode45求微分方程的数值解问题[T,Y]=ode45('function1',[0 3000],[2 0]);其中[0 3000]是起始时间和结束时间.那步长是在哪里设置的,怎么没有步长呢? matlab用ode45 解微分方程如何用ODE45解微分方程 其中初值给出方程如下,请帮忙写出方程在程序中如何表达,我用function dy=func(x,y,r)dy = zeros(size(y));dy(1) = y(1)dy(2) = sqrt((1+r^2)/r^2*log((1+r^2)/(1+y(1)^2))+y(1 matlab中ode45 及odeset的使用.在蝴蝶效应的解方程组中,方程是%x'=a*(X(2)-X(1))%y'=X(1)*(c-X(3))-X(2)%z'=X(1)*X(2)-b*X(3)使用ode45解方程.语句为lorenz=@(t,X) [a*(X(2)-X(1));X(1)*(c-X(3))-X(2);X(1)*X(2)-b*X(3)];[t X]=ode45(lorenz MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告:写在补充里面>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducin matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y) matlab中ODE45函数该如何使用?在ODE45中,总是有[t,y],t是说明ODE45只能用于时间相关求解,还是t只是表示任意参数的符号? 用matlab求y=f(x),y''=y^2的数值解,待解方程为:y=f(x),y''=y^2,即y的二阶导等于y的平方边界条件:y(0)=1,y'(100)=0用matlab函数ode45求解function dy=test(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=[y(1)]^1;end[X,Y]=ode45('test',[0 100 如何用ode45解这个微分方程 这个方程的matlab程序又该如何写?ay-b(y'^2)+c*siny=0 自变量为x a,b,c均为正实数 matlab中的ode45怎样设置定步长解微分方程?matlab中的ode45是变步长解法,定步长在simulink里有ode4和ode5,但是在matlab里面好像没有现成的函数,听说ode45可以实现定步长的输出,如何在matlab中(不是simu matlab中,ode45能解决含有(dy/dt)^2和d2y/d2t的方程么 matlab ode45求解微分方程 D2y-0.01*Dy.^2+2*y=sin(t),初始值0 matlab ode45求解微分方程后 怎样输出dy/dt,即如何plot(t,y') matlab解一阶微分方程方程 dy/dt+i=240*sin(2*pi*50*t)其中i=(0.0006455*exp(0.4182*(10*y/0.001+21)) + (-3.971e+004)*exp(-0.4854*(10*y/0.001+21)))*0.1/100这个i的方程是我用CFtoolbox得到的.我想用ODE45来得到这个微分方程的 matlab ode45的问题%定义函数function dy = For_y( t,y )dy=zeros(2,1);Dy(1)=y(2);Dy(2)=2;end %求解微分方程[Tget,Yget]=ode45(@For_y,[0 10],[0,0]);plot(Tget,Yget,'-*');%得出的图%同样的方程用dsolve函数就可以求y=dsolve( 方程组如下:X1(t) '=- x1(t)x22(t)+0.999+0.42cos(1.75t)X2(t) '=x1(t)x22(t)-x2(t)Y(t)=sin[x1(t)+x2(t)]用matlab中的ode45函数解方程 t∈[0,20],x1(0)=1.0,x2(0)=1.0它的程序是怎么实现的? matlab ode45 求解常微分方程组dx/dt = c-b*x-w*x*y;dy/dt = w*x*y-(b+m)*y;其中c,b,w,m为已知参数,x,y为未知的,分别对t求导,请问这个常微分方程组要怎么用ode45解呢? 应用Matlab的ode45解含有约束条件的多参数微分方程