programmer un waterfall

Bonjour

Je récupère les traces d'un analyseur de spectre et là, j'essaie faire un waterfall
J'aimerais que ce waterfall soit cyclique.
Pour le moment, j'ai réussi à en faire un waterfall avec un nombre de points en Z fixe par exemple 100 ou 500 en modifiant simplement la valeur finale de la boucle « for » qui gère l’acquisition.
X= ne nombre de point d’acquisition (751 pts) il sera transformé plus tard pour afficher la fréquence
Y= l’amplitude du signal
Z= le nombre d’itérations, par exemple 100

J'aimerais, qu'arriver à la valeur finale de Z (par exemple 100) la valeur continue à s'incrémenter, mais que la longueur du graphique n'augmente pas.
Un peu comme si le graphique avançait avec un nombre de points en Z toujours constant.

Je suppose d'arriver à la valeur finale de 100 par exemple, on insère la nouvelle valeur, on supprime la première valeur, et on décale tout le tableau d’une ligne.
Ça me semble assez lourd…
Est-ce que vous avez une idée du principe que je peux utiliser pour réaliser cette opération efficacement avec Matlab ?
Merci

l'extrait utile du code :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

%% Acquisition des traces actives for i=1:500 acquisition.i=i; disp(i) if parametres.canalActif(1)==1 write(ssa,'TRACe:DATA? 1'); trace= readline(ssa); data=textscan(trace,'%f ','Delimiter',',' ); %acquisition ASCII à lire comme double (%f) le séparateur entre 2 valeurs et la , acquisition.CH1(i,:)=cell2mat(data); parametres.n=size(acquisition.CH1(i),1); %détermine la taille des données normalement 751 end pause(0.2) end %% axe fréquences parametres.fstep=(parametres.ffin-parametres.fdeb)/parametres.n; acquisition.axef=( parametres.fdeb:parametres.fstep:parametres.ffin); acquisition.axef=acquisition.axef(1:(parametres.n)); %% Affiche la figure parametres.screen_size = get(0,'ScreenSize'); parametres.largeur = parametres.screen_size(3); parametres.hauteur = parametres.screen_size(4); figure('Name','Analyseur SIGLENT','Position',[100, (parametres.hauteur-(600+77+5)), 1024, 600]); %% tracé de 1 à 4 traces   waterfall(acquisition.CH1)

---

Pièce jointe 631334

je pense avoir trouvé

Code:

---

ch = circshift(ch,1)

---

:D

Vu le code de la fonction circshift, il serait plus efficace d'utiliser simplement les indices.

Par exemple :

Code:

---

ch = ch([2:end 1]);

---

Salut ;)

Ça ne fonctionne pas.
ch a une structure de 100x751 après passage dans ch = ch([2:end 1]); il se transforme en 1x75100

La commande ch = circshift(ch,1) avec tic toc donne une durée variant de 500µs à 1.5ms, ça me va très bien.

Je n'ai jamais trouvé de bonnes solutions pour arrêter un programme avec Matlab:
CTRL C, c'est brutal et termine mal le programme.
J'ai tenté un truc comme cela, ça fonctionne, mais il faut insister sur le bouton-STOP pour que ça le prenne en compte:?

Une idée pour améliorer la chose ?

Merci

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29

%% Acquisition de la trace A %:SWEep:TIME? parametres.n=751; acquisition.CH1 = zeros(300,751); global st st=0; f = figure; c = uicontrol; c.String = 'STOP'; c.Callback = @plotButtonPushed; colormap jet   while (st==0) write(ssa,'TRACe:DATA? 1'); trace= readline(ssa); data=textscan(trace,'%f ','Delimiter',',' ); acquisition.CH1(1,:)=cell2mat(data); meshc(acquisition.CH1); acquisition.CH1 = circshift(acquisition.CH1,1); pause(0.05) end       function plotButtonPushed(src,event) global st st=1; end

---

Pour le moment j'ai laissé tomber le waterfall en continu.

Je me concentre sur l'acquisition avec un nombre fixe d'acquisitions en synchronisant le balayage (sweep) de l'analyseur.

J'ai essayé de recycler le traitement des axes et "datatips" de la fonction plot à la fonction mesh mais ca ne fonctionne pas bien.

Je joins à ce message un workspace avec les données pour tester ;)

la partie du programme qui affiche le mesh :
Merci par avance pour le coup de main

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

%% tracé de la traces   ax=meshc(acquisition.CH1); %affichage du tracé colormap jet colorbar hold on %ax.XData=ax.XData/parametres.ech; %echelle en fonction du parametre echelle % dtt = ax.DataTipTemplate; %dtt.DataTipRows(1).Label = 'Freq ='; %dtt.DataTipRows(1).Format = parametres.unite;%https://fr.mathworks.com/help/matlab/ref/matlab.graphics.datatip.datatiptextrow.html %dtt.DataTipRows(2).Label = 'Amp ='; %dtt.DataTipRows(2).Format = parametres.unitedatatip;   xlabel(parametres.axet); zlabel(parametres.uniteY); ylabel('Durée'); axis([parametres.fdeb/parametres.ech parametres.ffin/parametres.ech -inf inf]) grid on %grid minor set(gca,'xtick',linspace(parametres.fdeb/parametres.ech,parametres.ffin/parametres.ech,11)) title('Analyse spectrale') hold off if parametres.plottools ==1 plottools end

---

J'ai bien avancé,

je récupère les données de l'analyseur dans acquisition.z
Durant l'acquisition en boucle, je récupère la durée que je mets dans :acquisition.t
Pour l'échelle des fréquences, j'ai bricolé un truc que je qualifierais de "pas bien" je duplique l'échelle des fréquences autant de fois qu'il y a d'acquisitions dans :acquisition.x

Je n'arrive pas à personnaliser le "dataTips" qui est X, Y, Z et qui devrait être f = ... HZ, t=...s, Ampl=...dB.

Avec mon tableau acquisition.x, j'ai un warning :
Warning: Number of values in 'XData' and 'YData' must equal the number of columns and rows in
'ZData', respectively.
> In matlab.graphics.chart.primitive.Contour/updateFill
In matlab.graphics.chart.primitive.Contour/doUpdate

le code :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

%% Affiche la figure parametres.screen_size = get(0,'ScreenSize'); parametres.largeur = parametres.screen_size(3); parametres.hauteur = parametres.screen_size(4); figure('Name','Analyseur SIGLENT','Position',[100, (parametres.hauteur-(600+77+5)), 1024, 600]); %% tracé en 3D ax=meshc(acquisition.x,acquisition.t,acquisition.z); %affichage du tracé colormap jet % des couleurs bien visibles colorbar % ajout de la barre de couleurs % hold on %% ne fonctionne pas avec ax=meshc %ax.XData=ax.XData/parametres.ech; %echelle en fonction du parametre echelle % dtt = ax.DataTipTemplate; %dtt.DataTipRows(1).Label = 'Freq ='; %dtt.DataTipRows(1).Format = parametres.unite;%https://fr.mathworks.com/help/matlab/ref/matlab.graphics.datatip.datatiptextrow.html %dtt.DataTipRows(2).Label = 'Amp ='; %dtt.DataTipRows(2).Format = parametres.unitedatatip; %% fonctionne xlabel(parametres.axet); zlabel(parametres.uniteY); %en 2D Z est Y ylabel('Durée en s'); grid on %grid minor set(gca,'xtick',linspace(parametres.fdeb/parametres.ech,parametres.ffin/parametres.ech,11)) % change l'interval des graduations title('Analyse spectrale')

---

le Workspace pour essai si une ame charitable veut bien jeter un coup d'oeil :)

j'ai trouvé pour "meshc"

Code:

---

1 2	[x,y] = meshgrid(acquisition.axef,acquisition.t); ax=meshc(x,y,acquisition.z); %affichage du tracé

---

la personalisation des "Data Tips" ne fonctionne toujours pas :)