此示例明确显示了如何计算由基波和两次谐波组成的信号的dBc中的总谐波失真。根据野．

创建一个以1 kHz采样的信号。该信号由一个振幅为2的100 Hz基波和两个振幅为0.01和0.005的200和300 Hz谐波组成。明确获得总谐波失真，并使用野．

t = 0:0.001:1 - 0.001;x = 2 * cos(2 *π* 100 * t) + 0.01 * cos(2 *π* 200 * t) + 0.005 * cos(2 *π* 300 * t);tharmdist = 10 * log10 ((0.01 ^ 2 + 0.005 ^ 2) / 2 ^ 2)

tharmdist = -45.0515

r = (thd (x)

r = -45.0515

创建一个以1khz采样的信号。信号由振幅为2的100hz基频和振幅为0.01、0.005和0.0025的200hz、300hz和400hz的三个谐波组成。

将谐波数设置为3。这包括基本的。因此，THD的计算采用100hz、200hz和300hz的功率。

t = 0:0.001:1 - 0.001;x = 2 * cos(2 *π* 100 * t) + 0.01 * cos(2 *π* 200 * t) +．..0.005*cos（2*pi*300*t）+0.0025*sin（2*pi*400*t）；r=thd（x，1000,3）

r = -45.0515

指定等于3的谐波数将忽略THD计算中400 Hz时的功率。

创建一个以1khz采样的信号。信号由振幅为2的100hz基频和振幅为0.01、0.005和0.0025的200hz、300hz和400hz的三个谐波组成。

获取信号的周期图PSD估计，并将PSD估计作为thd的输入。将谐波数设置为3。这包括基本的。因此，THD的计算采用100hz、200hz和300hz的功率。

t = 0:0.001:1 - 0.001;fs = 1000;x = 2 * cos(2 *π* 100 * t) + 0.01 * cos(2 *π* 200 * t) +．..0.005*cos（2*pi*300*t）+0.0025*sin（2*pi*400*t）；[pxx，f]=周期图（x，rectwin（长度（x）），长度（x），fs）；r=thd（pxx，f，3，psd的）

r = -45.0515

通过输入汉明窗获得的功率谱和窗的分辨率带宽来确定THD。

创建一个采样频率为10khz的信号。该信号由一个振幅为2的100hz基频和三个振幅为0.01、0.005和0.0025的300、500和700hz的奇数谐波组成。将谐波数指定为7。确定螺纹。

fs = 10000;t = 0:1 / fs: 1 - 1 / f;x = 2 * cos(2 *π* 100 * t) + 0.01 * cos(2 *π* 300 * t) +．..0.005*cos（2*pi*500*t）+0.0025*sin（2*pi*700*t）；[sxx，f]=周期图（x，汉明（长度（x）），长度（x），fs，“权力”）;rbw = enbw(汉明(长度(x)), fs);r = (thd (sxx f rbw 7“权力”）

r = -44.8396

产生一个信号，类似于输出弱非线性放大器与2.1 kHz音调作为输入。信号以10khz采样1秒。计算并绘制信号的功率谱。使用Kaiser窗口β= 38进行计算。

Fs = 10000;f = 2100;t = 0:1 / Fs: 1;x =双曲正切(罪(2 *π* f * t) + 0.1 + 0.001 * randn(1、长度(t));周期图(x, kaiser(长度(x), 38), [], Fs,“权力”）

在噪音中，4.2 kHz、6.3 kHz、8.4 kHz、10.5 kHz、12.6 kHz、14.7 kHz等频率的谐波明显突出。除第一个频率外，所有频率均大于奈奎斯特频率。谐波分别被混叠为3.7 kHz、1.6 kHz、0.5 kHz、2.6 kHz和4.7 kHz。

计算信号的总谐波失真。默认情况下，野将混叠谐波视为噪声的一部分。

(thd (x, Fs, 7);

重复计算，但现在将混叠谐波视为信号的一部分。

thd（x，Fs，7，“别名”）;

创建一个采样频率为10khz的信号。该信号由一个振幅为2的100hz基频和三个振幅为0.01、0.005和0.0025的300、500和700hz的奇数谐波组成。将谐波数指定为7。确定THD，谐波功率，以及相应的频率。

fs = 10000;t = 0:1 / fs: 1 - 1 / f;x = 2 * cos(2 *π* 100 * t) + 0.01 * cos(2 *π* 300 * t) +．..0.005 * cos(2 *π* 500 * t) + 0.0025 *罪(2 *π* 700 * t);[r, harmpow, harmfreq] = (thd (x 10000 7);[harmfreq harmpow]

ans =7×2300.0000 -43.0103 399.0000 -281.9259 500.0000 -49.0309 599.0000 -282.1066 700.0000 -55.0515

偶数次谐波的幂是 $$-3.00$$dB，对应的振幅为 $$10^{-15}$$．

生成频率为2.5 kHz的正弦波，在50 kHz处采样。向信号中添加标准偏差为0.00005的高斯白噪声。将结果通过弱非线性放大器。绘制THD。

fs = 5 e4;f0 = 2.5 e3;N = 1024;t = (0: n - 1) / fs;ct = cos(2 *π* f0 * t);CD = ct + 0.00005*randn(size(ct));Amp = [1e-5 5e-6 -1e-3 6e-5 1 25e-3];胡志明市= polyval (amp、cd);(thd(胡志明市,fs);

图中显示了用于计算比率的频谱和作为噪声处理的区域。直流电平不在计算范围内。基频和谐波被标记。