why doubling is done from 2 to end-1 only in plotting single sided spectrum: P1(2:end-1) = 2*P1(2:end-1); ?
3 次查看(过去 30 天)
显示 更早的评论
The full code is:
Y = fft(X);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
plot(f,P1);
title('Single-Sided Amplitude Spectrum of X(t)');
8 个评论
采纳的回答
John D'Errico
2018-2-27
编辑:John D'Errico
2018-2-27
end (when used as a subscript) is a shorthand to indicate the last position for the corresponding index.
So V(end) refers to the last element of a vector V. V(2:end) takes all elements but the first in a vector. And V(2:end-1) is the 2nd through penultimate element. This also applies to array indexing.
Therefore
P1(2:end-1) = 2*P1(2:end-1);
merely doubles the indicated elements in the vector P1.
I suppose use of the end keyword might have been confusing to some. They might have created some new keyword name for use in subscripting. They might, for example, have given us a lastIndex keyword, which might be clearer to some. But that in turn is longer to write, and establishes an entirely new keyword, possibly disabling some older code. Since end already existed as a keyword, but one that had no valid context in a subscript, it was deemed safe to use for indexing.
1 个评论
P Sun
2021-6-10
Hi John, thank you for your answer.
Do you have any idea about why don't use the first and last values?
Many thanks.
另请参阅
类别
在 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!