MATLAB音频找频谱中心和带宽

要找到音频的频谱中心和带宽，需要进行频谱分析。MATLAB可以使用fft函数对音频进行频谱分析。具体步骤如下：

1. 导入音频文件并将其转换为数字信号。
2. 对信号进行窗函数处理，以减小频谱泄漏的影响。
3. 对窗函数处理后的信号进行FFT变换，得到频谱数据。
4. 计算频率和幅度。
5. 计算频谱中心和带宽。

以下是MATLAB代码示例：

% 导入音频文件并将其转换为数字信号
[y, Fs] = audioread('audio.wav');

% 窗函数处理
win = hann(length(y));
y = y .* win;

% FFT变换
nfft = 2^nextpow2(length(y));
Y = fft(y, nfft)/length(y);
f = Fs/2*linspace(0,1,nfft/2+1);

% 计算频率和幅度
Y = 2*abs(Y(1:nfft/2+1));
f = f';

% 计算频谱中心和带宽
fc = sum(f .* Y) / sum(Y);
bw = sqrt(sum((f - fc).^2 .* Y) / sum(Y));

disp(['频谱中心：', num2str(fc)]);
disp(['频谱带宽：', num2str(bw)]);

其中，`audio.wav`为音频文件名，需要将其放在MATLAB当前目录下。`y`为数字信号，`Fs`为采样率。在上述代码中，使用了汉宁窗函数进行处理，`nfft`为FFT变换的点数。`f`为频率，`Y`为幅值。最后，计算频谱中心和带宽并输出结果。