Untitled BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH S K C 0 0 3 9 5 9 CÔNG TRÌNH MÃ SỐ S KC0 0 7 3 2 7 Tp Hồ Chí Minh, i BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM[.]
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH CƠNG TRÌNH 1&.+&3751*751*,0 0é3+1*9ơ3+ặ17ậ&+&ẩ&7+é1*6&$ %207/397/ộ7521*01*758𝛼2 𝛼1 +𝛼2 =1, có nghĩa tín hiệu 𝑥1 phân công với công suất phát lớn tín hiệu 𝑥2 Trong nghiên cứu này, hai liệu 𝑥1 𝑥2 hai liệu khác gởi đến nguồn đích Mảng anten Mảng anten hệ thống gồm nhiều anten đặt gần nhau, cách khoảng cách nhằm tạo chùm xạ cho mảng khác vị trí khác [14] Mảng anten có nhiều loại khác hai chiều, chiều, tuyến tính, khơng tuyến tính…[14] Trong nghiên cứu này, mảng anten xem xét mảng tuyến tính [14]: Khi đó, hệ số xếp (AF: array factor) mảng là: AF ( N , ) sin N (kdcos ) kdcos sin( ) Hình 2a Đồ thị xạ chiều mảng gồm anten (2) Trong đó: N số anten nút nguồn S k=2π/λ: thừa số bước sóng (với λ bước sóng tín hiệu) d: khoảng cách anten(d= λ/2(1-1/N) β: hệ số pha dòng cung cấp (trong nghiên cứu β=-k*d) ϕ: hệ số pha tính theo đơn vị rad mặt phẳng tọa độ cực Theo công thức (2) giá trị AF ( N , ) phụ thuộc vào biến: ϕ(0≤ϕ≤π), số lượng anten N Như vậy, N ϕ thay đổi giá trị AF ( N , ) thay đổi theo Hình 2a Hình 2b đồ thị AF ( N , ) với N=5 tọa độ cực tọa độ cầu Hình 2b Đồ thị xạ chiều mảng gồm anten Sau tín hiệu 𝑥1 𝑥2 ghép theo cơng thức (1), tín hiệu 𝑥𝑆 đưa đến mảng anten có hệ số AF công thức (2) để truyền đi, tín hiệu ngõ nút S có dạng công thức (3) xS _ AF AF ( N , SD ) 1PS x1 2 PS x2 (3) Với AF ( N , SD ) hệ số AF ( N , ) theo hướng nút D Tín hiệu xS-AF truyền môi trường vô tuyến đến nút D Bảo mật lớp vật lý mảng NOMA kết hợp mảng anten Lý thuyết bảo mật lớp vật lý mạng NOMA Tín hiệu từ nguồn phát (nút S) sau qua môi trường vô tuyến đến thu Trong nghiên cứu này, nút đích xem nằm vị trí trường vùng xa mảng anten (khi ta khơng xem xét đến ảnh hưởng anten đến tín hiệu thu nút đích mà xét ảnh hưởng chung với hệ số kênh truyền) [14] Khi đó, tín hiệu nút D nhận là: yD 1 AF ( N , SD ) AF ( N , SD ) Ps x1hSD Ps hSD x2 nSD (4) Trong đó, ℎ𝑆𝐷 hệ số kênh truyền fading Rayleigh S D, AF ( N ,SD ) giá trị mảng anten theo hướng nút D, 𝑛𝑆𝐷 nhiễu cộng D (là biến ngẫu nhiên có phân phối Gauss với giá trị trung bình phương sai 𝑁0 ) Trong kỹ thuật NOMA, nút đích D giải mã liệu 𝑥1 trước (vì 𝑥1 phân bổ công suất cao 𝑥2 ) Sau TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ ………… giải mã thành công 𝑥1 , kỹ thuật khử nhiễu (SIC) sử dụng [15] để tách lấy tín hiệu 𝑥2 (nút đích loại bỏ thành phần 1 AF ( N , SD ) Ps x1hSD khỏi tín hiệu nhận được) Tín hiệu cịn lại liệu 𝑥2 cộng với nhiễu: yD _ x 2 AF ( N , SD ) Ps hSD x2 nSD (5) 𝑆𝑁𝑅𝐸,𝑥2 = 𝐺𝑆𝐸 𝑃𝑠 𝛼2 g𝑆𝐸 𝑁0 (16) Xét kênh truyền mơ hình hệ thống kênh truyền fading Rayleigh, có CDF PDF độ lợi kênh g 𝑖 [16]: 𝐹g𝑖 (𝑥) = − 𝑒 −𝜆𝑖𝑥 Tương tự, nút E thu tín hiệu sau: yE 1 AF ( N , SE ) AF ( N , SE ) Với 𝜆𝑖 = Ps x1hSE Ps hSE x2 nSE (6) Trong đó, ℎ𝑆𝐸 hệ số kênh truyền fading Rayleigh S E, AF ( N ,SE ) giá trị mảng anten theo phương ϕSE (theo phương nút nghe E), 𝑛𝑆𝐸 nhiễu cộng E (là biến ngẫu nhiên có phân phối Gauss với giá trị trung bình phương sai 𝑁0 ) Nút nghe E giải mã tín hiệu 𝑥1 trước Sau giải mã thành công 𝑥1 , nút E sử dụng kỹ thuật SIC để giải mã tín hiệu 𝑥2 sau: yE x AF ( N , SE ) Ps hSE x2 nSE g SD hSD g SE hSE GSD AF ( N , SD ) (10) GSE AF ( N , SE ) (11) Từ công thức (4), (5), (8), (9), (10), (11), ta có tỷ số SNR đạt D 𝑥1 𝑥2 là: 𝑆𝑁𝑅𝐷,𝑥1 = 𝐺𝑆𝐷 𝑃𝑆 𝛼1 g𝑆𝐷 /𝑁0 𝐺𝑆𝐷 𝑃𝑆 𝛼2 gSD /𝑁0 +1 (12) 𝑆𝑁𝑅𝐷,𝑥2 = 𝐺𝑆𝐷 𝑃𝑠 𝛼2 g𝑆𝐷 𝑁0 (13) = 𝐺𝑆𝐷 𝑄𝛼1 g𝑆𝐷 𝐺𝑆𝐷 𝑄𝛼2 g𝑆𝐷 +1 𝑥 𝐶𝐷 𝑖 = log (1 + 𝑆𝑁𝑅𝐷,𝑥𝑖 ) (19) 𝑥 𝐶𝐸 𝑖 = log (1 + 𝑆𝑁𝑅𝐸,𝑥𝑖 ) (20) Với i = 1, Dung lượng kênh truyền bảo mật [16]: 𝑥 Trong 𝑁0 𝐺𝑆𝐸 𝑃𝑆 𝛼1 g𝑆𝐸 /𝑁0 𝐺𝑆𝐸 𝑃𝑆 𝛼2 g𝑆𝐸 /𝑁0 +1 (15) = 𝐺𝑆𝐸 𝑄𝛼1 g𝑆𝐸 𝐺𝑆𝐸 𝑄𝛼2 g𝑆𝐸 +1 (21) Xác suất dừng bảo mật lớp vật lý mạng NOMA kết hợp mảng anten Xác suất dừng bảo mật xác suất mà dung lượng bảo mật 𝐶𝑖 (trong công thức (21)) nhỏ ngưỡng tốc độ 𝐶𝑡ℎ cho trước [17] Các biểu thức đưa nhằm đánh giá hiệu xác suất dừng bảo mật mạng thứ cấp đánh giá khả nghe nút nghe E Như vậy, xác suất dừng bảo mật lớp vật lý liệu 𝑥1 là: 𝑥 𝑥 Thay công thức (16), (17) vào (19) công thức 𝑥1 cho 𝑃𝑜𝑢𝑡 sau: 𝑥 𝑃𝑜𝑢𝑡 = Pr[(log (1 + 𝑆𝑁𝑅𝐷,𝑥1 ) − log (1 + 𝑆𝑁𝑅𝐸,𝑥1 )) < 𝐶𝑡ℎ ] 1+𝑆𝑁𝑅𝐷,𝑥1 = Pr [log ( (14) Trong mơ hình này, ta xem xét trường hợp tệ nút E áp dụng nguyên lý SIC tốt Từ công thức (6), (7), (8) (9), (10), (11), ta có tỷ số SNR đạt nút nghe E để lấy 𝑥1 𝑥2 là: 𝑆𝑁𝑅𝐸,𝑥1 = + = max{𝑥, 0} 𝑥 = 𝐺1 𝑄𝛼2 g 𝑆𝐷 𝑃𝑠 [𝑥]+ 𝑥 𝐶𝑖 = [𝐶𝐷 𝑖 − 𝐶𝐸 𝑖 ] 𝑃𝑜𝑢𝑡 = Pr[𝐶1 < 𝐶𝑡ℎ ] = Pr[(𝐶𝐷 − 𝐶𝐸 ) < 𝐶𝑡ℎ ] (22) Trong công thức (12) (13), Q tỷ số mức cơng suất tín hiệu công suất nhiễu xác định sau: 𝑄= (18) Dung lượng kênh truyền nút đích nút nghe lén: (7) (9) (17) Trong đó: 𝛽 hệ số suy hao kênh truyền (8) 𝛽 𝑑𝑖 𝑓g𝑖 (𝑥) = 𝜆𝑖 𝑒 −𝜆𝑖𝑥 Các đại lượng gSD, gSE, GSD, GSE, cho công thức bên =𝐺𝑆𝐸 𝑄𝛼2 g 𝑆𝐸 1+𝑆𝑁𝑅𝐸,𝑥1 1+𝑆𝑁𝑅𝐷,𝑥1 sau: Đặt 𝜃 = = Pr [( 𝐶𝑡ℎ 1+𝑆𝑁𝑅𝐸,𝑥1 ) < 𝐶𝑡ℎ ] ) < 2𝐶𝑡ℎ ] (23) − Công thức (20) viết lại G Q g GSE Q1 g SE Poutx1 Pr SD SD 1 GSE Q g SE 1 GSDQ g SD (24) TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ ………… Với điều kiện xấu nhất, nút nghe E giải mã tín hiệu 𝑥1 cách tốt (tức nút nghe có loại bỏ nhiễu đồng kênh) Khi đó, xác suất dừng từ công thức (24) viết lại sau: G Q g x1 Pr SD SD 1 GSE Q1 g SE Pout G Q g SD SD G g Q Pr g SE S D SD ( 1)GS E Q G xQ f g SD ( x) Pr g SE S D dx ( 1)GS E Q f g SD ( x) Pr GS D xQ dx G xQ , GS D xQ dx f gSD ( x) Pr g SE S D ( 1)GS E Q GSD g SD g SE 1 GSE GSDQ g SD 1 f gSD ( x) Pr , dx G Q SE GSD x 1 GSE GSD Q g SD 1 1 GSE Q1 SD e SD x Pr GS D xQ dx (25) ∞ SD e (26) ∞ Poutx1 f gSD ( x) Pr x dx 1 GSDQ1 GSD g SD g SE 1 G G Q g 1 1 G Q , SD SE SD SE dx f gSD ( x) Pr x 1 GSDQ1 e SE GS D xQ ( 1) GS E g SE Q (27) 𝜃 (28) 𝛼1 − 𝜃𝛼2 𝛼1 𝛼2 (𝜃 + 1)𝑄𝐺𝑆𝐸 𝛼2 (𝜃 + 1)𝐺𝑆𝐸 𝑄(𝑄𝐺𝑆𝐷 𝛼2 𝑥 + 1) Tương tự công thức (22), xác suất dừng bảo mật dx n = − 𝑒 −𝜆𝑆𝐷 𝑛 + 𝜆𝑆𝐷 Với: 𝑛 = 𝐺𝑆𝐷 𝑄(𝛼1 − 𝜃𝛼2 ) ∞ −𝜆𝑆𝐷 𝑚 𝑒 + ∫𝑚 𝜆𝑆𝐷 𝑒 −𝜆𝑆𝐷𝑥 𝑒 −𝜆𝑆𝐸𝐴 𝑑𝑥 − SD x (30) Nếu 𝛼1 > 𝜃𝛼2 : 𝐴= SD Nên: 𝑃𝑜𝑢𝑡 = ∫0 𝑓g𝑆𝐷 (𝑥) 1𝑑𝑥 + ∫0 𝑓g𝑆𝐷 (𝑥) 0𝑑𝑥 =1+0=1 Với G xQ SD e x Pr g SE S D dx ( 1)GS E Q GSD g SD g SE 1 G G Q g 1 1 G Q , SE SD SE SD 0 Pr GSD x 1 GSE GSD Q g SD 1 1 GSE Q1 𝑥 SD e SD x 1 SD e SD n dx Pr[(𝛼1 − 𝜃𝛼2 )𝑥𝐺𝑆𝐷 𝑄 ≤ 𝜃] = 𝑃𝑜𝑢𝑡 =1− Q Nếu 𝛼1 ≤ 𝜃𝛼2 thì: 𝑉ớ 𝑚 = Thay cơng thức (19) (20) vào (29), được: Pr GS D g SDQ ( 1)GS E g SE Q GSD x 1 G G Q g 1 SD SE SD dx f gSD ( x) Pr 1 GSE Q1 𝑥 𝑥 P Pr 1 GS D g SDQ 2Ct h (1 GS E g SE Q ) GSD x g SE 1 G G Q g 1 SD SE SD dx f g SD ( x) Pr G Q SE 𝑥 𝑥 = Pr[𝐶𝐷2 − 𝐶𝐸2 < 𝐶𝑡ℎ ] 𝑃𝑜𝑢𝑡 (29) x2 out GSD g SD g SE 1 G G Q g 1 SE SD SD dx Pr G Q SE lớp vật lý liệu 𝑥2 là: 𝜃 𝑄𝐺SD𝛼2 𝑒 −𝜆𝑆𝐷 𝑛 𝐺 𝜆 𝑆𝐷 𝑆𝐸 𝜆𝑆𝐷 +(𝜃+1)𝐺 𝑆𝐸 Ta biết rằng: xác suất dừng nằm khoảng từ đến Giá trị Pout nhỏ (tiến gần 0) tính bảo mật hệ 𝑥 𝑥 thống tốt (xác suất để xảy (𝐶𝐷 − 𝐶𝐸 < 𝐶𝑡ℎ ) nhỏ, 𝑥1 𝑥1 tức xác suất (𝐶𝐷 − 𝐶𝐸 > 𝐶𝑡ℎ ) lớn) Ngược lại, Pout tiến đến tính bảo mật hệ thống (xác 𝑥 𝑥 suất để xảy (𝐶𝐷 − 𝐶𝐸 < 𝐶𝑡ℎ ) lớn, tức xác suất 𝑥1 𝑥1 (𝐶𝐷 − 𝐶𝐸 > 𝐶𝑡ℎ ) nhỏ) Công thức (28) cho thấy: xác suất tỷ lệ nghịch với hệ số 𝛼1 , 𝐺𝑆𝐷 , 𝑄; ngược lại xác suất dừng tỷ lệ thuận với hệ số 𝐺𝑆𝐸 Mặt khác, hệ số 𝐺𝑆𝐷 tỷ lệ thuận với bình phương số anten nút S Do số anten nút nguồn tăng lên xác suất dừng giảm Trong hệ số 𝐺𝑆𝐸 phụ thuộc vào hai thơng số vị trí nút nghe số lượng anten nút nguồn Trong mơ hình này, khoảng cách từ nút đích đến nút nguồn từ nút E đến nút nguồn đền nên hệ số 𝜆𝑆𝐷 1, nên không ảnh hưởng đến giá trị xác suất dừng 𝑥2 Công thức (30) cho ta thấy: xác suất dừng 𝑃𝑜𝑢𝑡 tỷ lệ nghịch TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ ………… 𝑥 tỷ lệ thuận với hệ số với hệ số 𝛼2 , 𝐺𝑆𝐷 , 𝑄 Hệ số 𝑃𝑜𝑢𝑡 𝐺𝑆𝐸 Kết mô đánh giá Các kết mô thực phần mềm MatLab theo phương pháp Monte-Carlo với số phép thử 10000 lần Công thức (28) công thức (30) xác suất dừng bảo mật lớp vật lý theo lý thuyết 𝑥1 𝑥2 Đánh giá theo hệ số anpha: Hình thể xác suất dừng bảo mật lớp vật lý theo hệ số phân chia công suất tín hiệu 𝑥1 𝑥2 Các thơng số mơ phỏng: số anten S N=5; anten D =1, E =1; Q=1và 5dB; Rt=0.2; GSE=1(hay ϕSE=800); λSD=1, λSE=1 Hình Xác suất dừng theo hệ số phân chia cơng suất 𝛼1 với GSE=1 Kết qủa Hình thể lý thuyết mô phù hợp với Với hình này, hệ số 𝛼1 tăng lên xác suất dừng 𝑥1 giảm 𝑥2 tăng lên Với tỷ lệ phân chia công suất 𝛼1 , Q tăng lên xác suất dừng 𝑥1 tăng 𝑥2 giảm Điều giải thích rằng, cơng suất tín tăng, nút nghe nghe tín hiệu tốt dẫn đến xác suất dừng tăng (bảo mật đi) Khi hệ số GSE giảm xuống (GSE =0.1 hay ϕSE =620), cịn thơng số khác giữ ngun, kết mơ cho Hình Với Hình 4, xác suất dừng 𝑥1 𝑥2 giảm 𝛼1 tăng (giống Hình 3) Khi Q tăng lên, xác suất dừng 𝑥2 giảm(giống Hình 3), cịn 𝑥1 giảm (khác so với Hình 3) Hình Xác suất dừng theo hệ số phân chia công suất 𝛼1 với GSE=0.1 Đánh giá theo số lượng anten Hình thể xác suất dừng 𝑥1 𝑥2 số lượng anten nút S thay đổi Các thơng số mơ cịn lại: Q=1dB 5dB; λSD=1, λSE=1; số anten D E 1; Rt=0.2; GSE=0.1 Hình Xác suất dừng theo số lượng anten nút S Kết Hình rằng, số lượng anten mảng tăng lên xác suất dừng bảo mật lớp vật lý 𝑥1 𝑥2 giảm Ngoài ra, tăng hệ số Q_dB (với số lượng anten) xác suất dừng 𝑥1 𝑥2 giảm Khi số lượng anten tăng lên đồng nghĩa với nhiều vấn đề phát sinh (nhưng yếu tố không xét nghiên cứu này) Đánh giá theo hệ số SNR (Q(dB)) Hình thể xác suất dừng 𝑥1 𝑥2 hệ số Q_dB (mức cơng suất tín hiệu) thay đổi S thay đổi Các thơng số mơ cịn lại: số anten S N=5; số anten D =1, E =1; GSE=0.1(hay ϕSE =620) 1; λSD=1, λSE=1; Rt=0.2 Kết Hình chứng tỏ xác suất dừng tín hiệu 𝑥2 có phù hợp lý thuyết mơ Với tín hiệu 𝑥1 : xác suất dừng lý thuyết mô có khác Q_dB lớn 5dB (với trường hợp Gse =0.1), trường hợp Gse=1 có khác rỏ Sự khác biệt từ công thức (24) xuống (25) giả sử mẫu số [GSEQα2gSE+1] xem gần Giả sử xem hệ số GSEQα2gSE nhỏ nhiều Tuy nhiên, TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ ………… Q tăng hay GSE lớn giả sử khơng cịn Khi đó, xác suất dừng bảo mật 𝑥1 theo lý thuyết lớn theo mơ Hình Xác suất dừng theo hệ số Q(dB) Kết từ Hình cịn giúp ta thấy rằng, GSE lớn (nút nghe E gần nút D) Q_dB tăng nút nghe giải mã tín hiệu từ nút S tốt dẫn đến việc bảo mật giảm (xác suất dừng bảo mật vật lý tăng lên) Đánh giá theo vị trí nút nghe E Trong mô đây, cho vị trí nút E thay đổi tồn mặt phẳng xung quanh điểm S Khi đó, giữ nút E, nút S nút D tạo góc 0≤ϕSE≤1800, kết cho Hình Các thơng số mơ cịn lại: số anten nút S =5; số anten D E 1; Q=5dB; λSD=1, λSE=1; Rt=0.2 Hình Xác suất dừng theo vị trí nút nghe E (đơn vị độ) Kết Hình chứng tỏ xác suất dừng 𝑥1 𝑥2 lý thuyết mô tương đối phù hợp Tuy nhiên, nút E thay đổi vị trí G SE lớn có sử khác biệt lý thuyết mơ tín hiệu 𝑥1 Khi đó, xác suất dừng bảo mật 𝑥1 theo lý thuyết lớn theo mô Kết Hình vị trí khác nút E, xác suất dừng bảo mật khác Đánh giá tỷ lệ lỗi bit BER hệ thống Trong nghiên cứu này, tác giả cịn mơ xác suất lỗi BER (Bit Error Ratio) tín hiệu Trong mơ tỷ lệ BER đây, giá trị GSE=0.1, α1 =0.6 Phương pháp Monte-Carlo sử dụng với tổng số bit truyền trường hợp 106 bit Bộ giải mã SIC sử dụng để giải mã tín hiệu nút D nút E Tín hiệu 𝑥2 phụ thuộc vào việc giải mã tín hiệu x1 Kết cho Hình Hình Hình Xác xuất BER hệ thống theo Q(dB) với tách tín hiệu lý tưởng nút D nút E (với anten nút S=2) Hình Xác xuất BER hệ thống theo Q(dB) với tách tín hiệu lý tưởng nút D nút E (với anten nút S=4) Trong Hình Hình 9, Q tăng lên BER 𝑥1 𝑥2 giảm (cả nút S nút E) BER Hình 𝑥1 𝑥2 nút D tốt so với BER Hình BER nút E 𝑥1 𝑥2 hai Hình thay đổ Điều chứng tỏ rằng, số anten nút S tăng lên tính bảo mật tốt Kết luận Trong nghiên cứu này, đưa mơ hình truyền thơng vơ tuyến với kỹ thuật NOMA kết hợp với mảng anten đánh giá vấn đề bảo mật lớp vật lý mơ hình Chúng tơi phân tích từ lý thuyết đến dạng công thức đơn giản để thể mối quan hệ đại lượng Cuối cùng, tiến hành mô phỏng, đánh giá kết đạt được, để từ thấy đắn lý thuyết đề chưa đạt đưa giả thuyết để đơn giản tính tốn TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ ………… Nghiên cứu tài liệu tham khảo cho hướng nghiên cứu vấn đề liên quan đến bảo mật lớp vật lý mảng anten Tài liệu tham khảo [16] J G Andrews et al., "What Will 5G Be?," IEEE Journal on Selected Areas in Communications, vol 32, no 6, pp 1065-1082, June 2014 [17] Y.-L Tseng, ‘‘LTE-advanced enhancement for vehicular communication,’’ IEEE Wireless Commun., vol 22, no 6, pp 4–7, Dec 2015 [18] O Maraqa, A S Rajasekaran, S Al-Ahmadi, H Yanikomeroglu and S M Sait, "A Survey of Rate-Optimal Power Domain NOMA With Enabling Technologies of Future Wireless Networks," in IEEE Communications Surveys & Tutorials, vol 22, no 4, pp 2192-2235, Fourthquarter 2020 [19] A S Marcano and H L Christiansen, “Performance of NonOrthogonal Multiple Access (NOMA) in mmWave wireless communications for 5G networks”, 2017 International Conference on Computing, Networking and Communications (ICNC), 2017, pp 969-974 [20] T E A Alharbi, K Z Shen and D K C So, "Full-Duplex Cooperative Non-Orthogonal Multiple Access System With Feasible Successive Interference Cancellation," 2020 IEEE 91st Vehicular Technology Conference (VTC2020-Spring), 2020, pp 16 [21] Y Liu, Z Qin, M Elkashlan, Y Gao and L Hanzo, "Enhancing the Physical Layer Security of Non-Orthogonal Multiple Access in Large-Scale Networks," in IEEE Transactions on Wireless Communications, vol 16, no 3, pp 1656-1672, March 2017 [22] Yuanwei Liu, Zhiguo Ding, Yue Gao, and Maged Elkashlan, Zhijin Qin, “Physical Layer Security for 5G Non-orthogonal Multiple Access in Large-scale Networks”, Queen Mary University of London, London, UK, 2016 [23] L Liu, R Zhang, and K.-C Chua, “Secrecy wireless information and power transfer with MISO beamforming,” IEEE Trans Signal Process., vol 62, pp 1850–1863, Apr 2014 [24] A Mukherjee and A L Swindlehurst, “Robust beamforming for security in mimo wiretap channels with imperfect CSI,” IEEE Trans Signal Process., vol 59, pp 351–361, Jan 2011 [25] Z Sheng, H D Tuan, T Q Duong and H V Poor, "Beamforming Optimization for Physical Layer Security in MISO Wireless Networks," in IEEE Transactions on Signal Processing, vol 66, no 14, pp 3710-3723, 15 July15, 2018 [26] E Yaacoub and M Al-Husseini, "Achieving physical layer security with massive MIMO beamforming," 2017 11th European Conference on Antennas and Propagation (EUCAP), 2017, pp 1753-1757 [27] Ghous, Mujtaba, Ziaul H Abbas, Ahmad K Hassan, Ghulam Abbas, Thar Baker, and Dhiya Al-Jumeily 2021 "Performance Analysis and Beamforming Design of a Secure Cooperative MISONOMA Network" Sensors 21, no 12: 4180 [28] N D Anh and P N Son, "Performance Analysis and Evaluation of Underlay Two-Way Cooperative Networks with NOMA," 2020 5th International Conference on Green Technology and Sustainable Development (GTSD), 2020, pp 103-108 [29] Constantine A Balanis, Antenna theory analysis and design, third edition, Published by John Wiley & Sons, Inc, 2005 [30] B C Jung and S.-W Jeon, S.-H Lee, “Successive Interference Cancellation With Feedback for Random Access Networks”, IEEE Communications Letters., vol 21, no 4, pp 825 - 828, 2017 [31] P N Son, V P Tuan, S Park and H Y Kong, "Closed-form Analysis of a Decode-and-Forward Scheme under Physical Layer Security over General Fading Channels," 2018 5th NAFOSTED Conference on Information and Computer Science (NICS), 2018, pp 1-5 [32] A Hyadi, Z Rezki and M Alouini, "An Overview of Physical Layer Security in Wireless Communication Systems With CSIT Uncertainty," in IEEE Access, vol 4, pp 6121-6132, 2016 (BBT nhận bài: …/…/201 , hoàn tất thủ tục phản biện: …/…/201 ) (The Board of Editors received the paper on …/…./201…, its review was completed on …/…/201…) S K L 0