BCB 3 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 THIẾT KẾ VÀ THỰC HIỆN BỘ LỌC IIR TRÊN MATLAB Giảng viên hướng dẫn: Huỳnh Văn Phận Họ tên sinh viên: Nguyễn Trần Tú Anh MSSV: 1910770 Lớp: L06 Câu Viết chương trình Matlab thực giải thuật xử lý mẫu dạng trực tiếp để thực lọc IIR tổng qt có bậc M có tín hiệu ngõ vào x Code: clc; close all; a = input('Nhap he so mau:'); b = input('Nhap he so tu:'); M = length(a); Lv= length(b); xn= input('Nhap gia tri ngo vao:'); w = zeros(1,M+1); L = length(xn); v = zeros(1,L); for i = : L v(1) = xn(i); for j = : M w(1) = w(1)-a(j).*w(j+1); end for r = : Lv w(1) = w(1)+b(r).*v(r); end y(i) = w(1); for k = L : -1 : v(k) = v(k-1); end for t = M+1 : -1 : w(t) = w(t-1); end end Kết mô phỏng: Câu Cho lọc IIR có hàm truyền Vẽ sơ đồ khối dạng trực tiếp lọc Vẽ đáp ứng tần số lọc (có thể sử dụng hàm freqz Matlab) Tìm ngõ y(n) ngõ vào x(n) = [1, 2, -1, 1] Code Matlab: clc; close all; a=[2 -3 4]; b=[1 0.2 -0.3 0.5]; [H, w] = freqz(a,b); H_magnitude = abs(H); plot(w/pi,H_magnitude); xlabel('w'); ylabel('Amplitude'); title('Dap ung tan so'); Kết mơ phỏng: Sử dụng chương trình Matlab viết câu hỏi chuẩn bị để xác định ngõ lọc ngõ vào ý x ( n ) = [ 1, 2, −1,1] Kiểm tra tính xác với kết tính Code Matlab: clc; close all; a = [1 0.2 -0.3 0.5]; b = [2 -3 4]; M = length(a); Lv= length(b); xn= [1 -1 1]; w = zeros(1,M+1); L = length(xn); v = zeros(1,L); for i = : L v(1) = xn(i); for j = : M w(1) = w(1)-a(j).*w(j+1); end for r= : Lv w(1) = w(1)+b(r).*v(r); end y(i) = w(1); for k = L : -1 : v(k) = v(k-1); end for t = M+1 : -1 : w(t) = w(t-1); end end Kết mô phỏng: Kết tính ý giống với ý Câu Viết chương trình Matlab thực giải thuật xử lý mẫu dạng tắc để thực lọc IIR tổng quát có bậc M có ngõ vào x Sau tính ngõ với thơng số lọc ngõ vào câu hỏi chuẩn bị Code Matlab: clc; close all; a = [1 0.2 -0.3 0.5]; b = [2 -3 4]; M = length(a); Lv= length(b); xn= [1 -1 1]; w = zeros(1,M+1); Lx= length(xn); y = zeros(1,Lx); for i = : Lx for j = : M w(1) = xn(i)-a(j).*w(j+1); end for r = : Lv y(i)=y(i)+b(r).*w(r); end for t = M : -1 : w(t) = w(t-1); end end Kết mô phỏng: Câu Viết chương trình Matlab thực giải thuật xử lý mẫu (dạng nối tiếp) để thực lọc IIR dạng nối tiếp với K = tầng lọc IIR bậc với ma trận hệ số lọc sos ma trận trạng thái trung gian W có ký hiệu (Để thuận tiện, sinh viên cho ma trận sos, tín hiệu ngõ vào x(n) để kiểm tra chương trình) Code Matlab: clc clear K=input('So tang cua bo loc:'); B=zeros(K,3); A=zeros(K,3); for i=1:K fprintf('cot:%d\n',i); for i1=1:3 b=input('Gia tri cua tu:'); B(i,i1)=b; end end for i=1:K fprintf('cot:%d\n',i); for i1=1:3 a=input('Gia tri cua mau:'); A(i,i1)=a; end end w=zeros(K,3); k=input('nhap cot:'); x=zeros(1,k); for i=1:k g=input('Nhap gia tri cho x:'); x(1,i)=g; end k1=input('So mau can tinh cua y:'); while k1