24小时热门版块排行榜    

查看: 1353  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhuwx1999

金虫 (小有名气)

[求助] 用matlab对微分方程作图,请高手帮忙!已有1人参与

matlab零基础,请高手帮忙写个程序,十分感谢!!!

方程为:dy/dt=a*y^0.266-b*y^3.532+c*y^2.266,其中a=0.01162, b=0.00135, c=0.0000051246.
             t=0, y=0.148.

想作出y-t曲线图(最好同时能导出y-t数据),请大家帮帮忙,再次感谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuwx1999

金虫 (小有名气)

引用回帖:
3楼: Originally posted by 月只蓝 at 2014-11-17 20:40:45
随着时间的推移,即增大计算时间,y会稳定某一数值:
clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, ,0.148   );

figure(1)
plot(t,y,'b-o')


...

谢谢你!
我还有两个小问题:
1、[t y]=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, [0:1000],0.148   )公式中的0.148是不是代表dy/dt的初值?
2、另外,想把公式中的0.266改为-0.266,用y^(-0.266)还是直接y^-0.266?
4楼2014-11-18 08:15:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuwx1999

金虫 (小有名气)

引用回帖:
3楼: Originally posted by 月只蓝 at 2014-11-17 20:40:45
随着时间的推移,即增大计算时间,y会稳定某一数值:
clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, ,0.148   );

figure(1)
plot(t,y,'b-o')


...

不好意思,再问一个问题:
根据实验,y-t曲线应该出现一个峰值(dy/dt小于等于0),不知是不是计算的时候有什么假设条件?
5楼2014-11-18 08:35:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuwx1999

金虫 (小有名气)

引用回帖:
6楼: Originally posted by 月只蓝 at 2014-11-18 08:52:10
问题1
答:是的,正是初值。
问题2
答:y^(-0.266)比较保险;此外注意要在英文状态下输入。...

可能我没有表达清楚。
问题1:公式中0.148表示的是y的初值还是dy/dt的初值?
8楼2014-11-18 13:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuwx1999

金虫 (小有名气)

引用回帖:
9楼: Originally posted by 月只蓝 at 2014-11-18 14:03:05
对于一阶常微分方程,初值专指t=0时刻,y的数值。...

谢谢你!

[ 发自小木虫客户端 ]
10楼2014-11-18 18:43:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuwx1999 的主题更新
信息提示
请填处理意见