Khảo sát hệ số PWV từ bóng thám không và hệ thống định vị vệ tinh GPS

HÀNH PH H HU KH KHÔNG VÀ H Chuyên ngành: K TP H , 11/2013 CƠNG TRÌNH ÀN THÀNH T – Cán b TP HCM PGS TS Cán b Cán b Lu Tr Chí Minh, ngày…… tháng…… Thành ph n H i (Ghi rõ h ng ên, h PGS TS Nguy ng i h c Bách Khoa, … nh giá lu n v n th c s g m: àm, h lu g …… … …… ……………… ình …………… ên ………… Xác nh n c a Ch t ch H i ng ngành sau lu n v n ã c s a ch a (n u có) CH HQG TP h giá LV Khoa qu lý chuyên C ÒA XÃ H -T -H NHI H ên h ên: Hu …………………… MSHV: 11224490 05/03/1984…………………………… Chuyên ngành: K I –B Bình Thu ……………………… Mã s ÀI: “Kh 60520503 àh v II NHI ÀN Thu th àd Tính h àt So sánh k Rút k ên ki III NGÀY GIAO NHI : 24/06/2013………………………………………… IV NGÀY HOÀN THÀNH NHI V CÁN B 22/11/2013…………………………… PGS TS NGUY ………………… Tp HCM, ngày CÁN B (H ên ch tháng ÀO T (H PGS TS NGUY ên ch PGS TS NGUY …………… (H ên ch L h ày t TS NGUY truy òng bi ã dành nhi ình h th àm lu Em xin dành s ãnh ịa t ãt Trong q trình h ên c àm lu th i em xin chân ên c h ài Ngun Mơi T trình h M l ên c ã giúp th ãr ,t àn thành àn thi ình n ình, nhiên không th nghiên c , Cô Em xin trân tr Tp H Hu TĨM T – PWV) khí quy c xác hi àv M dùng tín hi h Trong Lu ti ành thu th ày, h àd T d ên i thám khơng tính h tính h ành ki ên d cung c tr phân tích m ngày So sánh k d l -0.9mm 1.2mm, ngồi giá tr h nh i th bóng thám không c ày, t l ng th d GPS ù CBHD: PGS.TS M M 1.1 GI 1.2 TÍNH C ÀI 1.2.1 Máy phóng x – Water vapor Radiometer (WVR) 1.2.2 H 1.2.3 H 1.2.4 H b ên m based remote sensers) 1.2.5 H 1.2.6 H 1.3 M À PH 1.3.1 M ên c 1.3.2 ên c 1.4 Ý NGH ÊN C ph À TH ên c 10 1.4.1 Ý ngh 10 1.4.2 Ý ngh 10 T 11 2.1 S 12 2.2 T 2.2.1 Tình hình nghiên c ÌNH HÌNH NGHIÊN C ên th 2.2.2 Tình hình nghiên c 13 13 17 i HV: HU CBHD: PGS.TS KI 20 3.1 GI 21 3.1.1 Khái ni 21 3.1.2 Bóng th 22 3.1.3 Thi 24 3.2 THU BĨNG THÁM KHƠNG VƠ TUY 28 3.2.1 Thu [10] 28 3.2.2 Thu òng D Qu 31 3.2.3 vô ên h 32 3.3 THU LI TÍCH T PWV T 34 3.3.1 3.3.1.1 Thi 37 39 3.3.1.2 L 3.3.1.3 S àm 40 40 3.3.1.4 Hi 41 3.3.1.5 B 42 a Gi b Cơng th c Tóm t 42 43 Kalman 49 3.3.2 Thu 50 3.3.3 Mơ hình 54 ii HV: HU CBHD: PGS.TS 3.3.4 56 , K 58 4.1 THU TH TUY VÔ T 59 4.1.1 Ph 62 a C 64 b C io 65 c C 65 d C 66 e C 67 4.1.2 Ph 68 a 68 b 69 c 69 4.2 TH BĨNG THÁM KHƠNG VƠ TUY 4.3 THU TH 70 71 4.4 MÔ T ÀD 73 4.5 X TR 4.5.1 Gi ÒA - TP.HCM 76 ình tính h 76 4.5.2 K 78 4.6 X 4.6.1 Gi 86 ình tính h iii 86 HV: HU CBHD: PGS.TS 4.6.2 Ki d 88 a V 88 b Gi 90 c K 93 4.6.3 K 96 4.7 103 À KI .109 5.1 CÁC V 110 5.2 .111 5.3 KI .112 113 iv HV: HU CBHD: PGS.TS C M Trang HV: HU 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); elseif vt==11 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); t11=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); ten_vt(11)=t11(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); tong_l3(11)=t11(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); sine(11)=t11(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_tong_l3(11)=t1(2)-t11(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_map1(11)=t1(4)-t11(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); delta_map2(11)=t1(5)-t11(6); elseif vt==12 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); t11=fscanf(fid,'%f',6); t12=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); ten_vt(11)=t11(2); ten_vt(12)=t12(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); tong_l3(11)=t11(3); tong_l3(12)=t12(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); sine(11)=t11(4); sine(12)=t12(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_tong_l3(11)=t1(2)-t11(3); delta_tong_l3(12)=t1(2)-t12(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_map1(11)=t1(4)-t11(5); delta_map1(12)=t1(4)-t12(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); delta_map2(11)=t1(5)-t11(6); delta_map2(12)=t1(5)-t12(6); end vt_base=ten_vt(1); %input('') %CHUONG TRINH TINH MA TRAN A: % 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 % -% A; %input('') %tinh ma tran khoi dau X0 % -% % -% 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; %input('') 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 % -% elseif ten_vt(1)==31 for i=1:so_vt if first(31,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(31,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==32 for i=1:so_vt if first(32,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(32,ten_vt(i))=0; unknown=unknown+1; end end % -% elseif ten_vt(1)==33 for i=1:so_vt if first(33,ten_vt(i)) X(ten_vt(i)+1)=delta_tong_l3(i)-A(i,1)*X(1); first(33,ten_vt(i))=0; unknown=unknown+1; end end % -% end % -% %input('') % % %Prediction step Qx(1,1)=Qx(1,1)+0.1^2*15/3600; unknown=unknown+1;n_obser=n_obser+1; %Compute residual % delta_tong_l3-sine-delta_map1 V=(delta_tong_l3)-A*X; %input('') %Update step %tinh ma tran QL0 % -% Ql=(1/(sine(1))^2)*eye(vt,vt) ; for e=1:vt-1 if sine(e)~=0 Ql(e,e)=Ql(e,e)+1/(sine(e)^2); end end Ql=Ql*xigma0; %(sovetinh;sovetinh) %Ql=xigma0^2*diag((0.01+1./A(:,1)).^2); %input('') Z=Qx*A'; %(33;4) Qv=Ql+A*Z; % -% %Check raw errors ii=1; % -% while ii3.28*mu*sqrt(Qv(ii,ii)) h=find(A(ii,2:33)==1); k=h+1; %Restart ambiguity Qx(:,k)=0;Qx(k,:)=0; Qx(k,k)=1.0;first(k-1)=1; Qx; %Delete bad observation Z(:,ii)=[]; V(ii,:)=[]; Qv(ii,:)=[];Qv(:,ii)=[]; so_vt=so_vt-1;n_bad=n_bad+1; else ii=ii+1; end end % -% if so_vt>1 Qv=inv(Qv); X=X+Z*Qv*V; Qx=Qx-Z*Qv*Z'; % input('') VtPV=VtPV+V'*Qv*V; % unknown=so_vt+1; thamso=thamso+so_vt; muy=sqrt(VtPV/(thamso-so_vt+1)); 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,'%.2d %.16f \n',((n_epoch*15)+4*3600+57*60+0) ,X(1)); %fprintf(fid1,'%d %d %.16f %.16f %.16f %.16f %.16f\n','numesat','namesat','L3','sine','delta_l3','delta_map(1)',' delta_map(2)'); end fclose(fid);%fclose(fid1); fclose('all'); % -% function tinhhesoPWV_GPS fid=importdata('TZD_16012013.txt'); fid1=importdata('file chua ZHD_16012013.txt'); time1=12.22690300;%thoi gian trung binh cua bong tham khong theo Liu PWV_btk1=32.92702943;%do lon cua he so PWV cua bong tham khong theo Liu time2=12.22690300;%thoi gian trung binh bong tham khong theo Noaa PWV_btk2=32.63645962;%do lon cua he so PWV cua bong tham khong theo Noaa % -% %plot((fid(:,1)/3600),fid(:,2),'r.');hold on PWV_GPS=(fid(:,2)-fid1(:,1))*0159 plot((fid(:,1)/3600),PWV_GPS,'b.');hold on plot(time1,PWV_btk1,'r^'); hold on plot(time2,PWV_btk2,'ms'); hold on grid on; xlabel('Thoi gian thu ve tinh(h)'); ylabel('He so PWV cua GPS 16-1(mm)');hold on % % legend('PWV-GPS','PWV-bong.Liu','PWV-bong.Noaa'); %PWVmean=mean(PWV_GPS); fclose('all'); ... WVR - Radiosonde 18 (3 PWV) 13 (2 PWV) VLBI - GPS (1.5 PWV) 10 (1.5 PWV) VLBI - WVR (0.5 PWV) (1.5 PWV) VLBI - Radiosonde 24 (PWV) 11 (2 PWV) 2.2: Giá tr ênh l ình bóng thám khơng, phóng x (water... (VLBI)), GPS ênh l GPS 2mm v ình gi 2mm PWV tính t D óng thám khơng k ên, h àt t B Ave Diff In Std Dev In Diff ZWD (mm) Of ZWD (mm) WVR - GPS (1 PWV) (1.5 PWV) GPS - Radiosonde 12 (2 PWV) 14 ( 2PWV) ... VƠ TUY 3.1.1 Khái ni Là thi vô thám không vô àt 30km thám không vô d theo lên cao b bóng th quy thám khơng vơ ên t ình bóng bay lên khí v thơng s nt t MHz ho Bóng thám khơng vơ l ã Trang 21 HV: