马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
请教一下各位高手,我照着书本键入了以下的程序,但是得到的结果却于书上的差距很大,请各位帮忙看看,谢谢了!!!
function dx=jisi(t,x)
a=20;b=40;c=15;
s=sqrt((c-x(1))^2+(a*t-x(2)^2));
dx=[b*(c-x(1))/s;b*(a*t-x(2))/s];
ts=0:0.1:2;
x0=[0,0];
opt=odeset('reltol',1e-2,'abstol',1e-3);
[t,x]=ode45(@jisi,ts,x0,opt);
[t,x],
plot(t,x),grid,
gtext('x(t)'),gtext('y(t)'),pause
plot(x(:,1),x(:,2)),grid,
gtext('x'),gtext('y')
作出的y(x)图形书本上最后是沿着x=15竖直向上的,没有回转,而且当t=1时,x(t)=3.9561,和实际的运算结果也不一样,是什么原因造成的呢?麻烦各位帮忙看看,不胜感激!!!
[ 本帖最后由 eight 于 2007-8-12 10:52 编辑 ] |