% FIGURE10-11.TXT % This text file contains the data generated in MATLAB to produce Figure 11 in the paper. Note that Figure % 10 was produced in a similar manner but without any response trajectory being plotted. % The time, total displacement and total velocity data from the Figure7-9.txt data set was used in % conjunction with a data set called switch information matrix, which is included below for the X=10 % and Omega = 50 rad/s case. first_cycle = 0; last_cycle = 5 ; freq=50; ymin=-1.5e-3; ymax = 1e-3; zmin=-0.1; zmax=0.1; hold on; for cycle_num = first_cycle:last_cycle hSurface = patch([(cycle_num*2*pi/freq),(cycle_num*2*pi/freq),(cycle_num*2*pi/freq),(cycle_num*2*pi/freq)], [ymin,ymax,ymax,ymin], [zmin,zmin,zmax,zmax],'g') set(hSurface,'FaceColor',[0.5 0 0],'FaceAlpha',0.1); end for cycle_num = first_cycle:last_cycle-1 hSurface = patch([(cycle_num*2*pi/freq),(cycle_num*2*pi/freq),((cycle_num+1)*2*pi/freq),((cycle_num+1)*2*pi/freq)], [0 0 0 0], [zmin,zmax,zmax,zmin],'g') set(hSurface,'FaceColor',[0.5 0.5 0.5],'FaceAlpha',0.5); end box on; grid on; ax = gca; ax.BoxStyle = 'full'; xlim([(first_cycle*2*pi)/freq (last_cycle*2*pi)/freq]); ylim([ymin ymax]); zlim([zmin zmax]); set(gca,'YTick',[-1.5e-3 0 1e-3]) plot3(time,total_disp_resp,total_vel_resp'k','LineWidth',2);hold on; for cycle_num = first_cycle:last_cycle plot3(time((100*cycle_num)+1),total_disp_resp((100*cycle_num)+1),total_vel_resp((100*cycle_num)+1),'bo','MarkerSize',6,'MarkerFaceColor','b') end [row_num,col_num]=size(switch_information_matrix); for row = 1 : 2 : row_num plot3(switch_information_matrix(row,2),0,switch_information_matrix(row,5),'go','MarkerSize',6,'MarkerFaceColor','g') end for row = 2 : 2 : row_num plot3(switch_information_matrix(row,2),0,switch_information_matrix(row,5),'ro','MarkerSize',6,'MarkerFaceColor','r') end xlabel('Time (s)'); ylabel('Displacement (m)'); zlabel('Velocity (m/s)') set(gca,'FontName','Times New Roman','FontSize',16) view(-37.5,30) pbaspect([5,1,1]) % Switch information data for 10 forcing cycles of X=10N, Omega=50rad/s example. % Column 1 - active system. Column 2 - Switch Time. Column 3 - Number of complete % cycles at switch point. Column 4 - Switch Time from start of forcing cycle. % Column 5 - System velocity at switch point. 2.0000 0.0339 0 0.0339 -0.0860 1.0000 0.0980 0 0.0980 0.0558 2.0000 0.1608 1.0000 0.0352 -0.0400 1.0000 0.2224 1.0000 0.0968 0.0808 2.0000 0.2878 2.0000 0.0364 -0.0545 1.0000 0.3481 2.0000 0.0968 0.0737 2.0000 0.4122 3.0000 0.0352 -0.0506 1.0000 0.4738 3.0000 0.0968 0.0750 2.0000 0.5391 4.0000 0.0364 -0.0516 1.0000 0.5994 4.0000 0.0968 0.0751 2.0000 0.6648 5.0000 0.0364 -0.0518 1.0000 0.7251 5.0000 0.0968 0.0750 2.0000 0.7904 6.0000 0.0364 -0.0517 1.0000 0.8507 6.0000 0.0968 0.0750 2.0000 0.9161 7.0000 0.0364 -0.0517 1.0000 0.9764 7.0000 0.0968 0.0750 2.0000 1.0418 8.0000 0.0364 -0.0517 1.0000 1.1021 8.0000 0.0968 0.0750 2.0000 1.1674 9.0000 0.0364 -0.0517 1.0000 1.2277 9.0000 0.0968 0.0750 2.0000 1.2931 10.0000 0.0364 -0.0517 set(gcf,'Position',[277.8000 252.6000 500 300])