clear?all
close?all
fs=1024
K=64;??%段數
nfft=512;
f=fs*（0:nfft/2-1）/nfft;?%歸一化頻率
?
load?E:\datadecompose\LMD\gear\gearlw\CQLB11c1.dat
?
LB111=[CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1];
%?LB11=[CQLB11c1]
?x=LB111;
%??xx=LB11
%??ft（xx1024）
%?%??
%??[pFXbl]=myblp（xx1024）
%?%?x=FXbl;
N=length（x）
??
M=fix（N/K）?%?每段的數據長度??每段數據長128?
for?i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
%計算累積量
?
for?t=-（nfft/2-1）:nfft/2?????
c1（t+nfft/2）=third_slice_cumulant（xkt）;??%三階累積量對角切片
end
?
%?FFT變換
%?subplot（211）?
C1=c1.*hann（length（c1））‘;??
Y=fft（C1nfft）;???%對三階累積量對角切片求FFT
P1=Y.*conj（Y）/nfft;???%%%??P1
%?f=fs*（0:nfft/2-1）/nfft;?%歸一化頻率
%?plot（fP1（1:nfft/2））;????%畫1.5維譜圖
%?xlabel（‘頻率‘）
%?ylabel（‘幅值‘）
?
?
load?E:\datadecompose\LMD\gear\gearlw\CQLB11c2.dat
LB112=[CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2];
%?LB12=[CQLB11c2];
?
x=LB112;
%?N=length（x）
%?K=64??????%段數
%?M=fix（N/K）?%?每段的數據長度??每段數據長128?
?
for?i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
?????????
%計算累積量
?
for?t=-（nfft/2-1）:nfft/2?????
c2（t+nfft/2）=third_slice_cumulant（xkt）;??%三階累積量對角切片
end
?
%?FFT變換
%?subplot（211）?
C2=c2.*hann（length（c2））‘;??
Y=fft（C2nfft）;??%對三階累積量對角切片求FFT
P2=Y.*conj（Y）/nfft;???%%%??P1
%?f=fs*（0:nfft/2-1）/nfft;?%歸一化頻率
%?plot（fP2（1:nfft/2））;????%畫1.5維譜圖
%?xlabel（‘頻率‘）
%?ylabel（‘幅值‘）
?
?
load?E:\datadecompose\LMD\gear\gearlw\CQLB11c3.dat
LB113=[CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3];
%?LB13=[CQLB11c3];
x=LB113;
%?N=length（x）
%?K=64??????%段數
%?M=fix（N/K）?%?每段的數據長度??每段數據長128?
?
for?i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
?????????
%計算累積量
for?t=-（nfft/2-1）:nfft/2?????
c3（t+nfft/2）=third_slice_cumulant（xkt）;??%三階累積量對角切片
end
?
%?FFT變換
%?subplot（211）?
C3=c3.*hann（length（c3））‘;??
Y=fft（C3nfft）;???%對三階累積量對角切片求FFT
P3=Y.*conj（Y）/nfft;???%%%??P1
%?f=fs*（0:nfft/2-1）/nfft;?%歸一化頻率
%?plot（fP3（1:nfft/2））;????%畫1.5維譜圖
%?xlabel（‘頻率‘）
%?ylabel（‘幅值‘）
?
?
?
load?E:\datadecompose\LMD\gear\gearlw\CQLB11c4.dat
LB114=[CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4];
%?LB14=[CQLB11c4];
?
x=LB114;
%?N=length（x）
%?K=64??????%段數
%?M=fix（N/K）?%?每段的數據長度??每段數據長128?
?
for?i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
?????????
%計算累積量
?
for?t=-（nfft/2-1）:nfft/2?????
c4（t+nfft/2）=third_slice_cumulant（xkt）;??%三階累積量對角切片
end
?
%?FFT變換
%?subplot（211）?
C4=c4.*hann（length（c4））‘;??
Y=fft（C4nfft）;???%對三階累積量對角切片求FFT
P4=Y.*conj（Y）/nfft;???%%%??P1
?
X=[c1;c2;c3;c4]
SLB11=svd（X）;


%?clear?all
%?close?all
%?fs=1024
%?K=64;??%段數
%?nfft=512;
%?f=fs*（0:nfft/2-1）/nfft;?%歸一化頻率
?
load?E:\datadecompose\LMD\gear\gearlw\CQLB12c1.dat
?
LB121=[CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1];
%?LB12=[CQLB12c1]
?x=LB121;
%??xx=LB12
%??ft（xx1024）
%?%??
%??[pFXbl]=myblp（xx1024）
%?%?x

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????42213??2013-04-11?11:10??hospeciPRG\QDqiepian.m

?????文件????????683??2013-04-11?11:39??hospeciPRG\SQD.mat

?????文件??????22174??2013-04-11?12:46??hospeciPRG\ZCqiepian.m

?????文件????????211??2013-03-29?21:45??hospeciPRG\third_slice_cumulant.m

?????文件????????619??2013-04-11?11:39??hospeciPRG\SLB.mat

?????文件????????415??2013-04-11?14:35??hospeciPRG\SZC.mat

?????文件??????37202??2013-04-11?11:39??hospeciPRG\LBqiepian.m

?????目錄??????????0??2013-04-11?13:24??hospeciPRG

-----------?---------??----------?-----??----

???????????????103517????????????????????8