假设你有FFt分解数据，用Matlab编写一段计算THD的程序

在MATLAB中，计算总谐波失真（Total Harmonic Distortion, THD）通常需要先对信号进行傅里叶变换（Fast Fourier Transform, FFT），然后分析非线性成分。以下是简单的步骤来编写一个计算THD的MATLAB程序：

% 假设data是一个包含离散时间信号的向量
data = ... % 输入你的实际数据

% 使用FFT函数进行频谱分析
N = length(data);
Y = fft(data);           % 对数据进行FFT
Pxx = abs(Y/N);          % 平方并取平均值，得到每个频率点的功率谱密度

% 确定基频索引（假设输入是单周期信号）
f_s = Fs;                 % 采样率
base_freq_idx = round(f_s/(2*N));   % 根据采样频率计算基频位置

% 提取基频分量和剩余频率的功率
P_base = Pxx(base_freq_idx);
P_other = sum(Pxx(1:base_freq_idx-1) + Pxx(base_freq_idx+1:end));

% 计算总谐波失真
THD = 10*log10(P_other/P_base);    % 转换为dB表示，10*log10表示对数转换

% 结果
fprintf('The Total Harmonic Distortion (THD) is %.2fdB.\n', THD);