comment transformer le domaine de fréquence dans le domaine temporel

J'avais créé un 3 trois différents signal de fréquence et de filtrer le signal que je ne veux pas. Mais quand je l'aide ifft dans matlab, il montre une mauvaise graphique.Comment transformer mon domaine de fréquence du spectre de nouveau dans mes 3 secondes domaine graphique? Ci-dessous mon code est comme ci-dessous:

clc
clear all
Fs = 8192;
T = 1/Fs;                                             
%create tones with different frequency
t=0:T:1;
t2=1:T:2;
t3=2:T:3;

y1 = sin(2*pi*220*t); 
y2 = sin(2*pi*300*t2);
y3 = sin(2*pi*440*t3);                          

at=y1+y2+y3;


figure;
plot(t,y1,t2,y2,t3,y3),title('Tones with noise');

[b,a]=butter(2,[2*290/Fs,2*350/Fs],'stop');
e=filter(b,a,at);
et=(ifft(abs(e)));
figure,
plot(et)

vous effacez les informations de phase lorsque vous ne abs(ffs(..
hmmm...je ne suis pas d'obtenir ce que l'on entend
Je vois que vous editted la question et il l'a...
ok..j'ai tort essuyer..XD

InformationsquelleAutor Ambbie | 2012-12-26

1

Comme il est maintenant, et est dans le domaine de fréquence, en raison de la fft. Vous n'avez pas besoin de la fft. juste plot(e) et vous aurez le temps de domaine filtré de la forme d'onde. Yo peut vérifier les performances du filtre dans la freq. domaine par fft si, juste
```
plot(abs(fftshift(fft(fftshift(e)))));
xlim([4000 5000])
```
Edit:

Votre code comme il est écrit sur la question a le bug suivant: at a exactement 1 seconde de l'info en elle (ou 8192 éléments). Si vous plot(at) vous verrez la somme des fréquences d'accord, mais ils ont tous happen dans le même temps. C'est la façon de résoudre le problème:
```
clear all
Fs = 8192; % or multiply by 3 if needed
T = 1/Fs;                                             
%create tones with different frequency
t=0:T:3;

y1 = sin(2*pi*220*t).*(t<1); 
y2 = sin(2*pi*300*t).*(t<2 & t>=1);
y3 = sin(2*pi*440*t).*(t>=2);                          

at=y1+y2+y3;

[b,a]=butter(2,[2*290/Fs,2*350/Fs],'stop');
e=filter(b,a,at);
figure,
plot(t,e)
```
```
dt=t(2)-t(1);
N=length(at);
df=1/(N*dt);        % the frequency resolution (df=1/max_T)
if mod(N,2)==0
    f_vector= df*((1:N)-1-N/2);     % frequency vector for EVEN length vectors: f =[-f_max,-f_max+df,...,0,...,f_max-df]
else
    f_vector= df*((1:N)-0.5-N/2);   % frequency vector for ODD length vectors f =[-f_max,-f_max+fw,...,0,...,f_max]
end
freq_vec=f_vector;
fft_vec=fftshift(fft(e));

plot(freq_vec,abs(fft_vec))
xlim([0 1000])
```
- mais comment ifft heure de domaine qui est de 0 seconde à 3 seconde que j'ai créé avant. Ma question est de filtrer les 300Hz..donc..je veux montre mon signal sans 300Hz dans le domaine temporel..savez-vous comment faire cela?
- si je veux pleinement filtre mes 300 HZ? est-il possible?
- Bien sûr, vous pouvez, mais n'ai-je répondu à votre question?
- yup...u m'aide beaucoup..^^ reali merci, merci beaucoup. Partout, la Façon de concevoir le filtre qui peut entièrement filtré à la fréquence de 300 Hz?
- Et que signifie, par tis code "y1 = sin(2*pi*220*t).*(t<1); "...pourquoi le t est nécessaire .*(t<1)?
- Si vous prenez juste en bas de la limite inférieure de votre filtre à partir de 290 à 250, et d'exécuter le code que vous aurez tuer presque tous les 300 Hz contribution. Vous pouvez jouer avec les limites (290 et 350 Hz) jusqu'à ce que vous optimiser votre filtre...
- Grâce natan...^^ Et que signifie, par tis code "y1 = sin(2*pi*220*t).*(t<1); "...pourquoi le t est nécessaire .*(t<1)?
- c'est une condition logique qui définit la condition true (valeur 1) pour tous les éléments de t dans la gamme de t<1 et false (ou 0) pour le reste, rendant ainsi y1 l'égalité des sin(2*pi*220*t) seulement à l'intervalle [0,1], et y1=0 pour le reste de la gamme, [1 3].
- oci...Merci beaucoup!!!
- hai natan, je demande à un qn? pourquoi, après la fft, mon 220 fréquence supposons apparaître à 220 partie, mais il apparaît à 660 Hz-il? tous les m fréquence de multiplier par trois?
- voir mon ajout à la réponse sur la façon de la fft d'un signal...
InformationsquelleAutor bla

Vous devez vous connecter pour publier un commentaire.