用ode45求解微分方程时为什么位移及速度值是0？
我用ode45解方程时位移ss和速度vv值结果均为零，不知为什么，附上方程和程序，请指教y(1,1)=0;y(1,2)=0;beta=0.4488;w=157.1;
N=ceil(20*pi/beta);
dim(1)=1;
for n=0:1:N %N＝140（相当于t=0.4s）
q(n+1)=n;
if mod(n,2)==0; %判断奇数偶数
tspan=;
initial=;
options = odeset('RelTol',1e-6,'AbsTol',1e-14);
=ode45(@even,tspan,initial,options,q(n+1));
else
tspan=;
initial=; %本次初值为上次计算末值
options = odeset('RelTol',1e-6,'AbsTol',1e-14);
=ode45(@odd,tspan,initial,options,q(n+1));
end
dim=size(y);
tt(:,n+1)=t; %提取时间
ss(:,n+1)=y(:,1); %提取位移
vv(:,n+1)=y(:,2); %提取速度
end
function dy=even(t,y,q)
dy=zeros(2,1);
cca=0.5;ch=81.2;kca=25392;tao6=1e+5;m=3.24;beta=0.4488;w=157.1;
dy(1)=y(2);
dy(2)=(-(ch*((cos(w*t-q*beta)).^2+(cos(w*t-q*beta+2*beta)).^2+(cos(w*t-q*beta+4*beta)).^2+(cos(w*t-q*beta+6*beta)).^2)+cca)*y(2)... -(tao6*((cos(w*t-q*beta)).^2+(cos(w*t-q*beta+2*beta)).^2+(cos(w*t-q*beta+4*beta)).^2+(cos(w*t-q*beta+6*beta)).^2)+kca)*y(1))/m;
function dy=odd(t,y,q)
dy=zeros(2,1);
cca=0.5;ch=81.2;kca=25392;tao6=1e+5;m=3.24;beta=0.4488;w=157.1;
dy(1)=y(2);
dy(2)=(-(ch*((cos(w*t-(q-1)*beta)).^2+(cos(w*t-(q-1)*beta+2*beta)).^2+(cos(w*t-(q-1)*beta+4*beta)).^2)+cca)*y(2)... -(tao6*((cos(w*t-(q-1)*beta)).^2+(cos(w*t-(q-1)*beta+2*beta)).^2+(cos(w*t-(q-1)*beta+4*beta)).^2)+kca)*y(1))/m; 发一张方程的图片
页:
[1]