белый шум в матлабе Опции

[vadimuzzz]

есть такая функция y=awgn(x,SNR_dB), вроде все с ней понятно: взял BPSK сигнал x, подмешал шум, декодировал, посчитал коэфф. ошибок. повторил на нескольких SNR по многу раз, результат усреднил, получил кривую ошибок. а теперь строю ту же кривую с пом. функции berawgn(snr_axis,'psk',2,'nondiff') - результат совершенно другой.
сгенерил шум по-другому:

Код

SNR=10.^(SNRdB/10);
No=Eb/SNR;
N=sqrt(No/2)*randn(n);

добавил к сигналу, построил кривую - результат совпал с berawgn. посмотрел по справочникам - графики совпадают с тем, что выдает berawgn. что же

[vadimuzzz]

понятно, но мне и нужен белый шум, отсюда и randn. кажется понял в чем дело - awgn первым параметром принимает отсчеты мощности, а я давал амплитуды. видимо, в этом причина

[samurad]

понятно, но мне и нужен белый шум, отсюда и randn. кажется понял в чем дело - awgn первым параметром принимает отсчеты мощности, а я давал амплитуды. видимо, в этом причина

Судя по хелпу, y=awgn(x,SNRdB) принимает в х амплитуды отчетов сигнала-вектора, и по умолчанию считает, что амплитуды нормированы, так что мощность вектора х есть 0 дБВт; иначе вызов этой функции другой: y=awgn(x,SNRdB,XPowdBW). Кроме того, eсли у вас несколько отчетов на бит, то их амплитуда должна быть в соответственное число раз меньше, чем амплитуда принятого бита, т.к. awgn() не понимает время.

[vadimuzzz]

Судя по хелпу, y=awgn(x,SNRdB) принимает в х амплитуды отчетов сигнала-вектора, и по умолчанию считает, что амплитуды нормированы, так что мощность вектора х есть 0 дБВт; иначе вызов этой функции другой: y=awgn(x,SNRdB,XPowdBW).

я тоже сначала так подумал, но там есть еще вариант y=awgn(x,SNRdB,'measured'), который сам меряет вх. мощность. результаты тоже не сошлись, тогда я проверил "в лоб":

Код

SNRdB=3;
SNR=10.^(SNRdB/10);
Eb=1;
No=Eb/SNR;
n=1000000;
N1=sqrt(No/2)*randn(n,1);
BPSK_sig=-2*randint(n,1)+1;
s1=BPSK_sig+N1;
s2=awgn(BPSK_sig,SNRdB,'measured');
N2=s2-BPSK_sig;

и, что характерно, std(N1)~0.5; std(N2)~0.707. 0.707 подозрительно похож на sqrt(2)/2, из чего я и сделал вывод, что awgn работает с мощностями, в то время как я давал амплитуды.

[samurad]

я тоже сначала так подумал, но там есть еще вариант y=awgn(x,SNRdB,'measured'), который сам меряет вх. мощность. результаты тоже не сошлись, тогда я проверил "в лоб":

Код

SNRdB=3;
SNR=10.^(SNRdB/10);
Eb=1;
No=Eb/SNR;
n=1000000;
N1=sqrt(No/2)*randn(n,1);
BPSK_sig=-2*randint(n,1)+1;
s1=BPSK_sig+N1;
s2=awgn(BPSK_sig,SNRdB,'measured');
N2=s2-BPSK_sig;

и, что характерно, std(N1)~0.5; std(N2)~0.707. 0.707 подозрительно похож на sqrt(2)/2, из чего я и сделал вывод, что awgn работает с мощностями, в то время как я давал амплитуды.

Если отбросить подозрительность и рассмотреть факты (продолжая ваш код), то получается:

Код

BPSK_sig_pow = sum(abs(BPSK_sig.^2))./n   % измеренная мощность BPSK_sig
BPSK_sig_pow =
   1.0
% Для заданного SNR = 3dB, мощность шума должна быть в 2 раза меньше, т.е. 1/2
% Проверяем:
N1_pow = sum(abs(N1.^2))./n   % измеренная мощность шума N1
N1_pow =
   0.25  % ошибка
N2_pow = sum(abs(N2.^2))./n   % измеренная мощность шума N2
N2_pow =
   0.5  % верно

% Для проверки awgn(BPSK_sig,SNRdB):
s21 = awgn(BPSK_sig,SNRdB);  % можно, т.к. Eb = 1 W
N21 = s21 - BPSK_sig;
N21_pow = sum(abs(N21.^2))./n   % измеренная мощность шума N21
N21_pow =
   0.5  % опять верно

Вы ошиблись в моделировании шума N1: для реального сигнала BPSK весь шум так же реальный, поэтому правильно моделировать N1 так

Код

N11 = sqrt(No)*randn(n,1);
N11_pow = sum(abs(N11.^2))./n   % измеренная мощность шума N11
N11_pow =
   0.5 % верно

[vadimuzzz]

Вы ошиблись в моделировании шума N1: для реального сигнала BPSK весь шум так же реальный, поэтому правильно моделировать N1 так

Код

N11 = sqrt(No)*randn(n,1);

тогда не сходятся результаты с berawgn(snr_axis,'psk',2,'nondiff')
обратите внимание, везде берется No/2
http://www.dsplog.com/2007/08/05/bit-error...psk-modulation/

http://www.dsplog.com/db-install/wp-conten...criptbpsk_ber.m

http://gaussianwaves.blogspot.com/2010/04/...ation-over.html