Cho hệ thống ttbb và nq được mô tả bởi phương trình sai phân

BỘ CÔNG THƯƠNG BÀI ĐÁNH GIÁ THƯỜNG XUYÊN TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP HÀ NỘI Tên học phần: XỬ LÝ TÍN HIỆU SỐ Khoa: Điện tử Mã học phần: FE6049 Hình thức đánh giá: Báo cáo thực hành PHIẾU GIAO BÁO CÁO THỰC HÀNH Tên lớp: 20223FE6049002 Nhóm: Họ tên: Đỗ Chí Hào Mã sinh viên .2018606461 Phạm Văn Huy Mã sinh viên 2020605434 Nguyễn Bá Thế Mã sinh viên 2020607405 Phạm Quốc Việt Mã sinh viên 2021606689 I Nội dung thực hành Khóa: 16 Bài 1: Cho hệ thống TTBB NQ mô tả phương trình sai phân: 3y(n) = 2y(n-2) + y(n-3) + 2x(n) – x(n-2) + 3x(n-3) a) Biểu diễn giản đồ điểm cực, điểm khơng b) Biểu diễn tín hiệu đáp ứng xung hệ thống với n [-20;20] c) Biểu diễn đồ thị phần thực, phần ảo, đáp ứng biên độ, đáp ứng pha đáp ứng tần số d) Biểu diễn tín hiệu vào tín hiệu với n { với n=−3 x ( n)= [-20;20] biết tín hiệu vào với n=0 n−1 với ≤ n≤ với n cònlại Bài 2: Cho tín hiệu : x1(t)=2cos(2 300t), x2(t)= 3cos(2 500t), x3(t)= 2sin(2 800t), x4(t)=2cos(2 1200t), x(t)= x1(t) + x2(t)+ x3(t)+ x4(t) - Biểu diễn x1(t), x2(t), x3(t), x4(t) miền thời gian figure1 - Biểu diễn tín hiệu x(t) miền thời gian miền tần số figure - Thiết kế lọc FIR dùng hàm fir1 với N=70 để lọc tín hiệu có tần số 500Hz Biểu diễn đồ thị mô phổ biên độ H(ej ) theo đơn vị dB, tín hiệu sau lọc miền thời gian miền tần số figure Bài 3: Cho tín hiệu : x1(t)=sin(2 700t), x2(t)= 4sin(2 1000t), x3(t)= 2cos(2 1200t), x4(t)=3cos(2 1500t), x(t)= x1(t) + x2(t)+ x3(t)+ x4(t) - Biểu diễn x1(t), x2(t), x3(t), x4(t) miền thời gian figure1 - Biểu diễn tín hiệu x(t) miền thời gian miền tần số figure - Thiết kế lọc FIR dùng cửa sổ chữ nhật để lọc tín hiệu có tần số 700Hz, 1000Hz Biểu diễn đồ thị mô phổ biên độ H(ej ) theo đơn vị dB, tín hiệu sau lọc miền thời gian miền tần số figure Bài 4: Cho tín hiệu : x1(t)=2sin(2 450t), x2(t)= 4sin(2 650t), x3(t)= 2cos(2 850t), x4(t)=3cos(2 1250t), x(t)= x1(t) + x2(t)+ x3(t)+ x4(t) - Biểu diễn x1(t), x2(t), x3(t), x4(t) miền thời gian figure1 - Biểu diễn tín hiệu x(t) miền thời gian miền tần số figure - Thiết kế lọc IIR Chebysev có Rp = 0,2dB, Rs = 70dB để lọc tín hiệu có tần số 450Hz, 650Hz, 1250Hz Biểu diễn đồ thị mô phổ biên độ H(ej ) theo đơn vị dB, tín hiệu sau lọc miền thời gian miền tần số figure II Yêu cầu sản phẩm + Bài làm trình bày vào giấy A4, đóng bìa ghi rõ họ tên thành viên nhóm, lớp kèm theo phiếu giao + Đối với tập phải trình bày lại đề bài, phần code (cả chương trình hàm, khơng để dạng ảnh) phần hình ảnh kết mô + Thời gian: tuần, bắt đầu thực từ ngày 15/12/2023 đến ngày 27/12/2023 Hà Nội, ngày 14 tháng 12 năm 2022 Giảng viên giảng dạy Phạm Thị Thanh Huyền BÁO CÁO THỰC HÀNH Bài 1: Cho hệ thống TTBB NQ mô tả phương trình sai phân: 3y(n) = 2y(n-2) + y(n-3) + 2x(n) – x(n-2) + 3x(n-3) a) Biểu diễn giản đồ điểm cực, điểm khơng b) Biểu diễn tín hiệu đáp ứng xung hệ thống với n [-20;20] c) Biểu diễn đồ thị phần thực, phần ảo, đáp ứng biên độ, đáp ứng pha đáp ứng tần số d) Biểu diễn tín hiệu vào tín hiệu với n { với n=−3 x ( n)= với n=0 [-20;20] biết tín hiệu vào n−1 với ≤ n≤ với n cònlại Bài làm n = -20:20; num = [3 -2 -1]; dem = [2 -1 3]; subplot(331); zplane(num,dem); title('Toa diem cuc khong'); xlabel('thuc'); ylabel('ao'); subplot(332); h = impz(num,dem,n+1); stem(n,h); title('Dap ung xung cua he thong'); xlabel('n'); ylabel('h(n)'); w = -4*pi:8*pi/511:4*pi; subplot(333); h = freqz(num,dem,w); plot(w/pi,real(h)); grid on; title('Phan thuc cua X(e^(j\omega))'); xlabel('omega/\pi'); ylabel('Bien do'); subplot(334); plot(w/pi,imag(h)); grid on; title('Phan ao cua X(e^(j\omega))'); xlabel('omega/\pi'); ylabel('Bien do'); subplot(335); plot(w/pi,abs(h)); grid on; title('Pho bien cua X(e^(j\omega))'); xlabel('omega/\pi'); ylabel('Bien do'); subplot(336); plot(w/pi,angle(h)); grid on; title('Pho pha cua X(e^(j\omega))'); xlabel('omega/\pi'); ylabel('Bien do'); x = 2.*impseq(-3,-20,20) + 3.*impseq(0,-20,20) + ((2*n-1).*(stepseq(2,-20,20) stepseq(6,-20,20))); subplot(337); stem(n,x); title('Dau vao'); xlabel('n'); ylabel('x(n)'); subplot(338); y = filter(num,dem,x); stem(n,y); title('Dau ra'); xlabel('n'); ylabel('y(n)'); function [x,n] = impseq(n0,n1,n2) if ((n0 < n1) | (n0 > n2) | (n1 > n2)) error('arguments must satisfy n1