%自己編的EOF程序
clear;?clc
fid=fopen（‘f:\ga\pt\wave\jijie\pt.dat‘‘r‘）;
hgtn=fread（fid‘float‘）;
aay=reshape（hgtn2461258）;???
%a=mean（aay3）;
bb=aay（:6:）+aay（:7:）+aay（:8:）;????????
ab=squeeze（bb）;

%資料的標準化
an=ab‘;
avn=mean（an）;
ssn=std（an1）;
n=58;
jp=an-repmat（avnn1）;?
bn=jp./repmat（ssnn1）;

a=bn;
%提取主分量！
%先進行時空轉換
s=a*a‘;
[vd]=eig（s）;

j=n;????????????????????????%將特征值按從大到小排序
for?i=1:n
????dd（ii）=d（jj）;
???j=j-1;
end

j=n;???????????????????