Mercurial > repos > jcb-mpl > features_extract
annotate src/mz_6bis_durbin_watson.sci @ 0:c0a2e5e78725 draft
Uploaded
| author | jcb-mpl |
|---|---|
| date | Tue, 27 Apr 2021 14:48:05 +0000 |
| parents | |
| children |
| rev | line source |
|---|---|
| 0 | 1 function x_dw=mz_6bis_durbin_watson(x_tr_mz) |
| 2 | |
| 3 // pas de calcul matriciel car pas assez de memoire | |
| 4 // x_dw doit varier entre 0 et 2 | |
| 5 // 0=ce qu on recherche | |
| 6 | |
| 7 [n,q]=size(x_tr_mz.d); | |
| 8 | |
| 9 // calcul de la norme des differences et des vect-colonne de x_tr_mz | |
| 10 x_diff=x_tr_mz.d(2:n,:)-x_tr_mz.d(1:n-1,:); | |
| 11 rmsce=zeros(q,1); | |
| 12 rmsce2=zeros(q,1); | |
| 13 for i=1:q; | |
| 14 rmsce(i)=sqrt((x_diff(:,i)'*x_diff(:,i))/(n-1)); | |
| 15 rmsce2(i)=sqrt((x_tr_mz.d(:,i)'*x_tr_mz.d(:,i))/(n)); | |
| 16 end | |
| 17 | |
| 18 x_dw=rmsce./rmsce2; | |
| 19 | |
| 20 | |
| 21 endfunction |