这个例子显式地说明了如何在dBc中计算由基频和两次谐波组成的信号的总谐波失真。根据返回的结果检查显式计算th．

创建一个以1khz采样的信号。该信号由振幅为2的100赫兹基频和振幅为0.01和0.005的200和300赫兹的两个谐波组成。得到总谐波失真显式和使用th．

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采样的信号。该信号由200,300和400 Hz的振幅2和三个谐波的100 Hz基础，具有幅度为0.01,0.005和0.0025。

将谐波数量设置为3.这包括基本的。因此，在THD计算中使​​用100,200和300Hz的功率。

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) + 0.0025 *罪(2 *π* 400 * t);r = (thd (x, 1000年,3)

r = -45.0515.

指定3次谐波数忽略THD计算中400hz的功率。

创建一个以1khz采样的信号。该信号由200,300和400 Hz的振幅2和三个谐波的100 Hz基础，具有幅度为0.01,0.005和0.0025。

获取信号的时期PSD估计，并使用PSD估计作为THD的输入。将谐波数量设置为3.这包括基本的。因此，在THD计算中使​​用100,200和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 *π* 300 * t) + 0.0025 *罪(2 *π* 400 * t);[pxx f] =周期图(x, rectwin(长度(x)),长度(x), fs);r = (thd (pxx, f, 3,'psd'）

r = -45.0515.

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

创建一个以10khz采样的信号。该信号由振幅为2的100赫兹基波和振幅为0.01、0.005和0.0025的3次奇数次谐波组成。将谐波数指定为7。确定螺纹。

FS = 10000;T = 0：1 / FS：1-1 / FS;x = 2 * cos（2 * pi * 100 * t）+ 0.01 * cos（2 * pi * 300 * t）+......0.005 * cos(2 *π* 500 * t) + 0.0025 *罪(2 *π* 700 * t);[sxx f] =周期图(x,汉明(长度(x)),长度(x), fs,'力量'）;rbw = enbw(汉明(长度(x)), fs);r = (thd (sxx f rbw 7'力量'）

r = -44.8396.

生成类似于2.1 kHz音调作为输入的弱非线性放大器的输出的信号。信号在10 kHz时采样1秒。计算并绘制信号的功率谱。使用kaiser窗口β= 38进行计算。

FS = 10000;f = 2100;t = 0：1 / FS：1;X = Tanh（SIN（2 * PI * 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。

计算信号的总谐波失真。默认情况下,th将锯齿谐波视为噪音的一部分。

THD（X，FS，7）;

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

(thd (x, Fs 7'别名'）;

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

FS = 10000;T = 0：1 / FS：1-1 / FS;x = 2 * cos（2 * pi * 100 * t）+ 0.01 * cos（2 * pi * 300 * t）+......0.005 * cos(2 *π* 500 * t) + 0.0025 *罪(2 *π* 700 * t);[r，harmpow，harmfreq] = thd（x，10000,7）;[almfreq harmpow]

ans =7×2100.0000 3.0103 201.0000 -321.0983 300.0000 -43.0103 399.0000 -281.9259 500.0000 -49.0309 599.0000 -282.1066 700.0000 -55.0515

偶数次谐波的幂是 $$-3.0.0.$$分贝，相当于 $$1{0.}^{-15.}$$．

在50千赫采样产生频率2.5千赫的正弦。在信号中加入标准差为0.00005的高斯白噪声。将结果通过一个弱非线性放大器传递。绘制螺纹。

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);

该曲线显示用于计算与视为噪声的比率和区域的频谱。DC级别被排除在计算之外。基本和谐波标明。