Tín hiệu và hệ thống
110
t = (0:1/100:10-1/100); % Vector thời gian
x = sin(2*pi*15*t) + sin(2*pi*40*t); % Tín hiệu x(t)
y = fft(x); % Biến đổi DFT của x
m = abs(y); p = unwrap(angle(y)); % Biên độ và pha
f = (0:length(y)-1)*99/length(y); % Vector tần số
subplot(1,2,1);
plot(f,m); title('Bien do');
set(gca,'XTick',[15 40 60 85]);
subplot(1,2,2); plot(f,p*180/pi); title('Pha');
set(gca,'XTick',[15 40 60 85]);
xr = real(ifft(y)); % Biến đổi FFT ngược
figure;
plot(t,x,'b-',t,xr,'r--');
legend('Tin hieu goc','Tin hieu khoi phuc');
Kết quả như sau:
Hình 9.14.
Biến đổi DFT
