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))

[shaiHulud]

Bonjour,

La fonction randn génère des variables aléatoires x(1)...x(n), tandis que trapz attends une fonction en argument. Il faut donc créer la fonction n|-> x(n), par exemple

Code :

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

1
2
3
4
5
6
7
8
9
10
 
function f= timeserie2fun(x)
function xfloorn= wrap(n)
% vaut x(n) si n<=t<n+1 ; 0 si t<1 et x(end) si t> end
if x<1; xfloorn=0
elseif x>length(x); xfloorn= x(end);
else; xfloorn= x(n);
end
f= @wrap
end

utilisation

Code :

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

1
2
3
 
x= randn(10,1);
f= timeserie2fun(x)

Cela étant, il doit pas être dur de faire le calcul à la main en divisant l'intégrale sur les intervalles [n,n+1[

NB: je n'arrive pas a lire ta PJ, donc je n'ai peut être pas tout compris

[shaiHulud]

Je vois enfin ta PJ ! Es-tu sur que f(t) est la dernière valeur d'une séquence de nombre aléatoire ? N'est ce pas plutôt la densité d'une variable aléatoire ?