为什么在命令窗口能运行，M文件却出错？

kougewang

代码如下：

1. 
   
2. function solveyyy
   
3. clc;clear;
   
4. [t,x]=ode45(@dync1,[0,20],[77;0.25;0.25]);
   
5. figure;plot(t,x)
   
6. % function dync=dync1(t,x)
   
7. % dync=[ 50/m-g*cos(x(1))/(x(2));...
   
8. %     -50/m-g*sin(x(1))];
   
9. function dync=dync1(t,x)
   
10. global m;m=546;
    
11. global g;g=9.8;
    
12. dync=[ -0.1-g*sin(x(2));
    
13.     0.1-g*cos(x(2))/x(1);
    
14.     -0.1/(x(1)*cos(x(2)))];

复制代码

请各位高手指点。

sigma665

文件名是不是数字？

kougewang

文件名字是dync1.m
如果是数字就会出错吗？

ChaChing

改solveyyy.m试试!

[ 本帖最后由 ChaChing 于 2009-2-12 16:19 编辑 ]

kougewang

可以运行了，多谢多谢楼上各位朋友。:@)
还有个问题，就是解出微分方程后，我想画x(2)关于x(3)的曲线，怎么实现？

kougewang

知道了,命令是plot(x(:,2),(:,3))