在50kHz上采样2.5 kHz的频率为2.5 kHz的正弦样品。添加白色高斯噪声,使得信噪比(SNR)为80 dB。
FS = 5E4;f0 = 2.5e3;n = 2048;t =(0:n-1)/ fs;SNR = 80;X = COS(2 * PI * F0 * T);x = x + randn(size(x))* std(x)/ db2mag(snr);
通过多项式表示的弱非线性放大器通过结果。放大器在谐波的频率下引入虚假音调。
amp = [1e-5 5e-6 -1e-3 6e-5 1 25e-3];x = polyval(amp,x);
绘制信号频谱并向SNR注释,验证它具有预期值。这SNR.
功能将基底的功率率计算为噪声底板,忽略DC分量和谐波。
SNR(x,fs);
绘制信号频谱并注释总谐波失真(THD)。这th
功能将谐波的功率比计算为基本并忽略DC分量和噪声底板。
THD(X,FS);
绘制信号频谱并向噪声和失真率注释信号谱并向信号进行注释(SINAD)。这西拉
功能计算谐波和噪声底部基础的功率比。它只忽略了DC组件。
Sinad(X,FS);
验证SNR,THD和SINAD遵守等式
LHS = 10 ^( - snr(x,fs)/ 10)+ 10 ^(thd(x,fs)/ 10)
LHS = 7.2203E-08
RHS = 10 ^( - SINAD(X,FS)/ 10)
RHS = 7.1997e-08
绘制信号频谱并注释无杂散的动态范围(SFDR)。SFDR是最强壮的杂散成分的基础(“刺激”)的功率比。在这种情况下,刺激对应于第三次谐波。
SFDR(X,FS);