Bonjour,
Pourriez vous m'aider à résoudre mon problème? J'ai une fonction qui me retourne des nombres complexes alors qu'elle censée me retourner des réels.

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
function [ lvce ] = mLogvceN( param ) 
 
global wind h
 
c  = param( 1 ) ;                   
 
a0 = param( 2 ) ;                  
 
a1 = param( 3 ) ;                  
 
b1 = param( 4 ) ;                  
 
n = length( wind ) ;               
 
variance = a0 / ( 1 - a1 - b1 ) ; 
 
h = ones( n , 1 ) * variance ;
 
for i = 2 : n                       
 
   h( i ) = a0 + a1 * ( wind( i - 1 ) - c ) ^ 2 + b1 * h( i - 1 ) ;
 
end  
 
lvce = - 0.5 * ( sum( log( h ) ) + n * log( 2 * pi ) + sum( (( wind - c ) .^ 2) ./ h  ) ) ;
Merci d'avance.