BÁO cáo THÍ NGHIỆM môn xử lý số tín HIỆU bài 1 mô phỏng hệ thống và tín hiệu rời rạc bằng MATLAB

TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI VIỆN ĐIỆN TỬ - VIỄN THƠNG  BÁO CÁO THÍ NGHIỆM MƠN: XỬ LÝ SỐ TÍN HIỆU Mã học phần Sinh viên MSSV Lớp Mã lớp thí nghiệm : : : : : Hà Nội, 2020 BÀI Mô hệ thống tín hiệu rời rạc MATLAB A Tín hiệu hệ thống rời rạc miền n 1.1 Viết chương trình tạo dãy thực ngẫu nhiên xuất phát từ n1 đến n2 có giá trị biên độ theo phân bố Gauss với trung bình 0, phương sai Yêu cầu chương trình có tham số đầu vào đầu nhập theo câu lệnh với cú pháp: [x,n] = randnseq(n1,n2); Điền câu lệnh vào phần trống đây: 1.2 Viết chương trình tạo hàm lượng dãy u cầu chương trình có tham số đầu vào đầu nhập theo câu lệnh với cú pháp: Ex = energy(x,n); Điền câu lệnh vào phần trống đây: 1.3 Cho x ( n )= {1 , 2, , ,5 , , , ,5 , , ,2 , }−2 ≤ n≤ 10 Viết chương ↑ trình thể đồ thị dãy sau đây: a x ( n ) =2 x ( n−5 )−3 x (n+ 4) b x ( n ) =x ( 3−n )−x ( n ) x( n−2) Điền câu lệnh vào phần trống đây: Vẽ phác hoạ đồ thị vào phần trống đây: 1.4 Cho hệ thống mô tả phương trình sai phân tuyến tính hệ số sau: y ( n )− y ( n−1 )+ 0.9 y ( n−2 )=x (n) Sử dụng hàm filter MATLAB, viết chương trình thực cơng việc sau: a Biểu diễn đồ thị hàm đáp ứng xung đơn vị hệ thống với -20 ≤n ≤100 b Biểu diễn đồ thị dãy đáp ứng hệ thống với -20 ≤n ≤100 dãy đầu vào dãy nhảy đơn vị Điền câu lệnh vào phần trống đây: Vẽ phác hoạ đồ thị vào phần trống đây: B Tín hiệu hệ thống rời rạc miền Z, miền tần số liên tục ω, miền tần số rời rạc k 1.5 Cho dãy x ( n )=0,5 n u(n) a Dựa định nghĩa biến đổi Z, tìm biến đổi Z dãy b Kiểm chứng lại kết câu a hàm ztrans c Từ kết trên, tìm biến đổi Fourier x(n) Dùng MATLAB thể đồ thị phổ X ( e jω ) 501 điểm rời rạc khoảng [0,π] Điền câu lệnh vào phần trống đây: Vẽ phác hoạ đồ thị vào phần trống đây: 1.6 Cho dãy x(n) có dạng sau: x ( n )= {… , , , 1, , , , , ,0 , … } ↑ Đây dãy số xác định khoảng hữu hạn từ -1 đến Dựa công thức định nghĩa biến đổi Fourier, viết chương trình tính thể phổ dãy x(n) 501 điểm rời rạc khoảng [0,π] Cho dãy x ( n )=rect (n) Điền câu lệnh vào phần trống đây: Vẽ phác hoạ đồ thị vào phần trống đây: 1.7 Một hàm miền Z cho với công thức sau đây: X ( z )= z z −4 z +1 Hàm số X(z) viết dạng tỷ số hai đa thức theo z−1 sau X ( Z )= z z−1 0+ z−1 = = z 2−4 z +1 3−4 z −1+ z −2 3−4 z−1 + z−2 a Sử dụng lệnh residuez MATLAB, tính điểm cực, thặng dư điểm cực b Từ kết câu trên, viết công thức khai triển X(z) thành tổng phân thức đơn giản, từ tìm biến đổi Z ngược X(z), cho biết x(n) dãy nhân c Kiểm chứng lại kết câu b hàm iztrans Điền câu lệnh vào phần trống đây: 1.8 Cho hàm X(z) với công thức sau: X ( z )= ( 1−0,9 z−1 ) (1+0,9 z−1 ) a Viết chương trình tính điểm cực, thặng dư điểm cực hàm X(z) (gợi ý: dùng hàm poly MATLAB để khôi phục lại đa thức mẫu số từ mảng nghiệm đa thức - mảng điểm cực X(z)) b Từ kết câu trên, viết công thức khai triển X(z) thành tổng phân thức đơn giản, từ tìm biến đổi Z ngược X(z) miền |z|>0,9 Điền câu lệnh vào phần trống đây: 1.9 Cho hệ thống nhân biểu diễn phương trình sau: y ( n )−0,9 y ( n−1 )=x (n) a Tìm hàm truyền đạt hệ thống Sau thực cơng việc sau: b Dùng lệnh zplane MATLAB biểu diễn đồ thị mặt phẳng Z phân bố điểm cực điểm khơng c Tính biểu diễn đồ thị hàm đáp ứng tần số H (e jω )của hệ thống (bao gồm đáp ứng biên độ- tần số đáp ứng pha - tần số) 200 điểm rời rạc đường tròn đơn vị Điền câu lệnh vào phần trống đây: Vẽ phác hoạ đồ thị vào phần trống đây: 1.10 Tạo hàm thực việc biến đổi Fourier rời rạc thuận (đặt tên hàm dft) Fourier rời rạc ngược (đặt tên hàm idft) Dựa hàm dft xây dựng trên, tìm biến đổi Fourier rời rạc dãy có chiều dài N=20: { x ( n )= 0≤ n ≤ n lại Điền câu lệnh vào phần trống đây: 10 Vẽ phác hoạ đồ thị vào phần trống đây: 2.5 Thiết kế lọc thông thấp theo phương pháp cửa số với tham số đầu vào sau: ω p=0,2 π , R p =0,25 dB ω s=0,3 π , A s=50 dB 17 Tính biểu diễn đồ thị: a Dãy đáp ứng xung lọc lý tưởng b Dãy hàm cửa sổ Hamming c Hàm độ lớn tuyệt đối đáp ứng tần số d Hàm độ lớn tương đối tính theo dB đáp ứng tần số Điền câu lệnh vào phần trống đây: 18 Vẽ phác hoạ đồ thị vào phần trống đây: 2.6 Thiết kế lọc thông thấp theo phương pháp lấy mẫu tần số với tham số đầu vào sau: ω p=0,2 π , R p =0,25 dB ω s=0,3 π , A s=50 dB Giả sử ta chọn đáp ứng xung có chiều dài 60 tương đương với lấy 60 mẫu tần số khoảng [0,2π) Dải thơng có độ rộng 0,2π tương đương với mẫu nhận giá trị Giả sử tiếp trình tối ưu hoá nên chọn dải chuyển tiếp mẫu nhận giá trị T1 = 0,5925 T2 = 0,1099 Vậy dãy mẫu tần số cho sau: { h ( n ) = 1,1,1,1,1,1,1 , T ,T , 0⏟ , … ,0 T ,T , 1,1,1,1,1,1,1 43 mẫu0 Tính biểu diễn đồ thị: a Dãy mẫu tần số 19 } b Dãy đáp ứng xung bộlọc thực tế c Hàm độ lớn tuyệt đối đáp ứng tần số d Hàm độ lớn tương đối tính theo dB đáp ứng tần số Điền câu lệnh vào phần trống đây: 20 Vẽ phác hoạ đồ thị vào phần trống đây: 2.7 Thiết kế lọc thơng thấp theo phương pháp lặp (thuật tốn Parks McClellan) với tham số đầu vào sau: ω p=0,2 π , R p =0,25 dB ω s=0,3 π , A s=50 dB Trước tiên xuất phát từ độ dài dãy đáp ứng M theo công thức M= ω −ω −20 log √ δ δ 2−13 , với ∆ f = s p 2π 14,6 ∆ f Lặp cơng việc tìm lọc tối ưu theo nghĩa Chebyshev (dùng lệnh firpm) tăng M sau lần lặp để tìm lọc thoả mãn yêu cầu thiết kế, sau tính biểu diễn đồ thị: a Dãy đáp ứng xung lọc thực tế b Hàm độ lớn tuyệt đối đáp ứng tần số c Hàm độ lớn tương đối tính theo dB đáp ứng tần số d Hàm sai số E( ω) 21 Điền câu lệnh vào phần trống đây: 22 Vẽ phác hoạ đồ thị vào phần trống đây: B Thiết kế lọc có đáp ứng xung chiều dài vô hạn (bộ lọc số IIR) Các bước thực hành 2.8 Thiết kế lọc thông thấp tương tự, định dạng Chebyshev-I, cửa số với tham số đầu vào sau: ω p=0,2 π , R p =1dB ω s=0,3 π , A s=16 dB Viết chương trình tính biểu diễn đồthị: a Độ lớn đáp ứng tần số b Hàm đáp ứng pha lọc c Hàm độ lớn tương đối tính theo dB đáp ứng tần số d Hàm đáp ứng xung lọc tương tự Điền câu lệnh vào phần trống đây: 23 24 Vẽ phác hoạ đồ thị vào phần trống đây: 2.9 Chuyển đổi lọc với tham số cho phần 2.8 sang lọc số phương pháp biến đổi song tuyến Hàm bilinear cho phép thực việc chuyển đổi Tính biểu diễn đồ thị: a Độ lớn đáp ứng tần số b Hàm đáp ứng pha lọc c Hàm độ lớn tương đối tính theo dB đáp ứng tần số d Trễ nhóm theo tần số Điền câu lệnh vào phần trống đây: 25 26 Vẽ phác hoạ đồ thị vào phần trống đây: 2.10 Thực yêu cầu câu 2.9 theo phương pháp bất biến xung, dùng hàm impinvar MATLAB So sánh kết thu với câu Điền câu lệnh vào phần trống đây: 27 Vẽ phác hoạ đồ thị vào phần trống đây: 2.11 Tạo hàm thực việc chuyển đổi băng tần số, trả hàm truyền đạt lọc với tham số đầu vào hàm truyền đạt lọc thông thấp, hàm đa thức thể phép đổi biến số độc lập, ghi lại theo tên tệp zmapping.m: 28 Điền câu lệnh vào phần trống đây: 29 2.12 Viết chương trình chuyển đổi từ lọc thông thấp theo thiết kế câu 1.9 sang lọc thơng cao có tần số cắt ωc=0,6π Tính biểu diễn đồ thị a Độ lớn đáp ứng tần số b Hàm đáp ứng pha lọc c Hàm độ lớn tương đối tính theo dB đáp ứng tần số d Trễ nhóm theo tần số Điền câu lệnh vào phần trống đây: 30 Vẽ phác hoạ đồ thị vào phần trống đây: 31 ... ( n ) = 1, 1 ,1, 1 ,1, 1 ,1 , T ,T , 0⏟ , … ,0 T ,T , 1, 1 ,1, 1 ,1, 1 ,1 43 mẫu0 Tính biểu diễn đồ thị: a Dãy mẫu tần số 19 } b Dãy đáp ứng xung bộlọc thực tế c Hàm độ lớn tuyệt đối đáp ứng tần số d Hàm...BÀI Mơ hệ thống tín hiệu rời rạc MATLAB A Tín hiệu hệ thống rời rạc miền n 1. 1 Viết chương trình tạo dãy thực ngẫu nhiên xuất phát từ n1 đến n2 có giá trị biên độ... vị hệ thống với -20 ≤n ? ?10 0 b Biểu diễn đồ thị dãy đáp ứng hệ thống với -20 ≤n ? ?10 0 dãy đầu vào dãy nhảy đơn vị Điền câu lệnh vào phần trống đây: Vẽ phác hoạ đồ thị vào phần trống đây: B Tín hiệu