Trong project này chúng ta sẽ tập trung xử lý tín hiệu có độ dài hữu hạn. Điều này sẽ sử dụng hàm DTFT cho ở dưới đây. Project làm việc với những tín hiệu dạng xung và phổ tần số DTFT của chúng vì đây là những ví dụ dễ nhất cho tính toán phổ tần số của tín hiệu. Chú ý: Chúng ta cần 2 hàm để tính DTFT. Hàm MATLAB freqz sẽ dùng để xác định trong trường hợp tín hiệu infinite-length, nhưng một hàm mới sẽ cần để tính DTFT của tín hiệu finite-length. Hàm được gọi là DTFT(h,N) function[H,W]=DTFT(h,N) % DTFT calculate DTFT at N equally spacedfrequencies % usage: H=DTFT(h,N) % h: finite-length input vector, whose length is L % N: number of frequencies for evaluation over [-pi,pi] % ==> constraint N>=L % H: DTFT values (complex) % W: (2nd output) vector of freqs where DTFT is computed % N=fix(N); L=length(h); h=h(:); if(N<L) error('DTFT: # data samles cannot exceed # freq samples') end W=(2*pi/N)*[0:(N-1)]'; mid=ceil(N/2)+1; W(mid:N)=W(mid:N)-2*pi; W=fftshift(W); H=fftshift(fft(h,N)); Chúng ta không cần phải đưa vào biến độ dài L, vì nó sẽ được xác định bằng tính độ dài vector h. Ngoài ra, khi tín DTFT là hàm tuần hoàn vùng từ là vùng tần số âm, nên phải sử dụng hàm fftshift để sắp xếp lại miền tần số theo thứ tự. Sử dụng DTFT tính vector H, vẽ đồ thị trên miền [- ] với chú ý là H(1) là tần số lấy mẫu