- 1、本文档共51页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
EOF局限性 与区域范围选择有关,这给进行物理解释带来困难 过分强调与分析区域尺度相当的现象,因而使区域性相关结构反映不出来。 4.4 旋转EOF(REOF) REOF 基本思想 如果EOF分析截取了前K个空间型,累计解释场的总方差已达一定要求,例如累计方差贡献率达80%,可否(及怎样)将这K个空间型再做调整,使得调整后的K个空间型累计解释原场总方差的百分率保持不变,而单个空间型尽量反映场的局部相关结构。 这种对K个EOF再做调整的过程在数学上是在它们之间再做线性变换,几何形象是做旋转,所以称为REOF。 EOF分解 分量形式 % [eof,V] = varimax(eof[,sign,tol]) Varimax rotation % % Performs varimax rotation (Kaiser 1958) on the column vectors contained % in eof. We follow Kaisers notation. % % In: % eof: DxL matrix (LD) of column vectors (possibly the loadings from % a factor analysis). % sign: preferred majoritary sign for the vector components (-1: % negative, 1: positive, 0: dont care). Default 0. % tol: minimum relative increase of the objective function to keep on % iterating (default 1e-5). % Out: % eof: DxL matrix with the rotated loadings. % V: training curve for the objective function. % Copyright (c) 1997 by Miguel A. Carreira-Perpinan function [eof,V] = varimax(eof,sign,tol) % Argument defaults if nargin==1 sign=0; end; if nargin=2 tol=1e-5; end; [D,L] = size(eof); h = sqrt(sum(eof.^2))+exp(-700); % Communalities temp=eof./(h*ones(1,L)); V = [sum(sum(temp.^4))-sum(sum(temp.^2).^2)/D]; % Objective function V_old = V*(1-2*tol); while abs(V(length(V))-V_old) tol*V(length(V)) V_old = V(length(V)); for i=1:L-1 for j=i+1:L % Optimal angle to rotate columns i, j x = eof(:,i)./h; y = eof(:,j)./h; u = x.*x - y.*y; v = 2*x.*y; t = atan2( 2*(D*u*v-sum(u)*sum(v)), D*(u*u-v*v)-sum(u)^2+sum(v)^2 )/4; % Anticlockwise rotation of angle t (t+pi is valid, too) temp = [eof(:,i) eof(:,j)]*[cos(t) -sin(t); sin(t) cos(t)]; eof(:,i) = temp(:,1); eof(:,j) = temp(:,2); end end % New value of the objective function h = sqrt(sum(eof.^2))+exp(-700); % Communalities temp=eof./(h*ones(1,L)); V = [V sum(sum(temp.^4))-sum(sum(temp.^2).^2)/D]; % Objective function end % Sign inversion so that each column vector of eof has
您可能关注的文档
- 第四讲 重庆维普数据库.ppt
- 第一章 饭店服务中的人际交往心理20160315.ppt
- 第一章 工程定额的概念.ppt
- 第四讲 做情绪的主人.ppt
- 第一章 集成电路制造工艺.ppt
- 第四讲(魏晋南北朝文学).ppt
- 第四讲.金融创新和监管.ppt
- 第一章 楷书的学习.ppt
- 第四讲:高能束加工.ppt
- 第一章 运动系统2.doc
- 2024年05月山东交通职业学院招考聘用博士研究生50人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月安徽芜湖市弋江区老年学校(大学)工作人员特设岗位公开招聘2人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月山东东营河口区教育类事业单位招考聘用22人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月山东交通职业学院招考聘用100人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月山东威海职业学院招考聘用高层次人才2人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月安徽石台县事业单位工作人员33人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月山东滨州市博兴县事业单位公开招聘考察笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月安徽蚌埠固镇县湖沟镇选聘村级后备干部7人笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月山东省安丘市教育和体育局所属事业单位学校公开2024年招考232名工作人员笔试历年典型题及考点剖析附带答案含详解.docx
- 2024年05月山东临沂临港经济开发区工作人员(5人)笔试历年典型题及考点剖析附带答案含详解.docx
文档评论(0)