[디지털 신호처리 및 실습] 잡음에 따른 sine wave의 SNR - matlab을 이용한 원하는 dB값을 가지는 SNR계산, 변형방법, 파형 보고서

소개글

[디지털 신호처리 및 실습] 잡음에 따른 sine wave의 SNR - matlab을 이용한 원하는 dB값을 가지는 SNR계산, 변형방법, 파형 보고서에 대한 보고서 자료입니다.

목차

1) 잡음신호 발생 방법 및 청취
1. Matlab에서 random number발생 관련 함수 및 사용법 조사(rand, randn)
2. Uniform 및 normal(Gauss) distribution 잡음 3초씩 발생프로그램 작성 후, 소리확인

2) SNR 계산하기
1. SNR 계산방법 조사
2. 440Hz sine wave(Fs = 22050Hz)를 3초 발생후, 1)-2 잡음을 각각 섞었을 경우 SNR계산
3. 임의의 mean, var 갖도록 rand, randn 함수 변형방법 연구

3) 원하는 SNR(dB)을 갖도록 신호에 잡음을 혼합하기
1. 원하는 SNR을 갖도록 신호에 잡음을 혼합하기
2. sine wave에 잡음을 mix하여 60, 40 , 30 , 20, 10, 0, -10dB 신호 발생, 청취하는 프로그램 작성(GUI)

본문내용

e == 1
Noise = sqrt(var(x)/(10^(SNR/10)))*randn(1, kiri);
end
n = Noise;
y = x+n;
dByn = 10 *log10(var(x)./var(n))
end
2. sine wave에 잡음을 mix하여 60, 40 , 30 , 20, 10, 0, -10dB 신호 발생, 청취하는 프로그램 작성(GUI)
잡음(DistType = 0-uniform) 일때
dB = 60
f = 440; % 주파수
Fs = 22050; % 샘플링주파수
t = 0:1/Fs:3; % 시간
x = sin(2*pi*f*t); % 임의의 신호
dbdb = [60 40 30 20 10 0 -10]; %dB값
for i = 1:7
figure(i);
[y, n, dByn] = addnoise(x, dbdb(i), 0);
subplot(311); plot(t,x); % 임의의 신호
xlabel('time[sec]'); % x축 ‘시간’
ylabel('Amplitude'); % y축 '진폭‘
title('Original signal(x)');
axis([0 0.01 -1.5 1.5]); % x축, y축 조절
subplot(312); plot(t,n); % 노이즈 신호
xlabel('time[sec]'); % x축 ‘시간’
ylabel('Amplitude'); % y축 '진폭‘
title('Noise signal(n)');
axis([0 0.01 -1.5 1.5]); % x축, y축 조절
subplot(313); plot(t,y); % y = 노이즈 + 원래신호
xlabel('time[sec]'); % x축 ‘시간’
ylabel('Amplitude'); % y축 ‘진폭’
title('Result signal(Original + Noise signal)');
axis([0 0.01 -1.5 1.5]); % x축, y축 조절
pause % 신호대기
end
dB = 40
dB = 30
dB = 20
dB = 10
dB = 0
dB = -10
잡음(DistType = 1-Gaussian) 일때
dB = 60
f = 440; % 주파수
Fs = 22050; % 샘플링주파수
t = 0:1/Fs:3; % 시간
x = sin(2*pi*f*t); % 임의의 신호
dbdb = [60 40 30 20 10 0 -10]; %db값
for i = 1:7
figure(i);
[y, n, dByn] = addnoise(x, dbdb(i), 1);
subplot(311); plot(t,x); % 임의의 신호
xlabel('time[sec]'); % x축 ‘시간’
ylabel('Amplitude'); % y축 '진폭‘
title('Original signal(x)');
axis([0 0.01 -1.5 1.5]); % x축, y축 조절
subplot(312); plot(t,n); % 노이즈 신호
xlabel('time[sec]'); % x축 ‘시간’
ylabel('Amplitude'); % y축 '진폭‘
title('Noise signal(n)');
axis([0 0.01 -1.5 1.5]); % x축, y축 조절
subplot(313); plot(t,y); % y = 노이즈 + 원래신호
xlabel('time[sec]'); % x축 ‘시간’
ylabel('Amplitude'); % y축 ‘진폭’
title('Result signal(Original + Noise signal)');
axis([0 0.01 -1.5 1.5]); % x축, y축 조절
pause % 신호대기
end
dB = 40
dB = 30
dB = 20
dB = 10
dB = 0
dB = -10

키워드

SNR,   matlab,   dB,   계산값 ,   평균,   분산,   addnoise,   SNR계산,   sine wave,   디지털 신호처리