说明:采样频率 1kHz,采样点数 5k,pwelch 函数中的参数依次为:信号向
量、分段长度 500 采样点、重叠长度 300 采样点、500 点 DFT(使 100Hz 落在频
率窗内)、采样速率。代码运行结果如下:
fftshift 函数通过将频谱中的零频分量平移至数组中心,重新排列信号的傅里
叶变换。使用示例如下:
fs = 100; % sampling frequency
t = 0:(1/fs):(10-1/fs); % time vector
S = cos(2*pi*15*t);
n = length(S);
X = fft(S);
f = (0:n-1)*(fs/n); % frequency range
power = abs(X).^2/n; % power
plot(f,power);
Y = fftshift(X);
fshift = (-n/2:n/2-1)*(fs/n); % zero-centered frequency range
powershift = abs(Y).^2/n; % zero-centered power
plot(fshift,powershift);