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'));
|
---|