为什么用ode45求解时位移ss和速度vv是零？
为什么用ode45求解时位移ss和速度vv是零？是不是计算精度问题？请赐教；方程
1
程序：
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;
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
q(n+1)=n;
if mod(n,2)==0; %判断奇数偶数
tspan=;
initial=;
options = odeset('RelTol',1e-6,'AbsTol',1e-14);
=ode45(@evens,tspan,initial,options,q(n+1));
else
tspan=;
initial=; %本次计算的初值是上一次计算的末值
options = odeset('RelTol',1e-6,'AbsTol',1e-14);
=ode45(@odds,tspan,initial,options,q(n+1));
end
tt(:,n+1)=t; %取出时间
dim=size(y);
ss(:,n+1)=y(:,1); %取出位移
vv(:,n+1)=y(:,2); %取出速度
end
回复 楼主 的帖子
简单看了一下，你没有定义odds函数吧？回复 2楼 的帖子
抱歉，odds和evens函数就是程序中的odd和even函数，我发程序时漏掉一个s了，不过把补上后结果仍然是零。
页:
[1]