Code covered by the BSD License  

Highlights from
PCAtool

image thumbnail
from PCAtool by Guillaume MAZE
Compute everything you need for EOF,EEOF,CEOF,SVD,lagged SVD

dispSVD(typg,CHP1,CHP2,N,expvar,LAGS,dt)
% dispSVDlag(typg,CHP1,CHP2,N,expvar,LAGS,dt) Display SVDs lagged
%
% => Display SVDs lagged
% typg is plot type you want (see plotm.m).
% CHP contains SVDs and is of the form: CHP(iSVD,X,Y)
% N is number of SVD to display
% expvar contained the explained variance of SVDs
%
% Rq: it''s using the plotm function
%
%================================================================

%  Guillaume MAZE - LPO/LMD - December 2004
%  gmaze@univ-brest.fr

   function [] = dispSVD(typg,CHP1,CHP2,N,expvar,LAGS,dt)

load mapanom2
global ffile1 ffile2
% --------------------------------------------------
% AFFICHAGE
% --------------------------------------------------
if (typg~=1)&(typg~=7)&(typg~=3)&(typg~=6)
   CBAROR='horiz';
else
   CBAROR='vert';
end

   CHP1=squeeze(CHP1(:,N,:,:));
   CHP2=squeeze(CHP2(:,N,:,:));

Nlag=size(CHP1,1);

iw=round(Nlag/dt);jw=2;


f1=figur;clf; hold on
ilag=0;
for lag=1:dt:Nlag
  ilag=ilag+1;

    subplot(iw,jw,2*ilag-1)
C = squeeze(real(CHP1(lag,:,:)));
plotm(C,typg);caxis([-abs(xtrm(C)) abs(xtrm(C))]);
if(ilag==1);cx=caxis;end;caxis(cx);
colormap(mycolormap(mapanom,21));
c=colorbar(CBAROR);
titre1=strcat('lag:',num2str(LAGS(lag)),'y: ',ffile1);
title(titre1);

     subplot(iw,jw,2*ilag)
C = squeeze(real(CHP2(lag,:,:)));
plotm(C,typg);caxis([-abs(xtrm(C)) abs(xtrm(C))]);
if(ilag==1);cx=caxis;end;caxis(cx);
colormap(mycolormap(mapanom,21));
c=colorbar(CBAROR);
titre1=strcat('lag:',num2str(LAGS(lag)),'y: ',ffile2);
title(titre1);

end %for iN

suptitle(strcat('Lagged SVD, mode:',num2str(N)));

Contact us at files@mathworks.com