| 查看: 674 | 回复: 0 | |||
[交流]
关于下面这个eof程序，数据的输入格式（txt文档）改是如何？求指教
|
|
我不懂程序代码，正在认真学习，请帮忙解读下这个程序的过程，谢谢了。 clear; clc fid=fopen('C:\Users\Administrator\Desktop\test.txt','r'); hgtn=fread(fid,'float'); aay=reshape(hgtn,246,12,58); %a=mean(aay,3); bb=aay(:,6,  +aay(:,7,+aay(:,8,; ab=squeeze(bb); %资料的标准化 an=ab'; avn=mean(an); ssn=std(an,1); n=58; jp=an-repmat(avn,n,1); bn=jp./repmat(ssn,n,1); a=bn; %提取主分量！ %先进行时空转换 s=a*a'; [v,d]=eig(s); j=n; %将特征值按从大到小排序 for i=1:n dd(i,i)=d(j,j); j=j-1; end j=n; %将特征向量从大到小排序 for i=1:n for k=1:n vv(k,i)=v(k,j); end j=j-1; end sum=0.; %求累计方差 for i=1:n sum=sum+dd(i,i); end for i=1:n %求方差贡献 c(i)=dd(i,i)/sum; end %求得eof的距阵 x1=a'*vv(:,1); %x1为第一特征向量空间场 tri=vv(:,1); %tri为时间系数 fid=fopen('f:\ga\pt\wave\eof\1000v.dat','w'); out=fwrite(fid,x1,'float'); fclose(fid); fid=fopen('f:\ga\pt\wave\eof\ri.dat','w'); out=fwrite(fid,tri,'float'); fclose(fid); |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有26人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
点击这里搜索更多相关资源