|
0
|
1 function x_tr_mz3=mz_7_dwcleaning(x_tr_mz2,durbwats)
|
|
|
2
|
|
|
3
|
|
|
4 // cette fonction nettoie la matrice des EICs avec DW trop élevé
|
|
|
5
|
|
|
6 // x_tr_mz2: matrice EIC issue de la fonction mz_6_rescan_mz2
|
|
|
7
|
|
|
8 // durbwats: ex: dw=1.2
|
|
|
9
|
|
|
10
|
|
|
11
|
|
|
12 // nettoyage avec Durbin-Watson: très efficace
|
|
|
13 res_dw=mz_6bis_durbin_watson(x_tr_mz2);
|
|
|
14
|
|
|
15
|
|
|
16 tri=find(res_dw<durbwats);
|
|
|
17
|
|
|
18 if tri~=[] then
|
|
|
19 x_tr_mz3=x_tr_mz2(:,tri); // 1000*= rapide de sélectionner les bonnes colonnes qu'enlever les mauvaises
|
|
|
20 else
|
|
|
21 x_tr_mz3=x_tr_mz2;
|
|
|
22 end
|
|
|
23
|
|
|
24 // note: nettoyage avec eic_signal testé avec valeur=10000 -> inutile!
|
|
|
25
|
|
|
26
|
|
|
27
|
|
|
28 endfunction
|
