dy(t)/dt = k y(t) [1 - y(t)/w] y(0) = y0 用MATLAB能编程能求出y(t)吗?
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/25 19:47:17
dy(t)/dt = k y(t) [1 - y(t)/w] y(0) = y0 用MATLAB能编程能求出y(t)吗?
dsolve('Dy = k*y*(1 - y/w)','y(0) = y0')
...
ans =
-w/(exp(w*(log(-(w - y0)/y0)/w - (k*t)/w)) - 1)
再问: 我试过了,07版本的可以得出答案,但是10版本的就出来其他不正确的答案!
再答: 我用两个版本计算了一下,结果的表达式不同,但应是等价的。 syms x w k y0 t; y1 = -w/(exp(w*(log(-(w - y0)/y0)/w - (k*t)/w)) - 1) y2 = w/(1+exp(-k*t)*(w-y0)/y0) z1 = y1-y2 z1 = simple(z1) 结果: z1 = 0 ... 验算: syms x w k y0 t; y = -w/(exp(w*(log(-(w - y0)/y0)/w - (k*t)/w)) - 1) z = diff(y,t) - k*y*(1 - y/w) simple(z) ... 验算: syms x w k y0 t; y = w/(1+exp(-k*t)*(w-y0)/y0) z = diff(y,t) - k*y*(1 - y/w) simple(z)
...
ans =
-w/(exp(w*(log(-(w - y0)/y0)/w - (k*t)/w)) - 1)
再问: 我试过了,07版本的可以得出答案,但是10版本的就出来其他不正确的答案!
再答: 我用两个版本计算了一下,结果的表达式不同,但应是等价的。 syms x w k y0 t; y1 = -w/(exp(w*(log(-(w - y0)/y0)/w - (k*t)/w)) - 1) y2 = w/(1+exp(-k*t)*(w-y0)/y0) z1 = y1-y2 z1 = simple(z1) 结果: z1 = 0 ... 验算: syms x w k y0 t; y = -w/(exp(w*(log(-(w - y0)/y0)/w - (k*t)/w)) - 1) z = diff(y,t) - k*y*(1 - y/w) simple(z) ... 验算: syms x w k y0 t; y = w/(1+exp(-k*t)*(w-y0)/y0) z = diff(y,t) - k*y*(1 - y/w) simple(z)
请帮忙用Matlab求解微分方程dy/dt=[k^(t-1)*y-d]y
用matlab ode45求微分方程组 dx/dt+x+y=0 dy/dt+x-y=0 x(0)=0 y(0)=1 t=
求解dx/(x+t)=dy/(-y+t)=dt
用Matlab编程求解 二阶微分方程:4*d^2y(t)/dt^2+y(t)=dx(t)/d(t)-0.5x(t)
y= ∫[0,x](t-1)^3(t-2)dt,dy/dx(x=0)
x=f(t) y=g(t) 为什么dy/dx=(dy/dt)*(dt/dx)
dx/dt=x+t,dy/dt=-y+t,求x,y(t为常数).
用matlab解:y=dsolve('D2y+exp(-t)*Dy-sin(sqrt(t))*y=exp(-t)*sin
如果DY/ DT=-2Y并且若y=1当t =0时,t是什么当Y= 1/2?
微分方程d^2y/dt^2-dy/dt+y=1,y0=0,y1=0.怎么用matlab求解啊
MATLAB中,t=0:.1:2*pi; y=sint(t); plot(t,y)
d(t(dy/dt))/dx为什么等于t² d²y/dt²+t dy/dt