关于子程序间参数的调用疑问

bryant

程序作用：求解二阶二元微分方程
主程序：
global e xi Omega omega;
e=0.1;xi=0.01;Omega=0.5;omega=1;
[t,u]=ode45('ode3',[0,12],[0 0 0 0]);
子程序1：
function du=ode3(t,u)
global e xi Omega omega
du=zeros(4,1);%u1=x,u2=x',u3=y,u4=y'
du(1)=u(2);
du(2)=fcrack*k/Omega^2*sin(2*omega*t)/2*u(3)-2*xi/Omega*u(2)-(1/Omega^2-fcrack*k/Omega^2*cos(omega*t)^2)*u(1)+1/Omega^2+e*cos(omega*t);
du(3)=u(4);
du(4)=fcrack*k/Omega^2*sin(2*omega*t)/2*u(1)-2*xi/Omega*u(4)-(1/Omega^2-fcrack*k/Omega^2*sin(omega*t)^2)*u(3)+e*sin(omega*t);
子程序2：
function f=fcrack(alpha,theta,phi)
global e xi Omega omega
y=u(:,3),x=u(;,1);
phi=atan(y/x);
beta=0;
theta=mod(omega*t+beta-phi,2*pi);
if -pi/2+alpha<=theta&theta<=pi/2-alpha
    f=1;
elseif pi/2-alpha<=theta&theta<=pi/2+alpha
    f=1/2+1/2*cos((theta-pi/2+alpha)*pi/(2*alpha));
elseif pi/2+alpha<=theta&theta<=3*pi/2-alpha;
    f=0
elseif 3*pi/2-alpha<=theta&theta<=3*pi/2+alpha;
    f=1/2+1/2*cos((theta-3*pi/2-alpha)*pi/(2*alpha));
end
子程序1需要调用子程序2的frack，而子程序2又要调用子程序1的x跟y，即u1和u3
按上述程序运行时出现：Expected a variable, function, or constant, found ";".
请问如何改动程序才能正确调用

sogooda

怎么fcrack既是变量名又是函数名呢？
可以把你需要的参数都传递给主函数，然后再由主函数传递给所需的子函数。

eight

太混乱了，建议楼主补补编程基础，不可能循环调用的

ch_j1985

原帖由 bryant 于 2008-5-5 20:17 发表

                               
登录/注册后可看大图

程序作用：求解二阶二元微分方程
主程序：
global e xi Omega omega;
e=0.1;xi=0.01;Omega=0.5;omega=1;
[t,u]=ode45('ode3',[0,12],[0 0 0 0]);
子程序1：
function du=ode3(t,u)
global e xi Omega omega
...

楼主能否把微分方程帖一下？

bryant

是有点混乱，现在就是需要大家指点指点

bryant

那请问大概要怎样修改才能正确调用，谢谢指教