Frequence of the signal
古いコメントを表示
Hi , How I can calculate frequence of the signal in the matlab ?? Is there any function do that in matlab ?? can you give me sample example to do that .. Help me plazzzzze :"(
4 件のコメント
Dr. Seis
2011 年 12 月 13 日
Are you trying to find the frequency spectrum of your signal? Or is your signal only one frequency and you are trying to find what that one frequency is? Or are you trying to figure out the instantaneous frequency at some location in your signal?
Sean de Wolski
2011 年 12 月 13 日
Or how often your signal occurs?
fatimah
2011 年 12 月 13 日
Jan
2011 年 12 月 13 日
Who is plazzzzzze ?
回答 (1 件)
Paulo Silva
2011 年 12 月 13 日
From the MATLAB documentation of the fft function:
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid
x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
y = x + 2*randn(size(t)); % Sinusoids plus noise
plot(Fs*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
plot(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
%I added the next lines to find the value
%find maximum value, it should be the fundamental frequency (approximated)
[C,I]=max(2*abs(Y(1:NFFT/2+1)))
f(I)
6 件のコメント
fatimah
2011 年 12 月 13 日
Dr. Seis
2011 年 12 月 13 日
There will be some differences. I have seen enough TV that the signal of a beating heart on an EKG is not a perfect sinusoidal wave. It is also a little more complicated in that the pulse from a single heart beat is composed of more than one valve opening and closing. So, the lower frequency component should be associated with the frequency of each heart beat and the higher frequency component should be associated with orchestra of opening/closing valves during each heart beat.
fatimah
2011 年 12 月 13 日
Dr. Seis
2011 年 12 月 13 日
The above code is a good example. You will have to change things like the sampling frequency (i.e., Fs) and the vector you store your data in (i.e., x). Have you tried to copy/paste the above code into an .m file to see what it does or how it works? If you have questions about what each line does, try running it a few times with different values to see how it changes things (this is how most people learn new programming languages).
Paulo Silva
2011 年 12 月 13 日
1- Find the peaks (findpeaks function from S.P. toolbox or similar code)
2- Find the time T between two of those peaks
3- F=1/T
fatimah
2011 年 12 月 15 日
カテゴリ
ヘルプ センター および File Exchange で Spectral Measurements についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!