松动与耦合故障

zhuchang1121

这是我编的关于碰摩与松动耦合故障的程序(还没有进行无量纲化,没有加油墨力)大家看看有什么问题?
function dx=dafen1(t,x)
m1=32.1;
m2=50;
c=1050;
k=0.3*10^6;
w=1400;
kr=0.65*10^8;
c(31)=0.15;
c(33)=0.15
c(32)=0.05;
k(31)=50;
k(33)=50
k(32)=1.0;
u=0.032;
g=9.8;
f=0.12;
n=0.001;
q=0.0005;
w=1400;
e=sqrt(x(1)^2+x(2)^2);
if  e>=n
    N=kr*(e-n);
    Fx=-N*(x(1)-f*x(2))/e;
    Fy=-N*(f*x(1)+x(2))/e;
else
    N=0;
    Fx=0;
    Fy=0;
end
if x(3)< 0
c3=c(31);
k3=k(31);
else if (x(3)>0)&&(x(3)<=q)
c3=c(32);
k3=k(32);
else
c3=c(33);
k3=k(33);
end
dx=zeros(6,1);
dx(1)=x(4);
dx(4)=Fx/m1+u*w^2cos(w*t)-k*x(1)/m1-c*x(4)/m1;
dx(2)=x(5);
dx(5)= Fy/m1+u*w^2sin(w*t)-k*(x(2)-x(3))/m1-c*(x(5)-x(6))/m1-g;        
dx(3)=x(6);
dx(6)= k*(x(2)-x(3))/m2+c*(x(5)-x(6))/m2-g-c3*x(6)/m2-k3*x(3)/m2;
出现了这样的错误
??? Error: File: C:\MATLAB6p5\work\dafen1.m Line: 42 Column: 18
Missing operator, comma, or semicolon.
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\private\odearguments.m
On line 76  ==>     if (nargin(ode) == 2)           
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\ode45.m
On line 155  ==> [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0,
大家看看什么问题?谢谢我的QQ40942395

sigma665

估计你调用ode45格式错了

zhuchang1121

x0=[1,1,1,1,1,1];
>> options=odeset('RelTol',1e-5);
>> T=2*pi/1400;
>> [t,x]=ode45('dafen9',[1000*T:T/512:1100*T],z0,options);
这是调用程序.看不出错误啊!

sigma665

if end不匹配

建议选择smart index

ctrl+I


还有elseif，不是else if