用MATLAB对音频进行分析代码

MATLAB是一种强大的数值计算和数据可视化工具，常用于信号处理，包括音频分析。以下是使用MATLAB对音频文件进行基本分析的一个简单示例：

% 加载音频文件
[y, Fs] = audioread('your_audio_file.wav'); % 'Fs'表示采样频率

% 检查音频信息
disp(['音频样本数：', num2str(length(y))]);
disp(['采样频率：', num2str(Fs), ' Hz']);

% 音频波形显示
t = (0:length(y)-1)/Fs; % 时间轴
plot(t, y);
xlabel('Time (seconds)');
ylabel('Amplitude');
title('Audio Signal');

% 频谱分析（傅立叶变换）
Y = abs(fft(y)); % 计算幅度谱
Pxx = Y(1:length(Y)/2+1) / length(y); % 对称化并取有效值
[freqs, Pxx_spec] = freqz(y, [], Fs); % 频率向量和功率谱密度
subplot(2, 1, 2);
plot(freqs, Pxx_spec);
xlabel('Frequency (Hz)');
ylabel('Power/Frequency');
title('Spectrum Analysis');

% 跟踪特定频率成分
target_freq = 440; % 歌曲中的某个音调
index = find(freqs == target_freq);
if ~isempty(index)
    index = index(1);
    disp(['目标频率 (' num2str(target_freq) ' Hz) 的功率：', num2str(Pxx_spec(index))]);
else
    disp('目标频率未在频谱中找到');
end

这个例子展示了如何读取音频、查看波形、执行快速傅立叶变换（FFT）获取频谱以及定位特定频率。