求功率谱matlab

以下是在MATLAB中求功率谱的方法：
1. 直接法：

x = randn(1,1000); % 生成随机信号
N = length(x); % 信号长度
Pxx = (1/(N*Fs))*abs(fft(x)).^2; % 直接法求功率谱
f = 0:Fs/N:Fs/2; % 频率向量
plot(f,Pxx(1:N/2+1)) % 画出功率谱

2. Welch法：

x = randn(1,1000); % 生成随机信号
N = length(x); % 信号长度
window = hamming(128); % 窗函数
noverlap = 64; % 重叠长度
nfft = 256; % FFT长度
[Pxx,f] = pwelch(x,window,noverlap,nfft,Fs); % Welch法求功率谱
plot(f,Pxx) % 画出功率谱

其中，第一种方法是直接对信号进行FFT变换，然后求出幅度的平方作为功率谱。第二种方法是采用Welch法，即将信号分段，对每一段进行FFT变换，然后对所有段的功率谱进行平均，得到最终的功率谱。