matlab计算THD

要计算一个信号的THD，可以按照以下步骤进行：

1. 对信号进行FFT变换，得到频域信息。
2. 找到信号的基波频率和幅值，以及前n个谐波的频率和幅值。
3. 计算所有谐波的总功率，即除去基波的所有频率分量的幅值平方之和。
4. 计算总畸变功率，即所有谐波的幅值平方之和。
5. 计算THD，即总畸变功率与基波功率之比的平方根，再乘以100%。

下面是一个MATLAB代码示例：

% 生成一个包含基波和前5个谐波的信号
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
f1 = 50; % 基波频率
f = f1*[1:5]; % 前5个谐波频率
A = [1 0.2 0.1 0.05 0.01]; % 每个频率分量的幅值
x = sum(A.*sin(2*pi*f'*t)); % 信号

% 计算THD
N = length(x); % 信号长度
X = fft(x)/N; % FFT变换
X_mag = abs(X(1:N/2+1)); % 取正频率部分
f_axis = (0:N/2)*fs/N; % 频率轴
f_base = f_axis(abs(X_mag)==max(abs(X_mag))); % 基波频率
f_harm = f_axis(ismember(f_axis,f1*[2:5])); % 谐波频率
A_base = X_mag(abs(f_axis-f_base)<1e-3); % 基波幅值
A_harm = X_mag(ismember(f_axis,_harm)); % 谐波幅值
P_harm = sum(A_harm.^2); % 谐波功率
P_base = A_base^2/2; % 基波功率
THD = sqrt(P_harm/P_base)*100; % THD
disp(['THD = ' num2str(THD) '%']);

输出结果为：

THD = 7.9379%