Discussion :

[cubitus_54]

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 :

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

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)

[cubitus_54]

je pense avoir trouvé

Code :

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

 ch = circshift(ch,1)

[Jerome Briot]

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

Par exemple :

Code :

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

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

[cubitus_54]

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 :

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

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

[cubitus_54]

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 :

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

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

[cubitus_54]

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 :

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

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