Discussion :

[GammaZ]

Bonjour,

Je suis grand débutant en Matlab. J'ai créé une fonction à base d'opérateurs vectoriels (comme les produits scalaires et vectoriels). Cette fonction prend pour arguments des scalaires et renvoie un scalaire.

Admettons que ce soit cette fonction (c'est en fait un enchaînement de plusieurs fonctions, dont le détail n'est pas intéressant car ce n'est pas de là que ça vient) :

Code :

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

1
2
3
function r = TestScalaire(x,y)
r = dot([x,y],[2,3])
end

Cette fonction marche correctement (heureusement, je pense !).

J'aimerais intégrer cette fonction en prenant pour paramètre d'intégration les coordonnées d'un vecteur.

Code :

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

1
2
 
dblquad(@(x,y)TestScalaire(x,y), 0, 1, 0, 1);

Là, j'ai des erreurs provenant de ce produit scalaire, alors que rien ne devrait gêner son exécution. Mise à part une erreur de taille d'éléments puisque Matlab a l'air d'interpréter [x,y] comme une matrice de valeurs pour l'intégration de TestScalaire.

Code :

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

1
2
3
 
Error using ==> dot at xx
A and B must be the same size

La même chose se produit lorsque j'utilise dot() directement dans la fenêtre de commande :

Code :

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

dblquad(@(x,y)dot([x,y],[2,3]), 0, 1, 0, 1);

Voyez-vous une solution pour venir à bout de ce calcul sans développer tous les calculs vectoriels effectués dans mes fonctions ?

Merci d'avance !

[GammaZ]

Merci Dut pour ma mise en forme, j'utiliserai les balises de code pour mes prochains messages.

[mr_samurai]

salut,

tu peux aussi explicité le produit scalaire

Code :

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

dblquad(@(x,y)(2*x + 3*y), 0, 1, 0, 1)

++

[GammaZ]

Merci mr_samurai ; c'est ce que je proposais moi-même à la fin de mon post ; j'ai fini par réécrire le produit scalaire, chaque composante du produit scalaire ainsi que la norme pour ne manipuler que des scalaires au final.
Mais c'est dommage de se passer des fonctions de Matlab pour avoir le même résultat (étant donné que mes fonctions utilisaient très souvent des opérateurs vectoriels).

[mr_samurai]

C'est parce que ce n'est pas tout à faire vectoriel que ça crash ...
Tu peux toujours ajouter toi même les opérateurs vectoriels (ce qui se fait dans la plupart des implémentations) :

Code :

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

1
2
3
4
% Produit Scalaire : X , Y des vecteurs colonne
myDot = @(X,Y) Y'*X;
 
myDot([1;2;3;4],[5;6;7;8])

++