How to determine phase of a sine wave using FFT
187 次查看(过去 30 天)
显示 更早的评论
Dear MATLAB experts,
For a small project that I want to do, I need to compute the phase of a sine wave. However, I was unable to obtain a correct result. I have posted my code with this message (please see attachment) for your review and comments. If anyone could help me with this, I would greatly appreciate it.
Thank you!
Fuh
采纳的回答
Star Strider
2016-3-25
编辑:Star Strider
2016-3-25
The code hasn’t posted. Use the angle function on the output of the fft function:
t = [...]; % Time Vector
s = [...]; % Signal Vector
Ts = mean(diff(t)); % Sampling Time
Fs = 1/Ts; % Sampling Frequency
Fn = Fs/2; % Nyquist Frequency
L = length(s);
fts = fft(s)/L; % Normalised Fourier Transform
Fv = linspace(0, 1, fix(L/2)+1)*Fn; % Frequency Vector
Iv = 1:length(Fv); % Index Vector
amp_fts = abs(fts(Iv))*2; % Spectrum Amplitude
phs_fts = angle(fts(Iv)); % Spectrum Phase
Note — This is obviously UNTESTED CODE but it should work.
6 个评论
H S
2018-1-6
Dear Star Strider,
I think for the zero frequency (DC), your code after "amp_fts = abs(fts(Iv))*2;" should be corrected as follows:
amp_fts(1)=amp_fts(1)/2;
更多回答(1 个)
Fuh-Cherng
2016-3-27
1 个评论
Star Strider
2016-3-27
As always, my pleasure.
I would appreciate it then if you would Accept my Answer.
You can do that by deleting your Answer and then Accepting mine. (No one gets any credit for Accepting their own Answer, anyway.)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!