clc;
clear;
close?all;
%?WaveData=load（‘E:\數據\美國軸承庫\外圈\OR007X130-1796-6.0-12k.mat‘）;???%107
%?WaveData=WaveData.X130_DE_time;
%??sampleFreq???=?12*1000;???
%?Data=load（‘E:\形態學\STREAM20090804-104326-859--右下側--1A-2H--空爐搖爐2.wfs-2-2.5-4.5s.txt‘）;
%?ReSampleFre=2;
%?WaveData=decimate（DataReSampleFre）;
%?sampleLength?=?length（WaveData）;
%?sampleFreq=500000/ReSampleFre;
%?===============================
%?WaveData=load（‘E:\數據\美國軸承庫\48K\OR007X151-1730-3.3-48k.mat‘）;?
%?WaveData=WaveData.X151_DE_time;
%?sampleFreq???=?48*1000;?????
%=================
WaveData=load（‘C:\Users\huhuan\Desktop\大脫硫風機的數據\H_Fault_f25600_n16384.txt‘）;
%?WaveData=[WaveData（:2）;WaveData（:4）;?WaveData（:6）;WaveData（:8）];
sampleFreq???=?10000*2.56;
%?%==================================================================
%?WaveData=WaveData;
sampleLength=length（WaveData）;
%========================仿真信號==============
%?sampleFreq=100000
%?t=0:1/sampleFreq:1/32;
%?s1=0.05*exp（-2*pi*80*t）.*sin（2*pi*1000*t+pi）;
%?%s2=zeros（1length（s1））;
%?for?i=1:3
%?????s1=[s1?s1];
%?end
%?sampleLength?=?length（s1）;
%?time_plot_s1?=?0:1/sampleFreq:（sampleLength-1）?/?sampleFreq;?
%?plot（time_plot_s1s1）;
%?
%?%疊加噪聲
%?SignalLength?=?length（s1）;
%?%計算信號的標準差
%?Signal_STD?=?std（s1）;
%?%白噪聲信號?
%?NoiseLevel?=0.2;
%?init=2055415866;
%?randn（‘seed‘init）;
%?Whitenoise_Signal=NoiseLevel?*?Signal_STD?*?randn（1sampleLength）;
%?
%?t=0:1/sampleFreq:（1/sampleFreq）*（length（s1）-1）;
%?s2=0.05*sin（2*pi*100*t）;
%?WaveData1=s1+Whitenoise_Signal;%+s2;
%?WaveData=s1;
%==========================================================================
%?%?%高通濾波
%?[ba]????????=butter（40.8‘high‘）?;?%尖峰能量5K---65K
%?filter_data?=?filter（baWaveData）;
%?%filter_data=abs（filter_data1）;
%?WaveData=filter_data;
%?sampleLength=length（WaveData）;
%?fft_result?=?abs（fft（WaveData））?*?2?/?sampleLength;
%?
%?%畫圖的坐標變換以此為準
time_plot_s?=?（0:（sampleLength-1））/sampleFreq;?
fft_plot_Hz?=?（1:sampleLength/2-1）*sampleFreq/sampleLength;?%注意0頻分量，不能從1開始，如果不要直流分量，可從1開始
%?????????????????????????????????????????????????????????????%為了和去掉直流分量后fft
%?????????????????????????????????????????????????????????????%結果（從2開始）維數一樣，此
%?????????????????????????????????????????????????????????????%處減1
%?????????????????????????????????????????????????????????????
%????????????????????????????????????????????????????????????
%?
%?figure（2）;
%?subplot（211）
%?plot（time_plot_sWaveData）
%?xlabel（‘s‘）;
%?title（‘時域波形‘）;
%?subplot（212）
%?plot（fft_plot_Hzfft_result（2:sampleLength/2））
%?title（‘頻域波形‘）;
%?xlabel（‘Hz‘）;
%?hold?on;
%?%做包絡分析
%==========================================================================
%?%?帶通濾波
%?[ba]????????=butter（4[0.0390.1172]）?;?%共振區
%?filter_data?=?filter（baWaveData）;
%?%filter_data=abs（filter_data1）;
%?WaveData=filter_data;
%?
%?%包絡
%?envelop_hil?=?hilbert（WaveData）;
%?envelop_abs?=?abs（envelop_hil）;??
%?%fft變換
%?envelop_fft?=?abs（fft（envelop_abs））*2?/sampl