matalb频谱分析代码

以下是使用Matlab进行频谱分析的示例代码：

% 生成一个信号
t = linspace(0, 10, 1000);
f = sin(2*pi*5*t) + sin(2*pi*10*t);

% 计算信号的功率谱密度
[Pxx, freqs] = pwelch(f, hamming(1024), 512, 1024, 1/(t(2)-t(1)));

% 绘制信号的频谱图
plot(freqs, 10*log10(Pxx));
xlabel('Frequency');
ylabel('Power Spectral Density (dB/Hz)');

该代码与上述Python代码类似，首先生成一个包含两个正弦波的信号，然后使用`pwelch`函数计算信号的功率谱密度。在该函数中，第一个参数是信号，第二个参数是窗函数，这里使用的是汉明窗，第三个参数是重叠样本数，第四个参数是FFT点数，第五个参数是采样频率。最后，使用Matlab绘图函数将信号的频谱图绘制出来，使用`10*log10`将功率转换为分贝形式。