animer parcelle / trajectoire dans matlab / octave

Je suis en train d'animer cette spirale à l'aide de matlab /octave, je veux qu'il en spirale vers le haut ou vers le bas

t = 0:0.1:10*pi;
r = linspace (0, 1, numel (t));
z = linspace (0, 1, numel (t));
plot3 (r.*sin(t), r.*cos(t), z);

animer parcelle /trajectoire dans matlab /octave

J'ai essayé d'utiliser une boucle for pour l'animer, mais qui me donne juste une forme de cône voir le code et l'image ci-dessous

clear all, clc,clf,tic
t = 0:0.1:10*pi;
r = linspace (0, 1, numel (t));
z = linspace (0, 1, numel (t));

for ii=1:length(r)
    ii
    plot3 (r.*sin(t(ii)), r.*cos(t(ii)), z);
    hold on
    %pause (.00001)
end

Image animer parcelle /trajectoire dans matlab /octave

Découvrez cette réponse: stackoverflow.com/questions/16426378/...

OriginalL'auteur Rick T | 2013-05-15

  1. 1

    Le message suivant apparaît à travailler dans l'Octave 3.6.2

    t = 0:0.1:10*pi;
r = linspace (0, 1, numel (t));
z = linspace (0, 1, numel (t));

figure
axis([-1 1 -1 1 0 1])
hold on

for ii=1:length(r)
    plot3 (r(ii)*sin(t(ii)), r(ii)*cos(t(ii)), z(ii),'*');
    pause (.001)
end
    C'est gentil, mais il ne fournit pas une trajectoire, une séquence de points. Voir stackoverflow.com/a/17773313/181638 pour une solution qui anime une trajectoire.

    OriginalL'auteur am304

  2. 2

    Vous pouvez également utiliser la comet3() paquet, qui anime la trajectoire dans l'intrigue:

    delay = 0.001  % seconds
figure
comet3(r.*sin(t), r.*cos(t), z, delay);

    Cette anime en continu d'une trajectoire dont je préfère le plus discret de la séquence de"*".

    Le seul inconvénient est que la version de comet et comet3 livré avec Octave 3.6.4 sont lents, quel que soit le retard que vous utilisez. Mais cela peut être surmonté en utilisant l'astuce suivante courtoisie de andyras dans cette SORTE de question: 

    % plot the first point to get started
h = plot3(x(1),y(1),z(1),"b");
axis([min(x), max(x), min(y), max(y), min(z), max(z)]);

% refresh the plot in a loop through the rest of the data
for k = 1:length(z);
   set(h, 'XData', x(1:k));
   set(h, 'YData', y(1:k));
   set(h, 'ZData', z(1:k));
   pause (0.001); % delay in seconds
       % alternatively could provide a velocity function
       % pause(sqrt(vx(k)^2+vy(k)^2+vz(k)^2));  
endfor

    Note mineure: une fois que vous avez modifié la fonction, vous aurez besoin de force Octave de la recharger puisqu'il n'aura pas le faire par défaut. Vous pouvez redémarrer, ou mieux encore utiliser clear comet et clear comet3. Alors la prochaine fois que ces fonctions sont appelées, leurs définitions seront actualisées.

    OriginalL'auteur Assad Ebrahim

  3. 1

    Certainement pas la plus jolie, mais ce sont les premières modifications que vous devez apporter à votre code pour qu'il fasse quelque chose de proche de ce que vous voulez.

    t = 0:0.1:10*pi;
z = linspace (1, 0, numel (t));
for ii=1:length(t)
    plot3 (z(ii)*sin(t(ii)),z(ii)*cos(t(ii)), z(ii));
    hold on
    pause (.00001)
end
    Davin, Cela me donne juste une erreur", la ligne 0: Tous bords indéfinis ou hors de gamme, donc pas de complot."
    jamais utilisé octave avant, je suis exécution sur Matlab, on peut supposer que c'est une compatibilité/la version en question. Vous pouvez essayer d'ajouter de la clear ... ligne, ou le remplacement de plot3 avec scatter3, peut-être, qui l'aide. Si non, je suis frais d'idées.

    OriginalL'auteur davin