matlab计算thd

在 MATLAB 中计算 THD（Total Harmonic Distortion）需要以下步骤：

1. 首先，通过读取音频信号的数据，将其加载到 MATLAB 中，可以使用`audioread`函数。

2. 对音频信号应用快速傅里叶变换（FFT）以获取频谱信息。可以使用`fft`函数。

3. 计算基频频率的幅度，基频频率是频谱中的第一个峰值。可以使用`abs`函数计算幅度。

4. 对除基频以外的频率分量幅度进行平方，并相加得到总谐波分量的总能量。

5. 通过计算谐波能量和基频能量的比值，得到 THD。THD 定义为非基频频率幅度的和对基频频率幅度的比值。

以下是一个简单的 MATLAB 代码示例，用于计算 THD（假设音频数据已加载到变量 `audio` 中）：

audio = audioread('audio.wav');
Fs = 44100; % 样本率，根据音频文件实际情况进行设置
N = length(audio);
fft_audio = fft(audio);
frequencies = (0:N-1)*(Fs/N);
amplitude = abs(fft_audio/N);
fundamental_frequency = frequencies(find(amplitude == max(amplitude), 1));
harmonic_amplitude = sum(amplitude(frequencies ~= fundamental_frequency).^2);
THD = sqrt(harmonic_amplitude) / max(amplitude);
disp(['THD: ', num2str(THD)]);

这样，你就可以使用 MATLAB 计算给定音频信号的 THD。