1个回答
展开全部
>> f=inline('1/6-y/30','t','y');
>> [t,y]=ode45(f,[0,5],[0]);
>> plot(t,y)
另外，由高数的知识，可以得到解析解为：y=5-5*e^(-t/30)
把上述的解析值也放到上图中。
>> hold on
>> plot(t,5-5*exp(-t/30),'r*')
可以看到二者相当吻合。
>> [t,y]=ode45(f,[0,5],[0]);
>> plot(t,y)
另外，由高数的知识，可以得到解析解为：y=5-5*e^(-t/30)
把上述的解析值也放到上图中。
>> hold on
>> plot(t,5-5*exp(-t/30),'r*')
可以看到二者相当吻合。
追问
话说，没看懂...可以用求微分方程的方法来说不啦/(ㄒoㄒ)/~~谢谢
参考资料： 百度一下
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
