1 đào tạo Bộ giáo dục TRờng Đại học Bách khoa Hà nội *** Ngun Xu©n Trêng Xư lý tÝn hiƯu ¸p dơng ph¸t sè Ln ¸n tiến sỹ kỹ thuật Hà nội-2009 Bộ giáo dục đào tạo TRờng Đại học Bách khoa Hà nội *** Nguyễn Xuân Trờng Xử lý tín hiệu áp dụng phát số Chuyên ngành: Kỹ thuật mật mà M· sè: 62.52.73.01 Ln ¸n tiÕn sü kü tht Híng dÉn khoa häc: 1- PGS.TS NguyÔn Quèc Trung 2- TS Hà nội-2009 Nguyễn Phơng Li cam oan Tụi xin cam đoan cơng trình nghiên cứu khoa học riêng tơi Các số liệu, hình vẽ, biểu bảng, kết tính tốn luận án trung thực, chưa công bố đâu cơng trình khác Luận án tham khảo thư viện trường Cao đẳng Phát thanh-Truyền hình I, thư viện Đại học Bách khoa Hà nội, thư viện Quốc gia Hà nội Website: http://cdptth1.edu.vn Nghiên cứu sinh Nguyễn Xuân Trường Lời cảm ơn Sau bốn năm học tập, nghiên cứu đến hồn thành luận án tiến sỹ Để có kết này, bên cạnh nỗ lực cố gắng không ngừng thân, cịn có quan tâm tạo điều kiện Giáo sư, phó Giáo sư, Tiến sĩ, thầy giáo cô giáo, nhà khoa học giúp đỡ thực luận án Tôi xin chân thành cảm ơn Viện Đào tạo sau Đại học, trường Đại học Bách khoa Hà nội tạo điều kiện giúp đỡ tơi hồn thành bước nghiên cứu thực luận án nghiên cứu sinh Tôi xin chân thành cảm ơn PGS.TS Nguyễn Quốc Trung TS Nguyễn Phương, người thầy giáo định hướng tận tình bảo, hỗ trợ tài liệu, giúp đỡ tơi q trình nghiên cứu triển khai đề tài Tôi xin chân thành cảm ơn thầy giáo, cô giáo khoa Điện tử-Viễn thông, trường Đại học Bách khoa Hà nội tạo điều kiện cho tham gia buổi hội thảo, đóng góp ý kiến quý giá cho hình thức nội dung luận án Tôi xin chân thành cảm ơn Lãnh đạo Đài Tiếng nói Việt nam, Ban giám hiệu trường Cao đẳng Phát thanh-Truyền hình I tạo điều kiện thuận lợi cho nghiên cứu thực luận án Tơi xin chân thành cảm ơn Giáo sư, Phó Giáo sư, Tiến sỹ phản biện báo tơi, đóng góp nhiều ý kiến giá trị nội dung cách đặt vấn đề khoa học Tôi xin chân thành cảm ơn Ban tổ chức Hội nghị khoa học Hội Vô tuyến-Điện tử Việt nam (REV), trường Đại học Bách khoa Hà nội, trường Đại học Tự nhiên-Đại học Quốc gia Hà nội, tạo điều kiện cho tơi trình bày báo cáo khoa học Tơi xin chân thành cảm ơn thầy giáo, cô giáo Học viện kỹ thuật Quân sự, Đại học Cơng nghệ, Học viện Bưu chính-Viễn thơng, trường Cao đẳng Phát thanhTruyền hình I, đồng chí cán kỹ thuật Đài Tiếng nói Việt nam giúp đỡ việc thực luận án tiến sỹ Xin trân trọng cảm ơn ! Nghiên cứu sinh Nguyễn Xuân Trường Mục lục Trang Lời cam đoan …………………………………………………………………………… i Lời cảm ơn……………………………………………………………………………… ii Mục lục…………………………………………………………………………………… iii Danh mục ký hiệu ………………………………………………………………… v Danh mục chữ viết tắt ……………………………………………………………… vii Danh mục bảng liệu…………… ……………………………………………… ix Danh mục hình vẽ, đồ thị…………………………………………………………… x Chương mở đầu 1- Giới thiệu …………………………………………………………………… 2- Tính cấp thiết đề tài……………………………………………………… 3- Mục đích nghiên cứu………………………………………………………… 4- Đối tượng phạm vi nghiên cứu…………………………………………… 5- Phương pháp nghiên cứu…………………………………………………… 6- Tình hình nghiên cứu nước giới…………………………… 7- Ý nghĩa khoa học thực tiễn đề tài…………………………………… 8- Nội dung luận án………………………………………………………… 1 5 Chương 1: Mã hoá dải 10 1.1- Tổng quan mã hoá dải con……………………………………………………… 10 1.1.1- Giới thiệu………………………………………………………………… 10 1.1.2- Nguyên lý mã hoá dải (SBC)………………………………………… 12 1.1.3- Các thông số kỹ thuật SBC………………………………… 15 1.1.4- Các điều kiện mã hoá dải con………………………………………… 18 1.1.5- Điều kiện để tăng hệ số nén tín hiệu SBC………………………… 21 1.2- Bank lọc theo phân tích Wavelet…………………………………………………… 25 1.2.1- Bank lọc nhiều nhịp đơn phân giải……………………………………… 25 1.2.2- Bank lọc nhiều nhịp đa phân giải………………………………………… 27 1.3- Mã hoá dải ứng dụng xử lý âm số……………………………… 30 1.3.1- Chuẩn âm số MPEG/audio……………………………… 30 1.3.2- Thuật toán ATRAC……………………………………………………… 32 1.3.3- Thuật toán PASC………………………………………………………… 34 1.3.4- Bộ mã hoá dải APT-X………………………………………………… 35 1.4- Kết luận……………………………………………………………………………… 35 Chương 2: Mã hoá giải đa phân giải tuyệt đối 37 2.1- Giới thiệu ………………………………………………………………… 2.2- Đề xuất phương pháp xác định tổ hợp phân chia…………………………………… 2.2.1- Phương pháp xác định hệ số phân chia…………………………… 2.2.2- Xác định tổ hợp phân chia…………………………………………… 2.2.3- Xác định tổ hợp phân chia tối ưu………………………………………… 2.3- Thiết kế lọc FIR phương pháp dải chuyển tiếp-cửa sổ…………………… 37 37 37 41 42 44 2.3.1- Giới thiệu………………………………………………………………… 2.3.2- Thiết kế lọc FIR phương pháp trực tiếp………………………… 2.3.3- Thiết kế lọc FIR phương pháp dải chuyển tiếp………………… 2.3.4- Thiết kế lọc FIR phương pháp dải chuyển tiếp-cửa sổ………… 2.4- Kết luận……………………………………………………………………………… 44 46 49 58 62 Chương 3: Mã hoá dải đa phân dải tuyệt đối kênh SBC(632) kênh SBC(12642) 3.1- Giới thiệu…………………………………………………………………………… 3.2- Mã hoá dải với tổ hợp phân chia [6, 3, 2]……………………………………… 3.2.1- Sơ đồ khối SBC(632)………………………………………………… 3.2.2- Thiết kế bank lọc SBC(632)………………………………………… 3.2.3- Giải pháp biến đổi hệ số……………………………………… 3.2.4- Cấp pháp bit……………………………………………………… …… 3.2.5- So sánh SBC(632) với SBC(442)………………………………………… 3.3- Mã hoá dải với tổ hợp phân chia [12, 6, 4, 2]………………………………… 3.3.1- Sơ đồ khối SBC(12642)……………………………………………… 3.3.2- Thiết kế lọc SBC(12642)…………………………………… 3.3.3- Giải pháp biến đổi hệ số………………………………………………… 3.3.4- Cấp phát bit………………………………………………………………… 3.3.5- So sánh SBC(12642) với SBC(8842)…………………………………… 3.4- Kết luận……………………………………………………………………………… 63 63 63 63 64 72 78 79 85 85 86 96 99 100 106 Chương 4: Ứng dụng SBC(632) SBC(12642) phát số 107 4.1- Giới thiệu hệ thống phát số nay………………………………………… 107 4.1.1- Hệ thống phát số dải kênh………………………… 108 4.1.2- Hệ thống phát số EUREKA-147…………………………………… 110 4.2- Ứng dụng SBC(632) SBC(12642) phát số………………… 113 4.2.1- Ứng dụng SBC(632) SBC(12642) phát số qua Internet……………………………………………………………… 113 4.2.2- Ứng dụng SBC(632) SBC(12642) hệ thống phát số IBOC… 114 4.2.3- Ứng dụng SBC(632) SBC(12642) hệ thống phát số EUREKA-147……………………………………………… 115 4.3- Định hướng phát số vô tuyến Việt nam……………………………….…… 117 4.4- Kết luận…………………………………………………………………… 118 Kết luận kiến nghị 119 1- Kết luận………………………………………………………………………… 119 2- Kiến nghị……………………………… ……………………………………… 120 Danh mục báo tác giả liên quan tới luận án công bố…… 121 Danh mục tài liệu tham khảo…………………………………………………………… 122 Phụ lục …………………………………………………………………………………… 124 -*** - Danh mục ký hiệu 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 N TS TSv TSr p BL BB BH Btr ∆f ∆f B δp δs BOFDM BAM-IBOC BFM-IBOC G GdB GL G LdB GB G BdB GH G HdB H L(ejF) H B(ejF) H H(ejF ) kn Lb n Hbn Lbkn Bbkn Hbkn D I ni β3 β4 ∆f OFDM ε εTB ε(632) ε(442) Bậc lọc Chu kỳ lấy mẫu (đơn vị giây) Chu kỳ lấy mẫu tín hiệu đầu vào Chu kỳ lấy mẫu tín hiệu đầu Cấp phân chia phân tích Wavelet Độ rộng dải thơng lọc thông thấp (đơn vị rad) Độ rộng dải thông lọc thông dải (đơn vị rad) Độ rộng dải thông lọc thông cao (đơn vị rad) Độ rộng dải chuyển tiếp lọc (đơn vị rad) Độ rộng dải tần tín hiệu (đơn vị Hz) Độ rộng dải tới hạn tai người (đơn vị Hz) Độ gợn sóng dải thơng lọc Độ gợn sóng dải chắn lọc Độ rộng dải tần tín hiệu OFDM Độ rộng dải tần tín hiệu âm số hệ thống AM-IBOC Độ rộng dải tần tín hiệu âm số hệ thống FM-IBOC Đáp ứng biên độ lọc Đáp ứng biên độ lọc tính theo dB Đáp ứng biên độ lọc thông thấp Đáp ứng biên độ lọc thơng thấp tính theo dB Đáp ứng biên độ lọc thông dải Đáp ứng biên độ lọc thông dải tính theo dB Đáp ứng biên độ lọc thông cao Đáp ứng biên độ lọc cao tính theo dB Đáp ứng tần số lọc thông thấp theo thang tần số F Đáp ứng tần số lọc thông dải theo thang tần số F Đáp ứng tần số lọc thông cao theo thang tần số F Hàm cửa sổ Hệ số lọc thông dải phụ thuộc vào n Hệ số lọc thông cao phụ thuộc vào n Hệ số lọc thông thấp nhân với cửa sổ kn phụ thuộc vào n Hệ số lọc thông dải nhân với cửa sổ kn phụ thuộc vào n Hệ số lọc thông cao nhân với cửa sổ kn phụ thuộc vào n Hệ số phân chia phân chia Hệ số nội suy nội suy Hệ số phân chia dải thứ i Hệ số so sánh lỗi mã hóa dải kênh Hệ số so sánh lỗi mã hóa dải kênh Khoảng cách tần số hai sóng mang liền kề kỹ thuật OFDM Lỗi khơi phục mã hóa dải Lỗi khơi phục trung bình mã hóa dải Lỗi khơi phục trung bình mã hóa dải SBC(632) Lỗi khơi phục trung bình mã hóa dải SBC(442) 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 ε(12642) ε(8842) σ i2 bTB bTB(632) bTB(442) bTB(12642) bTB(8842) bi bTBi M N OFDM m f fmax Fmax F w ω fNy fS fc ω cH ω cB1 ω cB2 ω sL ω sH ω sB1 ω sB2 fSv fSr fB R(632) R(12642) RTB ηTB η(632) η(442) η(12642) η(8842) x(n) X(z) X(ejω) TOFDM Lỗi khơi phục trung bình mã hóa dải SBC(12642) Lỗi khơi phục trung bình mã hóa dải SBC(8842) Phương sai tín hiệu dải thứ i Số bít trung bình/mẫu mã hố dải Số bít trung bình/mẫu mã hố dải SBC(632) Số bít trung bình/mẫu mã hố dải SBC(442) Số bít trung bình/mẫu mã hố dải SBC(12642) Số bít trung bình/mẫu mã hố dải SBC(8842) Số bít cấp cho dải thứ i Số bít trung bình/mẫu cấp cho tín hiệu dải thứ i Số kênh mã hóa dải Số sóng mang kỹ thuật OFDM Số trạng thái điều chế số m-PSK Tần số (đơn vị Hz) Tần số cực đại tín hiệu âm (đơn vị Hz) Tần số cực đại tín hiệu âm chuẩn hố ω/π (khơn thứ ngun) Tần số chuẩn hố ω/π (khơng có thứ ngun) Tần số góc w=2πf (đơn vị rad/s) Tần số góc chuẩn hóa (đơn vị rad) Tần số Nyquist (đơn vị Hz) Tần số lấy mẫu (đơn vị Hz) Tần số cắt lọc (đơn vị Hz) Tần số cắt lọc thông cao (đơn vị rad) Tần số cắt lọc thông dải (đơn vị rad) Tần số cắt lọc thông dải Tần số giới hạn lọc thông thấp mức -20dB (đơn vị rad) Tần số giới hạn lọc thông cao mức -20dB (đơn vị rad) Tần số giới hạn lọc thông dải mức -20dB (đơn vị rad) Tần số giới hạn lọc thông dải mức -20dB (đơn vị rad) Tần số lấy mẫu tín hiệu đầu vào Tần số lấy mẫu tín hiệu đầu Tần số cảm giác (còn gọi tần số Bark) tai người (đơn vị Bark) Tốc độ bít SBC(632) Tốc độ bít SBC(12642) Tốc độ bít trung bình (đơn vị bps) Tỷ lệ nén liệu trung bình Tỷ lệ nén liệu mã hố dải SBC(632) Tỷ lệ nén liệu mã hoá dải SBC(442) Tỷ lệ nén liệu mã hoá dải SBC(12642) Tỷ lệ nén liệu mã hố dải SBC(8842) Tín hiệu miền thời gian rời rạc chuẩn hóa theo tần số lấy mẫu Tín hiệu miền z Tín hiệu miền tần số liên tục Thời gian tồn sóng mang kỹ thuật OFDM -*** - Danh mục chữ viết tắt 1- ADM 2- ADPCM 3- AM-IBOC 4- APT-X 5- ASPEC 6- ATRAC 7- BD-ROM 8- BFU 9- BPF 10- CD 11- CD-ROM 12- CELP 13- CVSD 14- DAB 15- DAB-T 16- DAT 17- DCC 18- DM 19- DVB-T 20- DVD 21- FFT 22- FM-IBOC 23- HPF 24- IBOC 25- IEC 26- ISDN 27- ISO 28- ITU 29- LM 30- LPAC 31- LPC 32- LPCM 33- LPF 34- LTAC 35- MD 36- MDCT 37- MPEG/audio 38- MPEG-1/audio 39- MPEG-4 Adaptive Delta Modulation Adaptive Differential Pulse Code Modulation Amplitude Modulation-In Band On Channel Audio Processing Technology-X Adaptive Spectral Perceptual Entropy Coding Adaptive Transform Acoustic Coding Blu-ray Read Only Memory Block Floating Unit Band Pass Filter Compact Disc Compact Disc-Read Only Memory Code Excited Linear Prediction Continuously Variable Slope Delta modulation Digital Audio Broadcasting Digital Audio Broadcasting-Terrestrial Digital Audio Tape Digital Compact Cassette Delta Modulation Digital Video Broadcasting-Terrestrial Digital Versatile Disc Fast Fourier Transform Frequency Modulation-In Band On Channel High Pass Filter In Band On Channel International Electrotechnical Commission Integrated Servises Digital Network International Organization for Standardization International Telecommunication Union Long Mode Lossless Predictive Audio Compression Linear Predictive Coding Linear Pulse Code Modulation Low Pass Filter Lossless Transform Audio Compression MiniDisc Modified Discrete Cosine Transform Motion Picture Experts Group/audio Motion Picture Experts Group-1/audio Motion Picture Experts Group-4 10 40- MUSICAM 41- OFDM 42- PAM 43- PASC 44- PC 45- PCM 46- PSD 47- PSK 48- QMF 49- QPSK 50- SBC 51- SM 52- SMR 53- TAK 54- AC-3 55- BST-OFDM-ISDB 56- EUREKA-147 57- HIPERLAN 58- IPL/VOA 59- Maskee 60- Masker 61- MP1 62- MP2 63- MP3 64- Non-tone 65- SBC(12642) 66- SBC(442) 67- SBC(632) 68- SBC(8842) 69- Skywave 2000 70- T2M 71- Tone Masking-pattern Universal Subband Integrated Coding And Multiplexing Orthogonal Frequency Division Multiplexing Pulse Amplitude Modulation Precision Adaptive Subband Coding Personal Coputer Pulse Code Modulation Power Spectral Density Phase Shift Key Quadrature Mirror Filter Quadrature Phase Shift Key Sub-Band Coding Short Mode Signal-to-Mask Ratio Tom’s verlustfreier AudioKompressor (một định dạng sound) AC3 (hay Dolby Digital) Hệ thống phát số Nhật Hệ thống phát số châu Âu Mạng tốc độ cao Hệ thống phát số Một âm đơn bị che thành phần âm đơn khác Một âm đơn có khả che thành phần âm đơn khác Gọi tắt chuẩn MPEG-1/audio lớp I Gọi tắt chuẩn MPEG-1/audio lớp II Gọi tắt chuẩn MPEG-1/audio lớp III Các thành phần âm đơn Mã hóa dải kênh dùng tổ hợp phân chia [12, 6, 4, 2] Mã hóa dải kênh dùng tổ hợp phân chia [4, 4, 2] Mã hóa dải kênh dùng tổ hợp phân chia [6, 3, 2] Mã hóa dải kênh dùng tổ hợp phân chia [8, 8, 4, 2] Hệ thống phát số Hệ thống phát số (Telefunken-Multicast) Một thành phần âm đơn -*** - 161 wsH=(wcH*tan(anpha)-1)/tan(anpha); Hbk0=(1+2*BH*tan(anpha))/(2*pi*tan(anpha)); Hb1=(tan(anpha)*(cos(n*wcH)-cos(n*wsH)))/pi; Hb2=Hb1./n; Hbn=Hb2./n; Hbkn=Hbn.*kn; HPF=[Hbk0 Hbkn]; %************************BANK LOC SBC(8842)*************************** %****************Thiet ke bo loc RLPF**************** RwcL=pi/8; RBL=RwcL; RwsL=(RwcL*tan(anpha)+1)/tan(anpha); RLbk0=(1+2*RBL*tan(anpha))/(2*pi*tan(anpha)); RLb1=(tan(anpha)*(cos(n*RwcL)-cos(n*RwsL)))/pi; RLb2=RLb1./n; RLbn=RLb2./n; RLbkn=RLbn.*kn; RLPF=[RLbk0 RLbkn]; %*************Thiet ke bo loc RBPF1******************* RwcB11=pi/8; RwcB21=pi/4; RBB1=RwcB21-RwcB11; RwsB11=((RwcB11*tan(anpha))-1)/tan(anpha); RwsB21=((RwcB21*tan(anpha))+1)/tan(anpha); RBbk01=(1+RBB1*tan(anpha))/(pi*tan(anpha)); RBb11=(tan(anpha)*(cos(n*RwcB21)+cos(n*RwcB11)-cos(n*RwsB21)-cos(n*RwsB11)))/pi; RBb21=RBb11./n; RBbn1=RBb21./n; RBbkn1=RBbn1.*kn; RBPF1=[RBbk01 RBbkn1]; %*************Thiet ke bo loc RBPF2******************* RwcB12=pi/4; RwcB22=pi/2; RBB2=RwcB22-RwcB12; RwsB12=((RwcB12*tan(anpha))-1)/tan(anpha); RwsB22=((RwcB22*tan(anpha))+1)/tan(anpha); RBbk02=(1+RBB2*tan(anpha))/(pi*tan(anpha)); RBb12=(tan(anpha)*(cos(n*RwcB22)+cos(n*RwcB12)-cos(n*RwsB22)-cos(n*RwsB12)))/pi; RBb22=RBb12./n; RBbn2=RBb22./n; RBbkn2=RBbn2.*kn; RBPF2=[RBbk02 RBbkn2]; %*************Thiet ke bo loc RHPF******************* RwcH=pi/2; RBH=pi-RwcH; RwsH=(RwcH*tan(anpha)-1)/tan(anpha); RHbk0=(1+2*RBH*tan(anpha))/(2*pi*tan(anpha)); RHb1=(tan(anpha)*(cos(n*RwcH)-cos(n*RwsH)))/pi; RHb2=RHb1./n; RHbn=RHb2./n; RHbkn=RHbn.*kn; RHPF=[RHbk0 RHbkn]; %**********************XU LY AM THANH QUA SBC(12642)******************** y0=filter(LPF,1,audio); y1=filter(BPF1,1,audio); y2=filter(BPF2,1,audio); y3=filter(HPF,1,audio); Py0=downsample(y0,12); Py1=downsample(y1,4); Py2=downsample(y2,2); 162 Py3=downsample(y3,2); varK0=var(Py0); varK1=var(Py1); varK2=var(Py2); varK3=var(Py3); %**********************XU LY AM THANH QUA SBC(8842)******************* Ry0=filter(RLPF,1,audio); Ry1=filter(RBPF1,1,audio); Ry2=filter(RBPF2,1,audio); Ry3=filter(RHPF,1,audio); RPy0=downsample(Ry0,8); RPy1=downsample(Ry1,8); RPy2=downsample(Ry2,4); RPy3=downsample(Ry3,2); RvarK0=var(RPy0); RvarK1=var(RPy1); RvarK2=var(RPy2); RvarK3=var(RPy3); %*****************TINH bTB8842 KHI bTB12642=8 bit/mau***************** bTB12642=8; E12642=4*0.5*pi*(sqrt(3))*(2^(2*bTB12642))*(varK0^(1/12))*(varK1^(1/4))*(varK2^(1/2))*(varK3^(1/2)) TG=E12642/(4*0.5*pi*(sqrt(3))*(RvarK0^(1/8))*(RvarK1^(1/8))*(RvarK2^(1/4))*(RvarK3^(1/2))); bTB8842=(log2(TG))/(-2) end; %========================KET THUC================================ Phụ lục P3.23 Bảng so sánh b TB SBC(12642) SBC(8842) lỗi khôi phục ε TT bTB(12642) b TB(8842) ε TT b TB(12642) bTB(8842) (x10 -6) 10.9704 0.000019 51 10.2342 10.1146 0.002134 52 10.5162 9.7926 0.010726 53 9.9797 10.1148 0.002138 54 9.9237 10.2211 0.001350 55 10.6345 10.4475 0.000379 56 10.2701 10.6560 0.000122 57 10.3283 8 10.3876 0.000498 58 10.7627 10.8320 0.000044 59 9.9281 10 10.0422 0.002584 60 10.2916 11 10.2930 0.000881 61 11.0886 12 13 14 15 16 17 18 19 20 8 8 8 8 11 2065 10.1918 10.9978 10.1167 10.3826 10.0629 9.8684 10.7509 10.7207 000006 0.000529 0.000019 0.002076 0.000563 0.002767 0.006547 0.000059 0.000076 62 63 64 65 66 67 68 69 70 8 8 8 8 10 0320 10.2309 9.9495 9.7479 9.8765 9.9621 10.1719 10.0509 9.9193 ε (x10-6) 0.001245 0.000274 0.005203 0.006313 0.000137 0.001005 0.000801 0.000064 0.005441 0.000958 0.000011 003399 0.001160 0.004803 0.017315 0.007564 0.005244 0.001622 0.001181 0.007403 163 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 Phụ lục P3.24 11.2955 10.1614 10.2752 10.7001 10.8824 10.2575 10.4446 10.4772 11.0120 10.9595 10.1415 10.6902 10.2304 10.6473 10.5926 9.7892 10.4456 10.0222 10.4668 10.0403 11.1260 11.2613 10.6500 10.0165 10.2447 10.5106 10.8973 10.0367 9.8961 10.1629 0.000003 0.001716 0.000929 0.000086 0.000024 0.001046 0.000244 0.000283 0.000016 0.000021 0.001551 0.000101 0.000993 0.000133 0.000166 0.010534 0.000406 0.003388 0.000337 0.003000 0.000009 0.000005 0.000125 0.003490 0.001225 0.000272 0.000024 0.003940 0.008119 0.001263 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 TB 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 10.6295 10.1126 10.2316 9.8847 9.8231 10.6910 10.4827 10.1294 10.6772 10.2738 9.5222 9.3763 10.6434 9.6873 10.6502 10.1772 10.0387 10.4867 11.2008 10.7517 10.9476 10.0546 9.7593 10.2010 10.5078 10.4476 10.5943 9.6255 9.7621 10.2656 10.340 0.000133 0.002342 0.001178 0.008465 0.010009 0.000098 0.000320 0.002200 0.000108 0.000906 0.046770 0.099740 0.000134 0.023783 0.000126 0.001515 0.002748 0.000260 0.000006 0.000074 0.000020 0.002979 0.013579 0.001479 0.000258 0.000339 0.000153 0.028426 0.016579 0.000799 0.004136 %============CHUONG TRINH SO SANH TY LE NEN DU LIEU================= %=========KHI CUNG LOI KHOI PHUC TRONG SBC(12642) VA SBC(8842)======== %***********************Tai nap du lieu am thanh************************* for i=1:100; y=wavread(['C:\MATLAB7\work\100file\' int2str(i) 'v.wav']); audio=y(5882:6763); %Doan am chieu dai 20ms %***************Bac bo loc, cua so Hamming va goc anpha**************** N=20; n=1:1:N; kn=0.54-0.46*cos((pi*(N-n))/N); anpha=(11/24)*pi; %************************BANK LOC SBC(12642)************************** %*************Thiet ke bo loc LPF******************* wcL=pi/12; BL=wcL; wsL=(wcL*tan(anpha)+1)/tan(anpha); Lbk0=(1+2*BL*tan(anpha))/(2*pi*tan(anpha)); 164 Lb1=(tan(anpha)*(cos(n*wcL)-cos(n*wsL)))/pi; Lb2=Lb1./n; Lbn=Lb2./n; Lbkn=Lbn.*kn; LPF=[Lbk0 Lbkn]; %*************Thiet ke bo loc BPF1******************* wcB11=pi/12; wcB21=pi/4; BB1=wcB21-wcB11; wsB11=((wcB11*tan(anpha))-1)/tan(anpha); wsB21=((wcB21*tan(anpha))+1)/tan(anpha); Bbk01=(1+BB1*tan(anpha))/(pi*tan(anpha)); Bb11=(tan(anpha)*(cos(n*wcB21)+cos(n*wcB11)-cos(n*wsB21)-cos(n*wsB11)))/pi; Bb21=Bb11./n; Bbn1=Bb21./n; Bbkn1=Bbn1.*kn; BPF1=[Bbk01 Bbkn1]; %*************Thiet ke bo loc BPF2******************* wcB12=pi/4; wcB22=pi/2; BB2=wcB22-wcB12; wsB12=((wcB12*tan(anpha))-1)/tan(anpha); wsB22=((wcB22*tan(anpha))+1)/tan(anpha); Bbk02=(1+BB2*tan(anpha))/(pi*tan(anpha)); Bb12=(tan(anpha)*(cos(n*wcB22)+cos(n*wcB12)-cos(n*wsB22)-cos(n*wsB12)))/pi; Bb22=Bb12./n; Bbn2=Bb22./n; Bbkn2=Bbn2.*kn; BPF2=[Bbk02 Bbkn2]; %*************Thiet ke bo loc HPF******************* wcH=pi/2; BH=pi-wcH; wsH=(wcH*tan(anpha)-1)/tan(anpha); Hbk0=(1+2*BH*tan(anpha))/(2*pi*tan(anpha)); Hb1=(tan(anpha)*(cos(n*wcH)-cos(n*wsH)))/pi; Hb2=Hb1./n; Hbn=Hb2./n; Hbkn=Hbn.*kn; HPF=[Hbk0 Hbkn]; %***********************BANK LOC SBC(8842)**************************** %*************Thiet ke bo loc RLPF******************* RwcL=pi/8; RBL=RwcL; RwsL=(RwcL*tan(anpha)+1)/tan(anpha); RLbk0=(1+2*RBL*tan(anpha))/(2*pi*tan(anpha)); RLb1=(tan(anpha)*(cos(n*RwcL)-cos(n*RwsL)))/pi; RLb2=RLb1./n; RLbn=RLb2./n; RLbkn=RLbn.*kn; RLPF=[RLbk0 RLbkn]; %*************Thiet ke bo loc RBPF1******************* RwcB11=pi/8; RwcB21=pi/4; RBB1=RwcB21-RwcB11; RwsB11=((RwcB11*tan(anpha))-1)/tan(anpha); RwsB21=((RwcB21*tan(anpha))+1)/tan(anpha); RBbk01=(1+RBB1*tan(anpha))/(pi*tan(anpha)); RBb11=(tan(anpha)*(cos(n*RwcB21)+cos(n*RwcB11)-cos(n*RwsB21)-cos(n*RwsB11)))/pi; RBb21=RBb11./n; RBbn1=RBb21./n; RBbkn1=RBbn1.*kn; RBPF1=[RBbk01 RBbkn1]; %*************Thiet ke bo loc RBPF2******************* 165 RwcB12=pi/4; RwcB22=pi/2; RBB2=RwcB22-RwcB12; RwsB12=((RwcB12*tan(anpha))-1)/tan(anpha); RwsB22=((RwcB22*tan(anpha))+1)/tan(anpha); RBbk02=(1+RBB2*tan(anpha))/(pi*tan(anpha)); RBb12=(tan(anpha)*(cos(n*RwcB22)+cos(n*RwcB12)-cos(n*RwsB22)-cos(n*RwsB12)))/pi; RBb22=RBb12./n; RBbn2=RBb22./n; RBbkn2=RBbn2.*kn; RBPF2=[RBbk02 RBbkn2]; %*************Thiet ke bo loc RHPF******************* RwcH=pi/2; RBH=pi-RwcH; RwsH=(RwcH*tan(anpha)-1)/tan(anpha); RHbk0=(1+2*RBH*tan(anpha))/(2*pi*tan(anpha)); RHb1=(tan(anpha)*(cos(n*RwcH)-cos(n*RwsH)))/pi; RHb2=RHb1./n; RHbn=RHb2./n; RHbkn=RHbn.*kn; RHPF=[RHbk0 RHbkn]; %**********************XU LY AM THANH QUA SBC(12642)******************* y0=filter(LPF,1,audio); y1=filter(BPF1,1,audio); y2=filter(BPF2,1,audio); y3=filter(HPF,1,audio); Py0=downsample(y0,12); Py1=downsample(y1,4); Py2=downsample(y2,2); Py3=downsample(y3,2); varK0=var(Py0); varK1=var(Py1); varK2=var(Py2); varK3=var(Py3); %**********************XU LY AM THANH QUA SBC(8842)****************** Ry0=filter(RLPF,1,audio); Ry1=filter(RBPF1,1,audio); Ry2=filter(RBPF2,1,audio); Ry3=filter(RHPF,1,audio); RPy0=downsample(Ry0,8); RPy1=downsample(Ry1,8); RPy2=downsample(Ry2,4); RPy3=downsample(Ry3,2); RvarK0=var(RPy0); RvarK1=var(RPy1); RvarK2=var(RPy2); RvarK3=var(RPy3); %*************************TINH TY LE NEN DU LIEU******************* RTB8842=256000; bTB8842=RTB8842/44100; E8842=4*0.5*pi*(sqrt(3))*(2^(2*bTB8842))*(RvarK0^(1/8))*(RvarK1^(1/8))*(RvarK2^(1/4))*(RvarK3^(1/2)); TG=E8842/(4*0.5*pi*(sqrt(3))*(varK0^(1/12))*(varK1^(1/4))*(varK2^(1/2))*(varK3^(1/2))); bTB12642=(log2(TG))/(-2) RTB12642=bTB12642*58800 TYLENENDULIEU12642=1411200/RTB12642 end; 166 %================================KET THUC========================= Phụ lục P3.25 Bảng kết tham số SBC(12642) so với SBC(8842), RTB(8842)=256kbps η4 TT b TB(12642) RTB(12642) TT bTB(12642) RTB (12642) 2.8346 166670 8.4668 51 3.5708 209960 3.6904 217000 6.5033 52 3.2888 193380 4.0124 235930 5.9815 53 3.8253 224930 3.6902 216980 6.5038 54 3.8813 228220 3.5839 210730 6.6967 55 3.1705 186430 3.3575 197420 7.1481 56 3.5349 207850 3.1490 185160 7.6216 57 3.4767 204430 3.4174 200940 7.0229 58 3.0422 178880 2.9730 174810 8.0727 59 3.8769 227960 10 3.7628 221250 6.3782 60 3.5133 206580 11 3.5120 206510 6.8337 61 2.7164 159730 12 2.5985 152790 9.2362 62 3.7730 221850 13 3.6132 212460 6.6422 63 3.5741 210160 14 2.8071 165060 8.5496 64 3.8555 226700 15 3.6883 216870 6.5071 65 4.0571 238550 16 3.4224 201240 7.0125 66 3.9285 230990 17 3.7421 220030 6.4136 67 3.8429 225960 18 3.9366 231470 6.0967 68 3.6331 213630 19 3.0541 179580 7.8584 69 3.7540 220740 20 3.0843 181360 7.7813 70 3.8856 228480 21 2.5095 147560 9.5636 71 3.1755 186720 22 3.6436 214240 6.5869 72 3.6924 217110 23 3.5298 207550 6.7993 73 3.5734 210120 24 3.1049 182570 7.7297 74 3.9203 230510 25 2.9226 171850 8.2120 75 3.9819 234140 26 3.5475 208590 6.7654 76 3.1140 183100 27 3.3604 197590 7.1421 77 3.3223 195350 28 3.3278 195670 7.2120 78 3.6756 216130 29 2.7930 164230 8.5929 79 3.1278 183910 30 2.8455 167320 8.4344 80 3.5312 207630 31 3.6635 215410 6.5511 81 4.2828 251830 32 3.1148 183150 7.7052 82 4.4287 260410 33 3.5746 210190 6.7140 83 3.1616 185900 η4 6.7213 7.2974 6.2740 6.1835 7.5697 6.7895 6.9030 7.8889 6.1906 6.8311 8.8351 6.3610 6.7150 6.2249 5.9156 6.1093 6.2452 6.6059 6.3931 6.1766 7.5579 6.4999 6.7163 6.1220 6.0273 7.7071 7.2239 6.5295 7.6732 6.7966 5.6038 5.4192 7.5912 35 3.2124 188890 7.4711 85 3.1548 185500 7.6075 36 37 38 39 40 41 42 43 44 4.0158 3.3594 3.7828 3.3382 3.7647 2.6790 2.5437 3.1550 3.7885 236130 197530 222430 196290 221360 157530 149570 185510 222760 5.9765 7.1442 6.3445 7.1894 6.3750 8.9584 9.4352 7.6070 6.3350 86 87 88 89 90 91 92 93 94 3.6278 3.7663 3.3183 2.6042 3.0533 2.8574 3.7503 4.0457 3.6040 213320 221460 195120 153120 179530 168010 220520 237890 211910 6.6155 6.3723 7.2326 9.2160 7.8603 8.3993 6.3994 5.9322 6.6594 167 45 46 47 48 49 50 3.5603 3.2944 2.9077 3.7682 3.9088 3.6420 209340 193710 170970 221570 229840 214150 6.7410 7.2850 8.2540 6.3690 6.1399 6.5897 95 96 97 98 99 100 TB 3.2972 3.3574 3.2107 4.1795 4.0429 3.5394 3.4646 193870 197420 188790 245750 237720 208120 203719 7.2790 7.1484 7.4750 5.7423 5.9364 6.7807 7.0333 Phụ lục chương Phụ lục P4.1 %============CHUONG TRINH SO SANH TY LE NEN DU LIEU================= %=========KHI CUNG LOI KHOI PHUC TRONG SBC(632) VA SBC(442)========== %***********************Tai nap du lieu am thanh************************* for i=1:100; y=wavread(['C:\MATLAB7\work\Viet\' int2str(i) 'v.wav']); audio=y(5882:6763); %Doan am chieu dai 20ms %***************Bac bo loc, cua so Hamming va goc anpha**************** N=20; n=1:1:N; kn=0.54-0.46*cos((pi*(N-n))/N); anpha=(11/24)*pi; %************************BANK LOC SBC(632)******************************** %***********************Thiet ke bo loc LPF******************************* wcL=pi/6; BL=wcL; wsL=(wcL*tan(anpha)+1)/tan(anpha); Lbk0=(1+2*BL*tan(anpha))/(2*pi*tan(anpha)); Lb1=(tan(anpha)*(cos(n*wcL)-cos(n*wsL)))/pi; Lb2=Lb1./n; Lbn=Lb2./n; Lbkn=Lbn.*kn; LPF=[Lbk0 Lbkn]; %***********************Thiet ke bo loc BPF****************************** wcB1=pi/6; wcB2=pi/2; BB=wcB2-wcB1; wsB1=((wcB1*tan(anpha))-1)/tan(anpha); wsB2=((wcB2*tan(anpha))+1)/tan(anpha); Bbk0=(1+BB*tan(anpha))/(pi*tan(anpha)); Bb1=(tan(anpha)*(cos(n*wcB2)+cos(n*wcB1)-cos(n*wsB2)-cos(n*wsB1)))/pi; Bb2=Bb1./n; Bbn=Bb2./n; Bbkn=Bbn.*kn; BPF=[Bbk0 Bbkn]; %************************Thiet ke bo loc HPF****************************** wcH=pi/2; BH=pi-wcH; wsH=(wcH*tan(anpha)-1)/tan(anpha); Hbk0=(1+2*BH*tan(anpha))/(2*pi*tan(anpha)); Hb1=(tan(anpha)*(cos(n*wcH)-cos(n*wsH)))/pi; Hb2=Hb1./n; Hbn=Hb2./n; 168 Hbkn=Hbn.*kn; HPF=[Hbk0 Hbkn]; %************************BANK LOC SBC(442)****************************** %***********************Thiet ke bo loc MLPF*************************** MwcL=pi/4; MBL=MwcL; MwsL=(MwcL*tan(anpha)+1)/tan(anpha); MLbk0=(1+2*MBL*tan(anpha))/(2*pi*tan(anpha)); MLb1=(tan(anpha)*(cos(n*MwcL)-cos(n*MwsL)))/pi; MLb2=MLb1./n; MLbn=MLb2./n; MLbkn=MLbn.*kn; MLPF=[MLbk0 MLbkn]; %***********************Thiet ke bo loc MBPF************************** MwcB1=pi/4; MwcB2=pi/2; MBB=MwcB2-MwcB1; MwsB1=((MwcB1*tan(anpha))-1)/tan(anpha); MwsB2=((MwcB2*tan(anpha))+1)/tan(anpha); MBbk0=(1+MBB*tan(anpha))/(pi*tan(anpha)); MBb1=(tan(anpha)*(cos(n*MwcB2)+cos(n*MwcB1)-cos(n*MwsB2)-cos(n*MwsB1)))/pi; MBb2=MBb1./n; MBbn=MBb2./n; MBbkn=MBbn.*kn; MBPF=[MBbk0 MBbkn]; %***********************Thiet ke bo loc MHPF************************** MwcH=pi/2; MBH=pi-MwcH; MwsH=(MwcH*tan(anpha)-1)/tan(anpha); MHbk0=(1+2*MBH*tan(anpha))/(2*pi*tan(anpha)); MHb1=(tan(anpha)*(cos(n*MwcH)-cos(n*MwsH)))/pi; MHb2=MHb1./n; MHbn=MHb2./n; MHbkn=MHbn.*kn; MHPF=[MHbk0 MHbkn]; %**********************XU LY AM THANH QUA SBC(632)******************** y0=filter(LPF,1,audio); y1=filter(BPF,1,audio); y2=filter(HPF,1,audio); Py0=downsample(y0,6); Py1=downsample(y1,2); Py2=downsample(y2,2); varK0=var(Py0); varK1=var(Py1); varK2=var(Py2); %**********************XU LY AM THANH QUA SBC(442)******************** My0=filter(MLPF,1,audio); My1=filter(MBPF,1,audio); My2=filter(MHPF,1,audio); MPy0=downsample(My0,4); MPy1=downsample(My1,4); MPy2=downsample(My2,2); MvarK0=var(MPy0); MvarK1=var(MPy1); MvarK2=var(MPy2); %*************************TINH TY LE NEN DU LIEU******************** RTB632=132000; 169 bTB632=RTB632/51450; E632=3*0.5*pi*(sqrt(3))*(2^(-2*bTB632))*(varK0^(1/6))*(varK1^(1/2))*(varK2^(1/2)) end; %================================KET THUC========================== Phụ lục P4.2 %============CHUONG TRINH SO SANH TY LE NEN DU LIEU================= %=========KHI CUNG LOI KHOI PHUC TRONG SBC(12642) VA SBC(8842)======= %***********************Tai nap du lieu am thanh************************* for i=1:100; y=wavread(['C:\MATLAB7\work\Viet\' int2str(i) 'v.wav']); audio=y(5882:6763); %Doan am chieu dai 20ms %***************Bac bo loc, cua so Hamming va goc anpha**************** N=20; n=1:1:N; kn=0.54-0.46*cos((pi*(N-n))/N); anpha=(11/24)*pi; %************************BANK LOC SBC(12642)************************** %*************Thiet ke bo loc LPF******************* wcL=pi/12; BL=wcL; wsL=(wcL*tan(anpha)+1)/tan(anpha); Lbk0=(1+2*BL*tan(anpha))/(2*pi*tan(anpha)); Lb1=(tan(anpha)*(cos(n*wcL)-cos(n*wsL)))/pi; Lb2=Lb1./n; Lbn=Lb2./n; Lbkn=Lbn.*kn; LPF=[Lbk0 Lbkn]; %*************Thiet ke bo loc BPF1******************* wcB11=pi/12; wcB21=pi/4; BB1=wcB21-wcB11; wsB11=((wcB11*tan(anpha))-1)/tan(anpha); wsB21=((wcB21*tan(anpha))+1)/tan(anpha); Bbk01=(1+BB1*tan(anpha))/(pi*tan(anpha)); Bb11=(tan(anpha)*(cos(n*wcB21)+cos(n*wcB11)-cos(n*wsB21)-cos(n*wsB11)))/pi; Bb21=Bb11./n; Bbn1=Bb21./n; Bbkn1=Bbn1.*kn; BPF1=[Bbk01 Bbkn1]; %*************Thiet ke bo loc BPF2******************* wcB12=pi/4; wcB22=pi/2; BB2=wcB22-wcB12; wsB12=((wcB12*tan(anpha))-1)/tan(anpha); wsB22=((wcB22*tan(anpha))+1)/tan(anpha); Bbk02=(1+BB2*tan(anpha))/(pi*tan(anpha)); Bb12=(tan(anpha)*(cos(n*wcB22)+cos(n*wcB12)-cos(n*wsB22)-cos(n*wsB12)))/pi; Bb22=Bb12./n; Bbn2=Bb22./n; Bbkn2=Bbn2.*kn; BPF2=[Bbk02 Bbkn2]; %*************Thiet ke bo loc HPF******************* wcH=pi/2; BH=pi-wcH; wsH=(wcH*tan(anpha)-1)/tan(anpha); Hbk0=(1+2*BH*tan(anpha))/(2*pi*tan(anpha)); Hb1=(tan(anpha)*(cos(n*wcH)-cos(n*wsH)))/pi; Hb2=Hb1./n; Hbn=Hb2./n; 170 Hbkn=Hbn.*kn; HPF=[Hbk0 Hbkn]; %***********************BANK LOC SBC(8842)*************************** %*************Thiet ke bo loc RLPF******************* RwcL=pi/8; RBL=RwcL; RwsL=(RwcL*tan(anpha)+1)/tan(anpha); RLbk0=(1+2*RBL*tan(anpha))/(2*pi*tan(anpha)); RLb1=(tan(anpha)*(cos(n*RwcL)-cos(n*RwsL)))/pi; RLb2=RLb1./n; RLbn=RLb2./n; RLbkn=RLbn.*kn; RLPF=[RLbk0 RLbkn]; %*************Thiet ke bo loc RBPF1******************* RwcB11=pi/8; RwcB21=pi/4; RBB1=RwcB21-RwcB11; RwsB11=((RwcB11*tan(anpha))-1)/tan(anpha); RwsB21=((RwcB21*tan(anpha))+1)/tan(anpha); RBbk01=(1+RBB1*tan(anpha))/(pi*tan(anpha)); RBb11=(tan(anpha)*(cos(n*RwcB21)+cos(n*RwcB11)-cos(n*RwsB21)-cos(n*RwsB11)))/pi; RBb21=RBb11./n; RBbn1=RBb21./n; RBbkn1=RBbn1.*kn; RBPF1=[RBbk01 RBbkn1]; %*************Thiet ke bo loc RBPF2******************* RwcB12=pi/4; RwcB22=pi/2; RBB2=RwcB22-RwcB12; RwsB12=((RwcB12*tan(anpha))-1)/tan(anpha); RwsB22=((RwcB22*tan(anpha))+1)/tan(anpha); RBbk02=(1+RBB2*tan(anpha))/(pi*tan(anpha)); RBb12=(tan(anpha)*(cos(n*RwcB22)+cos(n*RwcB12)-cos(n*RwsB22)-cos(n*RwsB12)))/pi; RBb22=RBb12./n; RBbn2=RBb22./n; RBbkn2=RBbn2.*kn; RBPF2=[RBbk02 RBbkn2]; %*************Thiet ke bo loc RHPF******************* RwcH=pi/2; RBH=pi-RwcH; RwsH=(RwcH*tan(anpha)-1)/tan(anpha); RHbk0=(1+2*RBH*tan(anpha))/(2*pi*tan(anpha)); RHb1=(tan(anpha)*(cos(n*RwcH)-cos(n*RwsH)))/pi; RHb2=RHb1./n; RHbn=RHb2./n; RHbkn=RHbn.*kn; RHPF=[RHbk0 RHbkn]; %**********************XU LY AM THANH QUA SBC(12642)****************** y0=filter(LPF,1,audio); y1=filter(BPF1,1,audio); y2=filter(BPF2,1,audio); y3=filter(HPF,1,audio); Py0=downsample(y0,12); Py1=downsample(y1,4); Py2=downsample(y2,2); Py3=downsample(y3,2); varK0=var(Py0); varK1=var(Py1); varK2=var(Py2); varK3=var(Py3); 171 %**********************XU LY AM THANH QUA SBC(8842)****************** Ry0=filter(RLPF,1,audio); Ry1=filter(RBPF1,1,audio); Ry2=filter(RBPF2,1,audio); Ry3=filter(RHPF,1,audio); RPy0=downsample(Ry0,8); RPy1=downsample(Ry1,8); RPy2=downsample(Ry2,4); RPy3=downsample(Ry3,2); RvarK0=var(RPy0); RvarK1=var(RPy1); RvarK2=var(RPy2); RvarK3=var(RPy3); %*************************TINH TY LE NEN DU LIEU********************** RTB8842=132000; bTB8842=RTB8842/44100; E8842=4*0.5*pi*(sqrt(3))*(2^(2*bTB8842))*(RvarK0^(1/8))*(RvarK1^(1/8))*(RvarK2^(1/4))*(RvarK3^(1/2)) end; %================================KET THUC======================= Phụ lục P4.3 Bảng lỗi khôi phục SBC(632) SBC(8842) tốc độ bít R=132kbps TT ε(632) (x10 -4) ε(8842) (x10 -4) TT ε(632) (x10 -4) ε(8842) (x10 -4) R=132kbps R=132kbps R=132kbps R=132kbps 0.00308450000 0.01187200000 51 0.12857000000 0.28495000000 0.19535000000 0.41375000000 52 0.03569600000 0.09284500000 0.77026000000 1.33080000000 53 0.46405000000 0.83661000000 0.19341000000 0.41456000000 54 0.54072000000 0.93940000000 0.14092000000 0.30328000000 55 0.01905800000 0.05479400000 0.04502900000 0.11657000000 56 0.10560000000 0.24181000000 0.01718000000 0.05000000000 57 0.09104200000 0.20889000000 0.05434800000 0.14085000000 58 0.00946220000 0.03046900000 0.00683570000 0.02308400000 59 0.43039000000 0.81453000000 10 0.23960000000 0.45302000000 60 0.10561000000 0.23733000000 11 0.09526800000 0.21869000000 61 0.00210690000 0.00845290000 12 0.00125550000 0.00537680000 62 0.28546000000 0.58766000000 13 0.05489700000 0.11417000000 63 0.12487000000 0.26413000000 14 0.00335640000 0.01241100000 64 0.39205000000 0.74061000000 15 0.19027000000 0.40367000000 65 1.23660000000 2.01920000000 16 0.06606900000 0.15823000000 66 0.56638000000 1.05430000000 18 0.54912000000 0.90232000000 68 0.15960000000 0.34053000000 19 20 21 22 23 24 25 26 27 00791410000 0.01070400000 0.00054465000 0.16816000000 0.09892500000 0.01230800000 0.00356970000 0.11056000000 0.02822400000 02765900000 0.03425000000 0.00265690000 0.35491000000 0.22499000000 0.03769300000 0.01363100000 0.24718000000 0.07470600000 69 70 71 72 73 74 75 76 77 10966000000 0.61802000000 0.01801000000 0.22292000000 0.12073000000 0.68917000000 0.70979000000 0.01422400000 0.04005500000 20960000000 1.09490000000 0.05259900000 0.45269000000 0.26859000000 1.19330000000 1.29530000000 0.04241600000 0.10331000000 172 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 0.03425200000 0.00279650000 0.00355780000 0.13966000000 0.01402900000 0.09149000000 0.01815200000 0.02205900000 0.75141000000 0.04896200000 0.30175000000 0.04059100000 0.25480000000 0.00169980000 0.00097284000 0.01748800000 0.29174000000 0.12769000000 0.03483900000 0.00430120000 0.35350000000 0.67104000000 0.12449000000 Phụ lục P4.4 0.09077600000 78 0.01058000000 79 0.01316900000 80 0.31220000000 81 0.04329100000 82 0.22600000000 83 0.05399700000 84 0.06246400000 85 1.30080000000 86 0.12465000000 87 0.57794000000 88 0.10656000000 89 0.52462000000 90 0.00700300000 91 0.00434650000 92 0.05078500000 93 0.59053000000 94 0.28449000000 95 0.09137000000 96 0.01376100000 97 0.68561000000 98 1.16270000000 99 0.26189000000 100 TB 0.21171000000 0.01540700000 0.09419700000 2.82890000000 6.00420000000 0.01821100000 1.72470000000 0.01731200000 0.15405000000 0.24480000000 0.02871100000 0.00117550000 0.01125700000 0.00314980000 0.26198000000 0.90780000000 0.15078000000 0.03228600000 0.03896800000 0.01991300000 1.96060000000 1.24640000000 0.07613800000 0.3039 0.43525000000 0.04580800000 0.21893000000 3.98860000000 6.94790000000 0.05396100000 2.54970000000 0.05141600000 0.32037000000 0.47961000000 0.08442800000 0.00518610000 0.03446900000 0.01209400000 0.53151000000 1.60860000000 0.32333000000 0.08643200000 0.10420000000 0.05760400000 2.79740000000 1.97180000000 0.19105000000 0.5038 %============CHUONG TRINH SO SANH TY LE NEN DU LIEU================= %=========KHI CUNG LOI KHOI PHUC TRONG SBC(12642) VA SBC(8842)======= %***********************Tai nap du lieu am thanh************************* for i=1:100; y=wavread(['C:\MATLAB7\work\Viet\' int2str(i) 'v.wav']); audio=y(5882:6763); %Doan am chieu dai 20ms %***************Bac bo loc, cua so Hamming va goc anpha**************** N=20; n=1:1:N; kn=0.54-0.46*cos((pi*(N-n))/N); anpha=(11/24)*pi; %************************BANK LOC SBC(12642)************************** %*************Thiet ke bo loc LPF******************* wcL=pi/12; BL=wcL; wsL=(wcL*tan(anpha)+1)/tan(anpha); Lbk0=(1+2*BL*tan(anpha))/(2*pi*tan(anpha)); Lb1=(tan(anpha)*(cos(n*wcL)-cos(n*wsL)))/pi; Lb2=Lb1./n; Lbn=Lb2./n; Lbkn=Lbn.*kn; LPF=[Lbk0 Lbkn]; %*************Thiet ke bo loc BPF1******************* wcB11=pi/12; wcB21=pi/4; BB1=wcB21-wcB11; wsB11=((wcB11*tan(anpha))-1)/tan(anpha); wsB21=((wcB21*tan(anpha))+1)/tan(anpha); 173 Bbk01=(1+BB1*tan(anpha))/(pi*tan(anpha)); Bb11=(tan(anpha)*(cos(n*wcB21)+cos(n*wcB11)-cos(n*wsB21)-cos(n*wsB11)))/pi; Bb21=Bb11./n; Bbn1=Bb21./n; Bbkn1=Bbn1.*kn; BPF1=[Bbk01 Bbkn1]; %*************Thiet ke bo loc BPF2******************* wcB12=pi/4; wcB22=pi/2; BB2=wcB22-wcB12; wsB12=((wcB12*tan(anpha))-1)/tan(anpha); wsB22=((wcB22*tan(anpha))+1)/tan(anpha); Bbk02=(1+BB2*tan(anpha))/(pi*tan(anpha)); Bb12=(tan(anpha)*(cos(n*wcB22)+cos(n*wcB12)-cos(n*wsB22)-cos(n*wsB12)))/pi; Bb22=Bb12./n; Bbn2=Bb22./n; Bbkn2=Bbn2.*kn; BPF2=[Bbk02 Bbkn2]; %*************Thiet ke bo loc HPF******************* wcH=pi/2; BH=pi-wcH; wsH=(wcH*tan(anpha)-1)/tan(anpha); Hbk0=(1+2*BH*tan(anpha))/(2*pi*tan(anpha)); Hb1=(tan(anpha)*(cos(n*wcH)-cos(n*wsH)))/pi; Hb2=Hb1./n; Hbn=Hb2./n; Hbkn=Hbn.*kn; HPF=[Hbk0 Hbkn]; %***********************BANK LOC SBC(8842)*************************** %*************Thiet ke bo loc RLPF******************* RwcL=pi/8; RBL=RwcL; RwsL=(RwcL*tan(anpha)+1)/tan(anpha); RLbk0=(1+2*RBL*tan(anpha))/(2*pi*tan(anpha)); RLb1=(tan(anpha)*(cos(n*RwcL)-cos(n*RwsL)))/pi; RLb2=RLb1./n; RLbn=RLb2./n; RLbkn=RLbn.*kn; RLPF=[RLbk0 RLbkn]; %*************Thiet ke bo loc RBPF1******************* RwcB11=pi/8; RwcB21=pi/4; RBB1=RwcB21-RwcB11; RwsB11=((RwcB11*tan(anpha))-1)/tan(anpha); RwsB21=((RwcB21*tan(anpha))+1)/tan(anpha); RBbk01=(1+RBB1*tan(anpha))/(pi*tan(anpha)); RBb11=(tan(anpha)*(cos(n*RwcB21)+cos(n*RwcB11)-cos(n*RwsB21)-cos(n*RwsB11)))/pi; RBb21=RBb11./n; RBbn1=RBb21./n; RBbkn1=RBbn1.*kn; RBPF1=[RBbk01 RBbkn1]; %*************Thiet ke bo loc RBPF2******************* RwcB12=pi/4; RwcB22=pi/2; RBB2=RwcB22-RwcB12; RwsB12=((RwcB12*tan(anpha))-1)/tan(anpha); RwsB22=((RwcB22*tan(anpha))+1)/tan(anpha); RBbk02=(1+RBB2*tan(anpha))/(pi*tan(anpha)); RBb12=(tan(anpha)*(cos(n*RwcB22)+cos(n*RwcB12)-cos(n*RwsB22)-cos(n*RwsB12)))/pi; RBb22=RBb12./n; RBbn2=RBb22./n; RBbkn2=RBbn2.*kn; RBPF2=[RBbk02 RBbkn2]; 174 %*************Thiet ke bo loc RHPF******************* RwcH=pi/2; RBH=pi-RwcH; RwsH=(RwcH*tan(anpha)-1)/tan(anpha); RHbk0=(1+2*RBH*tan(anpha))/(2*pi*tan(anpha)); RHb1=(tan(anpha)*(cos(n*RwcH)-cos(n*RwsH)))/pi; RHb2=RHb1./n; RHbn=RHb2./n; RHbkn=RHbn.*kn; RHPF=[RHbk0 RHbkn]; %**********************XU LY AM THANH QUA SBC(12642)****************** y0=filter(LPF,1,audio); y1=filter(BPF1,1,audio); y2=filter(BPF2,1,audio); y3=filter(HPF,1,audio); Py0=downsample(y0,12); Py1=downsample(y1,4); Py2=downsample(y2,2); Py3=downsample(y3,2); varK0=var(Py0); varK1=var(Py1); varK2=var(Py2); varK3=var(Py3); %**********************XU LY AM THANH QUA SBC(8842)****************** Ry0=filter(RLPF,1,audio); Ry1=filter(RBPF1,1,audio); Ry2=filter(RBPF2,1,audio); Ry3=filter(RHPF,1,audio); RPy0=downsample(Ry0,8); RPy1=downsample(Ry1,8); RPy2=downsample(Ry2,4); RPy3=downsample(Ry3,2); RvarK0=var(RPy0); RvarK1=var(RPy1); RvarK2=var(RPy2); RvarK3=var(RPy3); %*************************TINH TY LE NEN DU LIEU******************* RTB12642=256000; bTB12642=RTB12642/58800; E12642=4*0.5*pi*(sqrt(3))*(2^(2*bTB12642))*(varK0^(1/12))*(varK1^(1/4))*(RvarK2^(1/2))*(varK3^(1/2)) TG=E12642/(4*0.5*pi*(sqrt(3))*(RvarK0^(1/8))*(RvarK1^(1/8))*(RvarK2^(1/4))*(RvarK3^(1/2))); bTB8842=(log2(TG))/(-2); RTB8842=bTB8842*44100 end; %================================KET THUC=========================== Phụ lục P4.5 Bảng giá trị R(8842) ε(8842)=ε(12642) R(12642)=256kbps TT R(8842) ε(12642) (x10 -4) TT R(8842) R=256kbps 323090 0.00002923 51 290390 285230 0.00334820 52 302930 ε(12642) (x10 -4) R=256kbps 0.00196070 0.00043074 175 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 271490 285250 289950 299640 308880 297250 316840 281850 293090 333320 288210 324250 285310 297050 282990 274450 312760 311960 337210 287260 292230 310950 318660 291550 299770 301220 324810 322280 286430 310630 290350 308320 306230 270890 299830 281130 300780 281980 329760 335780 0.01658600 0.00335250 0.00211620 0.00059972 0.00019237 0.00078109 0.00006917 0.00407670 0.00138220 0.00000959 0.00084118 0.00002945 0.00325810 0.00088295 0.00433500 0.01024900 0.00009420 0.00011962 0.00000420 0.00269450 0.00146100 0.00013589 0.00003857 0.00163990 0.00038275 0.00044433 0.00002467 0.00003325 0.00243260 0.00015768 0.00155680 0.00021145 0.00026126 0.01652300 0.00063753 0.00531980 0.00052897 0.00470240 0.00001398 0.00000718 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 279260 276810 308090 292030 294670 313770 276960 293020 328170 281500 290340 278000 268950 274740 278500 287440 282860 276570 307970 285130 290330 275080 272350 310650 301430 285770 310050 292280 259110 252690 308610 266250 308850 288020 281930 301160 333140 313270 321920 282600 0.00816670 0.00990670 0.00021614 0.00158020 0.00125660 0.00010052 0.00854790 0.00150360 0.00001773 0.00534680 0.00182000 0.00752210 0.02726300 0.01186400 0.00822800 0.00257100 0.00182760 0.01163400 0.00020829 0.00367450 0.00185170 0.01328600 0.01571300 0.00015436 0.00050240 0.00346270 0.00016992 0.00141950 0.07335700 0.15640000 0.00020938 0.03746900 0.00019804 0.00237480 0.00430570 0.00041408 0.00000931 0.00011553 0.00003089 0.00467210 44 280920 0.00547310 94 288990 0.00232470 45 46 47 48 49 50 290930 302700 319740 281750 275580 287280 00192450 95 0.00042692 96 0.00003763 97 0.00619000 98 0.01274400 99 0.00198730 100 TB 302490 299820 306410 263710 269580 291900 295147 00040662 0.00053308 0.00023952 0.04453300 0.02609500 0.00125360 0.00648999 -HẾT - ... nhiều nhịp để phân tích tổng hợp tín hiệu Trong bank lọc nhiều nhịp, tần số lấy mẫu tín hiệu thay đổi q trình xử lý tín hiệu (qua phân chia nội suy) Do vậy, tín hiệu qua phân chia nội suy bị biến... Dạng tín hiệu âm gốc tín hiệu đầu bank lọc phân tích SBC(12642) Hình 3.25: Phổ tần tín hiệu vào x(n) fS=fNy Hình 3.26: Phổ tần tín hiệu dải đầu bank lọc phân tích Hình 3.27: Phổ tần tín hiệu. .. Dạng tín hiệu âm gốc tín hiệu đầu bank lọc phân tích SBC(632) Hình 3.9: Phổ tần tín hiệu vào x(n) f S=fNy Hình 3.10: Phổ tần ba tín hiệu dải đầu bank lọc phân tích Hình 3.11: Phổ tần tín hiệu