|
noise = sqrt(-2*100*log(uniform1)) *cos(2*3.141592* uniform2);
buffer[i][j]+=noise;
}
}
//영상처리된 이미지를 파일로 출력
for(i=0;i<height;i++)
{
fwrite(buffer[i], width, 1, output);
}
fcloseall();
}
2. 원본영상 3. AWGN이 첨가된 영상
노이즈에 의해 영상이 변형되었음을 알 수
|
|
|
noise_gen()
{
sigma = sqrt(energy_sym / (2 * pow(10, snr / 10)));
awgn_sig = sigma * gaussian1();
//gaussian noise=sigma*G(0,1)
return awgn_sig;
}
이론 참고 : [정보통신기술용어해설]
http://www.ktword.co.kr/abbr_view.php?m_temp1=4300 BPSK 기본적인 이론
C Simulation을 하기 위한 코드
|
|
|
noise_gen()
{
double gaussian;
gaussian = gaussian1();
sigma = sqrt(energy_sym / (2 * pow(10, snr / 10)));//
awgn_sig = sigma * gaussian;
return awgn_sig;
}
이론 참고 : [정보통신기술용어해설]
http://www.ktword.co.kr/abbr_view.php?m_temp1=4300 C 이용한 BPSK 시스템 Rayleigh 채널에서의 BER
|
|
|
noise;
->AWGN 채널을 거친 신호를 표현하였다.
utput = (receive+1)./2;
->수신신호를 비트 넘버로 변환하였다.
compare= xor(output , bit);
->수신신호와 송신신호비교(같으면0 다르면1을 의미한다)
error= sum(compare)/1000000;
->compare의 합(다른bit의 개수
|
|
|
noise_c=randn(1000000,1)*sqrt(var_n); %in-phase AWGN 표현
noise_s=randn(1000000,1)*sqrt(var_n); %quadri-phase AWGN 표현
SNR =[0 1 2 3 4 5 6 7 8 9 10]; %SNR을 0부터 10까지 나타냄
for a= 1:11 %for문을 11번 돌림
Eb(a)=2*var_n*10.^(SNR(a)/10); %SNR에 따른 bit에너지
E(a)=2*Eb(a); %symbol 에너지
si
|
|
메뉴펼치기
