whos

              Name        Size       Elements     Bytes    Density   Complex

                 a       1 by 6             6        48       Full      No 
             agsh1    3284 by 45       147780   1182240       Full      No 
                 b       1 by 6             6        48       Full      No 
             ensh1    3284 by 1          3284     26272       Full      No 
             hesh1    3284 by 10        32840    262720       Full      No 
                 j       1 by 1             1         8       Full      No 
              text       1 by 38           38       304       Full      No 
             vesh1    3284 by 45       147780   1182240       Full      No 
                vi    3284 by 44       144496   2311936       Full      Yes
               vif    3284 by 44       144496   2311936       Full      Yes
             vnsh1    3284 by 45       147780   1182240       Full      No 

Grand total is 768507 elements using 8459992 bytes

clear agsh1 ansh1 hesh1 vesh1 vnsh1

whos

              Name        Size       Elements     Bytes    Density   Complex

                 a       1 by 6             6        48       Full      No 
                 b       1 by 6             6        48       Full      No 
             ensh1    3284 by 1          3284     26272       Full      No 
                 j       1 by 1             1         8       Full      No 
              text       1 by 38           38       304       Full      No 
                vi    3284 by 44       144496   2311936       Full      Yes
               vif    3284 by 44       144496   2311936       Full      Yes

Grand total is 292327 elements using 4650552 bytes

%filtrage a frequence inertielle effective

[b,a]=fir1(600,[1/18 1/17]/15,kaiser(601,0));

[h,w]=freqz(b,a,2048,30);

subplot(211);plot(w,abs(h));axis([0 0.1 0 1.1])

max(abs(h))


ans =

    1.0039

vifin=filtfilt(b,a,vif(:,1))    %regardons d'abord le bin 1


 
vifin=filtfilt(b,a,vif(:,1));    %regardons d'abord le bin 1

subplot(212);plot(vifin)

subplot(212);plot(abs(vifin))

title('module de vif filtre frequence inertielle')

subplot(211);title('module de reponse en frequence du filtre fir')

xlabel('frequence en hr-1')

diary off

eiwt=exp((sqrt(-1)*2*pi*(1:3284)/494);

??? 2*pi*(1:3284)/494);
                      |
Improper function reference. A "," or ")" is expected.

eiwt=exp((sqrt(-1))*2*pi*(1:3284)/494);

eiwtmoins=exp((sqrt(-1))*(-2)*pi*(1:3284)/494);

moymob=ones(1,494)/494;

vifinplus=vifin.*eiwt';

vifinmoins=vifin.*eiwtmoins';

aplus=filtfilt(moymob,1,vifinplus);

amoins=filtfilt(moymob,1,vifinmoins);

figure(2)

subplot(211);plot(abs(aplus))

subplot(212);plot(abs(amoins))

 title('module, mm/s, composante circulaire anticyclo feffectif')

subplot(211);

 title('module, mm/s, composante circulaire cyclo feffectif')

figure(3)

subplot(212);plot(angle(amoins))

subplot(211);plot(angle(aplus))

%ici ce sont les graphes de la phase de vifin*eiwt

title('phase vifinmoins (bas)et vifinplus(haut)')



%normalement donc je peux reconstruire l'amplitude,aplus et amoins

who


Your variables are:

a              eiwt           moymob         vifinmoins     
amoins         eiwtmoins      text           vifinplus      
ans            ensh1          vi             w              
aplus          h              vif            
b              j              vifin          

whos

              Name        Size       Elements     Bytes    Density   Complex

                 a       1 by 1             1         8       Full      No 
            amoins    3284 by 1          3284     52544       Full      Yes
               ans       1 by 1             1         8       Full      No 
             aplus    3284 by 1          3284     52544       Full      Yes
                 b       1 by 601         601      4808       Full      No 
              eiwt       1 by 3284       3284     52544       Full      Yes
         eiwtmoins       1 by 3284       3284     52544       Full      Yes
             ensh1    3284 by 1          3284     26272       Full      No 
                 h    2048 by 1          2048     32768       Full      Yes
                 j       1 by 1             1         8       Full      No 
            moymob       1 by 494         494      3952       Full      No 
              text       1 by 38           38       304       Full      No 
                vi    3284 by 44       144496   2311936       Full      Yes
               vif    3284 by 44       144496   2311936       Full      Yes
             vifin    3284 by 1          3284     52544       Full      Yes
        vifinmoins    3284 by 1          3284     52544       Full      Yes
         vifinplus    3284 by 1          3284     52544       Full      Yes
                 w    2048 by 1          2048     16384       Full      No 

Grand total is 320496 elements using 5076192 bytes

vinertp=aplus.*eiwt';

vinertm=amoins.*eiwtmoins';

plot(angle(vinertm))

plot(angle(vinertp))

vifinplus=vifin.*eiwt.';

vifinmoins=vifin.*eiwtmoins.';

amoins=filtfilt(moymob,1,vifinmoins);

aplus=filtfilt(moymob,1,vifinplus);

subplot(212);plot(abs(amoins))

 title('module, mm/s, composante circulaire anticyclo feffectif')

subplot(211);plot(abs(aplus))

 title('module, mm/s, composante circulaire cyclo feffectif')

vinertp=amoins.*eiwt';

vinertm=aplus.*eiwtmoins.';

vinertp=amoins.*eiwt.';

%ces deux derniers ordes sont les bons car plus referent en fait a moins

plot(angle(vinertm))

plot(angle(vinertp))

%il faudrait peutre prendre abs(plus)*eiwtmoins' pour onde inertielle



plot(abs(vif))

(:,1)

??? (:
     |
Missing variable or function.

plot(abs(vif(:,1)))

hold on

plot(abs(vif(:,1)-vinertm)),'w')

??? :,1)-vinertm)),'w')
                      |
Missing operator, comma, or semi-colon.

plot(abs(vif(:,1)-vinertm),'w')



plot(abs(vinertm)),'r')

??? (abs(vinertm)),'r')
                      |
Missing operator, comma, or semi-colon.

plot(abs(vinertm),'r')

subplot(212)

plot(angle(vinertm),'r')

hold on

plot(angle(vif(:,1)-vinertm),'w')

title('phase des courants inertiels et des residus')

subplot(211)

title('modules du courant filtre 2 heures,du courant inertie et du residu')

subplot(212)

xlabel('temps en pas de deux minutes dynaproc shet1')

figure

plot(real(vif(:,1))

??? plot(real(vif(:,1))
                       |
Improper function reference. A "," or ")" is expected.

plot(real(vif(:,1)))

subplot(211)

plot(real(vif(:,1)))

hold on

plot(real(vinertm),'w')

title('composante Est de courant filtre2 heure et partie inertielle')

plot(real(vif(:,1)-real(vinertm),'g')

??? -real(vinertm),'g')
                       |
Improper function reference. A "," or ")" is expected.

plot(real(vif(:,1))-real(vinertm),'g')

subplot(212)

plot(imag(vif(:,1)))

plot(image(vif(:,1))-imag(vinertm),'g')

plot(image(vinertm),'w')

plot(imag(vif(:,1)))

hold on

plot(imag(vinertm),'w')

plot(imag(vif(:,1))-imag(vinertm),'g')

axis([0 3500 -200 200])



title('composante Nord de courant filtre2 heure et partie inertielle')

 save etuinert

diary off

