[21589] | 1 | % dispSVDnb(typg,CHP1,CHP2,N,expvar) Display SVDs
|
---|
| 2 | %
|
---|
| 3 | % => Display SVDs in white and black
|
---|
| 4 | % typg is plot type you want (see plotm.m).
|
---|
| 5 | % CHP contains SVDs and is of the form: CHP(iSVD,X,Y)
|
---|
| 6 | % N is number of SVD to display
|
---|
| 7 | % expvar contained the explained variance of SVDs
|
---|
| 8 | %
|
---|
| 9 | % Rq: it''s using the plotmNB function
|
---|
| 10 | %
|
---|
| 11 | %================================================================
|
---|
| 12 |
|
---|
| 13 | % Guillaume MAZE - LPO/LMD - December 2004
|
---|
| 14 | % gmaze@univ-brest.fr
|
---|
| 15 |
|
---|
| 16 | function [] = dispSVDnb(typg,CHP1,CHP2,N,expvar)
|
---|
| 17 |
|
---|
| 18 | global ffile1 ffile2
|
---|
| 19 | % --------------------------------------------------
|
---|
| 20 | % AFFICHAGE
|
---|
| 21 | % --------------------------------------------------
|
---|
| 22 | if (typg~=1)&(typg~=7)&(typg~=3)&(typg~=6)
|
---|
| 23 | CBAROR='horiz';
|
---|
| 24 | iw=N;jw=2;
|
---|
| 25 | else
|
---|
| 26 | CBAROR='vert';
|
---|
| 27 | iw=N;jw=2;
|
---|
| 28 | end
|
---|
| 29 |
|
---|
| 30 |
|
---|
| 31 | f1=figur;clf; hold on
|
---|
| 32 | for iN=1:N
|
---|
| 33 | subplot(iw,jw,2*iN-1)
|
---|
| 34 | C = squeeze(real(CHP1(iN,:,:)));
|
---|
| 35 | %[cs,h]=plotmNB(C,typg); whos cs h
|
---|
| 36 | [cs,h]=plotmNB(C,typg);
|
---|
| 37 | %clabel(cs,h,'labelspacing',100);
|
---|
| 38 | titre1=strcat('SVD',num2str(iN),':',ffile1);
|
---|
| 39 | titre2=strcat(' (',num2str(expvar(iN)),'%)');
|
---|
| 40 | title(strcat(titre1,titre2));
|
---|
| 41 |
|
---|
| 42 | subplot(iw,jw,2*iN)
|
---|
| 43 | C = squeeze(real(CHP2(iN,:,:)));
|
---|
| 44 | [cs,h]=plotmNB(C,typg);
|
---|
| 45 | %clabel(cs,h,'labelspacing',100);
|
---|
| 46 | titre1=strcat('SVD',num2str(iN),':',ffile2);
|
---|
| 47 | titre2=strcat(' (',num2str(expvar(iN)),'%)');
|
---|
| 48 | title(strcat(titre1,titre2));
|
---|
| 49 |
|
---|
| 50 | end %for iN
|
---|
| 51 |
|
---|
| 52 |
|
---|
| 53 |
|
---|
| 54 |
|
---|
| 55 | %titre1=strcat(num2str(lat(il1)),'/',num2str(lat(il2)));
|
---|
| 56 | titre2=strcat(' Total explained variance:',num2str(sum(expvar)),'%');
|
---|
| 57 | %titre3=strcat('(Analysed files: ',ffile1,' vs ',ffile2,')');
|
---|
| 58 | %titre=strvcat(titre2,titre3);
|
---|
| 59 | suptitle(titre2);
|
---|
| 60 |
|
---|
| 61 | set(f1,'Position',[378 39 313 647]); % Laptop screen
|
---|
| 62 | %set(f1,'Position',[369 55 316 899]); % Desktop screen
|
---|
| 63 | %set(f1,'Name',strcat(ffile,'<>',titre1,'<> REAL PART'));
|
---|