%~~~~~~~~~~~~~導入數據部分~~~~~~~~~~~~~%
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;


%?%~