.WAV 파일에서 FFT 적용

Question

.wav 파일에서 오는 특정 신호의 고속 푸리에 변환을 얻으려고합니다.

.wav 파일은 신호를 0.6 초마다 반복하는 것입니다.

저는 10 번 반복을 얻는 방법을 알아 내려고 결과가 비슷하다는 것을 보여주기 위해 비교하려고합니다.

이

지금은 내가 가지고있는 코드입니다 : 그래서

Fs = 44100; 
cj = sqrt(-1); 

[test,fs]= wavread('3b healthy2.wav'); % File data name 


dt = 1/Fs; 
time = 45.6; 
N = time/dt; 

left=test(:,1); 
right=test(:,2); 


I = left; 
Q = right; 

t = 0:dt:(time-dt); 

n = length(t); 

f = -Fs/2:Fs/n:Fs/2-Fs/n; 

s = I+cj.*Q; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% Smooth the signal 
ss = smooth(s,201); 

sf = fftshift(fft(ss(1:N))); % taking fft 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

figure(2) 

plot(f,(abs(sf))./max(abs(sf))) 

나는 시간 영역을 찾을 수 있습니다 거기에서 나는 하나의 반복은 45.8-45.2에 있음을 발견 프로그램을 분리하면. 나는이 프로그램을 실행 한 후

나는 그래프를 얻을 수 있지만 내가 볼 수있는, 그것은 하나의 반복하지만 45.2

편집 아래

까지 결합 된 모든 반복을 포함하지 않는에서의 그래프 무엇인지 신호는 45 초에 보입니다. 신호가 3 분 동안 계속 반복됩니다. 지금은 당신이 설정에 신호를 분석 할 수있는 창을 원하는 것 같습니다

Answer 1

그 45초 반복의 FFT이 필요합니다.

그렇다면 시간 영역에서 벡터를 자르면됩니다.

left=test(:,1); 
right=test(:,2); 
time = 45; 
interval =1; // per second, for 0.1 second use interval = 0.1; 
w_range = time*Fs: (time+interval)*Fs-1; 
I = left(w_range) // truncate here 
Q = right(w_range) // and here 
n = interval * Fs; 
f = -Fs/2:Fs/n:Fs/2-Fs/n; 
// continue FFT analysis ... 

의미가 있습니까? 원하지 않는 경우 댓글을 남겨주세요.