Comment générer une fonction de deux variables sans utiliser de boucle?

Supposons que j'ai une fonction y(t,x) = exp(-t)*sin(x)

Dans Matlab, j'définir

t = [0: 0.5: 5];
x = [0: 0.1: 10*2*pi];
y = zeros(length(t), length(x)); % empty matrix init

Maintenant, comment puis-je définir la matrice y sans utiliser de boucle, de telle sorte que chaque élément y(i,j) contient la valeur de la fonction désirée y à (t(i), x(j))? Ci-dessous est la façon dont je l'ai fait à l'aide d'une boucle for.

for i = 1:length(t)
    y(i,:) =  exp(-t(i)) .* sin(x);
end
Pas sûr qu'on puisse - est de la boucle de prendre de temps à s'exécuter, ou pourquoi voulez-vous?
La boucle n'est pas un problème ici, dans cet exemple simple. Mais je veux apprendre la façon de le faire sans boucle pour vraiment de grands vecteurs!

OriginalL'auteur Aamir | 2009-12-25

  1. 12

    Vos vecteurs en entrée x est 1xN et t est 1xM, matrice de sortie y est MxN. Pour vectoriser le code x et t doivent avoir la même dimension que y. 

    [x_,t_] = meshgrid(x,t);
y_ =  exp(-t_) .* sin(x_);

    Votre exemple est un simple 2D cas. Fonction meshgrid() travaille également en 3D. Parfois, vous ne pouvez pas éviter de la boucle, dans de tels cas, lorsque votre boucle peut aller de 1:N ou 1:M, choisir le plus court. Une autre fonction que j'utilise pour préparer vecteur pour vectorisé équation (vecteur x de multiplication de matrice) est diag().

    voila! merci pour une telle solution élégante.

    OriginalL'auteur Mikhail

  2. 6

    il n'est pas nécessaire pour meshgrid; utilisez simplement:

    y = exp(-t(:)) * sin(x(:)');    %multiplies a column vector times a row vector.
    +1 pour une simple et élégant à la matrice de multiplication de la solution.

    OriginalL'auteur shabbychef