其中:M=2 C=200 K=2X10^5 e=0.04,k=0,1,2...1000,A=14500,w=2.4
上面的方程按照下面这样表达错在哪里了?
function dx=mybit(t,x);
for k=0:10*2*pi
if 2*k*pi/2.4<=t&t<(2*k*pi+pi)/2.4
dx=[x(2);-150*x(2)-10^5*(x(1)-0.3).*(x(1)>0.3)-0.*(abs(x(1))<0.3)-10^5*(x(1)+0.3).*(x(1)<-0.3)+14500*sin(2.4*t)];
else
dx=[x(2);-150*x(2)-10^5*(x(1)-0.3).*(x(1)>0.3)-0.*(abs(x(1))<0.3)-10^5*(x(1)+0.3).*(x(1)<-0.3)];
end
function dx=mybit(t,x);
for k=0:10*2*pi
%这里的k取值应该是k=0:1000吧?
if 2*k*pi/2.4<=t&t<(2*k*pi+pi)/2.4
dx=[x(2);-150*x(2)-10^5*(x(1)-0.3).*(x(1)>0.3)-0.*(abs(x(1))<0.3)-10^5*(x(1)+0.3).*(x(1)<-0.3)+14500*sin(2.4*t)];
else
dx=[x(2);-150*x(2)-10^5*(x(1)-0.3).*(x(1)>0.3)-0.*(abs(x(1))<0.3)-10^5*(x(1)+0.3).*(x(1)<-0.3)];
end
%另外方程与你给的公式对不上
个人拙见,不一定正确