Discussion :

[Hellios]

Bonjour à tous !
Voila, je viens vous voir car pour un projet je dois réaliser une intégration (voir PJ).

Sachant que f(t) est un bruit blanc gaussien normalisé de moyenne nulle et d'intensité unité...

Tout d'abord savez-vous générer une telle fonction sous Matlab ?
(Je n'ai pas réussi à le faire et pour l'instant j'utilise la commande randn).

Voila mon code:

(lambda est un réel)

Code :

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

1
2
3
4
5
6
7
8
9
 
 
t=linspace(0,100,101);
 
for i = 1:t_max
u=0:1:t(i);
a=exp(lambda.*u).*randn(1,t(i)+1);
integral(i)=trapz(u,a);
end

Comme vous vous en doutez il ne marche pas et Matlab m'affiche:

??? Error using ==> permute
ORDER must have at least N elements for an N-D array

Error in ==> trapz at 44
  y = permute(y,perm);

Error in ==> Script at 65
integral(i)=trapz(u,a);

Si vous avez une quelconque idée pour réaliser cette intégration je suis fortement preneur !
Merci d'avance et bonne soirée.


(PS:Je possède la version 7.12.0.635 (R2011a))