Discussion :

[dido555555]

Salut tout le monde
je dois tranformer un algorithme en code matlab afin de calculer le coefficient de holder de chaque point d'une série temporelle,mais malheureusement je n'y arrive pas,
pourriez vous m'aider à résoudre mon problème?voici mon code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
x=-10:0.1:10; y=sin(x);% cest  un exemple  que j'ai pris au pif
N=1000;%idem
for i=0:N-1
       for k=-i:N-1-i %k doit être différent de zéro mais comment l'exprimer?
 
           R(i,k)=log(abs(y(i+k)-y(i)))/log((abs(k)/N));
            for s=1:N-1+i
                h(i,s)=min(R(i,s));
            end
       end
    for n=1:N-1-i 
 H(i)= sum(alpha^(n-1)*h(i,n))*(1-alpha)/(1-alpha^(N-1-i )) 
    end
end

Merci!!

[Caro-Line]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

k doit être différent de zéro mais comment l'exprimer?

=>

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

if k~=0

C'est la seule question ?

[dido555555]

Sinon c'etait pour savoir sil existait un code matlab pour calculer le coeffcient de holder pour un signal temporelle en utilisant la méthode des maxima des modules de la transformée d'ondelettes?
et je voudrais également savoir comment on peut définir sur matlab un intervalle ouvert: exemple faire la somme des (n+1)^2 pour n appartenant à l'intervalle [-2,0[ U ]0,2[
merci d'avance!

[Jean Dumoncel]

Salut,

pour ta première question, je ne sais pas (si elle n'existe pas, serait-elle compliqué à coder à ton avis?)

Pour la somme des (n+1)^2, une solution pourrait être :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
>> n=2;
>> sum(([-n:-1,1:n]+1).^2)

[dido555555]

salut,merci pour votre réponse.
En ce qui concerne le coefficient de Holder,j'ai trouvé un algorithme pour le calculer mais je l'ai trouvé dur à coder sur matlab même s'il ne fait que quelque lignes.
J'ai tout de même essayé de le coder mais sans grand résultat(voir premier message)
j'ai mis mis en pièce jointe l'algorithme,si quelqu'un arrive à trouver mes erreurs,ça m'aiderai beaucoup!
Voici mon code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
for i=0:N-1
       for k=-i:N-1-i 
 
           R(i,k)=log(abs(y(i+k)-y(i)))/log((abs(k)/N));
            for s=1:N-1+i
                h(i,s)=min(R(i,s));
            end
       end
    for n=1:N-1-i 
 H(i)= sum(alpha^(n-1)*h(i,n))*(1-alpha)/(1-alpha^(N-1-i )) 
    end
end