Giảm sai số hướng truyền cho hệ thống optical beamforming báo cáo tổng kết kết quả đề tài khcn cấp trường msđt t đđt 2013 17

ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA FOG BÁO CÁO TỔNG KẾT KẾT QUẢ ĐỀ TÀI KHCN CẤP TRƯỜNG Tên đề tài: Giảm sai số hướng truyền cho hệ thống optical beamforming Mã số đề tài: T-ĐĐT-2013-17 Thời gian thực hiện: 05/2013-05/2014 Chủ nhiệm đề tài: TS Phạm Quang Thái Cán tham gia đề tài: TS Phạm Quang Thái Thành phố Hồ Chí Minh – Tháng 03/2014 Danh sách cán tham gia thực đề tài 1-TS Phạm Quang Thái, đơn vị Khoa Điện – Điện tử, môn Viễn thông Mục lục Tổng quan 1.1 Mục đính nghiên cứu 1.2 Nội dung nghiên cứu Các kết nghiên cứu 2.1 Phân tích cấu trúc hệ thống beamforming quang học 2.2 Phân tích factorial Kết luận kết 13 Tài liệu tham khảo 14 Tổng quan 1.1 Mục đính nghiên cứu Kỹ thuật beamforming nghiên cứu ứng dụng rộng rãi nhiều ứng dụng thuộc lĩnh vực anten truyền sóng Bằng cách thay đổi pha biên độ tín hiệu cung cấp cho hệ thống mảng anten, đồ thị phát xạ (radiation pattern) mảng anten điều khiển theo mong muốn [1] Những lợi mảng anten so với anten thông thường độ lợi lớn, độ định hướng tốt, khả điều khiển hướng thu/phát khiến lĩnh vực tập trung nghiên cứu phát triển mạnh mẽ Để chuyển hướng phát anten đến góc θ, độ lệch pha ϕ hai anten kế cận mảng anten tính theo công thức [1]: ϕ= 2πd cos θ λ (1) Trong d khoảng cách hai anten liên tiếp mảng anten λ bước sóng tín hiệu Những mảng anten ban đầu thiết kế để hoạt động với băng thông hẹp Một mảng anten thơng thường có băng thơng khoảng 1% tần số tín hiệu [1] Với điều kiện này, mảng anten hoạt động theo yêu cầu Tuy nhiên, tần số tín hiệu băng thơng tín hiệu tăng dần, công thức (1) dẫn đến tượng trôi hướng phát chính, cịn gọi tượng beam-squint Trong cơng thức (1), khoảng cách d thường không đổi, λ thay đổi với băng thông lớn, độ lệch pha ϕ kéo theo hướng phát thay đổi Để khắc phục vấn đề này, thay làm trễ pha, ta làm trễ tín hiệu hai anten kế cận lượng thời gian τ theo công thức [2]: τ = d cos θ c (2) Với c số vận tốc ánh sáng Trong (2), độ trễ thời gian không phụ thuộc vào bước sóng tránh tượng beam-squint Tuy nhiên, việc làm trễ thời gian hệ thống mảng anten điều khiển miền điện cần dùng đường dây truyền sóng vốn phức tạp thiết kế, với tín hiệu tần số cao băng thơng rộng Cách ly tín hiệu kênh triệt nhiễu EMI vấn đề cần giải Chính vậy, việc thưc beamforming miền quang coi giải pháp khả thi cho tín hiệu tần số cao băng thông rộng [3] Trong miền quang, tín hiệu làm trễ thiết bị quang có băng thơng lớn, nhỏ gọn khơng nhiễu điện từ Có hai phương pháp chủ yếu dùng để làm thay đổi thời gian truyền tín hiệu quang [3]: • Giữ ngun bước sóng, sử dụng thiết bị không tán sắc (non-dispersive) thay đổi quãng đường di chuyển ánh sáng • Giữ nguyên chiều dài đường truyền, dùng thiết bị có tán sắc cao (dispersive) thay đổi bước sóng truyền qua Trong phương pháp thứ nhất, nguồn phát xạ rẻ tiền hơn, hệ thống làm trễ có suy hao cao kích thước lớn [4-6] Ngược lại, phương pháp thứ hai dùng nguồn phát xạ chuyển đổi bước sóng đắt tiền (tunable laser), thiết bị làm trễ lại nhỏ gọn [7-9] Một loại thiết bị sử dụng nhiều cho phương pháp thứ hai fiber grating [10] Trong[7-9], hệ thống sử dụng fiber grating đạt chiều ưu điểm giảm số lượng thiết bị kích thước Trong hệ thống sử dụng phương pháp thứ hai, đặc tuyến group delay-wavelegnth thiết bị tán sắc liên quan trực tiếp đến độ chênh lệch pha hai phần tử kế cận mảng anten Thông thường, đặc tuyến cần có tính tuyến tính cao Tuy nhiên, trình chế tạo fiber grating, sai số gây nhiệt độ, rung động cấu trúc không lý tưởng vật liệu khiến đặc tuyến group delay-wavelegnth thiết bị ln có sai số phân bố gợn sóng ngẫu nhiên dọc theo băng thông thiết bị Các sai số khiến cho giá trị tần số khác băng thơng tín hiệu RF, giá trị trễ pha Như phân tích [11], việc sử dụng fiber grating dẫn đến sai số thời gian trễ biên độ tín hiệu Cụ thể, sai số hướng tuyền hệ thống tăng theo băng thông, tần số tín hiệu điều chế, góc xoay độ gợn sóng đặc tuyến group delaywavelegnth grating Sử dụng kết nghiên cứu [11], hệ thống beamforming kết hợp [12] bước đầu giảm sai số hướng truyền Tuy nhiên, mức độ ảnh hưởng thông số tác động lẫn thơng số chưa phân tích Do đó, điều kiện thông số hoạt động khác nhau, hệ thống [12] cần đo đạc mô lại Việc tốn thời gian, công sức số trường hợp thực hạn chế mặt trang thiết bị Một phân tích mang tính định lượng thơng qua mơ hình tán học sai số hướng truyền hàm theo yếu tố phân tích [11] cho mơ hình [12] đóng góp nhiều mặt thiết kế tối ưu hóa hệ thống [12] 1.2 Nội dung nghiên cứu Đề tài tập trung nghiên cứu cách thức giảm sai số hướng truyền cho hệ thống optical beamforming Cụ thể: • Phân tích mối liên hệ mức độ ảnh hưởng lên sai số hướng truyền yếu tố băng thơng tín hiệu điều chế, tần số tín hiệu điều chế, góc xoay đặc tuyến group delay-wavelength grating Phương pháp phân tích thống kê factorial sử dụng • Từ phân tích trên, mơ hình tốn học mối liên quan yếu tố đến sai số hướng truyền thiết lập • Các yếu tố ảnh hưởng mạnh đến sai số hướng truyền cải thiện thông qua chỉnh sửa cải tiến tương ứng hệ thống optical beamforming Qua đó, giảm sai số hướng truyền cho hệ thống Các kết nghiên cứu 2.1 Phân tích cấu trúc hệ thống beamforming quang học Một hệ thống beamforming quang học sử dụng grating điển hình, thể hình 1, tách thành bốn khối sau: Nguồn laser điều chỉnh bước sóng để cung cấp tin hiệu sóng mang miền quang Thiết bị điều chế điện quang (electro-optic modulator - EOM) để điều chế tín hiệu RF cần đưa đến mảng anten với tín hiệu sóng mang quang từ nguồn laser Sợi grating để gây trễ miền thời gian Các photo-diode để chuyển đổi tín hiệu quang miền điện Do thiết bị photo-diode thiết bị phát thay đổi cường độ ánh sáng băng thông điều chế thiết bị nhỏ so với tần số sóng mang, nên photo-diode đồng thời đóng vai trị giải điều chế biên độ Hình 1: sơ đồ hệ thống beamforming quang tổng quát Thay đổi hướng phát mảng anten thực cách thêm khoảng thời gian trễ khác vào tín hiệu phần từ mảng anten Trong hệ thống xét, tín hiệu quang sau điều chế với tín hiệu RF truyền qua sợi grating Bộ circulator có chức định hướng ánh sáng vào-ra theo hướng mũi tên Khi ánh sáng truyền vào sợi grating, thay đổi chiết suất theo điều kiện Braggs sợi làm ánh sáng bước sóng khác phản xạ ngược lại vị trí khác sợi Ánh sáng phản xạ theo circulator truyền tiếp đến phía photo-diode Do quãng đường di chuyển khác sợi grating, bước sóng có mức độ trễ miền thời gian khác biệt Để điều khiển mức trễ theo ý muốn, ta cần điều chỉnh ảnh sáng laser tương ứng phần tử mảng anten đến bước sóng khác Hình thể kết đo đạc đặc tuyến group delay-wavelegnth sợi grating sử dụng [12] Trong hình này, phần tơ đậm bước sóng tín hiệu quang điều chế chiếu vào sợi grating Có bốn tín hiệu cho bốn phần tử anten Hình vẽ thể mức độ trễ 0, nghĩa hướng phát anten broadside Giả sử cần chuyển hướng búp chính, ta cần chỉnh bước sóng bốn nguồn laser để đạt độ chênh lệch thời gian cần thiết Hình 2: kết đo đạc đặc tuyến group delay-wavelegnth từ kênh 27 đến kênh 30 sợi grating (đường chấm chấm) khu vực có tín hiệu điều chế (đường liền nét) Thực tế, trình chế tạo sợi graitng ln tạo giá trị lỗi ngẫu nhiên đặc tuyến group delay-wavelegnth Như vậy, xét thành phần tần số băng thông tín hiệu điều chế cho phần tử anten, độ trễ thời gian có khác biệt giá trị lỗi ngẫu nhiên Khi vị trí tín hiệu điều chế thay đổi sợi grating, chẳng băng thơng điều chế thay đổi, chuyển bước sóng nguồn laser, thay thay đổi tần số RF, giá trị khác biệt độ trễ miền thời gian thành phần tần số thay đổi khác Kết hệ thống beamforming quang sử dụng grating bị giới hạn nhiều mặt phân tích [11] Trong đó, ảnh hưởng xấu hướng búp thay đổi đến vị trí khác thành phần tần số RF Điều may mắn giá trị sai số sợi grating gây trình chế tạo, nên không thay đổi sau chế tạo xong Do đó, ảnh hưởng sai số đến hệ thống sai số hướng truyền tính tốn cho sợi grating Phân tích [11] chứng minh sai số hướng truyền có liên quan đến yếu tố bao gồm: băng thơng tín hiệu RF, tần số RF, góc xoay độ chênh lệch thời gian trễ cặp khu vực sợi grating có sóng mang Tuy nhiên, tính chất phức tạp cấu trúc hệ thống, mơ hình tốn học yếu tố dựa vật lý phức tạp Ngược lại, mơ hình tốn thống kê xây dựng phân tích factorial Để ngắn gọn, băng thơng tín hiệu RF, tần số RF góc xoay ký BW, RF θ Độ chênh lệch thời gian trễ cặp khu vực sợi grating có sóng mang thể giá trị χ2 test cho khác biệt phân bố giá trị độ trễ Phương pháp phân tích factorial dùng để đánh giá mức độ ảnh hưởng yếu tố đến phương sai sai số hướng truyền σ2 Lý thuyết chi tiết phân tích factorial tham khảo [13] Nói ngắn gọn, yếu tố BW, RF, θ χ2 cần hai giá trị nhỏ lớn biên giới hạn miền giá tị cần phân tích Giá trị σ2 có đo đạc mơ biên dùng để phân tích thống kê factorial tìm σ2 cho vùng khơng gian giá trị bên biên Để tính tốn thống kê cách xác, thực tế điều kiện tổ hợp giá trị biên đo đạc nhiều lần Tuy nhiên, trường hợp mà đo đạc khơng thể lặp lại được, tổ hợp đo đạc lần mức độ quan trọng yếu tố Tuy nhiên, phân tích factorial sử dụng số trường hợp cách dùng thêm phép biến đổi giá trị đo Cụ thể trường hợp này, giá trị sợi gratng đo đạc lần hạn chế thiết bị thời gian đo Tính tốn độ quan trọng yếu tố facatorial khơng có kết với giá trị σ2 tìm thể mục phân tích factorial Vì vậy, phép tán biến đổi σ2 áp dụng Các cách biến đổi thông dụng σ , ln(σ ) σ Trong đó, cách biến đổi σ dùng Do ta cần tối thiểu hóa giá trị σ2 sau để có sai số hướng truyền nhỏ Tối thiểu hóa σ2 đồng nghĩa với tìm σ có tiệm cận gần Mơ hình tốn học σ có dạng n σ = β + ∑ β j xij + ε i (3) j =1 Trong tổng số yếu tố cần phân tích n xj yếu tố BW, RF, θ, χ2 mối liên hệ yếu tố i giá trị biên yếu tố xét xj εi sai số ngẫu nhiên đo đạc βj giải từ hệ phương trình sau: n n n n i =1 i =1 nβ + β1 ∑ xi1 + β ∑ xi + +β k ∑ xik = ∑ yi i =1 i =1 n n n n n i =1 i =1 i =1 i =1 i =1 n n n n n i =1 i =1 i =1 i =1 i =1 β ∑ xi1 + β1 ∑ x 2i1 + β ∑ xi1 xi + +β k ∑ xi1 xik = ∑ xi1 yi (4) β ∑ xik + β1 ∑ xi1 xi1 + β ∑ xi1 xi + +β k ∑ x 2ik = ∑ xik yi Sử dụng phương trình (3), giá trị σ tương ứng với giá trị BW, RF, θ χ2 không gian giới hạn biên định trước ln tính tốn Do đó, (3) sử dụng để tối ưu hóa σ 2.2 Phân tích factorial Sợi grating sử dụng sản phẩm công ty Teraxion Đặc tuyến group delaywavelegnth đo đạc hình Thiết bị đo sử dụng model Photonic Dispersion and Loss Analyzer 86038B hãng Agilent Có tổng cộng 800 mẫu đo đạc cho đoạn chiều dài 0.015nm, với sai số đo dạc ±100 fs Hình 3: Giá trị đo đạc đặc tuyến group delay-wavelegnth sợi grating sử dụng Sợi grating có tổng cộng 48 kênh với đặc tuyến kênh gần tương tự Tuy nhiên, thể hình 4, giá trị sai số trễ nhóm trung bình kênh khác biệt dao động tầm rộng từ -25ps đến 10ps Như trình bày phần trước, giá trị sai biệt cặp kênh thể qua thơng số χ2 Hình 4: giá trị sai số độ trễ nhóm trung bình kênh Phần mềm tính tốn xác suất quản lý quy trình Minitab dùng để tính tốn phân tích factorial ước lượng ảnh hưởng BW, RF, θ χ2 lên phương sai sai số hướng truyền σ2 Giá trị biên BW 2GHz 4GHz Giá trị biên RF 8GHz đến 10GHz Các giá trị tham khảo từ điều kiện hoạt động [12] cho hệ thống beamfomring cho ứng dụng ultra-wide band Giá tị biên θ 0° 20° tầm xoay mảng anten Cuối cùng, giá trị biên χ2 giá trị χ2 nhỏ lớn cặp kênh sợi grating, 1651.998116 29722.96121 Do hạn chế điều kiện đo đạc, giá trị dặc tuyến hình đo lần Vì vậy, tập điều kiện biên cho phân tích factorial có giá trị Sử dụng cách thiết lập tập điều kiện Yates, 16 lần mô tiến hành để tính tốn σ2 cho 16 tập điều kiện bảng Bảng 1: tập điều kiện mô Factor Run No BW (GHz) RF (GHz) Θ (°) χ2 1651,998 1651,998 10 1651,998 4 10 1651,998 20 1651,998 20 1651,998 10 20 1651,998 10 20 1651,998 29722,961 10 29722,961 11 10 29722,961 12 10 29722,961 13 20 29722,961 14 20 29722,961 15 10 20 29722,961 16 10 20 29722,961 Mức độ ảnh hưởng yếu tố đến σ2 thể hình Yếu tố θ χ2 có ảnh hưởng mạnh đến σ2, thể độ dốc lớn đường biểu diễn Trong đó, RF có ảnh hưởng thấp đến σ2 Do θ có liên hệ trực tiếp đến giá trị độ trễ nhóm đặc tuyến group delay-wavelegnth χ2 liên hệ trực tiếp đến vùng chứa tín hiệu điều chế bên sợi grating, kết hình cho biết sợi grating có mức độ sai lệch giá trị độ trễ nhóm kênh đóng vai trị quan trọng việc giảm sai lệch búp 10 Main Effects Plot for Error variance Data Means BW (GHz) 3,0 RF (GHz) 2,5 2,0 Mean 1,5 1,0 10 θ(°) 3,0 χ2 2,5 2,0 1,5 1,0 Hình 5: ảnh hưởng BW, RF, θ χ2 lên σ2 Sự ảnh hưởng lẫn yếu tố thể hình Mỗi hình nhỏ mối liên hệ qua lại cặp yếu tố Sự cắt đường thể phụ thuộc lẫn Như vậy, thấy từ hình BW RF phụ thuộc vào θ χ2, θ χ2 không phụ thuộc lẫn Kết qua đưa đến hai ý Thứ nhất, kết từ hình khẳng đinh lại tính quan trọng việc chọn grating có mức sai số chênh lệch thấp kết luận từ hình Thứ hai, kết từ hình cho biết khơng thể tách riêng yếu tố tính tốn σ2 Interaction Plot for Error variance Data Means 10 20 1652,0 29723,0 BW (GHz) RF (GHz) θ(°) BW (GHz) RF (GHz) 10 θ(°) 20 χ2 Hình 6: ảnh hưởng liên hệ lân yếu tố lên σ2 Do bị giới hạn số lần đo đạc mô phỏng, yếu tố cần thiết không cần thiết mơ hình tốn σ2 theo BW, RF, θ χ2 khơng phân biệt hình Tuy nhiên, trình bày mục trước, việc tối thiểu hóa phương sai sai số σ2 xem tìm điểm gần độ lệch chuẩn σ Các giá trị mô biến đổi tương ứng mơ hình tốn σ theo BW, RF, θ χ2 bao gồm yếu tố quan trọng thể hình 11 Normal Plot of the Standardized Effects (response is Error variance, Alpha = 0,05) 99 Effect Type Not Significant Significant 95 F actor A B C D 90 Percent 80 70 60 50 40 30 Name BW (G H z) RF (G Hz) θ(°) χ2 20 10 -3 -2 -1 Hình 7: đồ thị phân bố xác suất Gauss mức độ ảnh hưởng đến σ2 Normal Plot of the Standardized Effects (response is Standard derivation, Alpha = 0,05) 99 95 D 90 Percent 70 60 50 40 30 AD ACD A B AB ABC ABD ABC D N ame BW (G Hz) RF (G Hz) θ(°) χ2 AC 20 BD 10 BCD F actor A B C D C 80 Effect Ty pe Not Significant Significant -50 -25 25 Standardized Effect 50 75 Hình 8: đồ thị phân bố xác suất Gauss mức độ ảnh hưởng đến σ Mơ hình tốn σ theo BW, RF, θ χ2 sau loại bỏ yếu tố không quan trọng mối liên hệ RF-θ θ-χ2 sau: σ = 1.2483 − 0.034 BW + 0.0514 RF + 0.2841θ + 0.3084 χ +0.1489 BW RF − 0.1306 BW θ − 0.084 BW χ (5) −0.1544 RF χ + 0.1427 BW RF θ + 0.1396 BW RF χ −0.1092 BW θ χ − 0.1712 RF θ χ + 0.1381BW RF θ χ Sử dụng phương trình (5), kết tối ưu hóa giá trị bên vế phải để đưa vế trái thể hình Giá trị hình thể đạt giá trị sai số khơng gian giá trị yếu tố có biên định Tuy nhiên, giá trị gần phương sai σ=0.5544 đạt với điều kiện BW=2 GHz, RF=8 GHz, θ=0° χ2=1651.9981 Điều kiện đưa vào chương trình mơ để tính tốn Kết tính tốn mơ σ2=0.31 σ=0.5567 Do có tương ứng kết mơ hình tốn mơ phỏng, mơ hình (5) chấp nhận 12 Kết luận kết Để đạt giá trị sai số hướng truyền nhỏ nhất, hệ thống [12] cần hoạt động BW=2 GHz, RF=8 GHz, θ=0° χ2=1651.9981 Chọn sợi grating có chênh lệch sai số cặp kênh truyền tối thiểu tối thiểu giá trị sai số hướng truyền Mơ hình tốn (5) sử dụng để ước lượng nhanh chóng giá trị độ lệch chuẩn sai số hướng truyền điểm thuộc miền giá trị quy định biên BW, RF, θ χ2 chọn ban đầu Các mô tiêu tốn thời gian bỏ qua Phương thức sử dụng đề tài ứng dụng cho hệ thống beamforming khác sử dụng gating để xác lập mơ hình tốn giảm thiểu việc đo đạc-mơ Tp.HCM, ngày 03 tháng 03 năm 2014 Chủ nhiệm đề tài (Ký ghi rõ họ tên) Tp.HCM, ngày tháng năm TL HIỆU TRƯỞNG 13 Tài liệu tham khảo [1] K Barton, Radar Technology Encyclopedia: Artech House Publishers, 1997 [2] R J Mailloux, “Phased-Array Theory And Technology,” Proceedings of the IEEE, vol 70, no 3, pp 246-291, 1982 [3] N A Riza, Selected papers on photonic control systems for phased array antennas, SPIE Optical Engineering Press, 1997 [4] H Subbaraman, M Yihong, and R T Chen, “Photonic Crystal Fiber-Based TrueTime-Delay Beamformer for Multiple RF Beam Transmission and Reception of an XBand Phased-Array Antenna,” Journal of Lightwave Technology, vol 26, no 13-16, pp 2803-2809, Jul-Aug, 2008 [5] V Kaman, X Z Zheng, R J Helkey, C Pusarla, and J E Bowers, “A 32-element 8bit photonic true-time-delay system based on a 288 x 288 3-D MEMS optical switch,” IEEE Photonics Technology Letters, vol 15, no 6, pp 849-851, Jun, 2003 [6] J Rabb, B L Anderson, W D Cowan, and O B Spahn, “Spherical Fourier Cell and Application for Optical True time delay,” Journal of Lightwave Technology, vol 27, no 5-8, pp 879-886, Mar-Apr, 2009 [7] M Rajabvand, and F Behnia, “Fiber Bragg grating-based band edge delay line for radio frequency (RF) photonics applications,” Aeu-International Journal of Electronics and Communications, vol 64, no 3, pp 267-272, 2010 [8] B Hunter, M E Parker, and J L Dexter, “Demonstration of a continuously variable true time delay beamformer using a multichannel chirped fiber grating,” IEEE Transactions on Microwave Theory and Techniques, vol 54, no 2, pp 861-867, Feb, 2006 [9] B M Jung, and H P Yao, “A Two-Dimensional Optical True Time-Delay Beamformer Consisting of a Fiber Bragg Grating Prism and Switch-Based FiberOptic Delay Lines,” IEEE Photonics Technology Letters, vol 21, no 9-12, pp 627629, May-Jun, 2009 [10] Y T Dai, and J P Yao, “Multi-channel Bragg gratings based on nonuniform amplitude-only sampling,” Optics Express, vol 16, no 15, pp 11216-11223, Jul, 2008 [11] Pham Q Thai, Arokiaswami Alphones, and Desmond R Lim, “Limitations by Group Delay Ripple on Optical Beam-forming with Chirped Fiber Grating”, IEEE Journal of Lightwave Technology, Vol.27, pp – 5619-5625, Dec-2009 [12] Pham Q Thai, Arokiaswami Alphones, and Desmond R Lim, “A Novel Simplified Dual Beam-former using Multi-Channel Chirped Fiber Grating and Tunable Optical Delay Lines”, IEEE Journal of Lightwave Technology: 2008 Microwave Photonics special issue, Vol.13-16, pp -2629-2634, Aug-2008 [13] Douglas C Montgomery, George C Runger, and Norma F Hubele, Engineering Stastictics, Wiley, 2010 14 ... đến sai số hướng truyền cải thiện thông qua chỉnh sửa cải tiến t? ?ơng ứng hệ thống optical beamforming Qua đó, giảm sai số hướng truyền cho hệ thống Các k? ?t nghiên cứu 2.1 Phân t? ?ch cấu trúc hệ thống. .. cách thức giảm sai số hướng truyền cho hệ thống optical beamforming Cụ thể: • Phân t? ?ch mối liên hệ mức độ ảnh hưởng lên sai số hướng truyền yếu t? ?? băng thơng t? ?n hiệu điều chế, t? ??n số t? ?n hiệu... lượng thơng qua mơ hình t? ?n học sai số hướng truyền hàm theo yếu t? ?? phân t? ?ch [11] cho mơ hình [12] đóng góp nhiều m? ?t thi? ?t kế t? ??i ưu hóa hệ thống [12] 1.2 Nội dung nghiên cứu Đề t? ?i t? ??p trung