Отношение сигнала к шуму и искажениям
[ возвращает общую шумовую и гармоническую степень искажения (в дБ) сигнала.r,totdistpow]
= sinad(___)
sinad(___) без выходных аргументов строит спектр сигнала в окне текущей фигуры и помечает его основной компонент. Это использует различные цвета, чтобы чертить основной компонент, значение DC и шум. SINAD появляется выше графика.
Создайте два сигнала. Оба сигнала имеют основную частоту рад/отсчет с амплитудой 1 и первая гармоника частоты рад/отсчет с амплитудой 0.025. Один из сигналов дополнительно имеет аддитивный белый Гауссов шум с отклонением .
Создайте два сигнала. Установите генератор случайных чисел на настройки по умолчанию для восстанавливаемых результатов. Определите 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 кГц и модульной амплитудой, произведенной на уровне 480 кГц. Сигнал дополнительно состоит из первой гармоники с амплитудой 0.02 и аддитивный белый Гауссов шум с отклонением .
Определите 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 кГц и модульной амплитудой, произведенной на уровне 480 кГц. Сигнал дополнительно состоит из первой гармоники с амплитудой 0.02 и аддитивный белый Гауссов шум со стандартным отклонением 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
Сгенерируйте синусоиду частоты 2,5 кГц, произведенные на уровне 50 кГц. Добавьте Гауссов белый шум со стандартным отклонением 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 и основной принцип исключены из шумового расчета. Основной принцип помечен.