Intégration runge kutta numérique

[fusionfroide]

Bonjour,

Je fais face à un problème d'intégration numérique si vous pouvez m'aider s'il vous plait.

J'un tableau de 2 colonnes : TIME / SIGNAL. Se sont des valeurs relevées sur un oscilloscope. Ces points décrivent une courbe que j'aimerais intégrer.

J'utilise pour l'intégration la méthode de Runge Kutta d'ordre 4. Mais le code me renvoie des erreurs qui me font douter sur l'utilité de cette méthode dans mon cas.

Je me demande si la méthode Runge Kutta puisse marcher sur des fonctions non continues ... Ici, la fonction c'est le vecteur "signal" qui prend des valeurs sur un nombre de point limité et si jamais l'indice de signal n'est pas entier le code beug. Dans mon cas, l'indice à peu de chance qu'il soit entier à cause des valeurs prisent par t(i) et h comme on peut le voir dans le code suivant :

 Cliquez pour afficher

N=50;
t=0:2.0000e-08:9.8000e-07;			
h = (t(50)-t(1))/N;  % the step size
a = t(1);

for i = 1:1:N
   k1 = signal(t(i));
   k2 = signal(t(i)+h/2);
   k3 = signal(t(i)+h/2);
   k4 = signal(t(i)+h);
   y(i+1) = y(i) + (k1 + 2*k2 + 2*k3 + k4)/6;
   t(i+1) = a + i*h;
end

Merci pour l'aide.
-----

[gg0]

Bonjour.

Tu sembles avoir identifié le problème. Il ne te reste plus qu'à changer de méthode, puisque ce que tu fais ne fonctionne pas.
Ton code utilise une méthode de Simpson qui nécessite d'avoir les valeurs intermédiaires, donc de travailler sur 4 valeurs successives. Du vecteur des valeurs, ... ce que tu ne fais pas. A part t(1), aucun des t(i) n'est un indice (entier) !!!

Tu pourrais commencer par une simple intégration par la méthode des rectangle centrés. Ou des trapèzes.

Cordialement.

[Tryss]

Ici, ça ne sert à rien d'avoir une méthode d'intégration d'ordre élevée, puisque tu as déjà une perte lors de l'échantillonnage.

Je ferrai tout simplement la méthode du point milieu sur chaque morceau :