Discussion :

[samo25]

Bonsoir a tous,
Je veux savoir comment résolu ce problème parce que je n'arrive pas à tracer la réponse de Q en fonction de P.

Code :

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

1
2
3
4
5
6
 
for P=-50:50:150
 
   Q(P)=dblquad(@(X,Y)bfun(X,Y,P),0,1,0,1);
end 
plot(Q);grid on

le message d'erreur est :

??? Subscript indices must either be real positive integers or logicals.

Error in ==> blas at 37
   Q(P)=dblquad(@(X,Y)bfun(X,Y,P),0,1,0,1);

Merci pour toute réponse.
Meilleures salutations.

[Jean Dumoncel]

Bonjour,

voir la :
Subscript indices must either be real positive integers or logicals.

[samo25]

salut magelan,
j'ai vue ce lien de FAQ Matlab, mais est ce qu'il y a une autre méthode pour résolu ce problème car je sais que variable P doit commencer par P=1 , j'ai schéma qui commencer par P=-50 comment il arrive a tracé !!!!!!!!

[samo25]

J'ai essayé a reprogrammer :

Code :

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

1
2
3
4
5
6
7
 
P=-50:50:150;
Q=[];
 
Q(P)=dblquad(@(X,Y)bfun(X,Y,P),0,1,0,1);
 
plot(Q);grid on

et la fonction bfun :

Code :

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

1
2
3
 
function z=bfun(X,Y,P)
z=1-P.*(1-(x/a).^2).*(1-(y/b).^2);

le message d'erreur est :

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> bfun at 2
z=1-P.*(1-(x/a).^2).*(1-(y/b).^2);
Error in ==> @(X,Y)bfun(X,Y,P)

Error in ==> blas at 4
   Q(P)=dblquad(@(X,Y)bfun(X,Y,P),0,1,0,1);

Meilleures salutations.

[Jean Dumoncel]

C'est parce que tu multiplie terme à terme des vecteurs de tailles différentes...

Pour ton premier code, il suffisait juste d'ajouter 51 à P pour que -50 devienne 1, -49 devienne 2, etc :

Code :

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

1
2
3
4
5
for P=-50:50:150
  
   Q(P+51)=dblquad(@(X,Y)bfun(X,Y,P),0,1,0,1);
end 
plot(Q);grid on

[Invité]

Bonjour,

Faisant par pas de 50, j'aurais plus vu:

Code :

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

1
2
3
4
5
P=-50:50:150;
for i = 1:length(P)  
   Q(i)=dblquad(@(X,Y)bfun(X,Y,P(i)),0,1,0,1);
end 
plot(Q);grid on

Ce qui posait problème était Q(P) avec

Code :

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

P = -50     0    50   100   150

Ne faudrait-il pas des majuscules pour être en accord avec les paramètres:

Code :

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

1
2
3
function z=bfun(X,Y,P)
z=1-P.*(1-(x/a).^2).*(1-(y/b).^2);

?

[samo25]

Je suis d'accord avec vous concernant le X et Y majuscule, ils doivent être le même.
J'ai utilisé le code suivant :

Code :

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

1
2
3
4
5
6
 
P=-50:50:150;
for i = 1:length(P)  
   Q(i)=dblquad(@(X,Y)bfun(X,Y,P(i)),0,1,0,1);
end 
plot(Q);grid on

Autre chose je veux savoir comment changer l'axe de X de i par les valeurs de P parce que dans le graphe l'axe de X est 1 à 4 c'est à dire le nombre i des valeurs de P, et merci (voir le graphe).

[Invité]

Bah simplement

Code :

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

plot(P,Q)

Pense à regarder la documentation pour des questions de ce genre.

[samo25]

Tellement j'ai perturbé, les bases de plot sont disparu 'plot(x,y)'. merci