Discussion :

[nlbmoi]

Bonjour

Je dois tracer le spectre en amplitude d'une fonction périodique, je ne suis pas un très grand expert de Scilab, aussi je suis bloqué sur quelques points :
- je souhaiterai que l'abscisse soit de la forme n/T avec n un entier positif et T la période du signal qui est quelconque.
- comment faire apparaître l'amplitude de chaque abscisse sous la forme d'un baton (en gros un histogramme très fin)
Une idée de ce que ça doit rendre pour éclairer mon questionnement :


Merci beaucoup pour votre aide

[Jerome Briot]

Un début de solution :

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
30
31
n = 0:0.05:15
E = 2
A = 400*E./((2*n+1).*n)
 
figure
scf()
 
width = 3
for k = 0:15
    loc(1,k+1) = find(n==k);
    if k==0 then
        lab(1,k+1) = msprintf("%d", k)
    else
        lab(1,k+1) = msprintf("%d/T", k)
 
        x = loc(1,k+1)
        y = A(loc(1,k+1))
        height = y
        xfrect(x,y,width,height)
        gce().background = color("gray");
 
    end
end
 
plot2d(A)
 
db = gca().data_bounds 
db($) = 300;
gca().data_bounds = db
 
gca().x_ticks = tlist(["ticks" "locations", "labels"], loc, lab)

[nlbmoi]

Bonsoir

Merci pour cette réponse.
Quelques demandes d'explication : je ne comprends pas l'instruction de la boucle for pour k

Est-il possible de également de faire le même raisonnement avec les valeurs d'ordonnées en utilisant une constante quelconque (E ici) : j'ai regardé dans ce que tu proposes mais je n'arrive pas à savoir comment tu intègres le n/T (mais comme ça apparait dans la boucle for et que je ne la comprends pas, peut-être que la réponse à ma première question me donnera la seconde) !