24小时热门版块排行榜    

查看: 1349  |  回复: 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

金虫 (小有名气)

引用回帖:
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的回帖
查看全部 10 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
CODE:
clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
[t y]=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, [0:100],0.148   );
[t y]
figure(1)
plot(t,y,'b-o')

用matlab对微分方程作图,请高手帮忙!
附图1.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-11-17 20:39:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
zhuwx1999: 金币+50, ★★★很有帮助 2014-11-18 08:18:52
dbb627: 金币+2, 谢谢指导 2014-11-18 09:07:00
随着时间的推移,即增大计算时间,y会稳定某一数值:
CODE:
clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
[t y]=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, [0:1000],0.148   );
[t y]
figure(1)
plot(t,y,'b-o')

用matlab对微分方程作图,请高手帮忙!-1
附图1.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2014-11-17 20:40:45
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见