DSP 매트랩이용하여 DFT 과정(소스코드 및 MATLAB 과정 포함)

목차

1. 개요
2. 실습결과 및 분석
3. 소스코드

본문내용

1. 개요
･ DFT 함수를 구현하고 DFT를 이용해 다양한 신호의 주파수를 해석한다.
･ 연속주기신호를 연속비주기신호로 주기를 확장하고 표본화를 통해 이산 비주기신호를 만든다.

<중 략>

x1 = myfun_dt_sinusoid(n,1,0.1,0);
// 진폭=1, 주파수=0.1, theta=0인 이산정현파
x2 = myfun_dt_sinusoid(n,0.5,0.2,0);
// 진폭=0.5, 주파수=0.2, theta=0인 이산정현파
x3 = myfun_dt_sinusoid(n,0.1,0.3,0);
// 진폭=0.1, 주파수=0.3, theta=0인 이산정현파
[f_hat, Xk] = myfun_dft(x1+x2+x3);
// 위 3가지 다른 정현파들을 하나로 하여 DFT함
figure(2)
subplot(211)
stem(n,x1+x2+x3); // 3개의 정현파를 합한 그래프
subplot(212)
stem(f_hat, abs(Xk)); // 3개의 정현파를 DFT한 그래프
clear;
N = 50;// N=50으로 표본화(50등분)
n = 0:N-1; //시간축 0~49 x1 = myfun_dt_rect_pulse(n,5);
// L이 5인 구형파
x2 = myfun_dt_rect_pulse(n,15);
// L이 15인 구형파
x3 = myfun_dt_rect_pulse(n,25);
// L이 25인 구형파
x4 = myfun_dt_rect_pulse(n,101);
// L이 101인 구형파
[f_hat, Xk1] = myfun_dft(x1); //x1을 DFT함
[f_hat, Xk2] = myfun_dft(x2);
[f_hat, Xk3] = myfun_dft(x3);
[f_hat, Xk4] = myfun_dft(x4);

참고 자료

없음