Discussion :

[Caroline_1]

Bonjour,

Voici mon programme :

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
32
33
34
35
36
37
38
39
clear all
close all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1998
load('fichier1.txt')
jour=fichier1(1:length(fichiers1),2);
mois=fichier1(1:length(fichier1),3);
annee=fichier1(1:length(fichier1),4);
heure=fichier1(1:length(fichier1),5);
 
conversion_mois % conversion du mois en seconde
t=((jour-1)*24*3600+heure*3600+mois);
y=fichier1(:,6);
 
c=zeros(1:2*(length(y)),1);
 
%Interpolation des données (pour pouvoir observer des fréquences allant jusque 1/10min sur le spectre)
for i=1:length(y)
    c(2*(i-1)+1)=y(i);
end
for i=1:length(c)/2
    c(2*i)=(c(2*i-1)+c(2*i+1))/2;
end
 
%%FFT
Y=fft((c));
N=length(c);
df=1/(max(t)-min(t));
f=(0:1:N/2)*df;
figure(1)
 
hold on
 
loglog(f(2:N/2+1),abs(Y(2:N/2+1)))
a=f(2:N/2+1);
b=abs(Y(2:N/2+1));
title('S(f) - 1950')
xlabel('f (Hz)')
ylabel('S(f)')
hold off

D'une part j'aimerai savoir si vous connaissez une commande permettant de faire une boucle pour charger les noms de fichiers (car j'ai beaucoup de fichiers à charger et à manipuler et je veux éviter les copier-collers)
D'autre part je veux filtrer mon signal avec un filtre passe bas pour n'observer le spectre que sur la bande de fréquence [1/5400 1/600].
J'ai essayé d'utiliser filter, fir1 mais je n'ai pas réussi et ne suis pas sûr de m'y être bien pris. Ce serait sympa si vous pouviez m'aider.

Merci !

[Jerome Briot]

D'une part j'aimerai savoir si vous connaissez une commande permettant de faire une boucle pour charger les noms de fichiers (car j'ai beaucoup de fichiers à charger et à manipuler et je veux éviter les copier-collers)

Voir la : Comment nommer une suite de fichiers : nom_1.ext, nom_2.ext... nom_N.ext ?

[phryte]

Bonjour.

un filtre passe bas pour n'observer le spectre que sur la bande de fréquence [1/5400 1/600].

Tu peux utiliser Butterworth.

Code :

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

[numz,denz]=butter(n,[f1,f2]/(fe/2))

avec n l'ordre du filtre (multiple de 2)
puis :

Code :

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

s=filter(numz,denz,signal)

[Caroline_1]

Tout d'abord merci. J'ai chargé mes fichiers comme ceci :

Code :

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

1
2
3
4
for n=1:1000
load(['fichier0' num2str(n) '.txt'])
A=[str2mat(['fichier0',num2str(n)])]
end

Mais le problème que j'ai maintenant c'est que A est une chaîne de caractère à savoir 'fichier01' et pas la matrice associée au fichier01 comme je le voudrais. Savez-vous comment résoudre cela ?

D'autre part j'ai utilisé le filtre suivant :

Code :

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

1
2
3
fs=1/600;
Wn=[1/5400 1/601]/fs;
b=fir1(length(y)-1,Wn,'bandpass')

et si je veux l'appliquer à mon signal défini par le vecteur y il suffit que je fasse y*1./b ? parce que je ne sais pas comment l'appliquer d'autant plus qu'avec cette solution il faut obligatoirement que l'ordre de mon filtre soit : length(y)-1.

Merci beaucoup à vous et bonne soirée !

[Caroline_1]

Désolé j'avais oublié ceci :

la solution ne serait pas :

Code :

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

filter=(b,1,y)

?
Dans ce cas comment est-ce que je choisis l'ordre de mon filtre ?

Encore bonne soirée !

[Jerome Briot]

Il faut faire comme ceci dans ce cas :

Code :

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

1
2
3
4
5
for n=1:1000
   str = ['fichier0' num2str(n)];
   S = load([str '.txt']);
   A = S.(str)
end