这将是它: Y'(T)= Y(t)的/ TT^2/Y^2 * T Y(1)= 1用matlab解这个微分方程?
我曾尝试: 功能hazi3b()
[T,Y] = ode45(@bfugveny,[112],1); 情节(T,Y, 'GX')
端
和:
功能DY = bfugveny(T,Y)
DY = Y(T)/吨 - 吨^ 2 /(Y^2 * T);
末
这将是它: Y'(T)= Y(t)的/ TT^2/Y^2 * T Y(1)= 1用matlab解这个微分方程?
我曾尝试: 功能hazi3b()
[T,Y] = ode45(@bfugveny,[112],1); 情节(T,Y, 'GX')
端
和:
功能DY = bfugveny(T,Y)
DY = Y(T)/吨 - 吨^ 2 /(Y^2 * T);
末
您不必在公式中写出y(t)
。
传入你的oracle的y
已经是y
-evaluated-at-time-t
的猜测。
所以尽量
dy = y/t - t^2/(y^2*t);
谢谢,它工作:) –
所以接受答案....感恩是好的,但接受对你和帮助你的人是好的。也投了票。 – duffymo
出了什么问题? –
试图访问y(1.22);索引必须是正整数或逻辑。 bfugveny中的错误(第3行) dy = y(t)/ t-t^2 /(y^2 * t); (1),y + f * hB(:1),odeArgs {:});在ode45中的错误(行262) f(:,2)= feval(odeFcn,t + hA hazi3b(第2行)中的错误 [T,Y] = ode45(@bfugveny,[112],1); –