BÁO cáo CUỐI kì môn THỰC HÀNH ANTEN GIẢN đồ SMITH CHART ANTEN nửa bước SÓNG

ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA ĐIỆN TỬ VIỄN THÔNG BÁO CÁO CUỐI KÌ Mơn: THỰC HÀNH ANTEN Giảng viên: Th.Ngơ Thanh Hãi Sinh viên thực hiện: MSSV: Trương Nhật Dữ 18200004 TP Hồ Chí Minh, 2021 Trương Nhật Dữ -18200004 PHỤ LỤC CHƯƠNG I: ĐƯỜNG DÂY TRUYỀN SÓNG Chương II: GIẢN ĐỒ SMITH CHART Câu 1: Hoàn thành bảng sau CHƯƠNG III: ANTEN NỬA BƯỚC SÓNG Câu 1: Khảo sát đặc tính xạ anten nửa bước sóng Matlab Câu 2: Khảo sát đặc tính xạ anten nửa bước sóng HFSS CHƯƠNG IV: HỆ THỐNG BỨC XẠ Câu 1: Cho hệ thống xạ sau Câu 2: Tương tự câu trục Oz,Ox CHƯƠNG 5: KÊNH TRUYỀN FADING DIỆN RỘNG VÀ DIỆN HẸP Câu 1: Khảo sát mơ hình suy hao đơn giản kết hợp shadowing Câu 2: Khảo sát mơ hình suy hao diện hẹp kênh fading Rayleigh 10 10 12 12 16 22 22 30 46 46 49 NHẬT XÉT 52 2|Page Trương Nhật Dữ -18200004 CHƯƠNG I: ĐƯỜNG DÂY TRUYỀN SĨNG Câu 1: Thiết kế đường dây truyền sóng HFSS với tham số thiết kế sau: a.Thao tác gán vật liệu: Kích song: 3|Page Trương Nhật Dữ -18200004 Ngắn mạch: Chọn tần số mô phỏng: Kiểm tra: 4|Page Trương Nhật Dữ -18200004 Tính tốn trở kháng đặc tính đường dây: Kết mơ ngắn mạch: perfect E Kết mô hở mạch: perfect H 5|Page Trương Nhật Dữ -18200004 Tham số Ngắn mạch 𝒁𝒊𝒏 L 1.22 × 10−9 0.0161 + j1.5378 𝑹𝟎 C Lý thuyết: 𝐿 𝑅0 = √𝑐 = 23.9 Hở mạch 0.52787 –j375.63 2.12 × 10−12 Mơ phỏng: 25.136 b.Tìm tần số mà cổng hưởng xảy đường dây cho trrường hợp hở mạch ngắn mạch Điều kiện để mạch cộng lưởng nối tiếp song song là: 6|Page Trương Nhật Dữ -18200004 - Trường hợp 1: 𝑍𝑖𝑛 = 𝑅0 𝑍𝐿 Nếu tải hở mạch (ZL =+∞) Zin =0, tương đương với ngắn mạch đầu vào Khi đó, trở kháng Zin tương đương với trở kháng mạch LC nối tiếp tần số cộng hưởng 𝑊0 = √𝐿𝐶 Ta gọi mạch cộng hưởng nối tiếp Nếu tải hở mạch (ZL =0) Zin =+∞, tương đương với hở mạch đầu vào Khi đó, trở kháng Zin tương đương với trở kháng mạch LC song song tần số cộng hưởng 𝑊0 = Ta gọi mạch cộng hưởng song song √𝐿𝐶 - Trường hợp 2: 𝑍𝑖𝑛 = 𝑍𝐿 Nếu tải hở mạch (ZL =0) Zin =0, tương đương với ngắn mạch đầu vào Khi đó, trở kháng Zin tương đương với trở kháng mạch LC nối tiếp tần số cộng hưởng 𝑊0 = √𝐿𝐶 Ta gọi mạch cộng hưởng nối tiếp Nếu tải hở mạch (ZL =+∞) Zin =+∞, tương đương với hở mạch đầu vào Khi đó, trở kháng Zin tương đương với trở kháng mạch LC song song tần số cộng hưởng 𝑊0 = √𝐿𝐶 Ta gọi mạch cộng hưởng song song Sơ đồ chung điều kiện cộng hưởng : 7|Page Trương Nhật Dữ -18200004 Lý thuyết: Với l=20mm, giá trị 𝜆 mà xảy cổng hưởng : Trường hợp 1: 𝑙 = (2𝑘 + 1) Trường hợp 2: 𝑙 = 𝑘 𝜆 = 𝜆 𝜆 = 𝜆 (K=0) => 𝜆 =0.08 m => f1 = 2.4GHz (K=0) => 𝜆 =0.04 m => f1 = 5GHz Mô phỏng: Trường hợp perf H, tức mạch hở, tương đương với tụ C có trường hợp ( theo sơ đồ điều kiện cổng hưởng) -𝑙= -𝑙 = 𝜆 𝜆 => f1 = 2.4 GHz => Zin = ,ZC =0 = > phần ảo Zin =0 => f2 =5GHz => Zin = +∞ , ZC=+∞ => phần ảo Zin = +∞ 8|Page Trương Nhật Dữ -18200004 Kiểm tra đồ thị : Tại f1 = 2.4GHz ta thấy Zin = => điều kiện cộng hưởng Tại f2 = GHz ta thấy Zin = +∞ => điều kiện cộng hưởng Trường hợp perf E, tức ngắn mạch, tương đương với tụ L có trường hợp ( theo sơ đồ điều kiện cổng hưởng) -𝑙= -𝑙 = 𝜆 𝜆 => f1 = 2.4 GHz => Zin = +∞,ZL =+∞= > phần ảo Zin =+∞ => f2 =5GHz => Zin = , ZL=0 => phần ảo Zin = Kiểm tra đồ thị: Tại f1 = 2.4GHz ta thấy Zin = +∞ => điều kiện cộng hưởng Tại f2 = GHz ta thấy Zin = => điều kiện cộng hưởng Vậy tần số cộng hưởng đường dây là: f1 = 2.4GHz, f2 = GHz 9|Page Trương Nhật Dữ -18200004 Chương II: GIẢN ĐỒ SMITH CHART Câu 1: Hoàn thành bảng sau: Nhập vào giản đồ Smith trở kháng tải: 𝑍𝐿 = 100 + j25 𝛺 Trở kháng đặc tính đường dây: 𝑍0 = 50 𝛺 Lý thuyết Kết giản đồ Smith Tính hệ số phản xạ tải 𝑍𝐿 − 𝑍0 13 𝛤𝑙 = = + 𝑗 𝑍𝐿 + 𝑍0 37 37 Tính hệ số sóng đứng VSWR 1+|𝛤| VSWR = 1−|𝛤| = 2.162 Tính hệ số suy hao phản hồi RL 𝑅𝐿 = −20 log|𝛤| = 8.692 Tìm trở kháng điểm cách tải 0.2𝜆 100+𝑗25+𝑗50tan(72) 𝑍(𝑥) = 5050+𝑗(100+𝑗25)tan(72) Tìm hệ số phản xạ điểm cách tải 0.2𝜆 2𝜋 𝛤(𝑥) = 𝛤𝑙 𝑒 −𝑗2 𝜆 0.2𝜆 = 0.368𝑒 𝑗17.1 𝑒 −𝑗0.8𝜋 = 0.368𝑒 −𝑗126.9 10 | P a g e Trương Nhật Dữ -18200004 +Trục Z-: b.Trục Ox: b.1.Độ lớn trường Anten tham chiếu mô hình xạ 3D, mặt cắt X-Y, Y-Z, X-Z Độ lớn trường Anten tham chiếu (biểu thức) 𝜋 𝜂 |sin ( 𝑐𝑜𝑠(𝜃)) tan(𝜃)| , 𝐹(𝜃, 𝜙) = {2𝜋 𝜂 , Đồ thị 3D clc; clear; close all; 𝜋 ; 𝜋 𝜃= 𝜃 ≠ 38 | P a g e Trương Nhật Dữ -18200004 eta = 120*pi; %============== 3D ================== theta1 = 0:pi/100:pi/2-pi/100; phi1 = 0:2*pi/100:2*pi; [phi1,theta1] = meshgrid(phi1,theta1); r1 = 1/(2*pi)*eta.*(sin((pi/2) *cos(theta1)).*tan(theta1)); [x1,y1,z1] = sph2cart(phi1,pi/2 - theta1,r1); theta2 = pi/2+pi/100:pi/100:2*pi; phi2 = 0:2*pi/100:2*pi; [phi2,theta2] = meshgrid(phi2,theta2); r2 = 1/(2*pi)*eta.*(sin((pi/2) *cos(theta2)).*tan(theta2)); [x2,y2,z2] = sph2cart(phi2,pi/2 - theta2,r2); theta3 = pi/2; phi3 = 0:2*pi/100:2*pi; [phi3,theta3] = meshgrid(phi3,theta3); r3 = eta/4; [x3,y3,z3] = sph2cart(phi3,pi/2 - theta3,r3); figure(1) surf(x1,y1,z1); hold on; surf(x2,y2,z2); plot(x3,y3,'r-','linewidth',1.5) xlabel('x'); ylabel('y'); zlabel('z') grid on Mặt cắt X-Y %== 2D, quét phi, theta = pi/2 == X-Y theta = pi/2; phi = 0:2*pi/100:2*pi; r = eta/4 * ones(1,length(phi)); figure(3) h = polarplot(phi,r,'r','linewidth',1.5); ax = ancestor(h, 'polaraxes'); ax.ThetaZeroLocation = 'top'; ax.ThetaDir = 'clockwise'; Mặt cắt Y-Z %== 2D, phi=pi/2, quét theta == Y-Z phi = pi/2; theta1 = 0:pi/100:pi; r1 = zeros(1,length(theta1)); for i = 1:length(theta1) if theta1(i) == pi/2 r1(i) = eta/4; else 39 | P a g e Trương Nhật Dữ -18200004 r1(i) = eta/(2*pi) * abs(sin(pi*cos(theta1(i))/2) * tan(theta1(i))); end end theta2 = pi:pi/100:2*pi; r2 = flip(r1); theta = [theta1 theta2]; r = [r1 r2]; figure(2) h = polarplot(theta,r,'r','linewidth',1.5); ax = ancestor(h, 'polaraxes'); ax.ThetaZeroLocation = 'top'; ax.ThetaDir = 'clockwise'; Mặt cắt X-Z %== 2D, phi=0,pi quét theta == X-Z phi1 = 0; theta1 = 0:pi/100:pi; r1 = zeros(1,length(theta1)); for i = 1:length(theta1) if theta1(i) == pi/2 r1(i) = eta/4; else r1(i) = eta/(2*pi) * abs(sin(pi*cos(theta1(i))/2) * tan(theta1(i))); end end phi2 = pi; theta2 = pi:pi/100:2*pi; r2 = zeros(1,length(theta2)); for i = 1:length(theta2) if theta2(i) == pi/2 r2(i) = eta/4; else r2(i) = eta/(2*pi) * abs(sin(pi*cos(theta2(i))/2) * tan(theta2(i))); end end theta = [theta1 theta2]; r = [r1 r2]; figure(4) h = polarplot(theta,r,'r','linewidth',1.5); ax = ancestor(h, 'polaraxes'); ax.ThetaZeroLocation = 'top'; ax.ThetaDir = 'clockwise'; Ta có : 𝐽𝑟𝑒𝑓 (𝑟̅ ) = 𝑍̂ 40 | P a g e Trương Nhật Dữ -18200004 𝑙 𝑒 −𝑗𝐾𝑛 ̅′ ) ̅̅̅̅̅̅ ̅ (𝑟̅ ) 𝑒 −𝑗𝐾(𝑟̂ 𝑟𝑟𝑒𝑓 𝐴𝑟𝑒𝑓 (𝑟) = ∫ 𝐽𝑟𝑒𝑓 𝑑𝑙 4𝜋𝑟 − 𝑙 ′ → 𝐾 𝑟̂ 𝑟̅𝑟𝑒𝑓 = 2𝜋 𝑟̂ 𝑧̂𝑧 ′ ≈ 0(𝑣ì 𝑙 ≪ 𝜆) 𝜆 𝑒 −𝑗𝐾𝑟 𝑙 𝑧̂ 4𝜋𝑟 𝐸̅𝑟𝑒𝑓 (𝑟̅ ) = −𝑗𝜔𝜇[𝐴̅𝑟𝑒𝑓 (𝑟̅ ) − (𝐴̅𝑟𝑒𝑓 (𝑟̅ ) 𝑟̂ ) 𝑟̂ ] → 𝐴̅𝑟𝑒𝑓 (𝑟̅ ) = 𝑒 −𝑗𝐾𝑟 (−𝑗𝜔𝜋𝑙)[𝑧̂ − (𝑧̂ 𝑟̂ ) 𝑟̂ ] = 4𝜋𝑟 𝑧̂ = cos 𝜃 𝑟̅ − sin𝜃 𝜃̂ 𝑧̂ 𝑟̂ = 𝑐𝑜𝑠𝜃 𝑧̂ − (𝑧̂ 𝑟̂ ) 𝑟̂ = −𝑠𝑖𝑛𝜃 𝜃̂ → 𝐸̅𝑟𝑒𝑓 (𝑟̅ ) = 𝑒 −𝑗𝐾𝑟 𝑗𝜔𝜇𝑙 𝑠𝑖𝑛𝜃 𝜃̂ 4𝜋𝑟 b.2.Phân tích tính tốn hệ số xếp cho mảng Anten trên, ARFAC? 𝑁−1 𝑁−1 𝐴𝑅𝐹𝐴𝐶 = ∑ 𝐼𝑖 𝑒 𝑗𝑘𝑟̂ 𝑟̅ = ∑ 𝐼0 𝑒 𝑗𝑖𝛽 𝑒 𝑗𝑘𝑖𝑑𝑠𝑖𝑛𝜃𝑠𝑖𝑛𝜙 𝑖=0 𝑖=0 (𝑟̂ 𝑟̅ = 𝑟̂ 𝑥̂ 𝑖𝑑 = (𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙 𝑥̂ + 𝑠𝑖𝑛𝜃𝑠𝑖𝑛𝜙 𝑦̂ + 𝑐𝑜𝑠𝜃 𝑧̂ ) 𝑥̂ 𝑖𝑑 = 𝑖𝑑𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙) 𝑁−1 𝐴𝑅𝐹𝐴𝐶 = 𝐼0 ∑[𝑒 𝑗(𝛽+𝑘𝑑𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙) ]𝑖 𝑖=0 Đặt : 𝛹 = 𝛽 + 𝑘𝑑𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙 𝑁−1 𝐴𝑅𝐹𝐴𝐶 = 𝐼0 ∑[𝑒 𝑗𝛹 ]𝑖 𝑖=0 = 𝐼0 1− 𝑒 𝑗𝛹𝑁 1− 𝑒 𝑗𝛹 41 | P a g e Trương Nhật Dữ -18200004 𝛹𝑁 = 𝛹 𝑠𝑖𝑛 𝐼0 𝑒 𝑗 (𝑁−1) 𝛹2 𝑠𝑖𝑛 b.3.Tìm điều kiện β để mảng anten xạ broadside endfire + Hướng Broadside : 𝑦+ ∶ 𝜃 = 𝑦− ∶ 𝜃 = 𝜋 𝜋 ,𝜙 = 𝜋 ,𝜙 = − 𝜋 ARFAC phải cực đại theo hướng Broadside  𝛹 = → 𝛽 + 𝑘𝑑𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙 = 𝛽=0 + Hướng Endfire : 𝑥+ ∶ 𝜃 = 𝜋 ,𝜙 = 𝜋 𝑥− ∶ 𝜃 = , 𝜙 = 𝜋 ARRFAC phải cực đại theo hướng Endfire  𝛹 = (y+) : 𝛽 + 𝑘𝑑𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙 = → 𝛽 = −𝑘 ∗ 𝑑 (y–) : 𝛽 + 𝑘𝑑𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜙 = → 𝛽 =𝑘∗𝑑 b.4 Khảo sát Matlab độ lớn hệ số ARFAC theo góc ϕ cho hai trường hợp broadsize endfire với giá trị d là: d = 0.1λ, d = 0.2λ, d = 0.25λ, d = 0.35λ, d = 0.5λ d = 0.8λ Chọn giá trị d phù hợp 𝐴𝑅𝐹𝐴𝐶 = 𝑁Ѱ ) , Ѱ = 𝛽 + 𝑘𝑑𝑠𝑖𝑛(𝜃)𝑐𝑜𝑠(𝜙) Ѱ sin( ) sin( 𝑗Ѱ 𝐼𝑂 𝑒 (𝑁−1) + Hướng Broadside: clc; clear; close all; N=4; I0=1; lamda = 0.3; d=0.35*lamda; k =2*pi/lamda; beta =0; phi = 0:2*pi/1000:2*pi; theta =pi/2; psi = beta + k*d.*sin(theta).*cos(phi); r = I0.*abs(sin(N.*(psi/2))./sin(psi/2)); figure(1) h = polarplot(phi,r,'r-','linewidth',1.5); 42 | P a g e Trương Nhật Dữ -18200004 ax = ancestor(h, 'polaraxes'); ax.ThetaZeroLocation = 'top'; ax.ThetaDir = 'clockwise'; + Hướng Endfire: Theo x-: clc; clear; close all; N = 4; I0 = 1; lamda = 0.3; d = 0.25*lamda; k = 2*pi/lamda; beta = k*d; phi = 0:2*pi/1000:2*pi; theta = pi/2; psi = beta + k*d.*sin(theta).*cos(phi); r = I0.*abs(sin(N.*(psi/2))./sin(psi/2)); figure(1) h = polarplot(phi,r,'r-','linewidth',1.5); ax = ancestor(h, 'polaraxes'); ax.ThetaZeroLocation = 'top'; ax.ThetaDir = 'clockwise'; Theo y+: clc; clear; close all; N = 4; I0 = 1; lamda = 0.3; d = 0.25*lamda; k = 2*pi/lamda; beta = -k*d; phi = 0:2*pi/1000:2*pi; theta = pi/2; psi = beta + k*d.*sin(theta).*cos(phi); r = I0.*abs(sin(N.*(psi/2))./sin(psi/2)); figure(1) h = polarplot(phi,r,'r-','linewidth',1.5); ax = ancestor(h, 'polaraxes'); ax.ThetaZeroLocation = 'top'; ax.ThetaDir = 'clockwise'; 43 | P a g e Trương Nhật Dữ -18200004 Khảo sát Matlab d=0.3𝜆 thích hợp để thiết kế hệ thống Broadside Endfire Vì d=0.3𝜆 khảo sát Matlab thể rõ tính định hướng Broadside Endfire bước sóng phụ thấp Broadside (x) 44 | P a g e Trương Nhật Dữ -18200004 Endfire(x) 45 | P a g e Trương Nhật Dữ -18200004 CHƯƠNG 5: KÊNH TRUYỀN FADING DIỆN RỘNG VÀ DIỆN HẸP Câu 1: Khảo sát mơ hình suy hao đơn giản kết hợp shadowing Kỹ thuật viên đo đạt suy hao đường truyền (Path-Loss) mơi trường truyền sóng vô tuyến theo khoảng cách cho bảng sau: Tần số phát trạm phát 𝑓 = 900𝑀𝐻𝑧 Cho phép khoảng cách tham chiếu 𝑑𝑜 = 1𝑚 Viết phương trình tổng qt cảu mơ hình suy hao đơn giản (miền tuyến tính miền dB) Trả lời: 𝑑 −𝛼 Tuyến tính: 𝑃𝑟 = 𝑃𝑡 𝐾 ( ) 𝑑 𝑜 𝑑 Giai dB: 10𝑙𝑜𝑔10 𝑃𝑟 = 10𝑙𝑜𝑔10 𝑃𝑡 + 10𝑙𝑜𝑔10 𝐾 − 10𝛼𝑙𝑜𝑔10 ( ) 𝑑 𝑃𝑟 (𝑑𝐵𝑚) = 𝑃𝑡 (𝑑𝐵𝑚) + 𝐾(𝑑𝐵) − 10𝛼𝑙𝑜𝑔10 ( 𝑑 ) 𝑑0 Tính hệ số K từ độ lợi đường môi trường không gian tự khoảng cách tham chiếu Trả lời: 46 | P a g e Trương Nhật Dữ -18200004 Hệ số K: 𝐾(𝑑𝐵) = 20𝑙𝑜𝑔10 ( 𝜆 4𝜋𝑑𝑜 3.108 ) = 20𝑙𝑜𝑔10 (4𝜋.900.106 ) = −31.526(𝑑𝐵) Phương trình mơ hình suy hao đơn giản viết lại là: 𝑃𝑟 (𝑑𝐵) = 𝑃𝑟 (𝑑𝐵𝑚) − 𝑃𝑡 (𝑑𝐵𝑚) = −31.526 − 10𝛼𝑙𝑜𝑔10 𝑑 𝑃𝑡 Khớp liệu đo với mơ hình suy hao đơn gairn phương pháp cực tiểu trung bình bình phương lỗi (MMSE) để tìm hệ số suy hao đường 𝛼 phương sai 𝜎𝜓,𝑑𝐵 cho mơ hình Shadowing Trả lời: Hàm trung bình bình phương lỗi: 𝐽 = [(−60 + 31.526 − 10𝛼𝑙𝑜𝑔10 5)2 + (−80 + 31.526 − 10𝛼𝑙𝑜𝑔10 25)2 + (−105 + 31.526 − 10𝛼𝑙𝑜𝑔10 65)2 + (−115 + 31.526 − 10𝛼𝑙𝑜𝑔10 110)2 + (−135 + 31.526 − 10𝛼𝑙𝑜𝑔10 400)2 + (−150 + 31.526 − 10𝛼𝑙𝑜𝑔10 1000)2 ] Tìm 𝛼 để J nhỏ nhất: 𝛼 ∗ = 𝐴𝑟𝑔𝛼 𝑀𝑖𝑛𝐽 𝜎𝜓𝑖𝑑𝐵 = 𝐸[(đ𝑜 đạ𝑡 − 𝑚ơ ℎì𝑛ℎ)2 ] Code Matlab: Phân tích xác suất dừng OP 47 | P a g e Trương Nhật Dữ -18200004 Code : clc; clear; close all; f = 900*10^6; % Tần số phát c = 3*10^8; lamda = c/f; % Bước sóng phát d0 = 1; % Khoảng cách tham chiếu K = 20*log10(lamda/(4*pi*d0));% Hệ số K alpha = 3.9581; sigmaPsiDb = 9.681; d = 100; Pt = 10; PtDbm = 10*log10(Pt); Pmin = -110:-100; % Ngưỡng dịch vụ % ============= Monte Carlo ============== Ntry = 10^5; OP_Monte = zeros(1,length(Pmin)); for i = 1:length(Pmin) PsiDb = sqrt(sigmaPsiDb)*randn(1,Ntry); % Tạo mẫu theo Gauss PrDbm = PtDbm + K 10.^alpha.*log10(d/d0)-PsiDb; % Công suất thu indexLess = find(PrDbm < Pmin(i)); OP_Monte(i) = length(indexLess)/Ntry; end % ============== Lý thuyết =============== OP_Theo = qfunc((PtDbm+K10*alpha*log10(d/d0)Pmin)/sqrt(sigmaPsiDb)) % =============== Đồ thị ================= figure(1) semilogy(Pmin,OP_Theo,'r','linewidth',1.4); hold on; semilogy(Pmin,OP_Theo,'ko','linewidth',1 4); xlabel('Pmin (dBm)'); ylabel('OP') 48 | P a g e Trương Nhật Dữ -18200004 OP_Theo = Câu 2: Khảo sát mô hình suy hao diện hẹp kênh fading Rayleigh Khảo sát đặc tính phân bố độ lớn hệ số kênh truyền, độ lớn bình phương hệ số kênh truyền, phần thực phần ảo hệ số kênh truyền tn theo phân phối Gaussian có trung bình phương sai 0.5 (code) Phần thực clc; clear; close all; Ntry = 10^6; sigma = 1/2; %OmegaP/2 (OmegaP = 2*sigma) % ======= Phan thuc, Phan ao ======= hI = sqrt(sigma)*randn(1,Ntry); hQ = sqrt(sigma)*randn(1,Ntry); % ============ Bien ============= abs_h = sqrt(hI.^2 + hQ.^2); % ======= Binh phuong bien ====== Phần ảo: abs_h_sq = hI.^2 + hQ.^2; % ==== Monte Carlo Ham mat PDF ==== [fhI,xI] = ksdensity(hI); [fhQ,xQ] = ksdensity(hQ); [fabs_h,x_abs] = ksdensity(abs_h); [fabs_h_sq,x_abs_sq] = ksdensity(abs_h_sq); % ===== Ly thuyet ham mat PDF ==== fhIL = 1/(sqrt(2*pi*sigma))*exp((1/2).*(xI.^2)./sigma); fhQL = 1/(sqrt(2*pi*sigma))*exp((1/2).*(xQ.^2)./sigma); fabs_hL = (x_abs./sigma).*exp((1/2).*(x_abs.^2)./sigma); : 49 | P a g e Trương Nhật Dữ -18200004 fabs_h_sqL = 1/(2.*sigma).*exp(x_abs_sq./(2.*sigma)); % ============ Do thi ============= figure(1) plot(xI,fhI,'ko','linewidth',1.4); hold on; plot(xI,fhIL,'r-','linewidth',1.4); xlabel('h_I'); ylabel('f(h_I)'); legend('Simulation','Theory') figure(2) plot(xQ,fhQ,'ko','linewidth',1.4); hold on; plot(xQ,fhQL,'b-','linewidth',1.4); xlabel('h_Q'); ylabel('f(h_Q)'); legend('Simulation','Theory') figure(3) plot(x_abs,fabs_h,'ko','linewidth',1.4); hold on; plot(x_abs,fabs_hL,'g-','linewidth',1.4); xlabel('|h|'); ylabel('f(|h|)'); legend('Simulation','Theory') figure(4) plot(x_abs_sq,fabs_h_sq,'ko','linewidth',1 4); hold on; plot(x_abs_sq,fabs_h_sqL,'m','linewidth',1.4); xlabel('|h|^2'); ylabel('f(|h|^2)'); legend('Simulation','Theory') Biên độ: Bình phương biên độ: Ngưỡng dịch vụ đầu thu Pmin =−10dBm:10dBm Đặc tính kênh truyền câu a, khảo sát xác suất dừng hệ thống theo lý thuyết mô (Code) clc; clear; close all; sigma = 1/2; PminDbm = -10:1:10; Pmin = 10.^(PminDbm/10)*10^-3; (Hình) 50 | P a g e Trương Nhật Dữ -18200004 % =========== Mo phong ============ Ntry = 10^6; OP_simul = zeros(1,length(Pmin)); for i = 1:length(OP_simul) % ====== Phan thuc, phan ao ======= hI = sqrt(sigma)*randn(1, Ntry); hQ = sqrt(sigma)*randn(1, Ntry); % ===== Cong suat thu tuc thoi ==== P_tt = hI.^2 + hQ.^2; indexLess = find(P_tt < Pmin(i)); OP_simul(i) = length(indexLess)/Ntry; end OP_simul; % ========== Ly thuyet ============ OP_theo = - exp(-1/(2.*sigma).*Pmin); figure(1) plot(PminDbm,OP_simul,'ko','linewidth',1.4) ; hold on; plot(PminDbm,OP_theo,'r-','linewidth',1.4); xlabel('Pmin'); ylabel('OP'); legend('Simulation', 'Theory') 51 | P a g e Trương Nhật Dữ -18200004 NHẬN XÉT Tp.HCM, Ngày……Tháng……Năm… Giảng viên Ngô Thanh Hãi 52 | P a g e ... TRUYỀN SÓNG Chương II: GIẢN ĐỒ SMITH CHART Câu 1: Hoàn thành bảng sau CHƯƠNG III: ANTEN NỬA BƯỚC SĨNG Câu 1: Khảo sát đặc tính xạ anten nửa bước sóng Matlab Câu 2: Khảo sát đặc tính xạ anten nửa bước. .. -18200004 CHƯƠNG III :ANTEN NỬA BƯỚC SÓNG Câu 1: Khảo sát đặc tính xạ anten nửa bước sóng Matlab Xác định vecto cường độ điện trường cường độ từ trường vùng xa anten nửa bước sóng Trả lời: Vecto... -18200004 Chương II: GIẢN ĐỒ SMITH CHART Câu 1: Hoàn thành bảng sau: Nhập vào giản đồ Smith trở kháng tải: