这个例子明确地展示了如何计算由基波和两个谐波组成的信号在dBc中的总谐波失真。返回的结果对显式计算进行检查野。

创建一个1khz采样的信号。该信号由振幅为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²+0.005²)/2²)

Tharmdist = -45.0515

R = thd(x)

R = -45.0515

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

设置谐波数为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 *π* 300 * t) + 0.0025 *罪(2 *π* 400 * t);R = thd(x,1000,3)

R = -45.0515

在THD计算中，指定等于3的谐波数会忽略400hz的功率。

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

得到信号的周期图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 *π* 300 * t) + 0.0025 *罪(2 *π* 400 * t);[pxx,f] = periodogram(x,rectwin(length(x))，length(x)，fs);R = thd(pxx,f,3，psd的）

R = -45.0515

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

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

Fs = 10000;T = 0:1/fs:1-1/fs;X = 2*cos(2* *100*t)+0.01*cos(2* *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的音调。信号以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千赫、6.3千赫、8.4千赫、10.5千赫、12.6千赫和14.7千赫的频率下，谐波从噪声中脱颖而出。除第一个频率外，所有频率都大于奈奎斯特频率。谐波分别被混叠为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的100 Hz基波和三个振幅为0.01、0.005和0.0025的300、500和700 Hz的奇数谐波组成。指定谐波数为7。确定THD，在谐波处的功率，以及相应的频率。

Fs = 10000;T = 0:1/fs:1-1/fs;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×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.00$$dB，对应的振幅为 $$10^{-15}$$。

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

Fs = 5e4;F0 = 2.5e3;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];SGN = polyval(amp,cd);(thd(胡志明市,fs);

该图显示了用于计算比率的频谱和被视为噪声的区域。DC级被排除在计算之外。基波和谐波被标记。