matlab计算谐波含量程序

以下是一个计算谐波含量的 Matlab 程序示例：

% 输入信号
x = [1 2 3 4 5 6 7 8 9 10];

% 采样频率
fs = 1000;

% 计算FFT
X = fft(x);

% 计算频率分量
f = (0:length(x)-1)*fs/length(x);

% 计算谐波的频率
harmonics = 2:5;

% 初始化谐波幅值
harm_ampl = zeros(size(harmonics));

% 循环计算每个谐波的幅值
for i = 1:length(harmonics)
    idx = find(round(f)==harmonics(i)*fs/length(x));
    harm_ampl(i) = abs(X(idx));
end

% 计算谐波总含量
THD = sqrt(sum(harm_ampl.^2))/abs(X(1))*100;

% 显示结果
disp(['Total Harmonic Distortion (THD) = ' num2str(THD) '%']);

这个程序使用 FFT 函数计算输入信号的频谱，并计算谐波的频率。然后使用循环计算每个谐波的幅值，并计算谐波总含量。最后，程序显示 THD 的百分比。