Discussion :

[Philou29200]

Bonjour,

Je suis débutant en MatLab et je dois calculer 2 intégrales doubles. Malheureusement quand je demande à Matlab de calculer la première il me renvoit les erreurs suivantes :

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
 
varIntegral = dblquad(MyEsperance,-inf,inf,-inf,inf)
 
??? Error using ==> mtimes
Inner matrix dimensions must agree.
 
Error in ==> quad at 62
y = f(x, varargin{:});
 
Error in ==> dblquad>innerintegral at 88
    Q(i) = feval(quadf, intfcn, xmin, xmax, tol, trace, y(i), varargin{:}); 
 
Error in ==> quad at 62
y = f(x, varargin{:});
 
Error in ==> dblquad at 64
Q = feval(quadf, @innerintegral, ymin, ymax, tol, trace, intfcn, ...

Et pour la deuxième intégrale double :

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
varIntegral = dblquad(MyEsperance2,-inf,inf,-inf,inf)
 
??? Error using ==> mpower
Matrix must be square.
 
Error in ==> quad at 62
y = f(x, varargin{:});
 
Error in ==> dblquad>innerintegral at 88
    Q(i) = feval(quadf, intfcn, xmin, xmax, tol, trace, y(i), varargin{:}); 
 
Error in ==> quad at 62
y = f(x, varargin{:});
 
Error in ==> dblquad at 64
Q = feval(quadf, @innerintegral, ymin, ymax, tol, trace, intfcn, ...

Je n'y connais rien du tout et j'ai besoin de calculer ces 2 intégrales c'est très important pour moi (c'est mon mémoire de fin d'études).

Voici mes deux intégrales :

Code :

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

1
2
3
MyEsperance = @(v,x)((A_0*exp(-phi*sigma_A/sigma_2*(x-theta_1*T_1-v)+x+v-theta_2*T_1)-P_0*exp(-sigma_P/sigma_2*x+v+theta_3*T_1))*heaviside(-(A_0*exp(-phi*sigma_A/sigma_2*(x-theta_1*T_1-v)+x+v-theta_2*T_1)-P_0*exp(-sigma_P/sigma_2*x+v+theta_3*T_1)-C))*1/sqrt(2*pi*T_1)*exp(-1/(2*T_1)*v^2)*1/sqrt(2*pi*T_1)*(exp(-(x-mu_B1*T_1)^2/(2*T_1))-exp(-(mu_B1*x-mu_B1^2*T_1/2-(abs(x+1/2*log(A_0/P_0))+1/2*log(A_0/P_0))^2)/(2*T_1))))
 
MyEsperance2 = @(v,x)(heaviside(-(A_0*exp(-phi*sigma_A/sigma_2*(x-theta_1*T_1-v)+x+v-theta_2*T_1)-P_0*exp(-sigma_P/sigma_2*x+v+theta_3*T_1)-C))*1/sqrt(2*pi*T_1)*exp(-1/(2*T_1)*v^2)*1/sqrt(2*pi*T_1)*(exp(-(x-mu_B1*T_1)^2/(2*T_1))-exp(-(mu_B1*x-mu_B1^2*T_1/2-(abs(x+1/2*log(A_0/P_0))+1/2*log(A_0/P_0))^2)/(2*T_1))))

Voici les données :

Code :

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

mu_A=0.1055; mu_P=0.1051; sigma_A=0.1295; sigma_P=0.1323; phi=0.9992; A_0 = 722.9; P_0 = 674; C = 40.5; T=1; mu = mu_A - mu_P + 1/2 * (sigma_P^2-sigma_A^2); sigma_1 = sigma_A * sqrt(1-phi^2); sigma_2 = -(phi*sigma_A - sigma_P); sigma = sqrt((sigma_1^2+sigma_2^2)/2); T_1 = sigma^2/2*T; mu_B1 = (mu_A - mu_P + 1/2 * (sigma_P^2-sigma_A^2))/sigma^2; theta_1=mu_B1; theta_2 = (-mu_A-mu_P+1/2*(sigma_1^2+sigma_2^2))/sigma^2; theta_3 = sigma_P/sigma_2*theta_1+(mu_P-1/2*sigma_P^2)*2/sigma^2;

Pourriez vous m'aider s'il vous plait. Je vous remercie par avance
Philippe

[salseropom]

Salut,
bienvenue sur ce forum.

Tes problèmes de dimensions me font penser à des vecteurs lignes et colonnes que tu cherches à additionner par exemple.

Pour ton erreur sur la 2e intégrale, ta matrice n'est pas carrée. Affiche là.

[Philou29200]

Salut,

Merci beaucoup pour tes conseils, je vais y réfléchir. C'est la premiere fois que j'utilise Matlab, et je ne comprends pas vraiment pourquoi mon calcul d'intégrales amène à considérer des matrices...

Je vais aussi essayer de trouver comment afficher ces matrices...

Bonne soirée