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 1.1- Tổng quan mã hoá dải con……………………………………………………… 1.1.1- Giới thiệu………………………………………………………………… 1.1.2- Nguyên lý mã hố dải (SBC)………………………………………… 1.1.3- Các thơng số kỹ thuật SBC………………………………… 1.1.4- Các điều kiện mã hoá dải con………………………………………… 1.1.5- Điều kiện để tăng hệ số nén tín hiệu SBC………………………… 1.2- Bank lọc theo phân tích Wavelet…………………………………………………… 1.2.1- Bank lọc nhiều nhịp đơn phân giải……………………………………… 1.2.2- Bank lọc nhiều nhịp đa phân giải………………………………………… 1.3- Mã hoá dải ứng dụng xử lý âm số……………………………… 1.3.1- Chuẩn âm số MPEG/audio……………………………… 1.3.2- Thuật toán ATRAC……………………………………………………… 1.3.3- Thuật toán PASC………………………………………………………… 1.3.4- Bộ mã hoá dải APT-X………………………………………………… 1.4- Kết luận……………………………………………………………………………… 10 10 10 12 15 18 21 25 25 27 30 30 32 34 35 35 Chương 2: Mã hoá giải đa phân giải tuyệt đối 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 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ố 4.1- Giới thiệu hệ thống phát số nay………………………………………… 4.1.1- Hệ thống phát số dải kênh………………………… 4.1.2- Hệ thống phát số EUREKA-147…………………………………… 4.2- Ứng dụng SBC(632) SBC(12642) phát số………………… 4.2.1- Ứng dụng SBC(632) SBC(12642) phát số qua Internet……………………………………………………………… 4.2.2- Ứng dụng SBC(632) SBC(12642) hệ thống phát số IBOC… 4.2.3- Ứng dụng SBC(632) SBC(12642) hệ thống phát số EUREKA-147……………………………………………… 4.3- Định hướng phát số vô tuyến Việt nam……………………………….…… 4.4- Kết luận…………………………………………………………………… 107 107 108 110 113 113 114 115 117 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 ∆fB δp δs BOFDM BAM-IBOC BFM-IBOC G GdB GL GLdB GB GBdB GH GHdB HL(ejF) HB(ejF) HH(ejF) kn Lbn Hbn Lbkn Bbkn Hbkn D I ni β3 β4 ∆fOFDM ε ε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 NOFDM 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ứ nguyên) 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ã hoá 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 khơng phải â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 bTB SBC(12642) SBC(8842) lỗi khôi phục ε TT bTB(12642) bTB(8842) ε TT bTB(12642) bTB(8842) (x10-6) 8 10.9704 0.000019 51 10.2342 10.1146 0.002134 52 10.5162 9.7926 0.010726 53 9.9797 8 10.1148 0.002138 54 9.9237 8 10.2211 0.001350 55 10.6345 8 10.4475 0.000379 56 10.2701 8 10.6560 0.000122 57 10.3283 8 10.3876 0.000498 58 10.7627 8 10.8320 0.000044 59 9.9281 8 10 10.0422 0.002584 60 10.2916 8 11 10.2930 0.000881 61 11.0886 8 12 11.2065 0.000006 62 10.0320 8 13 10.1918 0.000529 63 10.2309 14 10.9978 0.000019 64 9.9495 15 10.1167 0.002076 65 9.7479 8 16 10.3826 0.000563 66 9.8765 8 17 10.0629 0.002767 67 9.9621 8 18 9.8684 0.006547 68 10.1719 8 19 10.7509 0.000059 69 10.0509 8 20 10.7207 0.000076 70 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 0.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 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 Phụ lục P3.24 %============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 TT bTB(12642) RTB(12642) η4 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 34 3.1577 185670 7.6004 84 4.1177 242120 35 3.2124 188890 7.4711 85 3.1548 185500 36 4.0158 236130 5.9765 86 3.6278 213320 37 3.3594 197530 7.1442 87 3.7663 221460 38 3.7828 222430 6.3445 88 3.3183 195120 39 3.3382 196290 7.1894 89 2.6042 153120 40 3.7647 221360 6.3750 90 3.0533 179530 41 2.6790 157530 8.9584 91 2.8574 168010 42 2.5437 149570 9.4352 92 3.7503 220520 43 3.1550 185510 7.6070 93 4.0457 237890 44 3.7885 222760 6.3350 94 3.6040 211910 η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 5.8285 7.6075 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 17 0.23016000000 0.49928000000 67 0.44688000000 0.82294000000 18 0.54912000000 0.90232000000 68 0.15960000000 0.34053000000 19 0.00791410000 0.02765900000 69 0.10966000000 0.20960000000 20 0.01070400000 0.03425000000 70 0.61802000000 1.09490000000 21 0.00054465000 0.00265690000 71 0.01801000000 0.05259900000 22 0.16816000000 0.35491000000 72 0.22292000000 0.45269000000 23 0.09892500000 0.22499000000 73 0.12073000000 0.26859000000 24 0.01230800000 0.03769300000 74 0.68917000000 1.19330000000 25 0.00356970000 0.01363100000 75 0.70979000000 1.29530000000 26 0.11056000000 0.24718000000 76 0.01422400000 0.04241600000 27 0.02822400000 0.07470600000 77 0.04005500000 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 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 Phụ lục P4.4 %============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 43 44 45 46 47 48 49 50 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 308870 280920 290930 302700 319740 281750 275580 287280 0.01658600 53 0.00335250 54 0.00211620 55 0.00059972 56 0.00019237 57 0.00078109 58 0.00006917 59 0.00407670 60 0.00138220 61 0.00000959 62 0.00084118 63 0.00002945 64 0.00325810 65 0.00088295 66 0.00433500 67 0.01024900 68 0.00009420 69 0.00011962 70 0.00000420 71 0.00269450 72 0.00146100 73 0.00013589 74 0.00003857 75 0.00163990 76 0.00038275 77 0.00044433 78 0.00002467 79 0.00003325 80 0.00243260 81 0.00015768 82 0.00155680 83 0.00021145 84 0.00026126 85 0.01652300 86 0.00063753 87 0.00531980 88 0.00052897 89 0.00470240 90 0.00001398 91 0.00000718 92 0.00019546 93 0.00547310 94 0.00192450 95 0.00042692 96 0.00003763 97 0.00619000 98 0.01274400 99 0.00198730 100 TB -HẾT - 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 269550 288990 302490 299820 306410 263710 269580 291900 295147 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 0.02130800 0.00232470 0.00040662 0.00053308 0.00023952 0.04453300 0.02609500 0.00125360 0.00648999 ... thống phát số coi chuẩn phát số, ứng dụng số nước châu Âu Một số nước phát triển hệ thống phát số riêng Việc chọn hệ thống phát số nào, để phù hợp với điều kiện, hoàn cảnh Việt nam thực trạng phát. .. 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. .. thuật toán áp dụng mã hố-nén tín hiệu âm Để nén tín hiệu âm thanh, thuật tốn dùng mơ hình tâm lý- thính giác để lượng tử hố thích nghi thành phần tín hiệu tai người nghe Những thành phần tín hiệu ngưỡng