Nghiên cứu các phương pháp xác định chỉ số pwv bằng công nghệ đo gps

CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học : TS NGUYỄN NGỌC LÂU Cán chấm nhận xét : Cán chấm nhận xét : Luận văn thạc só bảo vệ HỘI ĐỒNG BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, Ngày tháng năm TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA KỸ THUẬT XÂY DỰNG CỘNG HOÀ XÃ HỘI CHỦ NGHIÃ VIỆT NAM Độc Lập - Tự Do - Hạnh Phúc -oOo Tp HCM, ngày tháng năm NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: TRANG CAO TUYÊN Phái: Nam Ngày, tháng, năm sinh: 14/10/1979 Nơi sinh: Vĩnh Long Chuyên ngành: Kỹ thuật trắc địa MSHV: 02207400 1- TÊN ĐỀ TÀI: NGHIÊN CỨU CÁC PHƯƠNG PHÁP XÁC ĐỊNH CHỈ SỐ PWV BẰNG CÔNG NGHỆ ĐO GPS 2- NHIỆM VỤ LUẬN VĂN:  Giới thiệu tổng quan sở định vị hệ thống GPS  Tìm hiểu ảnh hưởng tầng đối lưu vào tín hiệu GPS phương pháp khắc phục  Xác định độ trễ đối lưu từ trị đo pha GPS đánh giá độ tin cậy giá trị vừa tính, so sánh với liệu cung cấp tổ chức IGS  Nghiên cứu phương pháp chuyển đổi từ độ trễ đối lưu thiên đỉnh (TZD) sang số nước (PWV)  Đề xuất giải pháp việc xác định số nước phục vụ công tác dự báo thời tiết xử lý trị đo GPS độ xác cao 3- NGÀY GIAO NHIỆM VỤ : 01/02/2009 4- NGÀY HOÀN THÀNH NHIỆM VỤ : 30/11/2009 5- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN (Ghi đầy đủ học hàm, học vị ): TS NGUYỄN NGỌC LÂU Nội dung đề cương Luận văn thạc sĩ Hội Đồng Chuyên Ngành thông qua CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) TS NGUYỄN NGỌC LÂU CHỦ NHIỆM BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH (Họ tên chữ ký) TS NGUYỄN NGỌC LÂU KHOA QL CHUYÊN NGÀNH (Họ tên chữ ký) LỜI CÁM ƠN Xin chân thành cám ơn tồn thể q thầy Bộ mơn Địa Tin học trường Đại học Bách khoa TP.HCM đặc biệt Thầy TS.NGUYỄN NGỌC LÂU dành nhiều tâm huyết, tận tình hướng dẫn, bảo, truyền đạt kiến thức cho suốt thời gian thực Luận văn Tôi xin dành biết ơn sâu sắc gửi đến Ban lãnh đạo Chi nhánh Công ty cổ phần Tư vấn xây dựng điện 3- Xí nghiệp Khảo sát xây dựng điện tạo điều kiện động viên thời gian thực Luận văn Xin cám ơn tất bạn học, đồng nghiệp giúp giải vướng mắc trình thực Cuối cùng, xin cám ơn Gia đình hỗ trợ tinh thần thời gian để hồn thành tốt Luận văn Tp Hồ Chí Minh, tháng 11 năm 2009 Trang Cao Tun TĨM TẮT LUẬN VĂN THẠC SĨ Chỉ số nước khí (precipitable water vapor – PWV) cần thiết công tác dự báo thời tiết, kỹ thuật Insar Viễn Thám… Việc xác định số cách xác vấn đề quan tâm Một kỹ thuật phát triển gần dùng tín hiệu GPS để đo PWV thơng qua độ trễ qua tầng đối lưu [11, 13, 15, 18] Trong Luận văn này, học viên tiến hành nghiên cứu phương pháp xác định độ trễ đối lưu thiên đỉnh, khảo sát hàm ánh xạ chuyển đổi từ độ trễ đối lưu theo hướng thiên đỉnh sang độ trễ đối lưu theo hướng bất kỳ, phương pháp chuyển đổi từ độ trễ đối lưu sang số PWV, sử dụng thuật toán Kalman để xây dựng mơ hình liên tục từ số rời rạc Sau tiến hành phân tích, so sánh kết tính tốn với liệu cung cấp miễn phí trạm IGS gần nước Việt Nam (trạm Côn Minh – Trung Quốc) Áp dụng tính tốn cho số trạm đo GPS lãnh thổ Việt Nam phân tích mối quan hệ với thời tiết diễn khu vực, đề xuất giải pháp xử lý trị đo GPS độ xác cao MỤC LỤC MỞ ĐẦU 1 TÍNH CẦN THIẾT CỦA ĐỀ TÀI MỤC TIÊU NGHIÊN CỨU NỘI DUNG NGHIÊN CỨU PHƯƠNG PHÁP NGHIÊN CỨU 5 ĐỐI TƯỢNG VÀ PHẠM VI NGHIÊN CỨU Ý NGHĨA KHOA HỌC VÀ THỰC TIỄN 6.1 Ý NGHĨA KHOA HỌC 6.2 Ý NGHĨA THỰC TIỂN CHƯƠNG 1: CƠ SỞ ĐỊNH VỊ GPS .7 1.1 KHÁI QUÁT VỀ GPS 1.1.1 Bộ phận không gian 1.1.2 Bộ phận trạm điều khiển 1.1.3 Bộ phận người sử dụng 1.1.4 Ứng dụng 1.2 TÍN HIỆU VÀ TRỊ ĐO GPS 10 1.2.1 Tín hiệu 10 1.2.2 Trị đo 11 1.2.2.1 Mã giả khoảng cách 12 1.2.2.2 Pha sóng tải 12 1.3 NGUYÊN TẮC ĐỊNH VỊ BẰNG VỆ TINH 13 1.3.1 Cơng thức tốn học 13 1.3.2 Các sách Bộ Quốc Phịng Mỹ 14 1.4 CÁC NGUỒN SAI SỐ 15 1.4.1 Sai số quỹ đạo vệ tinh 16 1.4.2 Sai số đồng hồ vệ tinh 17 1.4.3 Sai số đồng hồ máy thu 18 1.4.4 Sai số tầng điện ly 18 1.4.5 Sai số tầng đối lưu 20 1.4.6 Sai số tượng đa đường 20 1.4.7 Sai số độ nhiễu tín hiệu 21 1.5 PHƯƠNG TRÌNH TRỊ ĐO VÀ CÁC PHƯƠNG PHÁP KẾT HỢP 21 1.5.1 Phương trình trị đo 21 1.5.2 Các phương pháp kết hợp 24 1.5.2.1 Hiệu hai máy thu vệ tinh 24 1.5.2.2 Hiệu hai vệ tinh máy thu 25 1.5.2.3 Hiệu đôi hai vệ tinh hai máy thu 26 1.5.2.4 Sử dụng L3 để loại bỏ ảnh hưởng tầng điện ly 27 1.5.2.5 Sử dụng trị đo hiệu L3 hai vệ tinh máy thu 28 1.5.2.6 Sử dụng sóng dải rộng dải hẹp để xác định trị nguyên đa trị kiểm soát thay đổi tầng điện ly 28 1.5.2.7 Thiết lập L4, L5 L6 để sửa chữa trượt chu kỳ phase 32 1.6 CÁC SẢN PHẨM CỦA TỔ CHỨC IGS 34 1.6.1 Tọa độ vệ tinh sai số đồng hồ vệ tinh 34 1.6.2 Tổng độ trễ tầng đối lưu theo hướng thiên đỉnh 37 1.7 TÓM TẮT VÀ KẾT LUẬN 39 CHƯƠNG 2: TÌM HIỂU VỀ TẦNG ĐỐI LƯU VÀ CÁC MƠ HÌNH HIỆU CHỈNH 40 2.1 SỰ HÌNH THÀNH TẦNG ĐỐI LƯU VÀ ẢNH HƯỞNG CỦA TÍN HIỆU GPS 40 2.2 MƠ HÌNH HÓA ẢNH HƯỞNG CỦA TẦNG ĐỐI LƯU VÀO TRỊ ĐO GPS 43 2.3 CÁC MƠ HÌNH HIỆU CHỈNH TẦNG ĐỐI LƯU 44 2.3.1 Mơ hình độ trễ đối lưu cho thành phần khơ 44 2.3.1.1 Mơ hình Saastamoinen 44 2.3.1.2 Mơ hình Mops 46 2.3.1.3 Mơ hình Hopfield 48 2.3.1.4 Mô hình Hopfield cải tiến 50 2.3.2 Mơ hình độ trễ đối lưu cho thành phần ướt 53 2.3.2.1 Mơ hình Ifadis: 53 2.3.2.2 Mơ hình Mendes: 53 2.3.2.3 Mơ hình Mops: 53 2.3.2.4 Mơ hình Hopfield: 54 2.3.2.5 Mơ hình Hopfield cải tiến: 55 2.3.3 Mơ hình khí chuẩn 55 2.4 HÀM ÁNH XẠ 57 2.4.1 Hàm ánh xạ CfA-2.2 58 2.4.2 Hàm ánh xạ Ifadis 59 2.4.3 Hàm ánh xạ MTT 59 2.4.4 Hàm ánh xạ NMF 60 2.4.5 Hàm ánh xạ GMF 61 2.5 TÓM TẮT VÀ KẾT LUẬN 64 CHƯƠNG 3: THUẬT TOÁN XÁC ĐỊNH ĐỘ TRỄ ĐỐI LƯU THIÊN ĐỈNH VÀ CHỈ SỐ HƠI NƯỚC TẠI CÁC TRẠM ĐO GPS 65 3.1 CƠ SỞ LÝ THUYẾT 65 3.1.1 Các phương pháp xử lý số liệu đo GPS 65 3.1.1.1 Phương pháp xử lý tương đối 65 3.1.1.2 Phương pháp xử lý tuyệt đối 66 3.1.2 Bộ lọc Kalman 66 3.1.2.1 Giới thiệu lọc Kalman 66 3.1.2.2 Công thức Kalman: 68 3.1.2.3 Thuật toán Kalman: 70 3.2 XÁC ĐỊNH ĐỘ TRỄ ĐỐI LƯU THIÊN ĐỈNH BẰNG PHƯƠNG PHÁP XỬ LÝ TUYỆT ĐỐI 72 3.2.1 Xác định ma trận số hạng tự L 72 3.2.2 Các mô hình chức 73 3.2.3 Mơ hình sai số 74 3.3 THUẬT TOÁN XÁC ĐỊNH ĐỘ TRỄ ĐỐI LƯU THIÊN ĐỈNH BẰNG PHƯƠNG PHÁP XỬ LÝ TUYỆT ĐỐI 76 3.4 XÁC ĐỊNH ĐỘ TRỄ ĐỐI LƯU THIÊN ĐỈNH TẠI TRẠM ĐO IGS 78 3.4.1 Số liệu sử dụng 78 3.4.2 Kết 79 3.5 PHƯƠNG PHÁP CHUYỂN ĐỔI ĐỘ TRỄ ĐỐI LƯU THIÊN ĐỈNH SANG CHỈ SỐ HƠI NƯỚC 81 3.5.1 Cơ sở lý thuyết 81 3.5.2 Kết 82 3.5.3 Kỹ thuật lọc nhiễu 82 3.6 TÓM TẮT VÀ KẾT LUẬN 85 CHƯƠNG 4: KẾT QUẢ XỬ LÝ TẠI MỘT SỐ TRẠM ĐO Ở VIỆT NAM 87 4.1 THIẾT KẾ LƯỚI 87 4.2 KẾT QUẢ 89 4.2.1 Giá trị TZD 89 4.2.2 Giá trị PWV 92 4.2.3 Nhận xét 95 4.3 SỰ LIÊN HỆ GIỮA GIÁ TRỊ PWV VÀ ĐIỀU KIỆN THỜI TIẾT TRÊN KHU VỰC 96 4.4 TÓM TẮT VÀ KÉT LUẬN 99 CHƯƠNG 5: KẾT LUẬN 100 5.1 CÁC VẤN ĐỀ ĐẠT ĐƯỢC 100 5.2 NHỮNG KHÓ KHĂN VÀ HẠN CHẾ 101 5.3 HƯỚNG PHÁT TRIỂN ĐỀ TÀI 102 TÀI LIỆU THAM KHẢO 103 Trang MỞ ĐẦU TÍNH CẦN THIẾT CỦA ĐỀ TÀI Chỉ số nước (precipitable water vapor – PWV) số lượng toàn nước điểm bề mặt trái đất, biểu diễn độ cao cột chất lõng tương đương Chỉ số nước độ trễ điện ly cần thiết công tác dự báo thời tiết, xử lý trị đo GPS xác cao, kỹ thuật INSAR (kỹ thuật dựa độ lệch pha tín hiệu radar để tính tốn xử lý ảnh)  Trong dự báo thời tiết, người ta muốn theo dõi tiến triển bão, mưa nên người ta thường dùng bóng thám khơng, sóng water vapor radiometer (WVR) [10] Nhưng hai phương pháp có trở ngại riêng: bóng thám khơng sử dụng ít, đắt tiền tuổi thọ khơng cao; cịn sóng WPR khơng dùng lúc trời mưa Vì dùng GPS công cụ để lấy số PWV phục vụ dự báo thời tiết + Những thập niên gần đây, thiên tai tàn phá dội trái đất chúng ta: sóng thần Thái Lan, lũ lụt Ấn Độ gần động đất Trung Quốc gây hậu nghiêm trọng: làm chết nhiều người, thiệt hại nhiều tiền của, tàn phá thiên nhiên, tài nguyên, môi trường Chúng ta làm giảm bớt hậu thiên tai dự báo trước khả xảy thiên tai để có cách đối phó, xử lý hợp lý + Có nhiều cách để dự báo thiên tai: theo kinh nghiệm thực tế: nhìn trời âm u đốn có mưa, sau mưa trời có nắng ; theo kinh nghiệm cha ơng: “chuồn chuồn bay thấp mưa, bay cao nắng, bay vừa râm”, bìm bịp kêu nước lớn ; theo khoa học: đo địa chấn đất, nghiên cứu quy luật di chuyển số loài động vật để dự báo động đất, lũ lụt ; ngành thiên văn học ngắm di chuyển trời để dự báo thời tiết; ngành khí tượng thủy văn sử Trang dụng bóng thám khơng để đo nhiệt độ, áp suất, độ ẩm tầng đối lưu để phục vụ dự báo thời tiết (trong thơng số độ ẩm quan trọng nhất); đặt trạm quan trắc số khu vực đặt biệt để theo dõi thơng số khí tượng, thủy văn; tạo vườn để nghiên cứu thay đổi theo thời tiết Ngồi ra, người ta cịn dùng ảnh vệ tinh để dự báo thời tiết + Việc dự báo thời tiết chiếm vai trò quan trọng hầu hết đời sống sinh hoạt hàng ngày người: nông dân cần biết diễn biến thời tiết ngày hay ngày hôm sau để bố trí cơng việc đồng thuận lợi; ngư dân cần biết diễn biến thời tiết để định thời điểm khơi đánh cá thời điểm để trở đất liền; nhà thầu xây dựng cần biết diễn biến thời tiết để có kế hoạch thực cơng việc tốt Vì hầu hết người xem dự báo thời tiết để xếp công việc thuận lợi + Một nhu cầu cần thiết cần có cơng cụ khác để phục vụ dự báo thời tiết, người ta nghiên cứu phương thức khác để thay bóng thám khơng cho có lợi kinh tế đảm bảo độ xác dự báo tương đương  Trong việc xử lý trị đo công nghệ GPS địi hỏi độ xác cao, xử lý đường đáy dài, độ trễ đối lưu chướng ngại lớn việc xác định độ cao chiếm vai trò quan trọng việc nâng cao độ xác kết xử lý Để cải thiện điều này, ta cần hiểu biết xác tầng đối lưu, mơ hình hiệu chỉnh độ trễ đối lưu có sẵn ưu nhược điểm + Hiện nay, người ta sử dụng tham số hiệu chỉnh tầng đối lưu chủ yếu dựa vào mơ hình hiệu chỉnh có sẵn, kết kết từ trung tâm IGS Những mơ hình có sẵn HOPFIELD, MENDES, MOPS cung cấp tham số với độ xác tương đối Cịn việc lấy kết từ trung tâm IGS cho kết xác vị trí đặt trạm đo, ten_vt(5)=t5(2) ; ten_vt(6)=t6(2); %input('') tong_l3(1)=t1(2); tong_l3(2)=t2(3); tong_l3(3)=t3(3); tong_l3(4)=t4(3); tong_l3(5)=t5(3); tong_l3(6)=t6(3); sine(1)=t1(3); sine(2)=t2(4); sine(3)=t3(4); sine(4)=t4(4); sine(5)=t5(4); sine(6)=t6(4); delta_tong_l3(1)=t1(2)-t1(2); delta_tong_l3(2)=t1(2)-t2(3); delta_tong_l3(3)=t1(2)-t3(3); delta_tong_l3(4)=t1(2)-t4(3); delta_tong_l3(5)=t1(2)-t5(3); delta_tong_l3(6)=t1(2)-t6(3); delta_map1(1)=t1(4)-t1(4); delta_map1(2)=t1(4)-t2(5) ; delta_map1(3)=t1(4)-t3(5) ; delta_map1(4)=t1(4)-t4(5) ; delta_map1(5)=t1(4)-t5(5) ; delta_map1(6)=t1(4)-t6(5) ; delta_map2(1)=t1(5)-t1(5); delta_map2(2)=t1(5)-t2(6); delta_map2(3)=t1(5)-t3(6); delta_map2(4)=t1(5)-t4(6); delta_map2(5)=t1(5)-t5(6); delta_map2(6)=t1(5)-t6(6); % -% elseif vt==7 t1=fscanf(fid,'%f',5); t2=fscanf(fid,'%f',6); t3=fscanf(fid,'%f',6); t4=fscanf(fid,'%f',6); t5=fscanf(fid,'%f',6); t6=fscanf(fid,'%f',6); t7=fscanf(fid,'%f',6); so_vt=vt; ten_vt(1)=t1(1); ten_vt(2)=t2(2) ; ten_vt(3)=t3(2) ; ten_vt(4)=t4(2) ; ten_vt(5)=t5(2) ; ten_vt(6)=t6(2) ; ten_vt(7)=t7(2) ; tong_l3(1)=t1(2); tong_l3(2)=t2(3); tong_l3(3)=t3(3); tong_l3(4)=t4(3); tong_l3(5)=t5(3); tong_l3(6)=t6(3); tong_l3(7)=t7(3); sine(1)=t1(3); sine(2)=t2(4); sine(3)=t3(4); sine(4)=t4(4); sine(5)=t5(4); sine(6)=t6(4); sine(7)=t7(4); delta_tong_l3(1)=t1(2)-t1(2); delta_tong_l3(2)=t1(2)-t2(3); delta_tong_l3(3)=t1(2)-t3(3); delta_tong_l3(4)=t1(2)-t4(3); delta_tong_l3(5)=t1(2)-t5(3); delta_tong_l3(6)=t1(2)-t6(3); delta_tong_l3(7)=t1(2)-t7(3); delta_map1(1)=t1(4)-t1(4); delta_map1(2)=t1(4)-t2(5) ; delta_map1(3)=t1(4)-t3(5) ; delta_map1(4)=t1(4)-t4(5) ; delta_map1(5)=t1(4)-t5(5) ; delta_map1(6)=t1(4)-t6(5) ; delta_map1(7)=t1(4)-t7(5) ; delta_map2(1)=t1(5)-t1(5); delta_map2(2)=t1(5)-t2(6); delta_map2(3)=t1(5)-t3(6); delta_map2(4)=t1(5)-t4(6); delta_map2(5)=t1(5)-t5(6); delta_map2(6)=t1(5)-t6(6); delta_map2(7)=t1(5)-t7(6); % -% elseif vt==8 t1=fscanf(fid,'%f',5); t2=fscanf(fid,'%f',6); t3=fscanf(fid,'%f',6); t4=fscanf(fid,'%f',6); t5=fscanf(fid,'%f',6); t6=fscanf(fid,'%f',6); t7=fscanf(fid,'%f',6); t8=fscanf(fid,'%f',6); so_vt=vt; ten_vt(1)=t1(1); ten_vt(2)=t2(2) ; ten_vt(3)=t3(2) ; ten_vt(4)=t4(2) ten_vt(5)=t5(2) ten_vt(6)=t6(2) ten_vt(7)=t7(2) ten_vt(8)=t8(2) ; ; ; ; ; tong_l3(1)=t1(2); tong_l3(2)=t2(3); tong_l3(3)=t3(3); tong_l3(4)=t4(3); tong_l3(5)=t5(3); tong_l3(6)=t6(3); tong_l3(7)=t7(3); tong_l3(8)=t8(3); sine(1)=t1(3); sine(2)=t2(4); sine(3)=t3(4); sine(4)=t4(4); sine(5)=t5(4); sine(6)=t6(4); sine(7)=t7(4); sine(8)=t8(4); delta_tong_l3(1)=t1(2)-t1(2); delta_tong_l3(2)=t1(2)-t2(3); delta_tong_l3(3)=t1(2)-t3(3); delta_tong_l3(4)=t1(2)-t4(3); delta_tong_l3(5)=t1(2)-t5(3); delta_tong_l3(6)=t1(2)-t6(3); delta_tong_l3(7)=t1(2)-t7(3); delta_tong_l3(8)=t1(2)-t8(3); delta_map1(1)=t1(4)-t1(4); delta_map1(2)=t1(4)-t2(5) ; delta_map1(3)=t1(4)-t3(5) ; delta_map1(4)=t1(4)-t4(5) ; delta_map1(5)=t1(4)-t5(5) ; delta_map1(6)=t1(4)-t6(5) ; delta_map1(7)=t1(4)-t7(5) ; delta_map1(8)=t1(4)-t8(5) ; delta_map2(1)=t1(5)-t1(5); delta_map2(2)=t1(5)-t2(6); delta_map2(3)=t1(5)-t3(6); delta_map2(4)=t1(5)-t4(6); delta_map2(5)=t1(5)-t5(6); delta_map2(6)=t1(5)-t6(6); delta_map2(7)=t1(5)-t7(6); delta_map2(8)=t1(5)-t8(6); % -% elseif vt==9 t1=fscanf(fid,'%f',5); t2=fscanf(fid,'%f',6); t3=fscanf(fid,'%f',6); t4=fscanf(fid,'%f',6); t5=fscanf(fid,'%f',6); t6=fscanf(fid,'%f',6); t7=fscanf(fid,'%f',6); t8=fscanf(fid,'%f',6); t9=fscanf(fid,'%f',6); so_vt=vt; ten_vt(1)=t1(1); ten_vt(2)=t2(2) ; ten_vt(3)=t3(2) ; ten_vt(4)=t4(2) ; ten_vt(5)=t5(2) ; ten_vt(6)=t6(2) ; ten_vt(7)=t7(2) ; ten_vt(8)=t8(2) ; ten_vt(9)=t9(2) ; tong_l3(1)=t1(2); tong_l3(2)=t2(3); tong_l3(3)=t3(3); tong_l3(4)=t4(3); tong_l3(5)=t5(3); tong_l3(6)=t6(3); tong_l3(7)=t7(3); tong_l3(8)=t8(3); tong_l3(9)=t9(3); sine(1)=t1(3); sine(2)=t2(4); sine(3)=t3(4); sine(4)=t4(4); sine(5)=t5(4); sine(6)=t6(4); sine(7)=t7(4); sine(8)=t8(4); sine(9)=t9(4); delta_tong_l3(1)=t1(2)-t1(2); delta_tong_l3(2)=t1(2)-t2(3); delta_tong_l3(3)=t1(2)-t3(3); delta_tong_l3(4)=t1(2)-t4(3); delta_tong_l3(5)=t1(2)-t5(3); delta_tong_l3(6)=t1(2)-t6(3); delta_tong_l3(7)=t1(2)-t7(3); delta_tong_l3(8)=t1(2)-t8(3); delta_tong_l3(9)=t1(2)-t9(3); delta_map1(1)=t1(4)-t1(4); delta_map1(2)=t1(4)-t2(5) ; delta_map1(3)=t1(4)-t3(5) ; delta_map1(4)=t1(4)-t4(5) ; delta_map1(5)=t1(4)-t5(5) ; delta_map1(6)=t1(4)-t6(5) ; delta_map1(7)=t1(4)-t7(5) ; delta_map1(8)=t1(4)-t8(5) ; delta_map1(9)=t1(4)-t9(5) ; delta_map2(1)=t1(5)-t1(5); delta_map2(2)=t1(5)-t2(6); delta_map2(3)=t1(5)-t3(6); delta_map2(4)=t1(5)-t4(6); delta_map2(5)=t1(5)-t5(6); delta_map2(6)=t1(5)-t6(6); delta_map2(7)=t1(5)-t7(6); delta_map2(8)=t1(5)-t8(6); delta_map2(9)=t1(5)-t9(6); % -% elseif vt==10 t1=fscanf(fid,'%f',5); t2=fscanf(fid,'%f',6); t3=fscanf(fid,'%f',6); t4=fscanf(fid,'%f',6); t5=fscanf(fid,'%f',6); t6=fscanf(fid,'%f',6); t7=fscanf(fid,'%f',6); t8=fscanf(fid,'%f',6); t9=fscanf(fid,'%f',6); t10=fscanf(fid,'%f',6); so_vt=vt; ten_vt(1)=t1(1); ten_vt(2)=t2(2) ; ten_vt(3)=t3(2) ; ten_vt(4)=t4(2) ; ten_vt(5)=t5(2) ; ten_vt(6)=t6(2) ; ten_vt(7)=t7(2) ; ten_vt(8)=t8(2) ; ten_vt(9)=t9(2) ; ten_vt(10)=t10(2) ; tong_l3(1)=t1(2); tong_l3(2)=t2(3); tong_l3(3)=t3(3); tong_l3(4)=t4(3); tong_l3(5)=t5(3); tong_l3(6)=t6(3); tong_l3(7)=t7(3); tong_l3(8)=t8(3); tong_l3(9)=t9(3); tong_l3(10)=t10(3); sine(1)=t1(3); sine(2)=t2(4); sine(3)=t3(4); sine(4)=t4(4); sine(5)=t5(4); sine(6)=t6(4); sine(7)=t7(4); sine(8)=t8(4); sine(9)=t9(4); sine(10)=t10(4); delta_tong_l3(1)=t1(2)-t1(2); delta_tong_l3(2)=t1(2)-t2(3); delta_tong_l3(3)=t1(2)-t3(3); delta_tong_l3(4)=t1(2)-t4(3); delta_tong_l3(5)=t1(2)-t5(3); delta_tong_l3(6)=t1(2)-t6(3); delta_tong_l3(7)=t1(2)-t7(3); delta_tong_l3(8)=t1(2)-t8(3); delta_tong_l3(9)=t1(2)-t9(3); delta_tong_l3(10)=t1(2)-t10(3); delta_map1(1)=t1(4)-t1(4); delta_map1(2)=t1(4)-t2(5) ; delta_map1(3)=t1(4)-t3(5) ; delta_map1(4)=t1(4)-t4(5) ; delta_map1(5)=t1(4)-t5(5) ; delta_map1(6)=t1(4)-t6(5) ; delta_map1(7)=t1(4)-t7(5) ; delta_map1(8)=t1(4)-t8(5) ; delta_map1(9)=t1(4)-t9(5) ; delta_map1(10)=t1(4)-t10(5) ; delta_map2(1)=t1(5)-t1(5); delta_map2(2)=t1(5)-t2(6); delta_map2(3)=t1(5)-t3(6); delta_map2(4)=t1(5)-t4(6); delta_map2(5)=t1(5)-t5(6); delta_map2(6)=t1(5)-t6(6); delta_map2(7)=t1(5)-t7(6); delta_map2(8)=t1(5)-t8(6); delta_map2(9)=t1(5)-t9(6); delta_map2(10)=t1(5)-t10(6); end vt_base=ten_vt(1); % tinh ma tran A % -% for j=1:vt A(j,1)=delta_map2(j); %if A(j,1)~=0 A(j,ten_vt(j)+1)=1; %end end % -% if ten_vt(1)==1 for i=1:so_vt if first(1,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(1,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==2 for i=1:so_vt if first(2,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(2,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==3 for i=1:so_vt if first(3,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(3,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==4 for i=1:so_vt if first(4,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(4,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==5 for i=1:so_vt if first(5,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(5,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==6 for i=1:so_vt if first(6,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(6,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==7 for i=1:so_vt if first(7,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(7,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==8 for i=1:so_vt if first(8,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(8,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==9 for i=1:so_vt if first(9,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(9,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==10 for i=1:so_vt if first(10,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(10,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==11 for i=1:so_vt if first(11,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(11,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==12 for i=1:so_vt if first(12,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(12,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==13 for i=1:so_vt if first(13,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(13,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==14 for i=1:so_vt if first(14,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(14,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==15 for i=1:so_vt if first(15,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(15,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==16 for i=1:so_vt if first(16,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(16,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==17 for i=1:so_vt if first(17,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(17,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==18 for i=1:so_vt if first(18,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(18,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==19 for i=1:so_vt if first(19,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(19,ten_vt(i))=0; unknown=unknown+1; end end % -% % -% elseif ten_vt(1)==20 for i=1:so_vt if first(20,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(20,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==21 for i=1:so_vt if first(21,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(21,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==22 for i=1:so_vt if first(22,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(22,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==23 for i=1:so_vt if first(23,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(23,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==24 for i=1:so_vt if first(24,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(24,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==25 for i=1:so_vt if first(25,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(25,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==26 for i=1:so_vt if first(26,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(26,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==27 for i=1:so_vt if first(27,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(27,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==28 for i=1:so_vt if first(28,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(28,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==29 for i=1:so_vt if first(29,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(29,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==30 for i=1:so_vt if first(30,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(30,ten_vt(i))=0; unknown=unknown+1; end end % -% end % -% % -% %Prediction step Qx(1,1)=Qx(1,1)+0.01^2*30/3600; %Qx(2,2)=Qx(2,2)+sigma_rclck^2; unknown=unknown+1;n_obser=n_obser+1; %Compute residual so_vt; ten_vt; delta_tong_l3; A; X ; vt_base; V=delta_tong_l3-A*X; X(1); %n_epoch*30 %input('') %Update step %tinh ma tran QL0 % -% Ql=(1/(sine(1))^2)*eye(vt,vt) ; for j=1:vt-1 if sine(j)~=0 Ql(j,j)=Ql(j,j)+1/(sine(j)^2); end end Ql=Ql*xigma0; % Ql=xigma0*diag((0.01+1./A(:,1)).^2); %input('') Z=Qx*A'; Qv=Ql+A*Z; % -% %Check raw errors i=1; % % -% while i5*mu*sqrt(Qv(i,i)) j=find(A(i,2:33)==1);k=j+1; %Restart ambiguity Qx(:,k)=0;Qx(k,:)=0; Qx(k,k)=1.0;first(k-1)=1; %Delete bad observation Z(:,i)=[]; V(i,:)=[]; Qv(i,:)=[];Qv(:,i)=[]; so_vt=so_vt-1;n_bad=n_bad+1; else i=i+1; end end % -% if so_vt>1 Qv=inv(Qv); X=X+Z*Qv*V; Qx=Qx-Z*Qv*Z'; VtPV=VtPV+V'*Qv*V; thamso=thamso+so_vt; muy=sqrt(VtPV/(thamso-so_vt+1)) X(1) %Plot % subplot(2,1,1); %plot(23+n_epoch*30/3600/24,X(1),'r.');hold on; grid on; % subplot(2,1,2); % plot(23+n_epoch*30/3600/24,X(1)*0.159*0.9,'r.');hold on; % xlabel('Day of September 2007','fontsize',13); % ylabel('PWV','fontsize',13); % ylabel('TZD','fontsize',13); n_obser=n_obser+so_vt; if n_obser>unknown mu=sqrt(VtPV/(n_obser-unknown));end end n_epoch=n_epoch+1; vt_base2=ten_vt(1); %input('') fprintf(fid4,'%.16f %.16f \n',(n_epoch*30),X(1)); fprintf(fid2,'%.16f %.16f \n',n_epoch*30/3600,X(1)*0.159); fprintf(fid3,'%.16f \n',muy); end %fprintf(fid1,'%f %f %.16f %.16f %.16f %.16f %.16f\n',so_vt,ten_vt,tong_l3,sine,delta_tong_l3,delta_map1,delta_map2); fclose(fid);fclose(fid1); fprintf('%d %d %f\n',n_obser,n_bad,mu); % Ve TZD o tram Kunminh fn= 'kunm2660.07zpd'; fid=fopen(fn,'r'); fid1=fopen('TZD_IGS2660.txt','w'); for j=1:61 fgets(fid);end; i=1; while ~feof(fid) s=fscanf(fid,'%*s%s',1); if s(1:2)~='07' break;end x1(i)=str2num(s(8:12))/3600; y1(i)=fscanf(fid,'%lf',1)/1000; fprintf(fid1,'%.16f %.16f \n',x1(i),y1(i)); fgets(fid); i=i+1; end %subplot(2,1,1); %plot(23+x1/24,y1,'b.'); %subplot(2,1,2); plot(23+x1/24,(y1-ZHD)*0.159,'b.'); hold on; fclose(fid); fclose(fid1); % n_epoch; xlabel('Day of September 2007','fontsize',13); ylabel('PWV','fontsize',13); %ylabel('TZD','fontsize',13); grid on; % -% Loc Nhieu -fn= 'TZD_kunm2660.txt'; fid=fopen(fn,'r'); fif1=fopen('locnhieu.txt','w'); i=1; while ~feof(fid) s=fscanf(fid,'%f',2); if isempty(s) break;end x(i)=s(1) ; y(i)=s(2); i=i+1; end x1=0:300:85550; y1=spline(x,y,x1); fprintf(fid1,'%.16f %.16f \n',x1,y1); %plot(23+x/3600/24,y,'o',23+x1/3600/24,y1); %plot(23+x1/3600/24,y1,'c.'); plot(23+x1/3600/24,(y1-ZHD)*0.159,'c.'); hold on; grid on; fclose(fid); ... TÀI: NGHIÊN CỨU CÁC PHƯƠNG PHÁP XÁC ĐỊNH CHỈ SỐ PWV BẰNG CÔNG NGHỆ ĐO GPS 2- NHIỆM VỤ LUẬN VĂN:  Giới thiệu tổng quan sở định vị hệ thống GPS  Tìm hiểu ảnh hưởng tầng đối lưu vào tín hiệu GPS phương. .. TIÊU NGHIÊN CỨU  Nghiên cứu phương pháp xác định độ trễ đối lưu từ trị đo GPS, từ chọn lựa phương pháp thích hợp điều kiện Việt Nam  Nghiên cứu phương pháp chuyển đổi từ độ trễ đối lưu sang số. .. TOÁN XÁC ĐỊNH ĐỘ TRỄ ĐỐI LƯU THIÊN ĐỈNH VÀ CHỈ SỐ HƠI NƯỚC TẠI CÁC TRẠM ĐO GPS 65 3.1 CƠ SỞ LÝ THUYẾT 65 3.1.1 Các phương pháp xử lý số liệu đo GPS 65 3.1.1.1 Phương pháp