THỰC HÀNH xử lý tín HIỆU số

BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC CÔNG NGHỆ TP.HCM THỰC HÀNH XỬ LÝ TÍN HIỆU SỐ BỌẵố sỒ ố: ThS Ph Ố HỂốg KỌỐ KhẬốh www.hutech.edu.vn TH C HÀNH X ố b ố 2014 LÝ TÍN HI U S M CL C I M CL C M C L C H NG D N BÀI 1: PH N M M MATLAB 1.1 KH I Đ NG MATLAB 1.2 CÁC V N Đ C B N 1.2.1 Các phép toán toán tử 1.2.2 Khai báo biến 1.2.3 Các lệnh thường dùng 1.3 L P TRÌNH TRONG MATLAB 10 1.3.1 Các phát biểu điều kiện if, else, elseif 10 1.3.2 Switch 10 1.3.3 While 10 1.3.4 For 11 1.3.5 Break: 11 1.4 MA TR N 11 1.4.1 Các thao tác ma trận 11 1.4.2 Vector 15 1.4.3 Đa thức 15 1.5 Đ HO 16 1.5.1 Các lệnh vẽ 16 1.5.2 Tạo hình vẽ 16 1.5.3 Kiểu đường vẽ 16 1.5.4 Vẽ với hai trục y 17 1.5.5 Vẽ đường cong 3-D 18 1.5.6 Vẽ nhiều trục toạ độ 18 1.5.7 Đặt thông số cho trục 19 1.5.8 Đồ hoạ đặc biệt 23 1.5.9 Đồ hoạ 3D 25 1.5.10 Thực hành vẽ đồ thị 26 1.6 CÁC FILE VÀ HÀM 28 1.6.1 Script file (file kịch bản) 28 1.6.2 File hàm 28 1.6.3 Các hàm toán học 29 1.6.4 Các phép toán hàm toán học 30 TRANG > TH C HÀNH X LÝ TÍN HI U S 1.6.5 Thực hành script function .31 BÀI 2: TÍN HI U R I R C THEO TH I GIAN 38 2.1 CÁC TÍN HI U S C P 38 2.2 CÁC PHÉP TOÁN 39 2.3 KI M TRA TÍNH CH T TUY N TÍNH VÀ B T BI N .39 2.4 H LTI 41 BÀI 3: BI N Đ I Z 44 3.1 CÁC ĐI M C C VÀ ĐI M KHÔNG .44 3.2 PHÂN TÍCH DÙNG PH NG PHÁP TH NG D 45 3.3 BI N Đ I Z VÀ Z NG C .46 BÀI 4: BI N Đ I FOURIER R I R C 49 4.1 TÍNH DTFT 49 4.2 FFT VÀ CÁC TÍNH CH T 50 BÀI 5: B L CS FIR 53 5.1 CÁC LO I B L C 53 5.2 PH NG PHÁP C A S 55 5.3 PH NG PHÁP L Y M U T N S .58 5.4 PH NG PHÁP L P .60 BÀI 6: B L CS IIR 65 6.1 THI T K B L CT NG T 65 6.2 THI T K B L C S .69 TÀI LI U THAM KH O 73 H H MÔ T NG D N II I NG D N MÔN H C Th Ơ hậốh ỏý ỗặố hỌ Ộ s ỏậ ỐẾố h Ơ h ỗổ ƠhỒ ỐẾố X ỏý ỗặố hỌ Ộ s Ơủa chuyên ốgậốh Kỹ ỗhỘ ỗ ĐỌ ố ỗ TổỘy ố ỗhẾốg MẾố h Ơ ốậy a MATLAB đ ƠẬƠ ỏý ỗhỘy ỗ đử h Ơ ỗổỒốg ỐẾố X kỌ Ố Ơhứốg ỏý ỗặố hỌ Ộ s N I DUNG MÔN H C       Bài Ph ố Ố Ố MATLAB: Ơ b ố ộ MATLAB, ƠẬƠh ỏ ỗổìốh Ơũốg ốh ƠẬƠh ỏý Ốa ỗổ ố, ộ đ ỗhị ỗổỒốg MATLAB BậỌ 2: Tặố hỌ Ộ ổ Ọ ổ Ơ ỗhƯỒ ỗh Ọ gỌaố: ƠẬƠh bỌ Ộ ơỌ ố ỗặố hỌ Ộ ộậ h ỗh ốg ổ Ọ ổ Ơ ỗhƯỒ ỗh Ọ gỌaố, ƠẬƠ ỗặốh Ơh ỗ ộậ đẬồ ứốg ớỘốg Ơủa h LTI Bài 3: BỌ ố đ Ọ ộậ ốg Ơ: ƠhỘy ố ỗặố hỌ Ộ ỗ ỐỌ ố ỗh Ọ gỌaố saốg ỐỌ ố ờ, tính Ơh ỗ Ơủa bỌ ố đ Ọ ộậ ƠhỘy ố ỗặố hỌ Ộ h Ộ ỗỷ ỗổẵố ỐỌ ố saốg ỐỌ ố ỗh Ọ gỌaố Bài 4: BỌ ố đ Ọ FỒỘổỌƯổ ổ Ọ ổ Ơ: ƠhỘy ố ỗặố hỌ Ộ ổ Ọ ổ Ơ ỗổẵố ỐỌ ố ỗh Ọ gỌaố saốg ỐỌ ố ỗ ố s , ơỂốg ƠẬƠ ỗhỘ ỗ ỗỒẬố FFT đ ớẬƠ địốh bỌ ố đ Ọ FỒỘổỌƯổ ổ Ọ ổ Ơ Bài 5: B ỏ Ơ s FIR: ỗhỌ ỗ k b ỏ Ơ FIR ỗhƯỒ yẵỘ Ơ Ộ ƠhỒ ỗổ Bài 6: B ỏ Ơ s IIR: ỗhỌ ỗ k b ỏ Ơ IIR ỗhƯỒ yẵỘ Ơ Ộ ƠhỒ ỗổ Ơ Ơ KI N TH C TI N Đ MẾố h Ơ Th Ơ hậốh X ỏý ỗặố hỌ Ộ s đẽỌ h Ọ sỌốh ộỌẵố Ơó ố ố ỗ ốg ộ X s YÊU C U MÔN H C Ng Ọ h Ơ ồh Ọ h Ơ đ y đủ ƠẬƠ bỘ Ọ ỏẵố ỏ ộậ ỏậỐ bậỌ ỗ đ y đủ CÁCH TI P NH N N I DUNG MÔN H C ỏý ỗặố hỌ Ộ TRANG > TH C HÀNH X LÝ TÍN HI U S Đ h Ơ ỗ ỗ ỐẾố ốậy, ốg Ọ h Ơ Ơ ố ỗh Ơ hậốh ỗhƯỒ h ỗổ Ơ bậỌ Ố Ọ ộậ ỗìỐ ỗhẵỐ ƠẬƠ ỗhẾốg ỗỌố ỏỌẵố ỔỘaố đ ố bậỌ h Ơ PH   ốg ố, ỏậỐ ƠẬƠ bậỌ ỗ ồ; đ Ơ NG PHÁP ĐÁNH GIÁ MÔN H C MẾố h Ơ đ Ơ đẬốh gỌẬ g Ố: ĐỌ Ố ỔỘẬ ỗổìốh: 30% Hìốh ỗhứƠ ộậ ố Ọ ơỘốg ơỒ gỌ ốg ộỌẵố ỔỘy ỗ địốh, ồhỂ h ộ Ọ ỔỘy Ơh đậỒ ỗ Ồ ộậ ỗìốh hìốh ỗh Ơ ỗ ỗ Ọ ố Ọ ỗ ƠhứƠ h Ơ ỗ ĐỌ Ố ỗhỌ: 70% Hìốh ỗhứƠ bậỌ ỗhỌ máy tính 60 phút BÀI 1: PH N M M MATLAB BÀI 1: PH N M M MATLAB    SaỘ khỌ h Ơ ớỒốg bậỌ ốậy, ốg Ọ h Ơ Ơó ỗh : Sử dụng phần mềm MATLAB Thực tạo script file hay function lưu trữ MATLAB Biết công cụ MATLAB 1.1 KH I Đ NG MATLAB MATLAB (MaỗổỌớ ỏabỒổaỗỒổy) ỏậ ồh ố Ố Ố ơỂốg đ gỌ Ọ ƠẬƠ bậỌ ỗỒẬố kỹ ỗhỘ ỗ, đ Ơ bỌ ỗ ỏậ ƠẬƠ bậỌ ỗỒẬố ỏỌẵố ỔỘaố đ ố Ốa ỗổ ố MATLAB ƠỘốg Ơ ƠẬƠ ỗỒỒỏbỒớƯs, ỗứƠ ƠẬƠ hậỐ Ố ổ ốg ỐẾỌ ỗổ ốg MATLAB đ gỌ Ọ ỔỘy ỗ ƠẬƠ ộ ố đ đ Ơ bỌ ỗ ốh ỏý ỗặố hỌ Ộ s , h ỗh ốg đỌ Ộ khỌ ố, Ố ốg ốƯỘổỒố, ưỘờờy ỏỒgỌƠ, ỐẾ ồh ốg ộ.ộ C a s bỌ Ộ ỗ ốg Ơủa Ơh ốg ỗổìốh MATLAB: Hình 1.1 - C a s kh Ọ đ ốg Ơủa MATLAB BÀI 1: PH N M M MATLAB C a s ỏậỐ ộỌ Ơ Ơủa MATLAB: Th Ố Ơ hỌ ố hậốh Workspace: Danh sách file bỌ ố ỗổỒốg b ốh hỌ ố Ơó ỗổỒốg ỗh Ố Ơ hỌ ố hậốh C a s ỏ ốh Nềỗ Sỗaổỗ: Ơhứa CẬƠ ỏ ốh đử toolbox ỗh Ơ hỌ ố Hình 1.2 – C a s ỏậỐ ộỌ Ơ Ơủa MATLAB  C as l ốh (Command window): Lậ Ơ a s gỌaỒ ỗỌ Ơhặốh Ơủa MATLAB b Ọ đừy ỏậ ố Ọ ốh gỌẬ ỗổị ƠẬƠ bỌ ố, hỌ ố ỗhị gỌẬ ỗổị, ỗặốh ỗỒẬố gỌẬ ỗổị Ơủa bỌ Ộ ỗhứƠ, ỗh Ơ ỗhỌ ƠẬƠ hậỐ Ơó s ố ỗổỒốg ỗh ơỒ ốg Ọ ơỂốg ỏ ỗổìốh ổa ỗổỒốg M-ưỌỏƯs CẬƠ ỏ ốh đ ỗhỌ ỏ ốh bằốg ồhặỐ EốỗƯổ Đ Ố Ơh ộỌ ố hỒ Ơ ƠẬƠ hậỐ Ơ ốh saỘ Ộ ốhắƠ ‘>>‘ th Ơ ốg ỗổìốh sỒ ố ỗh Ồ ỗổỒốg MATLAB, gỀ ỏ ốh: >>edit Hình 1.3 – C a s ƯơỌỗ đ sỒ ố sƠổỌồỗ ưỌỏƯ hay ưỘốƠỗỌỒố BÀI 1: PH N M M MATLAB SaỘ ốh ộậỒ đỒ ố Ơh ốg ỗổìốh saỘ: % Chuong trinh M-file x= 0:pi/6:2*pi; y=sin(x); plot(x,y); L Ộ Ơh ỗ ốg ỗổìốh ộ Ọ ỗẵố ưỌỏƯ plot_sin.m bằốg ƠẬƠh ốh ố Cỗổỏ+S hay ốh ố ộậỒ bỌ Ộ ốg SaộƯ GỌ Ọ ỗhặƠh đỒ ố Ơh % s khẾốg đ ốg ỗổìốh ỗổẵố: Dẽốg ỏậ ơẽốg Ơhề ỗhặƠh, ƠẬƠ ƠhỘ Ọ Ơ ơịƠh Dẽốg địốh ốghĩa ộƯƠỗỒổ ỗổỒốg khỒ ốg ỗ ồhặa saỘ Ộ đ ố 2 vậ Ố Ọ gỌẬ ỗổị ƠẬƠh ốhaỘ Ố ỗ khỒ ốg /6 Dẽốg gẬố bỌ ố y = sỌố(ớ) ộậ ơẽốg ộ đ ỗhị ỗổỒốg ỏậ ỗổ Ơ hỒậốh ộậ y ỏậ ỗổ Ơ ỗỘốg Hình 1.4 – L Ộ ưỌỏƯ ỗổỒốg Ơ a s EơỌỗ Th Ơ ỗhỌ Ơh ốg ỗổìốh ỗổẵố ỗổỒốg CỒỐỐaốơ ỚỌốơỒỚ bằốg ơẽốg ỏ ốh saỘ: >>plot_sin  C as CỒỐỐaốơ HỌstỒổỜ: CẬƠ ơẽốg đử ốh ỗổỒốg CỒỐỐaốơ ỚỌốơỒỚ (ƠẬƠ ơẽốg ốậy Ơó ỗh ỏậ ơẽốg ốh bỌ ố, Ơó ỗh ỏậ ơẽốg ỏ ốh) đ ỗa s Ơ gỌ ỏ Ọ ỗổỒốg Ơ a s CỒỐỐaốơ HỌsỗỒổy ộậ Ơ a s ốậy ƠhỒ ồhỰồ ốg ỏ Ọ ốh ốg ỏ ốh bằốg ƠẬƠh ốh đềồ ƠhỘ ỗ ỏẵố ƠẬƠ ỏ ốh hay bỌ ố  C as WỒổkspace: BÀI 5: B Bài 5.13 ThỌ ỗ k b ỏ Ơ ỗhẾốg ỗh ỗhƯỒ ồh ốh saỘ: L CS FIR 59 ốg ồhẬồ ỏ y Ố Ộ ỗ ố s ộ Ọ ƠẬƠ ỗhaỐ s p = 0.2; s = 0.3; Rp = 0.25 dB; As = 50 dB Ch ố đẬồ ứốg ớỘốg Ơó ƠhỌ Ộ ơậỌ 60 ứốg ộ Ọ 60 Ố Ộ ỗ ố s ỗhẾốg Ơó đ ổ ốg ỏậ 0.2 ỗ ốg đ ỗổỒốg khỒ ốg [0,2) D Ọ ốg ộ Ọ Ố Ộ ốh ố gỌẬ ỗổị GỌ s ỔỘẬ ỗổìốh ỗ Ọ Ộ hỒẬ Ơhỉ ổa ốẵố Ơh ố Ọ ƠhỘy ố ỗỌ Ố Ộ ốh ố ƠẬƠ gỌẬ ỗổị T1 = 0.5925 T2 = 0.1099 M Ộ ƠẬƠ ỗ ố s đ Ơ ƠhỒ ốh saỘ: H(k) = [1,1,1,1,1,1,1,T1,T2,0,0,…,0,0,T2,T1,1,1,1,1,1,1] (43 s 0) >>M = 60; alpha = (M-1)/2; L = 0:M-1; wl = (2*pi/M)*L; >>Hk = [ones(1,7),0.5925,0.1099,zeros(1,43), 0.1099, ones(1,6)]; % Dap ung tan so mau ly tuong >>Hdr = [1,1,0,0]; wdl = [0,0.2,0.3,1]; % Dap ung tan so ly tuong de bieu dien thi >>k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; >>angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)]; >>H = Hk.*exp(j*angH); >>h = real(ifft(H,M)); >>[db,mag,pha,grd,w] = freqz_m(h,1); >>[a,ww,L,Hr] = FIR_t2(h); a ChỘ Ọ Ố Ộ ỗ ố s : >>subplot(221); plot(wl(1:31)/pi,Hk(1:31),'o',wdl,Hdr); >>axis([0,1,-0.1,1.1]); >>title('Cac mau tan so: M=60, T2 = 0.5925, T1 = 0.1099'); >>xlabel(f[*\pi]'); ylabel('Hr(k)'); b ĐẬồ ứốg ớỘốg Ơủa b ỏ Ơ ỗh Ơ ỗ : >>subplot(222); stem(L,h); 0.5925, 60 BÀI 5: B L CS FIR >>axis([-1,M,-0.1,0.3]); >>title('Dap ung xung'); >>xlabel('n'); ylabel('h(n)'); c BỌẵố đ Ơủa đẬồ ứốg ỗ ố s : >>subplot(223); plot(ww/pi,Hr,wl(1:31)/pi,Hk(1:31),'o'); >>axis([0,1,-0.2,1.2]); >>title('Bien cua dap ung tan so'); >>xlabel('f[*\pi]'); ylabel('Hr(w)'); d BỌẵố đ Ơủa đẬồ ứốg ỗ ố s theo dB: >>subplot(224); plot(w/pi,db); >>axis([0,1,-100,10]); grid >>title('Bien cua dap ung tan so '); >>xlabel('f[*\pi]'); ylabel('dB'); Bài 5.14 ThỌ ỗ k b ốh saỘ: ỏ Ơ ỗhẾốg Ọ ỗhƯỒ ồh ốg ồhẬồ ỏ y Ố Ộ ỗ ố s ộ Ọ ƠẬƠ ỗhaỐ s s1 = 0.2; p1 = 0.35; p2 = 0.65; s2 = 0.8; Rp = dB; As = 60 dB Ch ố đẬồ ứốg ớỘốg Ơó ƠhỌ Ộ ơậỌ 40 ứốg ộ Ọ 40 Ố Ộ ỗ ố s ỗhẾốg Ơó đ ổ ốg ỏậ 0.3 ỗ ốg đ ỗổỒốg khỒ ốg [0,2) D Ọ ốg ộ Ọ Ố Ộ ốh ố gỌẬ ỗổị GỌ s ỔỘẬ ỗổìốh ỗ Ọ Ộ hỒẬ Ơhỉ ổa ốẵố Ơh ố Ọ ƠhỘy ố ỗỌ Ố Ộ ốh ố ƠẬƠ gỌẬ ỗổị T1 = 109021 T2 = 0.59417456 M Ộ ƠẬƠ ỗ ố s đ Ơ ƠhỒ ốh saỘ: H(k) = [0,0,0,0,0,T1,T2,1,1,1,1,1,1,1,T2,T1,0,…,0,T1,T2,1,1,1,1,1,1,1,T2,T1,0,0,0,0] (9 s 0) QỘẬ ỗổìốh ỗh Ơ hỌ ố ỗ 5.4 PH ốg ỗ bậỌ 5.13 NG PHÁP L P ĐẬồ ứốg ỗ ố s Ơủa ỏỒ Ọ b ỏ Ơ FIR: H() = P()Q() HậỐ saỌ s gỌ a b ỏ Ơ ỗh Ơ ỗ ộậ b ỏ Ơ ỏý ỗ ốg: BÀI 5: B Ta địốh ốghĩa: � = � [ ̂ � = KhỌ đó: � ]= � − � � � = ̂ � [̂ � − Paổks ộậ MƠCỏƯỏỏaố đử đ a ổa gỌ Ọ ồhẬồ s ớỘốg Ơủa b ỏ Ơ ỗ Ọ FIR 61 � ] � � ̂ � = � � − � � [ L CS � ] ốg ỗhỘ ỗ ỗỒẬố RƯỐƯờ đ ỗìỐ ổa đẬồ ứốg Ộ, ỗứƠ ỏậ g ố đềốg ỗhƯỒ ốghĩa ChƯbyshƯộ đ Ọ ộ Ọ Ố ỗ b ỏ Ơ ỏý ỗ ƠhỒ gỌẬ ỗổị M ỏậ ƠhỌ Ộ ơậỌ Ơủa ƠhỘ Ọ đẬồ ứốg ớỘốg ộ Ọ ƠẬƠ đỌ Ộ kỌ ố ổậốg bỘ Ơ ộ sóốg Ọ ỗhẾốg ộậ Ọ Ơhắố ốh đ ốg, g ố saỘ: XẬƠ địốh ỏỒ Ọ b ỏ Ơ, ỗặốh gỌẬ ỗổị R ộậ ớừy ốg ƠẬƠ hậỐ W(ω), Q(ω) LỒ Ọ b l Ơ R − FIR ỏỒ Ọ FIR ỏỒ Ọ FIR ỏỒ Ọ FIR ỏỒ Ọ P() − � ∑̅ � ∑ ̅ � ∑ ̅ � = − − ∑̅ = − Q() = = cos(/2) sin sin(/2) Xừy ốg bậỌ ỗỒẬố g ố đềốg bằốg ƠẬƠh ớẬƠ địốh ƠẬƠ hậỐ ̂ � , S ốg ỗhỘ ỗ ỗỒẬố ỗổaỒ đ Ọ RƯỐƯờ đ ỗìỐ ổa hậỐ ỗ Ọ  Ch ố ỏ y R+2 đỌ Ố ổ Ọ ổ Ơ, gỌ s  Tổẵố Ơ s Ộ P(ω) � ỏậ ƠẬƠ Ơ Ơ ỗổị Ơủa hậỐ saỌ s ỗ Ọ R+2 đỌ Ố ổ Ọ ổ Ơ ốóỌ ỗổẵố, hậỐ E(ω) ỏỘừố ồhỌẵố đ Ọ Ộ ộậ Ơó ỗổị ỗỘy ỗ đ Ọ bằốg Ố ỗ gỌẬ ỗổị δ ốậỒ đó, ỗặốh ố Ọ sỘy ỏ Ọ gỌẬ ỗổị δ ộậ hậỐ P(ω), ỗ hậỐ saỌ s E(ω), ỗặốh đ Ơ Ơ Ơ ỗổị ỗh Ơ Ơủa hậỐ saỌ s ỗặốh ổa 62 BÀI 5: B L CS FIR  XƯỐ ớỰỗ ớƯỐ ƠẬƠ gỌẬ ỗổị ổ Ọ ổ Ơ đ Ơ Ơh ố baố đ Ộ Ơó ỗh Ơ s ỏậ ƠẬƠ đỌ Ố Ốậ hậỐ saỌ s E(ω) đ ỗ Ơ Ơ ỗổị ộậ Ơó ỗổị ỗỘy ỗ đ Ọ bằốg ốhaỘ hay khẾốg N Ộ khẾốg, ỗìỐ ƠẬƠ đỌ Ố ỗ Ọ E(ω) đ ỗ Ơ Ơ ỗổị  TổỒốg ƠẬƠ đỌ Ố Ơ Ơ ỗổị Ơủa E(ω) ỏ y ổa R+2 đỌ Ố ộậ ỔỘay ộ ỏ ỏ Ọ ỗ  L ỏ Ọ ƠẬƠ b  T b Ơ Ơ 2, 3, ộậ ƠhỒ đ ố khỌ ỗ h ƠẬƠ đỌ Ố ổ Ọ ổ Ơ h Ọ ỗ ỗ ƠẬƠ đỌ Ố ổ Ọ ổ Ơ ƠỘ Ọ ƠỂốg, ỗặốh ổa hậỐ P(ω), ỗ ỗặốh ổa ƠẬƠ h s Ơủa P(ω) Tặốh ƠẬƠ gỌẬ ỗổị Ơủa ƠhỘ i đẬồ ứng xung h(n) KhỌ Ơh ố gỌẬ ỗổị M Ơậốg ƠhỘ ố ỗhì k ỗ ỔỘ ỏậ ỗhỘ đ g ố ộ Ọ yẵỘ Ơ Ộ bậỌ ỗỒẬố N Ộ ốh Ơ b ỏ Ơ Ơó hậỐ đẬồ ứốg ỗ ố s Ơậốg ộ Ọ gỌẬ ỗổị M ốậỒ Ốậ Ơh a ỗhỒ Ốửố đ Ơ yẵỘ Ơ Ộ ỗhì ồh Ọ ỗăốg gỌẬ ỗổị M đ ố khỌ ốậỒ ỗhỒ Ốửố ƠẬƠ đỌ Ộ kỌ ố ổậốg bỘ Ơ ƠhỒ δp ộậ δs (hay As Rp) M ỗ ƠẾốg ỗhứƠ ỏ a Ơh ố baố đ Ộ ƠhỒ ƠhỌ Ộ ơậỌ M Ơủa đẬồ ứốg ớỘốg ỏậ: = − √� � − Δ ớ� Δ = Trong MATLAB, ỗìỐ đẬồ ứốg ớỘốg Ơủa b ỏ Ơ ỗ Ọ ỏý ỗ ốg ƠhỒ ỗổ Ơđ � − �� � Ộ ộ Ọ gỌẬ ỗổị M ộậ hậỐ đẬồ ứốg ỗ ố s Ơ ỗh Ơ hỌ ố b Ọ hậỐ firpm Bài 5.15 T Ồ sƠổỌồỗ ưỌỏƯ saỘ đ bỌ Ộ ơỌ ố b ỏ Ơ ỗổẵố đ ỗhị: wp = 0.2*pi; ws =0.3*pi; Rp = 0.25; As = 50; delta_w = 2*pi/1000; wsi = ws/delta_w+1; delta1 = (10^(Rp/20)-1)/(10^(Rp/20)+1); delta2 = (1+delta1)*(10^(-As/20)); deltaH = max(delta1,delta2); deltaL = min(delta1,delta2); weights = [delta2/delta1 1]; deltaf = (ws-wp)/(2*pi); M = ceil((-20*log10(sqrt(delta1*delta2))-13)/(14.6*deltaf)+1) f = [0 wp/pi ws/pi 1]; m = [1 0]; BÀI 5: B h = firpm(M-1,f,m,weights); [db,mag,pha,grd,w] = freqz_m(h,[1]); Asd = -max(db(wsi:1:501)) while Asd log( + √ − ) =⌈ ⌉ log � + √� − � − =√ � Bài 6.2 T Ồ hậỐ u_chb1ap ốhằỐ ỗổ ộ � = h � �� s Ơủa ƠẬƠ đa ỗhứƠ ỗ s ộậ đa ỗhứƠ Ố Ộ s Ơủa hậỐ ỗổỘy ố đ ỗ ƠhỒ ỗhỌ ỗ k b ỏ Ơ ốg ChƯbyshƯộ I Ơó ỗ ố s Ơắỗ ỗỘỳ ý: function [b,a] = u_chb1ap(N,Rp,Omegac) % Bo loc thong thap dang Chebyshev-1 % tan so cat khong duoc chuan hoa 68 BÀI 6: B L CS IIR % [b,a] = u_chb1ap(N,Rp,Omegac) % b = cac he so da thuc tu so cua Ha(s) % a = cac he so da thuc mau so cua Ha(s) % N = Bac cua bo loc Chebyshev-I % Rp = Do gon dai thong theo don vi dB; Rp > % Omegac = tan so cat theo don vi radians/sec [z,p,k] = cheb1ap(N,Rp); a = real(poly(p)); aNn = a(N+1); p = p*Omegac; a = real(poly(p)); aNu = a(N+1); k = k*aNu/aNn; B = real(poly(z)); b0 = k; b = k*B; HậỐ s ỐẾ ỗ ỗổẵố ỗổ ộ hậỐ ỗổỘy ố ộ Ọ b Ơ N ƠhỒ ỗổ Ơh ố ƠhỒ ồhỂ h ỗ Ọ Ộ ộ Ọ ƠẬƠ Ơhỉ ỗỌẵỘ kỹ ỗhỘ ỗ yẵỘ Ơ Ộ đ Ộ ộậỒ Bài 6.3 T Ồ hậỐ afd_chb1 ỗổ ChƯbyshƯộ Ơó b Ơ ỗ Ọ Ơ B Ơ Ơủa b ỏ Ơ Ơó ỗh ỏ a ộ ỗhỌ ỗ k b ỏ Ơ ỗhẾốg ỗh ỗ ốg ỗ , địốh ốg Ộ: function [b,a] = afd_chb1(Wp,Ws,Rp,As) % Analog Lowpass Filter Design: Chebyshev-1 % [b,a] = afd_chb1(Wp,Ws,Rp,As) % b = cac he so da thuc tu so cua Ha(s) % a = cac he so da thuc mau so cua Ha(s) % Wp = tan so cat dai thong theo don vi rad/sec; Wp >0 % Ws = tan so cat dai chan theo don vi rad/sec; Ws>Wp >0 % Rp = Do gon dai thong theo don vi dB; (Rp > 0) % As = Do suy giam dai chan theo don vi +dB; (Ap > 0) if Wp