求各位高手修改修改此程序，不知道是那里出问题了

熊猫酒仙 · 发表于 2012-12-2 19:59

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

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

x

clear;clc;
function dy=dfun(t,y)
dy=zeros(4,1);
r=0.1+0.01*cos(2*t);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(9.8*sin(y(1))+r*y(4).^2*sin(y(1))*cos(y(1))+0.04*y(2)*sin(2*t))/r;
dy(4)=(0.04*y(4)*sin(2*t)*sin(y(1))-2*r*y(1)*y(4)*cos(y(1)))/(r*sin(y(1)));
Y0=[4*pi/3,0,0,0.01];
[T,Y]=ode45('dfun',[16 21],Y0);
plot(T,Y)

不知道哪里出了问题，程序运行不出来，求高人相助。

happy · 发表于 2012-12-3 09:56

脚本文件中不能定义函数
可以将函数独立写成函数文件，脚本文件调用
或者把上述脚本文件改写成函数文件

liuyong · 发表于 2012-12-5 05:52

好像function的变量不行吧

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

[编程技巧] 求各位高手修改修改此程序，不知道是那里出问题了

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