EOF分析matlab代码

%自己编的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;