帮我生成一段MATLAB频谱分析代码

以下是一个简单的MATLAB频谱分析代码，用于分析音频信号的频谱：

% 读取音频文件
[y, Fs] = audioread('example.wav');

% 计算信号长度
len = length(y);

% 计算FFT点数
nfft = 2^nextpow2(len);

% 计算频率轴
freq = Fs/2*linspace(0,1,nfft/2+1);

% 计算FFT
Y = fft(y, nfft)/len;
Y = Y(1:nfft/2+1);

% 计算功率谱密度
P = 2*abs(Y).^2;

% 绘制频谱图
plot(freq, 10*log10(P));
xlabel('Frequency (Hz)');
ylabel('Power/Frequency (dB/Hz)');
title('Power Spectrum');

说明：

1. 首先使用 audioread 函数读取音频文件，获取音频信号 y 和采样率 Fs。

2. 计算信号长度 len。

3. 计算FFT点数 nfft，这里使用 nextpow2 函数计算出最接近 len 的2的次幂。

4. 计算频率轴 freq，其中 linspace 函数用于生成等间隔的频率点。

5. 计算FFT，使用 fft 函数进行快速傅里叶变换，将信号 y 转换为频域信号 Y。

6. 计算功率谱密度 P，这里使用 abs 函数计算幅值，然后平方，再乘以2，得到电平在每个频率点的功率谱密度。

7. 绘制频谱图，使用 plot 函数绘制频谱曲线，使用 xlabel 和 ylabel 函数设置坐标轴标签，使用 title 函数设置图表标题。其中，使用 log10 函数将功率谱密度转换为对数单位。