[MATLAB数学相关] 帮忙看看一个常微分方程的求解程序 - MATLAB爱好者论坛-LabFans.com

lxcwym · 2009-08-08, 13:16

麻烦大家帮忙看看问题出在哪儿了？我怎么试都不行,thanks
Function dy=myfun(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=2*y(2)-2*y(1)+5*exp(2*x)*sin(x);
[x, y]=ode23(myfun,[0,1],[-2,-3])；

在水一方257 · 2009-08-09, 23:17

最后一句话后面的分号改为英文状态下的分号

lxcwym · 2009-08-10, 12:38

还是不行啊？麻烦各位帮我再看看吧，比较急！

haozhou · 2009-08-12, 20:50

function dy=myfun(x,y)

dy=[y(2);2*y(2)-2*y(1)+5*exp(2*x)*sin(x);]
tspan=[0,1];
y0=[-2,-3];
[x, y]=ode23(@myfun,tspan,y0);
y1=y(:,1);
y2=y(:,2);
plot(x,y1,x,y2);就可以了

lxcwym · 2009-08-13, 12:41

非常感谢！可以了，呵呵

2009-08-08, 13:16	#1
lxcwym 初级会员注册日期: 2009-08-08 年龄: 45 帖子: 8 声望力: 0	帮忙看看一个常微分方程的求解程序麻烦大家帮忙看看问题出在哪儿了？我怎么试都不行,thanks Function dy=myfun(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=2y(2)-2y(1)+5exp(2x)*sin(x); [x, y]=ode23(myfun,[0,1],[-2,-3])；

2009-08-09, 23:17	#2
在水一方257 初级会员注册日期: 2009-05-10 年龄: 35 帖子: 1 声望力: 0	回复: 帮忙看看一个常微分方程的求解程序最后一句话后面的分号改为英文状态下的分号

2009-08-10, 12:38	#3
lxcwym 初级会员注册日期: 2009-08-08 年龄: 45 帖子: 8 声望力: 0	回复: 帮忙看看一个常微分方程的求解程序还是不行啊？麻烦各位帮我再看看吧，比较急！

2009-08-12, 20:50	#4
haozhou 初级会员注册日期: 2009-08-12 年龄: 40 帖子: 3 声望力: 0	回复: 帮忙看看一个常微分方程的求解程序 function dy=myfun(x,y) dy=[y(2);2y(2)-2y(1)+5exp(2x)*sin(x);] tspan=[0,1]; y0=[-2,-3]; [x, y]=ode23(@myfun,tspan,y0); y1=y(:,1); y2=y(:,2); plot(x,y1,x,y2);就可以了

2009-08-13, 12:41	#5
lxcwym 初级会员注册日期: 2009-08-08 年龄: 45 帖子: 8 声望力: 0	回复: 帮忙看看一个常微分方程的求解程序非常感谢！可以了，呵呵