Discussion :

[mariellemm]

Bonjour à tous!

Je suis débutante sur matlab et je dois utiliser ce logiciel afin de récupérer la deuxieme colone d'un fichier texte et en faire la transformée de fourier rapide.
Puis je voudrais afficher à l'écran.

pour l'instant voici mon 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
 
clear all;
fiche=fopen('a47.txt','r') ;
donnee=fscanf(fiche,'%g %g',[2 inf]);
fclose(fid);
 
[s,t]=textread('a47.txt');
% calcul de fft
Y=fft(s,t);
 
% dessiner la fft
plot(f,2*abs(Y(1:NFFT/2+1))) 
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|'

Je vous avoue que je n'y comprend pas grand chose... J'ai simplement regarder différent code ressemblant à ceux que je voulais faire...
Quelqu'un pourrait-il m'expliquer comment récupérer ces données et les afficher?

Merci d'avance!

[Newenda]

pour l'instant voici mon 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
clear all;
fiche=fopen('a47.txt','r') ;
donnee=fscanf(fiche,'%g %g',[2 inf]); <---- un seul '%g' suffit
fclose(fid);

[s,t]=textread('a47.txt'); <---- pourquoi lire deux fois?
% calcul de fft
Y=fft(s,t);

% dessiner la fft
plot(f,2*abs(Y(1:NFFT/2+1)))  <------ NFFT doit être défini ainsi que f
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|'

Je te propose un truc du genre :

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
close all ;
clear all ;

fiche=fopen('a47.txt','r') ;
donnee=fscanf(fiche,'%g',[2 inf]); 
s = donnee(:,1);
t = donnee(:,2);
fclose(fid);

N = lenght(t) ;
NFFT = 2^nextpow2(N); 

Y = fft(t,NFFT)/N;

spectre_t = 2*abs(Y(1:NFFT/2+1)) ;
freq = freq_ech/2*linspace(0,1,NFFT/2+1); freq_ech = fréquence échantillonage

puis plot (freq , spectre_t)..

Mais ce n'est surement pas ça.. à quoi correspond t et s ? pourquoi dis tu que tu dois récupérer uniquement la 2ème colonne alors que tu te sers de 's' et 't' ?

bon courage

[mariellemm]

Bonjour newenda,

Tout d'abord merci de ton aide!
Par conter je ne comprend pas très bien, comment ça fonctionne...
Que signifie '%g',[2 inf]? je l'avais mis sans savoir!
Et puis je ne comprend ce qu'est NFFT...

En fait dans mon fichier texte, j'ai la première colonne qui est le temps et la deuxième les valeurs du signal. Est ce que pour la transformée je ne devrais pas avoir finalement une fonction dépendant du temps aussi?

Enfin les 2 dernières lignes de ton code, c'est pour faire le tracé?

Encore merci pour ton aide!

[duf42]

Bonjour,

Que signifie '%g',[2 inf]? je l'avais mis sans savoir!

Code :

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

1
2
 
help fscanf

Et puis je ne comprend ce qu'est NFFT...

La encore je pense qu'un tour sur la doc serait intéressant:

Code :

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

1
2
 
help fft

Bonne journée,
Duf

[Newenda]

'%g',[2 inf] > ouvrir des variables de type g (voir opérateur matlab) sur [2 inf]..
fait help fscanf pour en savoir plus.

Bon en gros pour tracer un fft tu doit faire un plot(freq, spi) où freq est le vecteur de tes coordonnées fréquentielles (attention ce n'est pas fréquence d'échantillonage) et spi est le module du signal au carré
donc :

Code :

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

1
2
 
freq = [0:1/tmax:f-1/tmax]

où f = fréquence d'échantillonnage, tmax = valeur de ton temps maximum (pour l'avoir suffit de faire un max(t) si t est ton vecteur de temps défini)
Tu peux choisir la fréquence d'échantillonnage f mais attention au théorème de Shannon ! voir Wikipédia, Le théorème de Nquyst-Shannon est le plus important en calcul du signal numérique, en gros il stipule que ta fréquence d'échantillonnage doit plus importante que celle de ton signal pour bien échantillonner

et

Code :

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

spi = abs(X).^2

où X = fft(x)

PS :

si tes vecteurs comportent N éléments alors :

N = f x tmax

cela te permet de connaitre f

[mariellemm]

oki, c'est déjà plus clair, je suis également aller consulter l'aide.
Ce pendant j'ai un problème avec le fscanf, il me met ça:

Code :

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

1
2
3
4
5
 
fiche=fopen('a47.txt','r') ;
donnee=fscanf(fiche,'%g',[2 inf]); 
??? Error using ==> fscanf
Invalid fid.

Je ne vois pas ou est l'erreur... fid normalement c'est le nom de fichier définis par le fopen il me semble.
Une idée de ce qui bloque?

[duf42]

A priori il ne doit pas etre dans le bon répertoire.

Tu devrais vérifier l'existence du fichier (EXIST) et/ou aussi utiliser le chemin complet pour l'ouvrir (FULLFILE).