matlab에서 convolution의 multiplication성질을 확인할 수 있는 소스코드 (한 신호는 sinc이고, 다른 한 신호는 임의의 신호)

소개글

고려대학교 전자공학과 3학년 학생입니다.
디지털신호처리실험 시간에 matlab으로 한 실험결과보고서입니다.

time domain에서의 신호 x(t)를 만들고 이의 푸리에변환은 X(f)라고 하고
sinc 함수인 h(t)를 만들어서 이의 푸리에변환은 H(f)라고 할 때

x(t)와 h(t)의 convolution한 결과의 푸리에변환과
X(f)*H(f)의 결과가
서로 같음을 보이는 소스코드입니다

소스코드 각 줄마다 주석을 자세히 달았습니다

1번은 소스코드이고
2번은 결과그래프입니다

목차

① 소스 코드
② 그래프

본문내용

① 소스 코드
clc
clear

nfft = 2048;
Fs = 20; % sampling frequency를 20Hz으로 한다.
t = -10:1/Fs:10; % t를 –10부터 10의 범위에서 1/Fs만큼의 간격으로 만든다.

x1 = cos(2*pi*1*t); % 시간 영역의 신호 x1(t)를 만든다.
x2 = 4*sin(2*pi*4*t); % 시간 영역의 신호 x2(t)를 만든다.
x3 = 2*cos(2*pi*7*t); % 시간 영역의 신호 x3(t)를 만든다.
x = x1 + x2 + x3; % 시간 영역의 신호 x1(t), x2(t), x3(t) 더하여 합성된 신호 x(t)를 만든다.

f = 5;
L = length(x) ; % x의 길이를 L이라고 한다.
X = fft(x,nfft)/L; % nfft-점 DFT를 만들고 x의 길이인 L로 나누어 X라고 정의한다.

참고 자료

없음