阶次分析matlab代码

%~~~~~~~~~~~~~导入数据部分~~~~~~~~~~~~~%
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（211）plot（tarray_time_amp）title（‘time dominant振动信号时域图‘）xlabel（‘时间time‘）ylabel（‘幅值amplitude‘）;
subplot（211）plot（tarray_time_amp）title（‘变转速数据‘）xlabel（‘时间time‘）ylabel（‘幅值amplitude‘）;
grid on;
% subplot（212）plot（tpluse）title（‘keyphasor键相脉冲仿真信号时域图‘）xlabel（‘时间time‘）ylabel（‘幅值amplitude‘）;
subplot（212）plot（tpluse）title（‘转速脉冲‘）xlabel（‘时间time‘）ylabel（‘幅值amplitude‘）;
grid on;


% %~