用matlab写一个求功率谱的代码

下面是一个简单的MATLAB代码，用于计算信号的功率谱密度。该代码使用快速傅里叶变换(FFT)计算信号的频率分量，并使用Matlab提供的内置函数计算功率谱密度。代码注释中有详细解释。

% 输入信号
x = randn(1, 1024); % 此处使用随机噪声信号，可以自行替换

% 计算功率谱
n = length(x);
xdft = fft(x); % 使用FFT计算信号的频率分量
psdx = (1./(n*2)).*abs(xdft).^2; % 计算功率谱密度

% 创建频率轴
freq = 0:1/n:1-1/n;
freq = (n/2)*freq;

% 画出功率谱图
figure;
plot(freq, 10*log10(psdx(1:n/2+1)))
grid on;
title('单边功率谱密度');
xlabel('频率 (Hz)');
ylabel('功率谱密度 (dB/Hz)');

注意，该代码计算的是单边功率谱密度，即只画出频率范围为0到$\frac{1}{2}$的功率谱图。如果要计算双边功率谱密度，需要改变频率轴的范围。