xu ly tin hieu so he thong

Xử Lý Tín Hiệu Số & Hệ Thống Xử Lý Tín Hiệu Số & Hệ Thống Nguyễn Duy Bình-Kỹ sư chất lượng cao K53 Xử Lý Tín Hiệu Số & Hệ Thống Project COMPUTING THE DTFT: FINITE-LENGTH SIGNALS Trong project tập trung xử lý tín hiệu có độ dài hữu hạn Điều sử dụng hàm DTFT cho Project làm việc với tín hiệu dạng xung phổ tần số DTFT chúng ví dụ dễ cho tính tốn phổ tần số tín hiệu Chú ý: Chúng ta cần hàm để tính DTFT Hàm MATLAB freqz dùng để xác định trường hợp tín hiệu infinite-length, hàm cần để tính DTFT tín hiệu finite-length Hàm gọi DTFT(h,N) function [H,W]=DTFT(h,N) % DTFT calculate DTFT at N equally spaced frequencies % 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=(pi-100*eps)); if(~isempty(jkl)) kk=[jkl 1:(jkl(1)-1)]; Wflipped(jkl)=Wflipped(jkl)-2*pi; Wflipped=Wflipped(kk); G=G(kk); end 2.1 Zero—Phase Signals, Tín hiệu pha Nguyễn Duy Bình-Kỹ sư chất lượng cao K53 Xử Lý Tín Hiệu Số & Hệ Thống Xử lý tín hiệu zero-phase điều khó khăn MATLAB hàm DTFT Project cho tín hiệu n=0 tín hiệu zero-phase phải có tính đối xứng xung quanh điểm n=0 Một cách để xử lý vấn đề ta thay đổi hàm DTFT cách thêm vào biến xác định thời điểm bắt đầu tín hiệu Thời điểm ban đầu làm thay đổi giá trị hàm DTFT theo hệ thức DTFT x[ n  n0 ]   e i n0 X (ei ) a Xây dựng hàm DTFT_n0(x,n0,N) để thực dịch chuyển thời gian tín hiệu biến n0 function [H,W]=DTFT_n0(h,n0,N) %#ok %-% usage [H,W]=DTFT_n0(h,n0,N) % h finite-length input vector % n0 xac dinh thoi diem bat dau tin hieu % N number of freuencies for evaluation over [-pi,pi] % H DTFT_n0 values % W (2nd output) vector of frequencies for evaluation over % [-pi,pi] % N=fix(N); [X,W]=DTFT(h,N); H=exp(i*n0*W).*X; %#ok b Kiểm tra DTFT_n0 cách tính DTFT xung 21 điểm bắt đầu n=-10 Kết hàm thực chẵn 1 with -1  n  0 elsewhere Ta sử dụng tín hiệu r[ n]   x=ones(1,21); [H,W]=DTFT_n0(x,10,25); H Kết tính với r[n] ,n0=10,N=25 sau : n=1 0.2492 + 0.0000i n=13 n=2 -0.6969 - 0.0000i n=14 n=3 1.0000 + 0.0000i n=15 n=4 -1.0856 - 0.0000i n=16 n=5 0.9126 + 0.0000i n=17 n=6 -0.4778 - 0.0000i n=18 n=7 -0.1831 - 0.0000i n=19 n=8 1.0000 + 0.0000i n=20 n=9 -1.8782 - 0.0000i n=21 n=10 2.7111 + 0.0000i n=22 n=11 -3.3951 - 0.0000i n=23 n=12 3.8438 + 0.0000i n=24 n=25 21.0000 3.8438 - 0.0000i -3.3951 + 0.0000i 2.7111 + 0.0000i -1.8782 + 0.0000i 1.0000 - 0.0000i -0.1831 - 0.0000i -0.4778 - 0.0000i 0.9126 - 0.0000i -1.0856 + 0.0000i 1.0000 - 0.0000i -0.6969 - 0.0000i 0.2492 + 0.0000i Theo bảng số liệu ta thấy phổ tần số hoàn toàn thực giá trị phần ảo Đồng thời giá trị thu đối xứng qua điểm n=0 nen hàm chẵn.Điều phù hợp với lý thuyết c Vẽ đồ thị DTFT Đồ thị phần thực Nguyễn Duy Bình-Kỹ sư chất lượng cao K53 Xử Lý Tín Hiệu Số & Hệ Thống Đồ thị biên độ Vậy từ đồ thị ta lại thấy DTFT tín hiệu r[n] chẵn đối xứng, phù hợp với lý thuyết Giá trị phần ảo cỡ khơng hồn tồn khơng sai số làm tròn q trình tính tốn Chúng ta hồn tồn bỏ qua mà khơng ảnh hưởng tới kết thu d Chứng minh phần ảo pha biến thien từ tới   L R (ei )   r[n]* e in   e in  L L L    (ein  e in )   2 cos(n ) 1 Từ chứng minh ta thấy phần ảo Pha cố định tới L R(ei ) biến thiên từ thành phần Mặt khác Nguyễn Duy Bình-Kỹ sư chất lượng cao K53 Xử Lý Tín Hiệu Số & Hệ Thống  * k N k  N  Do dó pha R(ei ) biến thiên từ tới Cũng cần phải ý đối xứng quan sát số xung lẻ, số xung chẵn phải kể tới trễ pha nửa tín hiệu 2 Triangular Pulse Cho tín hiệu đối xứng khác dạng tam giác sau L - n with  n