|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
其中的无穷大,编程时我取了1000
function dx=eq9(t,x)
dx=zeros(2,1);
dx(1)=x(2);
a=0;
for n=1:1:1000
b=(sin(1.88*n)/n)*cos(0.0175*n*t);
a=a+b;
end
c=1+0.254*a;
dx(2)=0.74-3.3*10^(-4)*sqrt(c)*x(2)-c*x(1);
tspan=linspace(700,1800,10000);
[t,x]=ode45('eq9',tspan,[0 0]);
plot(t,x(:,1),'-');
得到的响应怎么会是这样的呢?求大家帮看看程序是不是有问题啊?谢谢。其中方程做过无量纲处理。啮合刚度是方波。
下面这张图是我把阻尼取得比较大的时候,得到的响应。
|
|