МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ

УО «ПОЛОЦКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ»

 

 

 

Кафедра радиоэлектроники

 

Лабораторная работа № 3

" Исследование шумов"

 

 

 

 

 

Выполнил: Макридин И.А.

Группа 15-ПЭЛ

 

Проверила: Цыро А.А.

 

 

 

Новополоцк

2016 г.

Теория: Шум квантования — ошибки, возникающие при оцифровке аналогового сигнала. В зависимости от типа аналого-цифрового преобразования могут возникать из-за округления (до определённого разряда) сигнала или усечения (отбрасывания младших разрядов) сигнала.

Практическая часть

Периодический сигнал

Код программы:

clc;

Fd=600;

N=1024;

t=(0:N-1)/Fd;

T=N/Fd

df=1/T;

Fs1=60*df;

Fs2=50*df;

q0=2;

y=20

x=sin(2*pi*Fs1*t+q0)+1.5;

figure(1)

plot(t,x,'r')

axis([0,T,min(x)-2,max(x)+2])

title('Signal');

xlabel('Time');

ylabel('Amplitude');

grid on;

 

 

 

 

 

 

Произвольный сигнал

Код программы:

clc;

Fd=600;

N=1024;

t=(0:N-1)/Fd;

T=N/Fd

df=1/T;

Fs1=60*df;

Fs2=50*df;

q0=2;

x=6*tan(2*pi*Fs1*t+q0)+1.5;

figure(1)

plot(t,x,'r')

axis([0,T,min(x)-2,max(x)+2])

title('Signal');

xlabel('Time');

ylabel('Amplitude');

grid on;

 

 

 

 

 

 

 

 

Сгенерировать шум

Код программы:

clc;

Fd=600;

N=1024;

t=(0:N-1)/Fd;

T=N/Fd

df=1/T;

Fs1=60*df;

Fs2=6*df;

q0=2;

x=(t);

y=2;

%x=8*sin(2*pi*Fs1*t+q0)>0;

S=awgn(x,y)

figure(1)

plot(x,S,'r')

axis([0,T,min(x)-2,max(x)+2])

title('signal');

gridon;

 

 

 

Сигнал+шум

Код программы:

clc;

Fd=600;

N=1024;

t=(0:N-1)/Fd;

T=N/Fd

df=1/T;

Fs1=60*df;

Fs2=50*df;

q0=0;

x=sin(2*pi*Fs1*t+q0)+sin(2*pi*Fs1*t+q0);

figure(1)

plot(t,x);

axis([0,T,min(x)-2,max(x)+ 2])

title('signal');

xlabel('time');

ylabel('amplitude');

grid on;

y=x+2*randn(size(t));

figure(2);

plot(t,y);

axis([0,T,min(x)-2,max(x)+2])

title('signal+noise');

xlabel('time');

ylabel('amplitude');

grid on;

Гистограмма шума

Код программы:

clc;

Fd=600;

N=1024;

t=(0:N-1)/Fd;

T=N/Fd

df=1/T;

Fs1=60*df;

Fs2=50*df;

q0=2;

y=20

x=(t);

S=awgn(x,y)

b=100

G=hist(S,b)

figure(1)

bar(G)

axis tight

title('Noise');

grid on;

 

Гистограмма шума+ сигнала

Код программы:

clc;

Fd=600;

N=1024;

t=(0:N-1)/Fd;

T=N/Fd

df=1/T;

Fs1=60*df;

Fs2=50*df;

q0=2;

y=20

x=sin(2*pi*Fs1*t+q0)+1.5;

S=awgn(x,y)

b=100

G=hist(S,b)

figure(1)

bar(G)

axis tight

title('signal');

grid on;

 

Вывод этой лабораторной работе мы научились генерировать шумы и создавать гистограммы шумов и сигналов.