Last change
on this file since 10 was
10,
checked in by pinsard, 16 years ago
|
first commit with original work of Sebastien Gervois
|
File size:
683 bytes
|
Line | |
---|
1 | %%%%%%%%%%%%%%% FICHIER OLR
|
---|
2 | % ouverture / lecture du fichier
|
---|
3 | clear
|
---|
4 | close all
|
---|
5 |
|
---|
6 | fid=fopen('olrf10-30_ctl.80.dat','r');
|
---|
7 | siz=24*15*3538;
|
---|
8 | v=fread(fid,siz,'float');
|
---|
9 | fclose(fid);
|
---|
10 | olr=reshape(v,360,3538); clear v;
|
---|
11 |
|
---|
12 |
|
---|
13 | % calcul de la matrice de covariance
|
---|
14 | vaudou=zeros(360,360);
|
---|
15 | for a=1:360;
|
---|
16 | for b=1:360;
|
---|
17 | clear e d covaria;
|
---|
18 | e=olr(b,:);
|
---|
19 | d=olr(a,:);
|
---|
20 | covaria=cov(e,d);
|
---|
21 | vaudou(a,b)=covaria(2,1);
|
---|
22 | end;
|
---|
23 | end;
|
---|
24 |
|
---|
25 | % on diagonalise la matrice de covariance
|
---|
26 | [C,L]=eig(vaudou);
|
---|
27 |
|
---|
28 | % spatialisation du premier mode
|
---|
29 |
|
---|
30 |
|
---|
31 | % serie temporelle de la premiere vP
|
---|
32 |
|
---|
33 | delete('eof_1030_ircaam_ctl.mat');
|
---|
34 | save('eof_1030_ircaam_ctl.mat','C','L','olr')
|
---|
Note: See
TracBrowser
for help on using the repository browser.