매트랩을 이용한 퓨리에 코딩

목차

1. 주기함수 스펙트럼

2. 비주기함수 스펙트럼

본문내용

s = 1/ts;
ws = 2*pi*fs;
T = 1;
x = (abs(t)< T/2) + (abs(t)==T/2)*0.5;
figure(31);
subplot;
plot(t,x);
rangex = max(x) - min(x);
axis([min(t), max(t), min(x)-rangex/4, max(x) + rangex/4]);
xlabel('t');
ylabel('f(t)')
이를 이용한 진폭 및 위상 스펙트럼
-------------------------------------------------------------------------------------------------------
t = linspace(-30, 30, 1024);
dt = t(2) - t(1);
ts = dt;
fs = 1/ts;
ws = 2*pi*fs;
T = 0.5;
x = (abs(t)< T/2) + (abs(t)==T/2)*0.5;
BW = ws/2;
w = linspace(-BW, BW, length(t));
dw = w(2) - w(1);
F = fftshift(fft(x))/fs;
FT = F;
figure(31);
subplot(2,1,1);
plot(w,real(F));
axis([min(w), max(w), min(real(F)), max(real(F))]);
title('진폭');
subplot(2,1,2);
plot(w, angle(F)*180/pi);
axis([min(w), max(w), min(angle(F)*180/pi), max(angle(F)*180/pi)]);
title('위상');

키워드

매트랩,   매트랩 퓨리에,   매트랩 푸리에,   코딩