请教两自由度碰撞系统分岔图

chenhu2031

刚刚看了几篇关于两自由度碰撞系统的一些文章。在画分岔图中不知道，碰撞条件怎样加进去，初始位移和速度是怎样设定的.如下面的一段程序，请各位多多指教！
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

chenhu2031

有谁做过这个方面的，请指教

天天狂奔

了解一下，

                               
登录/注册后可看大图

学习一把。

chenhu2031

没有人吗

chenhu2031

大神在哪

chenhu2031

人勒?....

张张

你做出来了吗？想请教一下