sinad
SINAD 比
構文
説明
[ は信号の全ノイズおよび全高調波歪みのパワー (dB) を返します。r,totdistpow] = sinad(___)
出力引数を設定せずに sinad(___) を使用すると、現在の Figure ウィンドウに信号のスペクトルがプロットされ、基本成分にラベルが付けられます。スペクトルの基本成分、DC 値およびノイズが、さまざまな色で描画されます。プロットの上に SINAD と表示されます。
例
2 つの信号を作成します。どちらの信号も基本周波数は ラジアン/サンプルで振幅が 1、1 次高調波の周波数が ラジアン/サンプルで振幅が 0.025 とします。一方の信号にはさらに、分散が の加法性ホワイト ガウス ノイズを設定します。
これらの 2 つの信号を作成します。再現性のある結果を得るために、乱数発生器を既定の状態に設定します。加法性ノイズがない信号の SINAD を決定し、その結果を理論上の SINAD と比較します。
n = 0:159;
x = cos(pi/4*n)+0.025*sin(pi/2*n);
rng default
y = cos(pi/4*n)+0.025*sin(pi/2*n)+0.05*randn(size(n));
r = sinad(x)r = 32.0412
powfund = 1; powharm = 0.025^2; thSINAD = 10*log10(powfund/powharm)
thSINAD = 32.0412
加法性ノイズを伴う正弦波信号の SINAD を決定します。加法性ノイズの理論上の分散が組み込まれることで SINAD が近似される様子を示します。
r = sinad(y)
r = 22.8085
varnoise = 0.05^2; thSINAD = 10*log10(powfund/(powharm+varnoise))
thSINAD = 25.0515
基本周波数が 1 kHz、単位振幅が 1 の、480 kHz でサンプリングされた信号を作成します。また、この信号は、振幅が 0.02 の 1 次高調波と、分散が の加法性ホワイト ガウス ノイズで構成されます。
SINAD を決定し、その結果を理論上の SINAD と比較します。
fs = 48e4;
t = 0:1/fs:1-1/fs;
rng default
x = cos(2*pi*1000*t)+0.02*sin(2*pi*2000*t)+0.01*randn(size(t));
r = sinad(x,fs)r = 32.2058
powfund = 1; powharm = 0.02^2; varnoise = 0.01^2; thSINAD = 10*log10(powfund/(powharm+varnoise*(1/fs)))
thSINAD = 33.9794
基本周波数が 1 kHz、単位振幅が 1 の、480 kHz でサンプリングされた信号を作成します。また、この信号は、振幅が 0.02 の 1 次高調波と、標準偏差が 0.01 の加法性ホワイト ガウス ノイズで構成されます。再現性のある結果を得るために、乱数発生器を既定の状態に設定します。
信号のピリオドグラムを取得し、このピリオドグラムを sinad の入力として使用します。
fs = 48e4; t = 0:1/fs:1-1/fs; rng default x = cos(2*pi*1000*t)+0.02*sin(2*pi*2000*t)+0.01*randn(size(t)); [pxx,f] = periodogram(x,rectwin(length(x)),length(x),fs); r = sinad(pxx,f,'psd')
r = 32.2109
50 kHz でサンプリングされる、周波数 2.5 kHz の正弦波を生成します。標準偏差 0.00005 のホワイト ガウス ノイズを信号に付加します。結果を弱非線形増幅器に通します。SINAD をプロットします。
fs = 5e4; f0 = 2.5e3; N = 1024; t = (0:N-1)/fs; ct = cos(2*pi*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); sinad(sgn,fs);
プロットには、比率の計算に使用されるスペクトルと、ノイズとして処理される領域が示されます。DC レベルと基本波は、ノイズの計算では考慮されません。基本波にはラベルが設定されます。
