매트랩 주파수 변환 코드, 결과

소개글

매트랩 주파수 변환 코드

목차

없음

본문내용

f0 = 4; % 신호의 주파수
fs = 100; % 샘플링 주파수
Ts = 1/fs; % 샘플링 주기
t = 0: Ts : 1-Ts;

n = length(t); % 샘플링 수

y = 2*sin(2*pi*f0*t); % 신호

k = 0:n-1 % 0~N-1
T = n/fs; % 주파수 간격
freq = k/T; % 주파수 범위

Y = (fft(y)/n)*2; % 노말라이즈
cutoff = ceil(n/2);

Y = Y(1:cutoff);
freq= freq(1:cutoff);
subplot(2,1,1);

stem(t, y);
title(`Sampling Sine Wave`);
xlabel(`time(seconds)`); ylabel(`Amplitude`);

subplot(2,1,2);
stem(freq, abs(Y));
axis([0 20 0 2.5]); % 범위제한
grid on;
title(`Fourier transform`);
xlabel(`sampling number`);
ylabel(`Amplitude`);

참고 자료

없음