Задание №6

Синтезировать смесь гармонического сигнала и гауссовскогошума с отношением сигнал/шум S/N=-40dB, S/N= 0 dBи S/N=40dB. Построить гистограммы распределений.

%%%---------------noisysignals’ generation------------------%%%

SN1=-40; SN2=0; SN3=20; SN4=40; % отношения сигнал/шум

dis=cov(noise_g); % дисперсияшума

A1=10^(SN1/20)*sqrt(2*dis);% амплитудасигнала

s1=A1*sin(2*pi*fs1*t);

sn1=s1+noise_g;% аддитивнаясмесь сигнала и шума

%--------------------------------------------------------------%

figure(3)

subplot(1,2,1); plot(t,sn1); title ('signal + noise (-40dB)'); xlabel('time'); ylabel('amplitude');

gridon; ylim([(min(sn1)-1) (max(sn1)+1)]); xlim([0 0.2]);

subplot(1,2,2); hist(sn1,100);

title ('histogram'); xlabel('range of values'); ylabel('number of elements');

%--------------------------------------------------------------%

Рис. 15. Смесь полезного сигнала и шума (S/N= - 40 dB) и гистограмма распределения.

%--------------------------------------------------------------%

A2=10^(SN2/20)*sqrt(2*dis);

s2=A2*sin(2*pi*fs1*t);

sn2=s2+noise_g;

figure(4)

subplot(1,2,1); plot(t,sn2); title ('signal + noise (0dB)'); xlabel('time'); ylabel('amplitude');

gridon; ylim([(min(sn2)-1) (max(sn2)+1)]); xlim([0 0.2]);

subplot(1,2,2); hist(sn2,100);

title ('histogram'); xlabel('range of values'); ylabel('number of elements');

%--------------------------------------------------------------%

Рис. 16. Смесь полезного сигнала и шума (S/N= 0 dB) и гистограмма распределения.

%--------------------------------------------------------------%

A3=10^(SN3/20)*sqrt(2*dis);

s3=A3*sin(2*pi*fs1*t);

sn3=s3+noise_g;

figure(5)

subplot(1,2,1); plot(t,sn3); title ('signal + noise (20dB)'); xlabel('time'); ylabel('amplitude');

gridon; ylim([(min(sn3)-1) (max(sn3)+1)]); xlim([0 0.2]);

subplot(1,2,2); hist(sn3,100);

title ('histogram'); xlabel('range of values'); ylabel('number of elements');

%--------------------------------------------------------------%

Рис. 17. Смесь полезного сигнала и шума (S/N= 20 dB) и гистограмма распределения.

%--------------------------------------------------------------%

A4=10^(SN4/20)*sqrt(2*dis);

s4=A4*sin(2*pi*fs1*t);

sn4=s4+noise_g;

figure(6)

subplot(1,2,1); plot(t,sn4); title ('signal + noise (40dB)'); xlabel('time'); ylabel('amplitude');

gridon; ylim([(min(sn4)-1) (max(sn4)+1)]); xlim([0 0.2]);

subplot(1,2,2); hist(sn4,100);

title ('histogram'); xlabel('range of values'); ylabel('number of elements');

%--------------------------------------------------------------%

Рис. 18. Смесь полезного сигнала и шума (S/N= 40 dB) и гистограмма распределения.

Список литературы

1. Рабинер Л., Гоулд Р. Теория и применение цифровой обработки сигналов. М.: Мир, 1978.

2. Опенгейм А., Шафер Р. Цифровая обработка сигналов. М.: Связь, 1979.

3. Голденберг Л. М., Матюшкин Б. Д., Поляк М. Н. Цифровая обработка сигналов: Учеб.пособие для вузов. М.: Радио и связь, 1990.

4. Куприянов М. С., Матюшкин Б. Д. Цифровая обработка сигналов. Процессоры, алгоритмы, средства проектирования. СПб.: Политехника, 1998.

5. Дьяконов В. П., Абраменкова И. В. МATLAB 5.0/5.3. Система символьной математики. М.: Нолидж, 1999.

6. Гоноровский И.С. Радиотехнические цепи и сигналы. М.: Радио и связь, 1986.

7. Баскаков С. И. Радиотехнические цепи и сигналы. М.: Высш. шк., 1988.

Поиск по сайту: