急！请问为什么这个lorenz方程的matlab程序一运行就出现如下问题

chenyland · 发表于 2006-5-22 13:59

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

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

x

程序如下： function xdot=lorenzeq(t,x) xdot=[-8/3*x(1)+x(2)*x(3);-8620*x(2)+8620*x(3);-x(1)*x(2)+144*10^7*x(2)-x(3)]; t_final=100;x0=[0;0;1e-10]; [t,x]=ode45('lorenzq',[o,t_final],x0);plot(t,x), figure; plot(x(:,1),x(:,2),x(:,3)); 一运行就出现如下提示： ??? Input argument 'x' is undefined.
Error in ==> C:\MATLAB6p5\work\lorenzeq1.m On line 2 ==> xdot=[-8/3*x(1)+x(2)*x(3);-8620*x(2)+8620*x(3);-x(1)*x(2)+144*10^9*x(2)-x(3)]; 感觉应该没有问题啊，怎么老是这样呢？？请高手赐教

[此贴子已经被VibInfo于2006-5-22 14:40:06编辑过]

happy · 发表于 2006-5-22 14:39

<DIV class=quote>以下是引用chenyland在2006-5-22 13:59:36的发言： 
程序如下： function xdot=lorenzeq(t,x) xdot=[-8/3*x(1)+x(2)*x(3);-8620*x(2)+8620*x(3);-x(1)*x(2)+144*10^7*x(2)-x(3)]; t_final=100;x0=[0;0;1e-10]; [t,x]=ode45('lorenzq',[o,t_final],x0);plot(t,x), figure; plot(x(:,1),x(:,2),x(:,3)); 一运行就出现如下提示： ??? Input argument 'x' is undefined.
Error in ==> C:\MATLAB6p5\work\lorenzeq1.m On line 2 ==> xdot=[-8/3*x(1)+x(2)*x(3);-8620*x(2)+8620*x(3);-x(1)*x(2)+144*10^9*x(2)-x(3)]; 感觉应该没有问题啊，怎么老是这样呢？？请高手赐教</DIV>
红色标记部分明显错误，这种错误自己看看就应该发现的 别出了点问题就问，第一时间应该自己检查 第二 function xdot=lorenzeq(t,x) xdot=[-8/3*x(1)+x(2)*x(3);-8620*x(2)+8620*x(3);-x(1)*x(2)+144*10^7*x(2)-x(3)]; 这一部分保存成lorenzeq.m t_final=100;x0=[0;0;1e-10]; [t,x]=ode45('lorenzeq',[0,t_final],x0);plot(t,x), figure; plot(x(:,1),x(:,2),x(:,3)); 在命令窗口运行或者也写成m文件并运行

chenyland · 发表于 2006-6-4 13:49

谢谢这位大侠，你说的先保存为m文件我也试过，可能还是lorenzq那里的问题，最近上不了振动论坛！今天才行

zhxllhy · 发表于 2008-3-11 23:08

clear
global a b c
a=16;b=4;c=45;
[T,X]=ode45('Lorenz',[0,30],[12;4;0]);
hold on
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');
hold off
function dx=Lorenz(t,x)
global a b c
dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];

??? Error: Function definitions are not permitted at the prompt or in scripts.
请问怎么修改?谢谢!!!!

sogooda · 发表于 2008-3-12 08:19

原帖由 zhxllhy 于 2008-3-11 23:08 发表

 登录/注册后可看大图

clear
global a b c
a=16;b=4;c=45;
[T,X]=ode45('Lorenz',[0,30],[12;4;0]);
hold on
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');
hold off
function dx=Lorenz(t,x ...

这个问题比较基础了，function 后面的部分应该另存为一个m文件，并且名字应该是Lorenz.m

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

[综合讨论] 急！请问为什么这个lorenz方程的matlab程序一运行就出现如下问题

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

回复：（chenyland）急！请问为什么这个lorenz方程的m...

请教高手

评分