phase_track_modify.m

该程序在MATLAB环境下编译

%% 主程序
clf;
%%%%%%%%%%%%%%% load in data file %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% data = load('database.dat');      % load in data file
f_name_1='s_17.txt'; %转速、进给速度一定时，进给量由小变大进行比较
f_name_2='s_14.txt';
f_name_3='s_15.txt'; 
f_name_4='s_16.txt';
f_name_5='s_26.txt'
start=[5000, 9000, 7000, 6000];
fname=[f_name_5];
file_number=length(fname(:,1));
%%%%%%%%%%%%%%%%%%%%%%%%%%%
for f_sequence=1:file_number
   begin=500         % 开始行
   % numline=2048;   % 总行数
   numline=512;     % 总行数
   Nocolume=2;       % 所选列,1、2、3分布代表x、y、z方向的加速度
  x=FLRead(fname(f_sequence,:),begin,numline,Nocolume)';%%%%
end
   [C,L]=wavedec(x,3,'db4');         %%%%%%%%%%%%小波分解              
   data=wrcoef('a',C,L,'db4');       %%%%%%%%%%%小波重构振动信号
   N = length(data);                 % calculate number of data points
nextdata_y=zeros(N,1);               % create another array for return map
nextdata_z=zeros(N,1);
tau=10;
for (i=1:N-2*tau)
  nextdata_y(i)=data(i+tau);         % 构造Y方向延迟点，延迟步长为tau
  nextdata_z(i)=data(i+2*tau);       % 构造Z方向延迟点，延迟步长为2*tau 
end

figure(1)                            % now plot it (except last point)
plot3(data(1:N-2*tau),nextdata_y(1:N-2*tau),nextdata_z(1:N-2*tau),'-.');
                                     % 绘制相轨迹图, blue
grid on
xlabel('x(i)','FontSize',18),ylabel('x(i+10)','FontSize',18),zlabel('x(i+20)','FontSize',18)
title('phase track map');            % give it a title
axis square; 
box on;
hold on;

pause(2);                             
for k = 1:(N-2*tau)
  if k> 1
     delete(hdk);
  end
  hdk = plot3(data(k),nextdata_y(k),nextdata_z(k),'ro',...
                    'MarkerEdgeColor','r',...
                    'MarkerFaceColor','r',...
                    'MarkerSize',5);
  pause(0.2)
end
hold off