Table des Matières
Publicité
Exemple 2 – Afin de produire le signal donné par le spectre, nous modifions le
programme GDATA pour inclure une valeur absolue, de telle sorte que le
programme s'écrive comme suit :
<<
m a b << '2^m' EVAL
n << '(b-a)/(n+1)' EVAL
Dx << 1 n FOR j
'a+(j-1)*Dx' EVAL f ABS NEXT n
ARRY >> >> >> >>
Enregistrez cette version du programme sous GSPEC (Generate SPECtrum).
Lancez le programme avec m = 6, a = 0, b = 100. En mode RPN, utilisez :
6#0#100@GSPEC!
Une fois fini, appuyez sur la commande ` , pour garder une copie
supplémentaire du spectre de la matrice. Convertissez ensuite ce vecteur ligne
en un vecteur colonne, puis enregistrez-le dans ΣDAT. En suivant la méthode de
l'histogramme, le spectre obtenu dans cet exemple est donné ci-dessous. Dans
le cas présent, l'échelle horizontale est de 0 à 64, alors que l'échelle verticale
est de –1 à 10 :
Pour obtenir le signal du spectre ci-dessus, utilisez la fonction IFFT. Puisque
nous avons gardé une copie du spectre dans la pile (un vecteur ligne), vous
avez seulement besoin d'utiliser la fonction IFFT du menu MTH/FFT ou le
catalogue de commande, ‚N. Autrement, vous pourriez aussi saisir le
nom de la fonction (par exemple!: ~~ifft`. Le signal est affiché
dans une matrice (vecteur ligne) avec des nombres complexes. Nous nous
interessons seulement aux parties réelles de ces nombres. Pour obtenir leurs
parties réelles, utilisez la fonction RE du menu CMPLX (voir Chapitre 4), c'est-à-
dire entrez ~~re`. Ce qui est produit est un autre vecteur ligne.
Convertissez-le en un vecteur colonne, puis enregistrez le dans SDAT et tracez
un histogramme pour afficher le signal. Le signal de cet exemple est indiqué ci-
dessous, en utilisant une échelle horizontale de 0 à 64 et une échelle verticale
de –1 à 1 :
Page. 16-55
Publicité
Table des Matières
