版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

yuefengwang

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 307.6
帖子: 30
在线: 12.3小时
虫号: 2447389

[交流] 关于下面这个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);

回复此楼

» 猜你喜欢

268求调剂已经有9人回复
301求调剂已经有4人回复
0703化学336分求调剂已经有5人回复
301求调剂已经有4人回复
材料专硕306英一数二已经有9人回复
304求调剂已经有8人回复
材料工程专硕调剂已经有5人回复
344求调剂已经有4人回复
326求调剂已经有6人回复
296求调剂已经有4人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2013-09-05 13:11:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

C语言新手求助，关于数据类型转换问题！求指教！已经有10人回复
matlab模拟菲涅尔直边衍射已经有9人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 yuefengwang 的主题更新