Manipulation d'une fonction à base d'opérations vectorielles pour une intégration

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:

---

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:

---

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:

---

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:

---

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 !

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

salut,

tu peux aussi explicité le produit scalaire ;)

Code:

---

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

---

++

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

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:

---

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

---

++