Định vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWBĐịnh vị trong vô tuyến băng siêu rộng UWB
HỌC VIỆN CƠNG NGHỆ BƯU CHÍNH VIỄN THƠNG - TRẦN MINH TÚ ĐỊNH VỊ TRONG VÔ TUYẾN BĂNG SIÊU RỘNG UWB LUẬN VĂN THẠC SĨ KỸ THUẬT (Theo định hướng ứng dụng) HÀ NỘI – 2019 HỌC VIỆN CÔNG NGHỆ BƯU CHÍNH VIỄN THƠNG - TRẦN MINH TÚ ĐỊNH VỊ TRONG VÔ TUYẾN BĂNG SIÊU RỘNG UWB Chuyên ngành: Kỹ thuật Viễn thông Mã số: 8.52.02.08 LUẬN VĂN THẠC SĨ KỸ THUẬT (Theo định hướng ứng dụng) NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS VŨ VĂN SAN HÀ NỘI – 2019 i LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu riêng Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác Tác giả luận văn ký ghi rõ họ tên Trần Minh Tú ii LỜI CÁM ƠN Tôi xin trân trọng cám ơn thầy cô giáo tạo điều kiện cho môi trường học tập tốt, đồng thời truyền đạt cho vốn kiến thức quý báu, tư khoa học để phục vụ cho q trình học tập cơng tác tơi Tơi xin gửi lời cảm ơn đến bạn lớp Cao học Kỹ thuật viễn thơng M17CQTE01-B khóa 2017- 2019 giúp đỡ suốt thời gian học tập vừa qua Đặc biệt, tơi xin bày tỏ lòng biết ơn sâu sắc đến PGS.TS VŨ VĂN SAN tận tình bảo cho tơi suốt q trình học tập nghiên cứu, giúp tơi có nhận thức đắn kiến thức khoa học, tác phong học tập làm việc, tạo điều kiện thuận lợi để tơi hồn thành luận văn Cuối cùng, tơi xin gửi lời cảm ơn tới gia đình, bạn bè, đồng nghiệp động viên, giúp đỡ trình hồn thành luận văn Hà Nội, ngày tháng năm 2018 Trần Minh Tú iii MỤC LỤC LỜI CAM ĐOAN .i LỜI CÁM ƠN ii MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT v DANH SÁCH BẢNG vii DANH SÁCH CÁC HÌNH viii MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN VỀ HỆ THỐNG VÔ TUYẾN BĂNG SIÊU RỘNG UWB 1.1 Giới thiệu .4 1.2 Định nghĩa tín hiệu hệ thống vô tuyến băng siêu rộng UWB 1.2.1 Độ rộng băng tần phân đoạn 1.2.2 Tín hiệu UWB 10 1.3.Mơ hình hệ thống truyền thơng vơ tuyến băng siêu rộng UWB 11 1.4 Đặc điểm tín hiệu hệ thống UWB 19 1.5 Các lĩnh vực ứng dụng vô tuyến UWB 25 1.6 Kết luận chương 27 CHƯƠNG II: ĐỊNH VỊ TRONG VÔ TUYẾN BĂNG SIÊU RỘNG UWB 28 2.1 Giới thiệu 28 2.2 Ước tính tham số định vị 29 2.3 Phương pháp ước lượng tham số định vị 31 2.3.1 Phương pháp cường độ trường tín hiệu thu RSS .31 2.3.2 Phương pháp góc đến AOA 34 2.3.3 Phương pháp thời gian đến TOA .36 2.3.4 Phương pháp vi sai thời gian đến TDOA 42 2.3.5 Một số phương pháp kết hợp điển hình 44 2.4 Phương pháp định vị trí .45 iv 2.4.1 Định vị trí cầu 46 2.4.2 Định vị trí hyperbolic 47 2.4.3 Ước tính vị trí theo thuật tốn sai số bình phương nhỏ LSE 50 2.5 Giao thức định vị 51 2.6 Kết luận chương 58 CHƯƠNG 3: PHÂN TÍCH, MƠ PHỎNG, ĐÁNH GIÁ HIỆU NĂNG MỘT SỐ PHƯƠNG PHÁP ĐỊNH VỊ ĐIỂN HÌNH TRONG VƠ TUYẾN BĂNG SIÊU RỘNG UWB 59 3.1 Giới thiệu 59 3.2 Mô hình kịch mơ 59 3.3 Phân tích, mơ phỏng, đánh giá hiệu số phương pháp định vị điển hình 61 3.3.1 Tính tốn, mơ biểu diễn tính xác việc ước tính khoảng cách 61 3.3.2 Tính tốn, mơ biểu diễn tính xác việc ước tính góc đến AOA 63 3.3.3 Tính tốn biểu diễn tính xác việc ước tính thời điểm đến TOA 65 3.4 Giải pháp điển hình để cải thiện tính xác định vị 68 3.5 Kết luận chương 73 KẾT LUẬN VÀ KIẾN NGHỊ 74 DANH MỤC TÀI LIỆU THAM KHẢO .75 v DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT Từ viết tắt AOA Nghĩa tiếng Anh Angle of Arrival Additive White Gaussian AWGN Noise BPSK Binary Phase Shift Keying Code Division Multiple CDMA Access CIR Channel Impulse Response CRLB Cramer-Rao Lower Bound Defense Advanced Research DARPA Project Agency DSA Dynamic Spectrum Access DSDirect Sequence Spread UWB Spectrum UWB Equivalent Isotropically EIRP Radiated Power ESD Energy Spectral Densty European telecommunication ETSI Standards Institue Federal Communications FCC Commission FEC Forward Error Correction FH Frequency Hoping FSK Frequency Shift Keying GPS Global Positioning System Global System for Mobile GSM Telecommunications Inverse Discrete Fourier IDFT Transform Institute of Elictrical and IEEE Elictronics Engineers IRImpluse Radio-Ultra Wide UWB Band IR Impulse Radio ISI Inter-Symbol Interference Nghĩa tiếng Việt Góc tới (góc đến) Tạp âm Gausơ trắng cộng Khóa chuyển pha pha hai trạng thái Đa truy nhập phân chia theo mã Đáp ứng xung kim kênh Giới hạn Cramer-Rao Cơ quan phụ trách dự án nghiên cứu cao cấp Quốc phòng (ở Mỹ) Truy cập phổ tần động UWB trải phổ chuỗi trực tiếp Công suất phát xạ đẳng hướng tương đương Mật độ phổ lượng Viện tiêu chuẩn Viễn thông châu Âu Ủy ban truyền thông Liên bang (ở Mỹ) Hiệu chỉnh lỗi trước Nhảy tần Khố dịch tần Hệ thống định vi tồn cầu Hệ thống thơng tin di động tồn cầu Biến đổi Fourier rời rạc nghịch Viện kỹ sư Điện Điện tử Vô tuyến xung kim – Băng siêu rộng Vô tuyến xung kim Giao thoa ký hiệu (nhiễu xuyên ký hiệu) vi LOS LSE MB MCCDMA ML International Telecommunications Union International Telecommunications UnionRadio Sector Line of Sight Least Square Error Multi Band Multi Carrier Code Division Multiple Access Maximum Likelihood LSE Least Square Error MUI NLOS Multi User Interference None Line of Sight Orthogonal Frequency Division Multipex Pulse Amplitude Modulation Power Delay Profile Pseudo Noise Pulse Position Modulation Pricise Positioning System Power Spectral Density Phase Shift Keying Primary User Received Signal Strength Self – Positioning Algorithm Signal to Noise Ratio Standard Positioning Systerm Secondary User Time difference of Arrival Time Hopping Time Hopping Ultra Wide Band Time of Arrival Ultra Wideband Wireless Local Area Network Wireless Personal Access Network ITU ITU-R OFDM PAM PDP PN PPM PPS PSD PSK PU RSS SPA SNR SPS SU TDOA TH THUWB TOA UWB WLAN WPAN Liên minh Viễn thông Quốc tế Liên minh Viễn thông quốc tế phận vô tuyến Trực xạ (đường truyền thẳng) Sai số quân phương nhỏ Đa băng Đa truy nhập phân chia theo mã đa sóng mang Khả giống cực đại (hợp lý cực đại) Sai số bình phương nhỏ (bình phương sai số nhỏ nhất) Nhiễu đa người dùng Không trực xạ (truyền không trực tiếp) Ghép kênh phân chia theo tần số trực giao Điều chế biên độ xung Lý lịch trễ công suất Giả tạp âm Điều chế vị trí xung Hệ thống định vị xác Mật độ phổ cơng suất Khóa dịch pha Người dùng sơ cấp Cường độ trường tín hiệu thu Thuật tốn tự định vị Tỷ số tín hiêu tạp âm Hệ thống định vị tiêu chuẩn Người dùng thứ cấp Vi sai thời gian đến Nhảy thời gian Vô tuyến băng siêu rộng nhảy thời gian Thời gian đến (góc tới) Băng siêu rộng Mạng nội không dây Mạng truy nhập cá nhân không dây vii DANH SÁCH BẢNG Bảng 1: a) Phân bổ kênh UWB – tín hiệu băng rộng; b) Quy hoạch kênh UWB – tín hiệu băng hẹp .14 Bảng 2: Các tham số đặc trưng số hệ thống vô tuyến điển hình .20 Bảng 1: Các tham số kênh mơi trường kênh truyền sóng UWB 61 viii DANH SÁCH CÁC HÌNH Hình 1: Độ rộng băng tần lượng Hình 2: Phổ cơng suất hệ thống UWB phổ công suất hệ thống vơ tuyến hành điển hình Hình 3: Các tới hạn để thiết kế tín hiệu UWB 12 Hình 4: a)Phân bổ băng tần cho kênh UWB (hình bên trên); b) Quy hoạch băng tần cho kênh UWB 500 MHz (hình bên dưới) 13 Hình 5: a) Minh họa mơ hình khái niệm hệ thống UWB dạng đơn giản; b) Minh họa máy phát xung tín hiệu UWB 18 Hình 6: Minh họa khả năng, đặc tính, tính cách hệ thống UWB 22 Hình 1: a) Định vị theo nút trung tâm; b) Định vị tương đối 29 Hình 2: Minh họa ước tính khoảng cách dựa vào RSS khi: a) khơng có sai số (chính xác); b) có sai số, khoảng cách nhận có độ bất định biểu thị vùng màu xám 33 Hình 3: Minh họa ước tính góc đến AOA: a) nút tham chiếu đo xác định góc ψ nút đích; b) quan hệ khác thời điểm đến ψ dàn anten .34 Hình 4: Xung chữ nhật có độ rộng T đầu lọc thích hợp .38 Hình 5: Ảnh hưởng trễ trước lên đầu đồng mở cổng sớm muộn (early-late gate synchronizer ) 40 Hình 6: Quá trình trao đổi bắt tay đo cự ly thiết bị định vị UWB (được đề xuất Fleming and Kushner, 1997) 41 Hình 7: a) Ước lượng TDOA xác định đường hyperbol mà nút đích thơng qua có tiêu cự nút tham chiếu; b) Minh họa ước lượng lai ghép TOA/AOA .43 Hình 8: Minh họa định vị trí cầu nút Ni khơng gian hai chiều với N1, N2, N3 nút tham chiếu .47 Hình 9: Minh họa định vị hyperbolic nút Ni không gian hai chiều với nút tham chiếu N1, N2 N3 .49 77 Phụ lục 1A TĨM TẮT CÁC CHƯƠNG TRÌNH MƠ PHỎNG Sim_UWB 3.01 Tính tốn, mơ biểu diễn tính xác việc ước tính khoảng cách theo biểu thức (2.7) kịch truyền sóng mơ hình kênh IEEE 802.15.4a cho chương trình Prog3.01 Prog 3.01 [result_CRLB] = CRLB_RSS_Distance(d, n1, n2, n3, n4, sigma_sh1,… sigma_sh2, sigma_sh3, sigma_sh4) Chức Đầu vào Thực tính tốn biểu diễn giới hạn độ lệch chuẩn ước tính khoảng cách dựa vào RSS theo biểu thức (2.7) Các đầu vào chương trình gồm ba tham số sh ; n; d Thiết trị cho tham số đặc trưng này: số mũ suy hao độ lệch chuẩn cho bảng 3.1 Đầu Đầu chương trình ước tính khoảng cách theo biểu thức (2.7) Sim_UWB 3.02 Tính tốn, mơ biểu diễn tính xác việc ước tính góc đến AOA theo biểu thức (2.12) khảo sát theo SNR dùng dàn anten loại ULA có phần tử khoảng cách phần tử cm Tín hiệu chọn phương trình (7.8) xung tín hiệu UWB có độ rộng thay đổi, tín hiệu đến máy thu có góc đến / hệ số kênh a = Prog 3.02 [beta] = uwb70201_effective_bandwidth(Tp, f) Chức Đầu vào Đầu Sim_UWB 3.03 Chương trình thực tính băng tần hiệu dụng tín hiệu Chương trinh nhận hai tham số đầu vào là: độ rộng xung tín hiệu, Tp; tần số, f Đầu chương trình băng tần hiệu dụng tín hiệu theo phương trình (2.13) Chương trình tính tốn biểu diễn giới hạn Cramer-Rao sai số ước tính trễ cho biểu thức (2.18), ta ước lượng giá trị tín hiệu UWB Chương trình thực tính tốn biểu diễn tính xác việc ước tính góc đến AOA theo biểu thức (2.12) khảo sát theo thực Sim_UWB 7.03 dàn anten ULA có phần tử khoảng cách phần tử cm Tín hiệu chọn phương trình (2.8) xung tín hiệu UWB có độ rộng thay đổi, tín hiệu đến máy thu có SNR = dB hệ số kênh a = 78 Sim_UWB 3.04 Ước lượng biểu diễn tính xác thuật tốn LSE cho toán định vị định nghĩa hệ phương trình (2.28), hàm cấp độ sai số qua trình đo cự ly số lượng nút tham chiếu Trong Sim_UWB 7.04 này, ta thực phân tích tốn khơng gian hai chiều phép đồ thị hóa kết Sai số mơ hình hóa giá trị xác rút q trình Gausơ trắng có trị trung bình Việc ước lượng thực ba chương trình: Prog3.03 tạo vị trí nút; Prog3.04 chọn nút đích nút tham chiếu; Prog3.05 thực ước lượng vị trí theo thuật toán LSE Prog 3.03 [positions, ranges]=uwb30401_create_network(N, area_side, G) Chức Chương trình tạo vị trí tập nút vùng diện tích hình vng tính tốn khoảng cách cặp nút Đầu vào Chương trình nhận đầu vào gồm: số lượng nút, N; độ dài side vùng diện tích hình vng, area_side; cờ kích hoạt đồ họa đầu ra, G Đầu Chương trình trả đầu ra: (1) ma trận kích thước Nx2 chứa vị trí (X, Y) nút, positons; (2) ma trận kích thước NxN chứa khoảng cách cặp nút, ranges [Nx, Ref] = uwb30402_select_nodes(N, k) Prog 3.04 Chức Đầu vào Đầu Chương trình lựa chọn nút đích k nút tham chiếu từ tập nút Chương trình nhận đầu vào là: tổng số nút, N; số lượng nút tham chiếu, k Chương trình trả đầu ra: (1) số nguyên biểu trưng cho nút đích, Nx; (2) vector độ dài k chứa ID nút tham chiếu, Ref [PosNx, ErrNx] = uwb30403_find_LSE_position(positions, ranges, Nx, Ref, sigma_2, G) Prog 3.05 Chức Chương trình tìm nghiệm LSE cho tốn định vị khơng gian hai chiều Chương trình thực nhiệm vụ bản: (i) đưa sai số vào ước tính khoảng cách với giá trị sai số lấy từ trình Gausơ trung bình khơng; (ii) ước tính vị trí nút đích cách giải hệ phương trình tuyến tính (2.28); (iii) xác định sai số định vị định nghĩa khoảng cách Euclide vị trí ước tính vị trí xác Đầu vào Chương trình nhận đầu vào gồm: (1) ma trận chứa vị trí nút, positions; (2) ma trận chứa khoảng cách cặp nút, ranges; (3) ID nút đích Nx; (4) vector độ dài k chứa ID nút tham chiếu, Ref; (5) giá trị của phương sai tạp âm, sigma_2; (6) cờ cho cho phép đồ thị hóa kết quả, G Đầu Chương trình trả hai đầu ra: (1) ma trận biểu trưng cho vị trí ước tính nút đích, PosNx; (2) giá trị sai số tương ứng với vị trí xác, ErrNx 79 Phụ lục 1B Mà CHƯƠNG TRÌNH MÔ PHỎNG %======================================================================== %============================ Sim_UWB_301================================ %======================================================================== % [result_CRLB] = CRLB_RSS_Distance(d,n1,n2,n3,n4, % sigma_sh1,sigma_sh2,sigma_sh3,sigma_sh4) %======================================================================== % clc; clear; close all; % -% According to Table 3.01 d = 0:0.3:10; % Residential, LOS n1 = 1.79; sigma_sh1 = 2.22; % Residential, NLOS n2 = 4.58; sigma_sh2 = 3.51; % Indoor office, LOS n3 = 1.63; sigma_sh3 = 1.9; % Indoor office, NLOS n4 = 3.07; sigma_sh4 = 3.9; %======================================================================== % minimum SD: Eq(3.7) MSD1 = log(10)*sigma_sh1*d/n1/10; MSD2 = log(10)*sigma_sh2*d/n2/10; MSD3 = log(10)*sigma_sh3*d/n3/10; MSD4 = log(10)*sigma_sh4*d/n4/10; %======================================================================== num_fig = 74; h = figure(num_fig); set(h,'name','Sim_UWB_701, MSD_RSS_Distance'); P1 = plot(d,MSD1,'-or',d,MSD2,'-*m',d,MSD3,'-sb',d,MSD4,'-^k'); set(P1,'LineWidth',2); xlabel('Khoảng cách (m)','fontname','.Vntime','fontsize',14,'color','b'); ylabel('Độ lệch chuẩn nhỏ nhất, (m)','fontname','.Vntime','fontsize',16,'color','b'); title('Giới hạn -ớc tính khoảng cách dựa vào RSS môi tr-ờng truyền sóng khác nhau', 'fontname','.Vntime','color','b','fontsize',16); LS = legend('Trong nhà riêng, LOS (n=1,79; \sigma=2,22)','Trong nhà riêng, NLOS (n=4,58; \sigma=3,51)', 'Trong văn phòng, LOS (n=1,63; \sigma=1,9)','Trong văn phòng, NLOS (n=3,07; \sigma=3,90)'); set(LS,'fontname','.Vntime','fontsize',16); grid on; 80 result_CRLB = [MSD1; MSD2; MSD3; MSD4]; AX = gca; set(AX,'fontsize',14); %======================================================================== %============================ Sim_UWB_302=============================== %======================================================================== % clc; clear; close all; % -% CRLB vs SNR clc; clear; close all; SNR_dB N = -10:0.5:10; = length(SNR_dB); for i=1:N SNR=exp(SNR_dB*log(10)/10); end ONES fc T = ones(1,N); = 1e8; = 10^-9*[0.5 2]; for i=1:3 beta(i) = uwb30201_effective_bandwidth(T(i),fc); % Prog 3.02 end c N_a l phi = = = = 3*10^8; 4; 5*10^-2; pi/4; % Minimum SD: Eq(3.12) MSD1=sqrt(3)*c/(sqrt(2)*pi*beta(1)*sqrt(N_a*(N_a^21))*l*cos(phi))*ONES./(sqrt(SNR)); MSD2=sqrt(3)*c/(sqrt(2)*pi*beta(2)*sqrt(N_a*(N_a^21))*l*cos(phi))*ONES./(sqrt(SNR)); MSD3=sqrt(3)*c/(sqrt(2)*pi*beta(3)*sqrt(N_a*(N_a^21))*l*cos(phi))*ONES./(sqrt(SNR)); % num_fig = 76; h = figure(num_fig); set(h,'name',' Sim_UWB_702, MDS_AOA_SNR'); P1=plot(SNR_dB,MSD1,'-*r',SNR_dB,MSD2,'-sk',SNR_dB,MSD3,'-^b'); set(P1,'LineWidth',1.5); xlabel('SNR (dB) ','fontname','.Vntime','fontsize',14,'color','b'); ylabel('§é lƯch chn nhá nhÊt (rad)','fontname','.Vntime','fontsize',14,'color','b'); title('BiĨu diƠn giíi h¹n cđa -ớc tính AOA theo SNR độ rộng xung cđa tÝn hiƯu UWB kh¸c nhau', 'fontname','.Vntime','color','b','fontsize',14); SL = legend('T = 0.5 ns','T = ns','T = ns'); 81 set(SL,'fontname','.vntime','fontsize',16); grid on; AX = gca; set(AX,'fontsize',14); %%======================================================================= clc; clear all; phi = -1.5:0.05:1.5; N = length(phi); ONES = ones(1,N); fc = 1e8; T = 10^-9*[0.5 2]; for i=1:3 beta(i) = uwb30201_effective_bandwidth(T(i),fc); end c = 3*10^8; N_a = 4; l = 5*10^-2; SNR_dB = 5; SNR = exp(SNR_dB*log(10)/10); % Minimum SD: MSD1=sqrt(3)*c/(sqrt(2)*pi*beta(1)*sqrt(N_a*(N_a^21))*l*sqrt(SNR))*ONES./(cos(phi)); MSD2=sqrt(3)*c/(sqrt(2)*pi*beta(2)*sqrt(N_a*(N_a^21))*l*sqrt(SNR))*ONES./(cos(phi)); MSD3=sqrt(3)*c/(sqrt(2)*pi*beta(3)*sqrt(N_a*(N_a^21))*l*sqrt(SNR))*ONES./(cos(phi)); % h = figure(77); set(h,'name','Sim_UWB_702, MSD_AOA_Angle'); P1 = plot(phi,MSD1,'-*r',phi,MSD2,'-sk',phi,MSD3,'-^b'); set(P1,'LineWidth',1.5); xlabel('Gãc tíi (rad)','fontname','.Vntime','fontsize',14,'color','b'); ylabel('§é lƯch chn nhá nhÊt (rad)','fontname','.Vntime','fontsize',14,'color','b'); title('BiĨu diƠn giíi h¹n d-íi cđa -íc tÝnh AOA theo góc đến SNR = 5dB độ rộng xung cđa tÝn hiƯu UWB kh¸c nhau', 'fontname','.Vntime','color','b','fontsize',16); SL = legend('T = 0.5 ns','T = ns','T = ns'); set(SL,'fontname','.vntime','fontsize',16); axis([ -1.5 1.5 0.5]) grid on AX = gca; set(AX,'fontsize',14); ========================== %======================================================================== %============================ Sim_UWB_303: MSD_TOA_SNR.m ================ %======================================================================== % clc; clear; close all; % -SNR_dB = -10:0.7:10; N = length(SNR_dB); 82 for i=1:N SNR = exp(SNR_dB*log(10)/10); end ONES = ones(1,N); T fc = 10^-9*[0.5 2]; = 1e8; for i=1:3 beta(i)=uwb30201_effective_bandwidth(T(i),fc); % Prog 7.02 end c = 3*10^8; % -MSD1 = c/(2*sqrt(2)*pi*beta(1))*ONES./(sqrt(SNR)); MSD2 = c/(2*sqrt(2)*pi*beta(2))*ONES./(sqrt(SNR)); MSD3 = c/(2*sqrt(2)*pi*beta(3))*ONES./(sqrt(SNR)); % num_fig = 711; h = figure(num_fig); set(h,'name',' Sim_UWB_703, MDS_AOA_SNR'); P1 = plot(SNR_dB,MSD1,'-*r',SNR_dB,MSD2,'-sk',SNR_dB,MSD3,'-^b'); set(P1,'LineWidth',2); xlabel('SNR (dB) ','fontname','.Vntime','fontsize',14,'color','b'); ylabel('§é lƯch chn nhá (m)','fontname','.Vntime','fontsize',14,'color','b'); title('Khảo sát giới hạn -ớc tính khoảng cách theo ph-ơng pháp TOA theo SNR với độ réng xung kh¸c nhau', 'fontname','.Vntime','color','b','fontsize',14); legend('T = 0.5 ns','T = ns','T = ns'); grid on AX = gca; set(AX,'fontsize',14); % -% -Prog 2.01: uwb20101_genrect -% -function [beta] = uwb30201_effective_bandwidth(Tp,f) A = 1; smp = 1000; p = sin(2.*pi.*f.*linspace(0,Tp,smp)); sinpulse = zeros(1,3*smp); % Time window sinpulse(1+smp:2*smp) = p; % % Signal in frequency domian % fs = smp / Tp; % sampling dt = / fs; N = length(sinpulse); T = N * dt; df = / T; % X = fft(sinpulse); X = X/N; ss_X = 2.*X(1:floor(N/2)); n = length(ss_X); f0 = fs/N; % - 83 % Effective_bandwidth s1 = 0; s2 = 0; for i=1:n s1 = s1+(i*df*abs(ss_X(i)))^2*df; s2 = s2+(abs(ss_X(i)))^2*df; end beta = sqrt(s1/s2); %======================================================================= % ============================= Sim_UWB_304_2D ========================= %======================================================================== %======================================================================== % clc; clear; close all; % -% -% Generates the positions of a set of nodes % -N = 10; area_side = 50; G = 1; NumFig = 7161; [positions,ranges] = uwb30401_create_network(N,area_side,G,NumFig); % -% Selects a target node and k reference nodes from a set of nodes % -k = 3; [Nx,Ref] = uwb30402_select_nodes(N,k); % -% Determines the LSE solution to a positioning problem in a bidimensional space % -sigma_2 = 0; G = 1; NumFig = 7162; [PosNx,ErrNx] = uwb30403_find_LSE_position(positions,ranges,Nx,Ref,sigma_2,G,NumFig); %======================================================================== h= figure(716); set(h,'name','LSE solution to a positioning problem in a Bidimensional space (2_D)'); subplot(1,2,1); clc; clear all; load Sim_UWB_2D_H716a.mat; PT = scatter(positions(:,1),positions(:,2),'filled'); axis([0 area_side area_side]); xlabel('X [m]','fontname','.Vntime','fontsize',16); ylabel('Y [m]','fontname','.Vntime','fontsize',16); title('Tạo nút ban đầu', 'fontname','.Vntime','color','b','fontsize',16); box on; AX = gca; set(AX,'fontsize',14); % 84 subplot(1,2,2); clc; clear all; load Sim_UWB_2D_H716b.mat; scatter(positions(:,1),positions(:,2)); xlabel('X [m]','fontname','.Vntime','fontsize',16); ylabel('Y [m]','fontname','.Vntime','fontsize',16); box on; hold on; scatter(PosNx(1), PosNx(2), 200, 'filled', 'k','p'); scatter(positions(Nx,1),positions(Nx,2),200, 'filled','^'); for i=1:k % Reference nodes: 'r','s' scatter(positions(Ref(i),1),positions(Ref(i),2), 'filled','r','s'); end title('¦íc tính vị trí cho nút đích','fontname','.Vntime','color','b','fontsize',16); LT = legend('Nút khác','Vị trí nút đích','Vị trí hiệu dụng nút đích','Nút tham chiÕu'); % set(LT,'fontname','.Vntime','fontsize',14); hold off; AX = gca; set(AX,'fontsize',14); % % -% Prog 3.03: % [positions, ranges]=uwb30401_create_network(N, area_side, G) % % -% function [positions, ranges] = uwb30401_create_network(N,area_side,G,NumFig) for i = 1:N positions(i,1)=rand*area_side; positions(i,2)=rand*area_side; j=1; for j=1:(i-1) ranges(i,j)= sqrt((positions(i,1) - positions(j,1))^2 + (positions(i,2) - positions(j,2))^2); while(ranges(i,j)==0) X(i)=rand*50; Y(i)=rand*50; ranges(i,j)=sqrt((positions(i,1)- positions(j,1))^2+(positions(i,2)- positions(j,2))^2); end ranges(j,i)=ranges(i,j); end end save Sim_UWB_2D_H716a.mat; if G==1 85 h = figure(NumFig); set(h,'name','a) Generates the positions of a set of nodes'); PT = scatter(positions(:,1),positions(:,2),'filled'); axis([0 area_side area_side]); xlabel('X [m]','fontname','.Vntime','fontsize',16); ylabel('Y [m]','fontname','.Vntime','fontsize',16); title('Tạo nút ban đầu', 'fontname','.Vntime','color','b','fontsize',16); box on; end AX = gca; set(AX,'fontsize',14); % % % % % % -Prog 3.04 [Nx, Ref] = uwb70402_select_nodes(N, k) This function selects a target node and k reference nodes from a set of nodes function [Nx,Ref] = uwb30402_select_nodes(N,k) % Extraction of the target node Nx = ceil(N*rand); Check=zeros(1,N); for i=1:k Ref(i) = Nx; % Check: is Ref(i) different from Nx? while((Ref(i)==Nx)||(Ref(i)==0)) Ref(i)= ceil(N*rand); % Check: is Ref(i) already selected? if(Check(Ref(i))) Ref(i) = 0; else Check(Ref(i))=1; end end end % % -% Prog 3.05 % [PosNx, ErrNx] = uwb70403_find_LSE_position(positions,ranges,Nx,Ref,sigma_2,G) function [PosNx, ErrNx] = uwb30403_find_LSE_position(positions, ranges, Nx, Ref, sigma_2, G,NumFig) % Adding errors to the range estimation N = size(ranges,1); err_ranges = ranges + sqrt(sigma_2)*randn(N); 86 % Defining the linear problem % Matrix A k = length(Ref); for i=1:(k-1) A(i,1) = positions(Ref(i),1) - positions(Ref(k),1); A(i,2) = positions(Ref(i),2) - positions(Ref(k),2); end A=-2*A; % Matrix b b=zeros(2,1); for i=1:(k-1) b(i) = err_ranges(Ref(i),Nx)^2 - err_ranges(Ref(k),Nx)^2 - positions(Ref(i),1)^2 + positions(Ref(k),1)^2 - positions(Ref(i),2)^2 + positions(Ref(k),2)^2; end % Solving the problem PosNx=A\b; %Computing the error ErrNx = sqrt((PosNx(1)-positions(Nx,1))^2+(PosNx(2)- positions(Nx,2))^2); save Sim_UWB_2D_H716b.mat; %Graphical output if G==1 h = figure(NumFig); set(h,'name','b) determines the LSE solution to a positioning problem in a bidimensional space'); scatter(positions(:,1),positions(:,2)); xlabel('X [m]','fontname','.Vntime','fontsize',16); ylabel('Y [m]','fontname','.Vntime','fontsize',16); box on; hold on; scatter(PosNx(1), PosNx(2), 200, 'filled', 'k','p'); scatter(positions(Nx,1),positions(Nx,2),200, 'filled','^'); for i=1:k scatter(positions(Ref(i),1),positions(Ref(i),2), 'filled','r','s'); end title('Ước tính vị trí cho nút đích','fontname','.Vntime','color','b','fontsize',16); hold off; end AX = gca; set(AX,'fontsize',14); %======================================================================== % ========================== Sim_UWB_304_3D ============================ %======================================================================== 87 % -% Prog_3D: uwb30403_3D_find_LSE_position % -% clc; clear; close all; % -N area_side G num_fig [positions, = 5; = 10; = 1; = 7171; ranges] = uwb30401_3D_create_network(N,area_side,G,num_fig); % -k = 4; [Nx,Ref] = uwb30402_3D_select_nodes(N,k); % -G = 1; num_fig = 7172; sigma_2 = 0; [PosNx, ErrNx] = uwb30403_3D_find_LSE_position(positions, ranges,Nx, Ref,sigma_2, G,num_fig); %======================================================================== h= figure(717); set(h,'name','LSE solution to a positioning problem in a 3_Dimensional space (3_D)'); subplot(1,2,1); clc; clear all; load Sim_UWB_3D_H717a.mat; scatter3(positions(:,1),positions(:,2),positions(:,3),'filled'); axis([0 area_side area_side area_side]); xlabel('X [m]','fontname','.Vntime','fontsize',16); ylabel('Y [m]','fontname','.Vntime','fontsize',16); zlabel('Z [m]','fontname','.Vntime','fontsize',16); box on; hold on AX = gca; set(AX,'fontsize',12); % -subplot(1,2,2); clc; clear all; load Sim_UWB_3D_H717b.mat; scatter3(positions(:,1),positions(:,2),positions(:,3)); xlabel('X [m]','fontname','.Vntime','fontsize',16); ylabel('Y [m]','fontname','.Vntime','fontsize',16); zlabel('Z [m]','fontname','.Vntime','fontsize',16); box on; hold on; scatter3(PosNx(1), PosNx(2), PosNx(2), 200, 'filled', 'k','p'); scatter3(positions(Nx,1),positions(Nx,2),positions(Nx,3),200, 'filled','^'); for i=1:k 88 scatter3(positions(Ref(i),1),positions(Ref(i),2),positions(Ref(i),3), 'filled','r','s'); end title('Ước tính vị trí cho nút đích','fontname','.Vntime','color','b','fontsize',16); LT = legend('Nút khác','Vị trí nút đích','Vị trí hiệu dụng cđa nót ®Ých','Nót tham chiÕu'); % set(LT,'fontname','.Vntime','fontsize',12); hold off; AX = gca; set(AX,'fontsize',12); % -% uwb30401_3D_create_network function [positions, ranges] = uwb30401_3D_create_network(N,area_side,G,num_fig) for i = 1:N positions(i,1)=rand*area_side; positions(i,2)=rand*area_side; positions(i,3)=rand*area_side; j=1; for j=1:(i-1) ranges(i,j)= sqrt((positions(i,1) - positions(j,1))^2 + (positions(i,2) - positions(j,2))^2 + (positions(i,3) - positions(j,3))^2); while(ranges(i,j)==0) X(i)=rand*50 Y(i)=rand*50 ranges(i,j)=sqrt((positions(i,1)- positions(j,1))^2+(positions(i,2)- positions(j,2))^2+(positions(i,3) - positions(j,3))^2); end ranges(j,i)=ranges(i,j); end end save Sim_UWB_3D_H717a.mat; if G h = figure(num_fig); set(h,'name','H7_18: Create_nodes_3D'); scatter3(positions(:,1),positions(:,2),positions(:,3),'filled'); axis([0 area_side area_side area_side]); xlabel('X [m]'); ylabel('Y [m]'); zlabel('Z [m]'); box on; hold on end AX = gca; set(AX,'fontsize',12); 89 % -% Prog: uwb30402_3D_select_nodes % -function [Nx,Ref] = uwb30402_3D_select_nodes(N,k) % Extraction of the target node Nx = ceil(N*rand); Check=zeros(1,N); for i=1:k Ref(i) = Nx; % Check: is Ref(i) different from Nx? while((Ref(i)==Nx)||(Ref(i)==0)) Ref(i)= ceil(N*rand); % Check: is Ref(i) already selected? if(Check(Ref(i))) Ref(i) = 0; else Check(Ref(i))=1; end end end % -% uwb30403_3D_find_LSE_position % -function [PosNx, ErrNx] = uwb30403_3D_find_LSE_position(positions, ranges, Nx, Ref,sigma_2, G,num_fig) % Adding errors to the range estimation N = size(ranges,1); err_ranges = ranges + sqrt(sigma_2)*randn(N); % Defining the linear problem % Matrix A k = length(Ref); for i=1:(k-1) A(i,1) = positions(Ref(i),1) - positions(Ref(k),1); A(i,2) = positions(Ref(i),2) - positions(Ref(k),2); A(i,3) = positions(Ref(i),3) - positions(Ref(k),3); end A=-2*A; % Matrix b b=zeros(3,1); for i=1:(k-1) b(i) = err_ranges(Ref(i),Nx)^2 - err_ranges(Ref(k),Nx)^2 - positions(Ref(i),1)^2 + positions(Ref(k),1)^2 - positions(Ref(i),2)^2 + positions(Ref(k),2)^2 - positions(Ref(i),3)^2 + positions(Ref(k),3)^2; 90 end % Solving the problem PosNx=A\b; %Computing the error ErrNx = sqrt((PosNx(1)-positions(Nx,1))^2+(PosNx(2)- positions(Nx,2))^2+(PosNx(3)- positions(Nx,3))^2); save Sim_UWB_3D_H717b.mat; %Graphical output if G h = figure(num_fig); set(h,'name','H7_20 LSE_position_3D'); scatter3(positions(:,1),positions(:,2),positions(:,3)); xlabel('X [m]','fontname','.Vntime','fontsize',14); ylabel('Y [m]','fontname','.Vntime','fontsize',14); zlabel('Z [m]','fontname','.Vntime','fontsize',18); box on; hold on; scatter3(PosNx(1), PosNx(2), PosNx(2), 200, 'filled', 'k','p'); scatter3(positions(Nx,1),positions(Nx,2),positions(Nx,3),200, 'filled','^'); for i=1:k scatter3(positions(Ref(i),1),positions(Ref(i),2),positions(Ref(i),3), 'filled','r','s'); end title('Ước tính vị trí cho nút đích','fontname','.Vntime','color','b','fontsize',16); LT = legend('Nút khác','Vị trí nút đích','Vị trí hiệu dơng cđa nót ®Ých','Nót tham chiÕu'); % set(LT,'fontname','.Vntime','fontsize',12); hold off; end AX = gca; set(AX,'fontsize',12) ... phân tích nguyên lý độ rộng băng tần lượng phân đoạn định nghĩa vô tuyến băng siêu rộng UWB 1.2 Định nghĩa tín hiệu hệ thống vô tuyến băng siêu rộng UWB 1.2.1 Độ rộng băng tần phân đoạn Việc... pháp định vị điển hình; Giải pháp điển hình để cải thiện tính xác định vị 4 CHƯƠNG 1: TỔNG QUAN VỀ HỆ THỐNG VÔ TUYẾN BĂNG SIÊU RỘNG UWB 1.1 Giới thiệu Công nghệ vô tuyến băng siêu rộng UWB (Ultra-Wide... thông vô tuyến băng siêu rộng UWB 11 1.4 Đặc điểm tín hiệu hệ thống UWB 19 1.5 Các lĩnh vực ứng dụng vô tuyến UWB 25 1.6 Kết luận chương 27 CHƯƠNG II: ĐỊNH VỊ TRONG VÔ TUYẾN BĂNG