[전자통신공학]디지털 신호처리-(matlab)Aliasing Phenomenon 또는 중첩현상

소개글

Matlab을 이용하여 합성된 signal을 Fast Fourier Transform 시켜보고, transform시에 사용하는 sampling frequency 따른 Aliasing 현상을 직접 구현하였습니다.
매트랩 소스와 출력파형, 결과 분석이 모두 들어있습니다.

목차

1. Purpose
2. How to simulation
3. Source
4.Plotting & Analysis

본문내용

2개의 sinusoidal signal(sine wave)을 generate하고, 그 신호를 합성한 signal을 fft취하여 plotting한다. 이 때 sampling frequency를 변화시켜 Aliasing이 발생하지 않는 경우와 발생하는 경우를 모두 plotting 해보고 각 경우에 대해 다시 inverse fft를 하여 원래의 파형이 detection되는지를 비교해본다. 또 이 simulation의 기본 원리가 되는 Nyquist Theorem (Sampling정리) 에 대해 논의해본다.

3. Source
% Generate signal & Aliasing 현상을 보기 위해 fs를 조정
fs=1500; %sampling frequency를 1500Hz로 함. 후에 290, 110Hz로 바꿈
Ts=1/fs;
t= 0 : Ts : 0.1-Ts; %시간 축은 1/1500초 간격으로 data를 나타내도록 함
f1 =60;
f2= 150;
a1=2;
a2=5;
sig1=a1*sin(2*pi*f1*t);
sig2=a2*sin(2*pi*f2*t);
totalsig = sig1+sig2; %합성한 신호

figure(1); %한번의 실험에 3개의 plot창이 출력,
subplot(211); %첫번째 창에는 time domain에서 생성한 신호를 plotting
plot(t,sig1,t,sig2); %각각의 sine wave를 plotting
axis([0 0.1 -10 10])
grid
title(`sig1 and sig2 in time domain`)
xlabel(`Time(s)`)
ylabel(`Amplitude`)
legend(`sig1(t)`, `sig2(t)`)
subplot(212);
plot(t,totalsig); %2개의 sine wave를 합성한 신호를 plotting
axis([0 0.1 -10 10])
grid
title(`totalsig = sig1 + sig2`)
xlabel(`Time(s)`)
ylabel(`Amplitude`)
legend(`totalsignal(t)`)

참고 자료

없음