看看这个代码错在哪里？

xsj3917 · 发表于 2006-10-18 13:53

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

％ FUN2.M 带有数据集合的时变微分方程的例子
function  xdot = fun2 ( t, x, beta, omega, T, P )
pt = interp1 ( T, P, t );
xdot (2 ) = beta * x (2 ) – omega^2 * x (1) + pt ;
xdot (1 ) = x (2 );
xdot = xdot ( : );  % 使得xdot是一列
％ fun2.m 结束

在MATLAB中采用如下的形式调用这个函数：
beta = .1 ;
omega = 2 ;
A = .1 ;
w0 = 1.3 ;
theta = pi /4 ;
X0 = [ 0 ;1];
t0 = 0 ;
tf = 20 ;
T = t0 – eps : 0.1 : tf + theta +eps ;
P = A.* sin(w0 .*T –theta );
[ t, y ] = ode23 ( @ fun2 , [ t0, tf ], X0, [ ], theta, omega, T, P );
plot ( t, y (:,1));
hold on;
plot(t,y(:,2),'g');

运行提示错误：
??? T = t0 – eps : .1 : tf + theta +eps ;
         |
Error: Missing variable or function.

无水1324 · 发表于 2006-10-18 22:19

T = t0 – eps : 0.1 : tf + theta +eps ;
xdot (2 ) = beta * x (2 ) – omega^2 * x (1) + pt ;

两个式子的减号不对！
我试了一下，盖了之后就好了

xsj3917 · 发表于 2006-10-19 09:38

非常感谢，我说怎么老找不到错误呢

账号		自动登录	找回密码
密码			我要加入

[数学理论] 看看这个代码错在哪里？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。