%%%%%%%%%%-----FDD?Programe----%%%%%%%%%%%%?
?
clear?all
close?all

load?sec1?data;%?load??data?
AY0=?data（:[1?2?7?8?3:6]）;???%列向量
%?[BA]?=?butter（40.99‘low‘）;?
%?AY????=?filter（BAAY0）;?
Wp=20/256;Ws=30/256;Rp=1;Rs=20;???%通帶20Hz，阻帶30Hz
[nWn]=buttord（WpWsRpRs）;?????%計算濾波器的階數n和濾波參數
[filt_numfilt_den]?=?butter（nWn）;??%濾波器設計
AY=?filter（filt_numfilt_denAY0）;???%濾波操作
AY（:[468]）=?AY（:[468]）.*（-1）;?%?取相反數
AY=?detrend（AY）;???%去除線性分量
fs=?512;????%?sampling?frequency?
NFFT=?4096;??%?number?of?FFT?
channel=[1:8];?
G0=?zeros（（NFFT/2+1）length（channel）^2）;?%----為相關預留空間
j2=?0;?
for?i=1:length（channel）????%計算各個相關系數
????for?j=1:length（channel）?
????????j2?=?j2+1;?
????????[Pff2f1]=?cpsd（AY（:channel（i））AY（:channel（j））NFFTNFFT/4NFFTfs）;?
????????G0（:j2）=?Pff2;?
?????end?
end?
G=?zeros（length（channel）length（channel））;?%奇異值分解，然后處理數據
for?i=1:（NFFT/2+1）??????????????????????
????for?i2=1:length（channel）?
????????G（i2:）=?G0（i（（i2-1）*length（channel）+1）:（（i2-1）*length（channel）+length（channel）））;?
????end?
????[USV]=?svd（G）;?
????UU（（（i-1）*length（channel）+1）:（（i-1）*length（channel）+length（channel））:）?=?U;?
????SS（（（i-1）*length（channel）+1）:（（i-1）*length（channel）+length（channel））:）?=?S;?
end?

for?i=1:length（channel）?
????for?j=1:（NFFT/2+1）?
???