滚动轴承信号仿真阶次跟踪法分析

%~~~~~~~~~~~~~导入数据~~~~~~~~~~~~~%
clc
clear
Fs=71680;
N=Fs*5;
t=(0:N-1)/Fs;
Adc = 1; %直流分量幅度
% S=sin(2*pi*t.^2- pi/6)+sin(4*pi*t.^2- pi/6)+sin(8*pi*t.^2- pi/6)
% S2=Adc+ sin(2*pi* t.^2- pi/6);%参考轴的转速为n(t)=60t r/min

S=sin(2*pi*t.^2)+sin(4*pi*t.^2)+sin(8*pi*t.^2);
S2=Adc+ sin(2*pi* t.^2);


array_time_amp=S;    %导入时域振动信号
pluse=S2;             %导入脉冲信号
figure(1);
% subplot(2,1,1),plot(t,array_time_amp),title('time dominant振动信号时域图'),xlabel('时间time'),ylabel('幅值amplitude');
subplot(2,1,1),plot(t,array_time_amp),title('变转速数据'),xlabel('时间time'),ylabel('幅值amplitude');
grid on;
% subplot(2,1,2),plot(t,pluse),title('keyphasor键相脉冲仿真信号时域图'),xlabel('时间time'),ylabel('幅值amplitude');
subplot(2,1,2),plot(t,pluse),title('转速脉冲'),xlabel('时间time'),ylabel('幅值amplitude');
grid on;


% %~~~~~~~~~~~~~~~计算脉冲发生时刻~~~~~~~~~~~~~~~%
% Num_pluse1&