BCB 2 TN XLSTH

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA ĐIỆN – ĐIỆN TỬ   BÀI CHUẨN BỊ THÍ NGHIỆM XỬ LÝ SỐ TÍN HIỆU BÀI 2: BIẾN ĐỔI DFT VÀ GIẢI THUẬT FFT CHUẨN BỊ LÝ THUYẾT THÍ NGHIỆM Cho tín hiệu rời rạc x(n) = [1; -2; -3; -4] a Viết biểu thức hàm truyền (biến đổi Z) DTFT tín hiệu x(n)? Vẽ phổ biên độ- tần số pha-tần số tín hiệu trên? - Biểu thức hàm truyền (biến đổi Z) x(n): - DTFT x(n): - Phổ biên độ-tần số pha-tần số tín hiệu: b Tính DTF-4 điểm tín hiệu x(n)? Vẽ biên độ pha DTF-4 điểm tín hiệu trên? Nhận xét? - DFT điểm: - Biên độ pha tín hiệu: - Nhận xét: Nếu so sánh đáp án với hình ta nhận thấy kết mơ giống với tính tốn c Tính DTF-8 điểm tín hiệu x(n)? Vẽ biên độ pha DTF-8 điểm tín hiệu trên? Nhận xét? - Biên độ pha tín hiệu: - Nhận xét: Với hình ta nhận thấy đồ thị phổ biên độ đối xứng qua F=4 Đó điểm N/2 với N điểm Đồ thị phổ pha đối xứng ngược qua F=4 d Vẽ sơ đồ thực tính FFT-8 điểm dùng giải thuật phân chia miền thời gian? Giải thuật phân chia miền thời gian tính theo sơ đồ sau: Với Sau tính tốn: e Vẽ sơ đồ thực tính FFT-8 điểm dùng giải thuật phân chia miền tần số? Giải thuật phân chia miền tần số tính tốn theo sơ đồ sau: Với Sau tính tốn: Khảo sát vẽ đồ thị (dạng rời rạc) hàm cửa sổ Gaussian, Hamming với thông số độ dài 64, 128 giá trị = 1, 2, cửa sổ Gaussian * Độ dài 64: - Hàm cửa sổ Gaussian +α=1 +α=2 +α=5 - Hàm cửa sổ Hamming: * Độ dài 128: - Hàm cửa sổ Gaussian: +α=1 +α=2 +α=5 - Hàm cửa sổ Hamming Viết file-M, lưu với tên Bai2_NHOMx_STFT_CB.m, để thực chương trình biến đổi x n f  500 STFT với yêu cầu sau: tín hiệu   lấy mẫu với tần số s Hz khoảng thời x t  cos  50 t  gian t = đến t = s, cho    t  0.1 , x  t   cos  50 t   cos  100 t  0.1  t  0.2 x t  cos  50 t   cos  100 t   cos  150 t  ,   0.2  t  0.3 Sử dụng cửa số Gaussian có độ dài 128 độ dài chồng chéo 30 Vẽ đồ thị biên độ ma trận STFT vừa tìm được? Nêu nhận xét? clc; close all; Fs = 500; Ts = 1/Fs; T = 3; t_1 = 0:Ts:0.1; t_2 = 0.1:Ts:0.2; t_3 = 0.2:Ts:0.3; x1 = cos(50*pi*t_1); x2 = cos(50*pi*t_2) + cos(100*pi*t_2); x3 = cos(50*pi*t_3) + cos(100*pi*t_3) + cos(150*pi*t_3); X = [t_1, t_2, t_3]; Y = [x1,x2,x3]; [S,F,T] = stft(Y,Fs,"Window",gausswin(128),"OverlapLength",30); waterfall(F,T,abs(S(:,:,1))'); zlabel('Amplitude'); ylabel('Sec'); xlabel('Hz'); - Nhận xét: Ta thấy tín hiệu nằm vùng window khoảng [0; 128], với kết cần thu Đối xứng điểm  ... 0:Ts:0.1; t _2 = 0.1:Ts:0 .2; t_3 = 0 .2: Ts:0.3; x1 = cos(50*pi*t_1); x2 = cos(50*pi*t _2) + cos(100*pi*t _2) ; x3 = cos(50*pi*t_3) + cos(100*pi*t_3) + cos(150*pi*t_3); X = [t_1, t _2, t_3]; Y = [x1,x2,x3];... thông số độ dài 64, 128 giá trị = 1, 2, cửa sổ Gaussian * Độ dài 64: - Hàm cửa sổ Gaussian +α=1 +α =2 +α=5 - Hàm cửa sổ Hamming: * Độ dài 128 : - Hàm cửa sổ Gaussian: +α=1 +α =2 +α=5 - Hàm cửa sổ...  50 t   cos  100 t  0.1  t  0 .2 x t  cos  50 t   cos  100 t   cos  150 t  ,   0 .2  t  0.3 Sử dụng cửa số Gaussian có độ dài 128 độ dài chồng chéo 30 Vẽ đồ thị biên độ