매트랩,퓨리에 코딩

본문내용

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(\'위상\');

키워드

매트랩,   퓨리에 코딩,   스펙트럼