刚刚看了几篇关于两自由度碰撞系统的一些文章。在画分岔图中不知道,碰撞条件怎样加进去,初始位移和速度是怎样设定的.如下面的一段程序,请各位多多指教!
clear;clc
wcanshu=0.01:0.01:4;
T=2*pi./wcanshu;
hold on
for k=1:length(wcanshu)
[t,y]=ode23('hdbfun',[0:T/100:70*T],[1;0;1;0],[],wcanshu(k));
plot(wcanshu(k),y(500:100:1400,2),'linewidth',5);
title('随参数变化的分岔图');
xlabel('随外界激励频率omega的变化');
ylabel('位移')
end