Đề tài Khảo sát phương pháp bình sai lưới trắc địa tự do và ứng dụng trong xử lý số liệu quan trắc lún công trình

vµ 0iH øng víi c¸c ®iÓm cña l­íi c¬ së (dïng lµm ®iÓm ®Þnh vÞ) cña chu kú thø i ®­îc chän lµ ®é cao ®· b×nh sai tõ chu kú thø (i-1). §Ó ý (3.1) sÏ thÊy c¸c hÖ sè mja cña ph­¬ng tr×nh sè hiÖu chØnh trÞ ®o m ®­îc x¸c ®Þnh theo nguyªn t¾c sau: -1, nÕu j=i; mja 1, nÕu j=k; 0, nÕu j kh¸c i vµ k. HÖ ph­¬ng tr×nh sè hiÖu chØnh viÕt d­íi d¹ng ma trËn: LXAV  . 3. Thµnh lËp hÖ ph­¬ng tr×nh chuÈn HÖ ph­¬ng tr×nh chuÈn Èn sè viÕt d­íi d¹ng ma trËn: 0bRX (3.2) Trong ®ã:    n i i T ii T aapPAAR 1 b 4. Chän ®iÒu kiÖn ®Þnh vÞ Nh­ ®· biÕt, hÖ ph­¬ng tr×nh ®iÒu kiÖn ®Þnh vÞ trong b×nh sai l­íi tù do cã d¹ng tæng qu¸t: 0XC T (3.3) §èi víi l­íi ®é cao tù do, sè khuyÕt d=1, do vËy C lµ ma trËn cét gåm n phÇn tö (n lµ tæng sè ®iÓm trong l­íi). C¸c phÇn tö cña ma trËn C ®­îc lùa chän nh­ ®· nãi ë môc 2.4    n i i T ii T lapPLA 1 T tHHHX )......( 21  Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 34 - Líp Tr¾c ®Þa B-K48 5. TÝnh T vµ TTT §Ó tÝnh nghiÖm vµ ®¸nh gi¸ ®é chÝnh x¸c theo ma trËn gi¶ nghÞch ®¶o, cÇn ph¶i x¸c ®Þnh thªm ma trËn T. Ma trËn T ®­îc x¸c ®Þnh theo c«ng thøc (2.7) 1)(  BCBT T Trong ®ã, B lµ ma trËn hÖ sè cña phÐp chuyÓn ®æi to¹ ®é Helmert, trong ®ã: TB )11...111( §èi víi l­íi ®é cao tù do, )( BC T lµ ma trËn mét phÇn tö, v× vËy dÔ dµng chøng minh r»ng: (3.4) Suy ra: (3.5) víi k lµ sè ®iÓm tham gia ®Þnh vÞ. 6. LËp vµ gi¶i hÖ ph­¬ng tr×nh chuÈn më réng HÖ ph­¬ng tr×nh chuÈn më réng cã d¹ng (2.4): 0 00        b K X C CR T TÝnh ma trËn gi¶ nghÞch ®¶o: TT TTCCRR  1~ )( (3.6) TÝnh ma trËn (R + CCT)-1 b×nn th­êng, sau ®ã trõ ®i ma trËn TTT. Sau khi tÝnh ma trËn gi¶ nghÞch ®¶o sÏ t×m ®­îc nghiÖm theo c«ng thøc: bRX ~ (3.7) Dùa vµo vector nghiÖm X tiÕn hµnh ph©n tÝch sù æn ®Þnh cña c¸c ®iÓm l­íi c¬ së, nÕu ph¸t hiÖn thÊy cã mét hay nhiÒu ®iÓm cã dÊu hiÖu kh«ng æn ®Þnh th× tiÕn hµnh qu¸ tr×nh tÝnh lÆp, tøc lµ chän l¹i ®iÒu kiÖn ®Þnh vÞ C (sÏ ®­îc nãi râ ë phÇn quy tr×nh xö lý l­íi). NÕu c¸c ®iÓm l­íi lµ æn ®Þnh th× tiÕn hµnh b×nh sai l­íi vµ ®¸nh gi¸ ®é chÝnh x¸c. TT BB k TT 2 1 B k T 1 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 35 - Líp Tr¾c ®Þa B-K48 7. B×nh sai Sau khi t×m ®­îc nghiÖm cña hÖ ph­¬ng tr×nh chuÈn më réng, sÏ tiÕn hµnh tÝnh: - Vector chªnh cao b×nh sai: LAXhh  0 (3.8) - Vector ®é cao b×nh sai: XHH  0 (3.9) víi h0 vµ H0 lµ vector trÞ ®o vµ vector ®é cao gÇn ®óng. 8. §¸nh gi¸ ®é chÝnh x¸c - Sai sè trung ph­¬ng träng sè ®¬n vÞ: dtn PVV T 0 (3.10) - Sai sè trung ph­¬ng cña c¸c Èn sè: ~ 0 XXX Rm  (3.11) - Sai sè trung ph­¬ng cña hµm c¸c Èn sè: fRfm TF ~0 (3.12) Víi f lµ vector hÖ sè cña hµm cÇn ®¸nh gi¸. 3.2 s¬ ®å khèi vµ quy tr×nh xö lý l­íi ®é cao ®o lón Nh­ ®· tr×nh bµy ë ch­¬ng 2, viÖc xö lý hÖ thèng l­íi ®é cao ®o lón bao gåm hai nhiÖm vô: - Ph©n tÝch hÖ thèng l­íi ®é cao c¬ së, t×m ra nh÷ng ®iÓm ®é cao gèc æn ®Þnh t¹i thêi ®iÓm xö lý l­íi. - B×nh sai tæng thÓ hÖ thèng l­íi, ®Þnh vÞ l­íi theo nh÷ng ®iÓm ®é cao gèc æn ®Þnh ®Ó x¸c ®Þnh ®é cao cña c¸c ®iÓm ®o lón g¾n trªn c«ng tr×nh. VÒ nguyªn t¾c, hai nhiÖm vô nµy cã thÓ gi¶i quyÕt t¸ch biÖt. Sau khi ph©n tÝch, t×m ra nh÷ng ®iÓm ®é cao c¬ së æn ®Þnh vµ hiÖu chØnh cho nh÷ng ®iÓm kÐm æn ®Þnh, sÏ b×nh sai l­íi quan tr¾c nh­ mét l­íi phô thuéc víi sè liÖu gèc lµ ®é cao cña c¸c ®iÓm c¬ së. Tuy nhiªn nh­ vËy, l­íi quan tr¾c Ýt nhiÒu sÏ bÞ Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 36 - Líp Tr¾c ®Þa B-K48 t¸c ®éng cña sai sè sè liÖu gèc. H¬n n÷a, nÕu sö dông m¸y tÝnh ®iÖn tö sÏ lµm gi¶m kh¶ n¨ng tù ®éng ho¸ khi xö lý sè liÖu. §Ó kh¾c phôc nh÷ng h¹n chÕ nªu trªn, gi¶i ph¸p cña chóng t«i lµ tiÕn hµnh xö lý ®ång thêi 2 bËc l­íi theo thuËt to¸n ®· tr×nh bµy ë tiÕt 3.1. Theo ®ã, hai nhiÖm vô nãi trªn khi xö lý hÖ thèng l­íi ®é cao ®o lón ®­îc chóng t«i gi¶i quyÕt nhê qu¸ tr×nh tÝnh lÆp sau ®©y: - B­íc 1: B×nh sai tæng thÓ hai cÊp l­íi theo thuËt to¸n b×nh sai l­íi tù do, ®Þnh vÞ l­íi trong hÖ thèng c¸c ®iÓm ®é cao c¬ së (®é cao gèc). - B­íc 2: Ph©n tÝch hÖ thèng l­íi c¬ së ®Ó t×m ra nh÷ng ®iÓm ®é cao gèc æn ®Þnh theo tiªu chuÈn (1.11). - B­íc 3: B×nh sai l¹i hÖ thèng l­íi, ®Þnh vÞ l¹i m¹ng l­íi theo c¸c ®iÓm ®é cao c¬ së æn ®Þnh ®· t×m ®­îc ë b­íc 2. M« h×nh tÝnh lÆp nµy ®­îc chóng t«i cô thÓ ho¸ b»ng mét quy tr×nh xö lý sè liÖu l­íi ®o lón theo thuËt to¸n b×nh sai l­íi tù do. Cô thÓ nh­ sau:  Chu kú 1 : LÊy mét ®iÓm c¬ së lµm ®iÓm khëi tÝnh ®Ó tÝnh trÞ gÇn ®óng cho ®é cao tÊt c¶ c¸c ®iÓm trong hÖ thèng l­íi quan tr¾c lón. Sau ®ã, b×nh sai chung c¶ l­íi c¬ së vµ l­íi quan tr¾c nh­ mét l­íi tù do (cã sè khuyÕt d = 0) cã l­u ý ®Õn träng sè cña c¸c trÞ ®o trong mçi bËc l­íi. Kh«ng ®Æt vÊn ®Ò ph©n tÝch ®é æn ®Þnh cña c¸c ®iÓm ®é cao c¬ së.  Chu kú 2 : B×nh sai hÖ thèng l­íi theo thuËt to¸n b×nh sai l­íi tù do víi c¸c ®iÓm ®Þnh vÞ lµ tÊt c¶ c¸c ®iÓm ®é cao gèc, kÕt hîp ph©n tÝch ®é æn ®Þnh c¸c mèc gèc theo tiªu chuÈn (1.11). Nã sÏ xÈy ra mét trong c¸c tr­êng hîp sau: - Cã mét sè ®iÓm nµo ®ã trong hÖ thèng c¸c ®iÓm gèc kh«ng æn ®Þnh : Lo¹i lÇn l­ît mét trong sè nh÷ng ®iÓm gèc kh«ng æn ®Þnh, b¾t ®Çu tõ ®iÓm kÐm æn ®Þnh nhÊt ra khái danh s¸ch c¸c ®iÓm ®Þnh vÞ ,tiÕn hµnh ®Þnh vÞ l­íi theo nh÷ng ®iÓm ®Þnh vÞ cßn l¹i. NÕu kÕt qu¶ cho thÊy nh÷ng ®iÓm ®Þnh vÞ míi nµy ®Òu æn ®Þnh th× dõng l¹i vµ chÊp nhËn kÕt qu¶ b×nh sai. - Cã mét ®iÓm gèc kh«ng æn ®Þnh: Lo¹i ®iÓm gèc nµy ra khái danh s¸ch c¸c ®iÓm gèc, tiÕn hµnh ®Þnh vÞ l­íi theo c¸c ®iÓm gèc cßn l¹i. - TÊt c¶ c¸c ®iÓm gèc ®Òu æn ®Þnh: Qu¸ tr×nh b×nh sai diÔn ra b×nh th­êng, l­íi quan tr¾c ®­îc ®Þnh vÞ theo tÊt c¶ c¸c ®iÓm gèc.  ë chu kú thø i, ®é cao cña c¸c ®iÓm gèc æn ®Þnh ®­îc gi÷ nguyªn, cßn Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 37 - Líp Tr¾c ®Þa B-K48 nh÷ng ®iÓm kh«ng æn ®Þnh sÏ nhËn gi¸ trÞ ®é cao míi ®Ó ®Þnh vÞ tiÕp cho chu kú thø (i+1). ViÖc ph©n tÝch ®é æn ®Þnh còng ®­îc thùc hiÖn nh­ ë chu kú thø 2. S¬ ®å khèi cña quy tr×nh nµy ®­îc biÓu diÔn nh­ h×nh 3.2. H×nh 3.2. S¬ ®å khèi ch­¬ng tr×nh b×nh sai l­íi ®é cao ®o lón c«ng tr×nh Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 38 - Líp Tr¾c ®Þa B-K48 3.3 lËp tr×nh b×nh sai l­íi quan tr¾c ®é lón Mét trong nh÷ng yªu cÇu c¬ b¶n cña c«ng t¸c xö lý sè liÖu tr¾c ®Þa lµ b¶o ®¶m tÝnh ®óng ®¾n cña c¸c kÕt qu¶ tÝnh to¸n vµ t¨ng nhanh tèc ®é tÝnh to¸n. TÝnh ®óng ®¾n cña qu¸ tr×nh tÝnh to¸n cã thÓ ®¶m b¶o ®­îc nhê viÖc lùa chän nh÷ng thuËt to¸n vµ quy tr×nh xö lý ®óng. Cßn ®Ó t¨ng nhanh tèc ®é tÝnh to¸n th× gi¶i ph¸p duy nhÊt lµ øng dông c¸c tiÕn bé cña c«ng nghÖ tin häc. H¬n n÷a, gi¶i ph¸p nµy còng ®ång thêi n©ng cao ®­îc chÊt l­îng cña c¸c phÐp to¸n do kh¾c phôc ®­îc nh÷ng h¹n chÕ cña c¸c ph­¬ng tiÖn tÝnh to¸n ®¬n gi¶n. §ång thêi víi viÖc ®Ò xuÊt c¸c thuËt to¸n vµ quy tr×nh xö lý phï hîp víi ®Æc ®iÓm vµ b¶n chÊt cña c¸c m¹ng l­íi quan tr¾c theo dâi ®é lón c«ng tr×nh, trong cuèn ®å ¸n nµy chóng t«i ®· tiÕn hµnh x©y dùng mét ch­¬ng tr×nh phÇn mÒm m¸y tÝnh theo c¸c thuËt to¸n vµ quy tr×nh ®· nãi ë trªn ®Ó cã thÓ tù ®éng ho¸ xö lý hÖ thèng l­íi ®é cao ®o lón. Ch­¬ng tr×nh ®­îc viÕt b»ng ng«n ng÷ lËp tr×nh VISUAL BASIC 6.0. Sau ®©y sÏ lÇn l­ît giíi thiÖu cÊu tróc cña toµn bé ch­¬ng tr×nh. 3.3.1 CÊu tróc ch­¬ng tr×nh chÝnh Phï hîp víi c¸c thuËt to¸n vµ quy tr×nh ®· x©y dùng, ch­¬ng tr×nh xö lý l­íi ®é cao ®o lón cña chóng t«i ®­îc lËp cã cÊu tróc nh­ sau: a. Thµnh lËp hÖ ph­¬ng tr×nh sè hiÖu chØnh Ph­¬ng tr×nh sè hiÖu chØnh cña c¸c chªnh cao ®o ®­îc biÓu diÔn d­íi d¹ng tuyÕn tÝnh nh­ c«ng thøc (3.1): mikm lHHv  , mp HÖ ph­¬ng tr×nh sè hiÖu chØnh viÕt d­íi d¹ng ma trËn: 111   nttnn LXAV C¸c gi¸ trÞ hÖ ph­¬ng tr×nh sè hiÖu chØnh nµy ®­îc l­u gi÷ trong c¸c m¶ng riªng cã kÝch th­íc t­¬ng øng. 1 2 1 1 2 1 222 111 1 2 1 .... .. ........ .. .. ..                                       nntttnnnnnn L L L H H H tba tba tba V V V Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 39 - Líp Tr¾c ®Þa B-K48 Trong ch­¬ng tr×nh cña chóng t«i c¸c ®iÓm ®Þnh vÞ ®­îc s¾p xÕp ë cuèi danh s¸ch, nªn Tkmmm HHHHHX )......( 121   m cét k cét Trong ®ã: m + k = t b. Thµnh lËp hÖ ph­¬ng tr×nh chuÈn Khi chuyÓn tõ hÖ ph­¬ng tr×nh sè hiÖu chØnh ®Õn hÖ ph­¬ng tr×nh chuÈn theo nguyªn lý sè b×nh ph­¬ng nhá nhÊt sÏ thu ®­îc: R . X + b = 0 Víi PLAbPAAR TtTtt   1; Thùc hiÖn phÐp nh©n 2 ma trËn chóng t«i ®· viÕt 1 modul : Nhan_MaTran (). HÖ ph­¬ng tr×nh chuÈn, ®­îc l­u gi÷ trong c¸c m¶ng 11 ;;  tttt bXR . c. Chän ®iÒu kiÖn ®Þnh vÞ Víi lÇn lÆp ®Çu tiªn, chóng t«i chän ®iÒu kiÖn ®Þnh vÞ theo tiªu chuÈn ®· nãi ë tiÕt 2.4: ic = 1 øng víi c¸c ®iÓm ®é cao c¬ së. ic = 0 øng víi c¸c ®iÓm cña l­íi quan tr¾c. Ma trËn ®Þnh vÞ C sÏ cã d¹ng nh­ sau: TC )1....1110......000( m cét k cét CÇn l­u ý r»ng ®iÒu kiÖn nµy kh«ng ph¶i lµ cè ®Þnh, nã cã thÓ sÏ ®­îc x¸c lËp l¹i nhê qu¸ tr×nh tÝnh lÆp nÕu biÕn so s¸nh cña ch­¬ng tr×nh ph¸t hiÖn thÊy cã ®iÓm ®Þnh vÞ nµo ®ã cña l­íi c¬ së kh«ng æn ®Þnh. 0 .... .. ........ .. .. 1 2 1 1 2 1 21 22221 11211                               tntttttttt t t b b b H H H RRR RRR RRR Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 40 - Líp Tr¾c ®Þa B-K48 d. TÝnh ma trËn gi¶ nghÞch ®¶o Ma trËn gi¶ nghÞch ®¶o ®­îc tÝnh theo c«ng thøc (3.6) TT TTCCRR  1~ )( Thùc hiÖn phÐp lÊy tæng 2 ma trËn chóng t«i ®· viÕt modul: Cong_matran () §Ó x¸c ®Þnh 1~ )(  TCCRR chóng t«i thùc hiÖn phÐp nghÞch ®¶o ma trËn b×nh th­êng, phÐp nghÞch ®¶o ma trËn ®­îc chóng t«i viÕt bëi modul : Nghichdao_matran (). Cuèi cïng, ma trËn gi¶ nghÞch ®¶o ®­îc x¸c ®Þnh nhê phÐp lÊy hiÖu 2 ma trËn theo (3.6), trong ®ã TTT lµ ma trËn cã c¸c phÇn tö b»ng nhau vµ b»ng 1/k2. Ma trËn TTT ®­îc l­u gi÷ trong 1 m¶ng cã kÝch th­íc t­¬ng øng. Thùc hiÖn phÐp lÊy hiÖu 2 ma trËn ë ®©y chóng t«i viÕt thªm modul DoiDau_MaTran (), C¸c phÇn tö cña ma trËn gi¶ nghÞch ®¶o ®­îc l­u gi÷ trong m¶ng riªng. Sö dông ma trËn gi¶ nghÞch ®¶o ®Ó tÝnh nghiÖm theo (3.7). B×nh sai vµ ®¸nh gi ¸®é chÝnh x¸c ®­îc thùc hiÖn theo c¸c c«ng thøc tõ (3.8) ®Õn (3.12). Nh×n chung sau khi viÕt c¸c modul vÒ ma trËn nh­: Nghichdao_matran(), Cong_matran(), Nhan_matran(), Doidau_matran(), Chuyenvi_matran(), chóng t«i thùc hiÖn c¸c b­íc cña bµi to¸n b×nh sai th«ng th­êng. 3.3.2 CÊu tróc d÷ liÖu cña ch­¬ng tr×nh §Ó ch­¬ng tr×nh cã thÓ ho¹t ®éng ®­îc víi nh÷ng sè liÖu cña l­íi ®· cho, chóng t«i ®· tæ chøc tÖp d÷ liÖu cho ch­¬ng tr×nh víi cÊu tróc nh­ sau:            ~~ 2 ~ 1 ~ 2 ~ 22 ~ 21 ` 1 ~ 12 ~ 11 ~ .. ........ .. .. tttt t t RRR RRR RRR R Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 41 - Líp Tr¾c ®Þa B-K48 a. X¸c ®Þnh c¸c tham sè kü thuËt cña l­íi - Tªn c«ng tr×nh. - Chu kú ®o, Thêi ®iÓm ®o - Sè l­îng ®iÓm ®Þnh vÞ. - Sè l­îng ®iÓm quan tr¾c. - Sè l­îng chªnh cao ®o. b. M· ho¸ th«ng tin l­íi - C¸c ®iÓm cña l­íi ®­îc ®¸nh sè liªn tôc tõ 1 ®Õn hÕt theo nguyªn t¾c: c¸c ®iÓm cña l­íi quan tr¾c ®­îc ®¸nh sè tr­íc, sau ®ã ®Õn c¸c ®iÓm cña l­íi c¬ së. - Tªn c¸c ®iÓm cña l­íi ®­îc l­u gi÷ trong m¶ng riªng còng theo nguyªn t¾c trªn. - Gi¸ trÞ ®é cao gèc æn ®Þnh ®­îc l­u trong mét m¶ng kh¸c theo thø tù t¨ng dÇn cña m· sè c¸c ®iÓm gèc. - C¸c chªnh cao ®o ®­îc m· ho¸ theo nguyªn t¾c: m· sè cña ®iÓm ®Çu, m· sè cña ®iÓm cuèi, gi¸ trÞ chªnh cao ®o, sè tr¹m ®o. 3.4 ch­¬ng tr×nh nguån vµ tÖp d÷ liÖu Trong tr¾c ®Þa, c¸c ng«n ng÷ lËp tr×nh gióp cho chóng ta gi¶i quyÕt nh÷ng bµi to¸n cã khèi l­îng tÝnh to¸n lín mét c¸ch ®¬n gi¶n vµ nhanh chãng. Víi ­u thÕ cã nhiÒu tÝnh n¨ng, c¸c ®iÒu khiÓn míi cho phÐp ta viÕt ch­¬ng tr×nh øng dông kÕt hîp víi giao diÖn, c¸ch xö lý vµ tÝnh n¨ng cña Office, nªn ng«n ng÷ lËp tr×nh VisualBasic 6.0 (VB 6.0) lµ c«ng cô h÷u Ých gióp chóng ta trong c«ng t¸c xö lý sè liÖu tr¾c ®Þa. §Æc biÖt trong phÇn mÒm chóng t«i viÕt, chóng t«i ®· khai th¸c ­u ®iÓm m¶ng ®éng cña VB 6.0. M¶ng nµy cã thÓ thay ®æi kÝch cì, lµ mét trong nh÷ng ­u ®iÓm cña VB 6.0, m¶ng ®éng gióp qu¶n lý bé nhí mét c¸ch hiÖu qu¶. Ta cã thÓ dïng mét m¶ng lín trong thêi gian ng¾n, sau ®ã xo¸ bá ®Ó tr¶ vïng nhí cho hÖ thèng. Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 42 - Líp Tr¾c ®Þa B-K48 PhÇn mÒm chóng t«i viÕt gåm 4 th­ môc: - B×nh sai l­íi tù do - B×nh sai l­íi truyÒn thèng. - TÝnh c¸c th«ng sè lón - Dù b¸o lón. 3.4.1 Giíi thiÖu tÖp d÷ liÖu cña ch­¬ng tr×nh (vÝ dô) - NhËp tªn c«ng tr×nh. - NhËp chu kú quan tr¾c. - NhËp ngµy quan tr¾c. 4 25 38 0.95 1 2 29.07 1 2 3 -70.18 1 3 4 17.62 1 4 5 63.59 1 5 6 -87.44 1 6 7 95.13 1 7 8 0.06 1 8 9 -100.65 1 9 10 -28.69 1 .......... 23 9 84.45 1 20 22 63.75 1 22 21 61.74 1 21 24 -114.7 2 26 29 -226.03 5 26 27 22.47 3 29 27 248.81 3 28 27 209.58 4 1 N1 2 N2 3 N3 . . . . 24 N24 25 N25 26 M1 8618.26 27 TC 8641.7 28 M2 8429.74 29 M3 8389.9 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 43 - Líp Tr¾c ®Þa B-K48 3.4.2 Giíi thiÖu ch­¬ng tr×nh nguån Sau ®©y chóng t«i xin giíi thiÖu mét phÇn cña ch­¬ng tr×nh nguån cña modul B×nh sai l­íi tù do: Private Sub mnuBinhsaitudo_Click() Dim pass As Boolean pass = False Dim DaLap As Boolean DaLap = False Dim tencongtrinh As String tencongtrinh = txtTencongtrinh Dim Chukydo As String Chukydo = txtChukydo Dim Ngaydo As String Ngaydo = txtngaydo Dim ncs As Long, nqt As Long, nh As Long ncs = VSGrid.TextMatrix(0, 0) nqt = VSGrid.TextMatrix(0, 1) nh = VSGrid.TextMatrix(0, 2) Dim nqtP As Long, ncsP As Long nqtP = nqt ncsP = ncs Dim lanlap As Long lanlap = 0 Dim tcod As Double tcod = VSGrid.TextMatrix(0, 3) Dim id() As Long Dim ic() As Long Dim h() As Double Dim nt() As Long Dim i, j As Long ReDim id(1 To nh), ic(1 To nh), h(1 To nh), nt(1 To nh) For i = 1 To nh id(i) = VSGrid.TextMatrix(i, 0) ic(i) = VSGrid.TextMatrix(i, 1) h(i) = VSGrid.TextMatrix(i, 2) nt(i) = VSGrid.TextMatrix(i, 3) Next Dim Tendiem() As String ReDim Tendiem(1 To nqt + ncs) For i = 1 To (nqt + ncs) Tendiem(i) = VSGrid.TextMatrix(nh + i, 1) Next Dim hgd() As Double ReDim hgd(1 To nqt + ncs) For i = 1 To nqt hgd(i) = 0 Next j = 0 For i = (nqt + 1) To (nqt + ncs) j = j + 1 hgd(i) = VSGrid.TextMatrix(nh + nqt + j, 2) Next Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 44 - Líp Tr¾c ®Þa B-K48 Dim c() As Double ReDim c(1 To nqt + ncs, 1 To 1) For i = 1 To nqt c(i, 1) = 0 Next For i = nqt + 1 To nqt + ncs c(i, 1) = 1 Next Dim solantinhlap As Long With CommonDialog1 .DialogTitle = "In file nao" .CancelError = True .Filter = "Text file *.txt|*.txt|All Files | *.*" .ShowSave End With str = CommonDialog1.FileName Dim FileNum As Long FileNum = FreeFile Open str For Output As FileNum Print #FileNum, String(16, " ") & "KET QUA TINH TOAN BINH SAI LUOI DO CAO DO LUN" Print #FileNum, String(14, " ") & "*************************************************" Print #FileNum, String(12, " ") & tencongtrinh Print #FileNum, String(19, " ") & "CHU KY: " & Chukydo & " - " & "NGAY DO: " & Ngaydo Print #FileNum, Print #FileNum, Print #FileNum, String(20, " ") & "CAC CHI TIEU KY THUAT CUA LUOI :" Print #FileNum, String(20, " ") & "==============================" Print #FileNum, String(20, " ") & "1." & "So luong diem quan trac:" & String(6, " ") & nqt Print #FileNum, String(20, " ") & "2." & "So luong diem DINH VI :" & String(6, " ") & ncs Print #FileNum, String(20, " ") & "3." & "So luong chenh cao do:" & String(6, " ") & nh solantinhlap = 1 Print #FileNum, Print #FileNum, Print #FileNum, String(7, " ") & "PHAN TICH DO ON DINH CUA CAC DIEM CO SO:" Print #FileNum, String(7, " ") & "===============================" 7: Dim mtA() As Double ReDim mtA(1 To nh, 1 To nqt + ncs) For i = 1 To nh For j = 1 To (nqt + ncs) mtA(i, j) = 0 Next mtA(i, id(i)) = -1 mtA(i, ic(i)) = 1 Next Dim l() As Double ReDim l(1 To nh, 1 To 1) For i = 1 To nh l(i, 1) = (hgd(ic(i)) - hgd(id(i))) - h(i) Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 45 - Líp Tr¾c ®Þa B-K48 Next Dim p() As Double ReDim p(1 To nh, 1 To nh) For i = 1 To nh p(i, i) = 1 / nt(i) Next Dim R() As Double ReDim R(1 To (nqt + ncs), 1 To (nqt + ncs)) R = Nhan_MaTran(chuyenvi_matran(mtA(), nh, (nqt + ncs)), _ Nhan_MaTran(p(), mtA(), nh, nh, (nqt + ncs)), (nqt + ncs), nh, (nqt + ncs)) Dim b() As Double ReDim b(1 To nh, 1 To 1) b = Nhan_MaTran(chuyenvi_matran(mtA(), nh, (nqt + ncs)), _ Nhan_MaTran(p(), l(), nh, nh, 1), (nqt + ncs), nh, 1) Dim CTC() As Double ReDim CTC(1 To nqt + ncs, 1 To nqt + ncs) CTC = Nhan_MaTran(c(), chuyenvi_matran(c(), nqt + ncs, 1), nqt + ncs, 1, nqt + ncs) Dim cc() As Double ReDim cc(1 To nqt + ncs, 1 To nqt + ncs) For i = 1 To nqt + ncs For j = 1 To nqt + ncs cc(i, j) = 1000000 * CTC(i, j) Next Next Dim r_c() As Double ReDim r_c(1 To ncs + nqt, 1 To ncs + nqt) r_c = cong_matran(R(), cc(), nqt + ncs, nqt + ncs) Dim Rn() As Double ReDim Rn(1 To nqt + ncs, 1 To nqt + ncs) Rn = nghichdao_matran(r_c(), nqt + ncs) Dim x() As Double ReDim x(1 To nqt + ncs, 1 To 1) x = DoiDau_MaTran(Nhan_MaTran(Rn(), b(), nqt + ncs, nqt + ncs, 1), nqt + ncs, 1) . . . . . . . . . Dim CTX() As Double ReDim CTX(1 To nqt + ncs, 1 To 1) CTX = Nhan_MaTran(chuyenvi_matran(c(), nqt + ncs, 1), x(), 1, nqt + ncs, 1) For i = nqt + 1 To nqt + ncs If Abs(x(i, 1)) > tcod Then pass = True End If Next If pass = False Then GoTo 6 End If pass = False 5: max = Abs(x(nqt + 1, 1)) imax = nqt + 1 For i = nqt + 1 To nqt + ncs For j = 1 To 1 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 46 - Líp Tr¾c ®Þa B-K48 If (c(i, 1) 0) And Abs(x(i, j)) > tcod And Abs(x(i, j)) > max Then max = Abs(x(i, j)) imax = i Else max = max End If Next Next Print #FileNum, String(7, " ") & "Loai diem " & Tendiem(imax) & " ra khoi danh sach cac diem DINH VI" 9: Dim ViTriDiem1 For i = 1 To (nqt + ncs) If c(i, 1) = 1 Then ViTriDiem1 = i Exit For End If Next c(imax, 1) = 0 Dim TG As Double Dim tg1 As String If ViTriDiem1 = imax Then GoTo Kodoivitri End If TG = c(imax, 1) c(imax, 1) = c(ViTriDiem1, 1) c(ViTriDiem1, 1) = TG tg1 = Tendiem(imax) Tendiem(imax) = Tendiem(ViTriDiem1) Tendiem(ViTriDiem1) = tg1 TG = hgd(imax) hgd(imax) = hgd(ViTriDiem1) hgd(ViTriDiem1) = TG For i = 1 To nh Select Case id(i) Case Is = imax id(i) = ViTriDiem1 Case Is = ViTriDiem1 id(i) = imax End Select Select Case ic(i) Case Is = imax ic(i) = ViTriDiem1 Case Is = ViTriDiem1 ic(i) = imax End Select Next . . . . . . . . . . . . Dim Hbs() As Double ReDim Hbs(1 To nqt + ncs, 1 To 1) Hbs = cong_matran(ho(), x(), nqt + ncs, 1) Dim V() As Double ReDim V(1 To nh, 1 To 1) V = cong_matran(Nhan_MaTran(mtA(), x(), nh, nqt + ncs, 1), l(), nh, 1) Dim hhbs() As Double Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 47 - Líp Tr¾c ®Þa B-K48 ReDim hhbs(1 To nh, 1 To 1) hhbs = cong_matran(hh(), V(), nh, 1) 'Sai so trung phuong trong so don vi Dim vpv() As Double Dim mo() As Double ReDim vpv(1 To 1, 1 To 1) ReDim mo(1 To 1, 1 To 1) vpv = Nhan_MaTran(Nhan_MaTran(chuyenvi_matran(V(), nh, 1), p(), 1, nh, nh), V(), 1, nh, 1) For i = 1 To 1 For j = 1 To 1 mo(i, j) = Sqr(vpv(i, j) / (nh - (nqt + ncs) + 1)) Next Next Dim Rxx() As Double ReDim Rxx(1 To nqt + ncs, 1 To 1) For i = 1 To nqt + ncs Rxx(i, 1) = Sqr(Rn(i, i)) Next Dim mx() As Double ReDim mx(1 To nqt + ncs, 1 To 1) mx = Nhan_MaTran(Rxx(), mo(), nqt + ncs, 1, 1) Dim fRf() As Double ReDim fRf(1 To nh, 1 To nh) fRf = Nhan_MaTran(Nhan_MaTran(mtA(), Rn(), nh, nqt + ncs, nqt + ncs), chuyenvi_matran(mtA(), nh, nqt + ncs), nh, nqt + ncs, nh) Dim cfRf() As Double ReDim cfRf(1 To nh, 1 To 1) For i = 1 To nh cfRf(i, 1) = Sqr(fRf(i, i)) Next Dim mf() As Double ReDim mf(1 To nh, 1 To 1) mf = Nhan_MaTran(cfRf(), mo(), nh, 1, 1) Print #FileNum, Print #FileNum, String(7, " ") & "Lay " & ncs & " diem o bang " & solantinhlap & " de dinh vi luoi!" Print #FileNum, String(7, " ") & "====================================" Print #FileNum, Print #FileNum, Print #FileNum, String(7, " ") & "BINH SAI VA VA DANH GIA DO CHINH XAC: " Print #FileNum, String(7, " ") & "===============================" Print #FileNum, Print #FileNum, Print #FileNum, String(7, " ") & "Tri do sau binh sai va do chinh xac cua tri do :" Print #FileNum, String(7, " ") & "================================================" i = solantinhlap + 1 Print #FileNum, String(60, " ") & "BANG: " & i Print #FileNum, String(7, " ") & "==============================================================" Print #FileNum, String(7, " ") & "| T. | TUYEN DO |CH. CAO DO| Vh | SO |CH. CAO SAU| Mh |" Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 48 - Líp Tr¾c ®Þa B-K48 Print #FileNum, String(7, " ") & "| T. | DAU-CUOI | (mm) | (mm) |TRAM | BINH SAI | (mm) |" Print #FileNum, String(7, " ") & "==============================================================" j = 1 3: For i = 1 To nh Print #FileNum, String(7, " ") & "|" & " " & Dinhdangso(i, 3, 0, True, False) & "|" & " " _ & Dinhdangso(Tendiem(id(j)), 6, 4, True, False) _ & Dinhdangso(Tendiem(ic(j)), 5, 4, True, False) _ & "|" & " " & Dinhdangso(h(i), 9, 2, True, True) & "|" & " " _ & Dinhdangso(V(i, 1), 5, 2, True, True) & "|" _ & String(1, " ") & String(1, " ") & nt(i) & String(2, " ") & "|" & " " & Dinhdangso(hhbs(i, 1), 10, 2, True, True) _ & "|" & String(2, " ") & Dinhdangso(mf(i, 1), 4, 2, True, False) & "|" j = j + 1 Next If j < nh + 1 Then GoTo 3 End If Print #FileNum, String(7, " ") & "==============================================================" Print #FileNum, Print #FileNum, Print #FileNum, String(7, " ") & "Do cao binh sai va sai so cua chung :" Print #FileNum, String(7, " ") & "===============================" j = solantinhlap + 2 Print #FileNum, String(60, " ") & "BANG: " & j Print #FileNum, String(7, " ") & "==============================================================" Print #FileNum, String(7, " ") & "| T. | TEN DIEM | DO CAO | Mx | GHI |" Print #FileNum, String(7, " ") & "| T. | | (mm) | (mm) | CHU |" Print #FileNum, String(7, " ") & "==============================================================" For i = 1 To nqt Print #FileNum, String(7, " ") & "|" & " " & Dinhdangso(i, 3, 0, True, False) & "|" & String(5, " ") & Dinhdangso(Tendiem(i), 9, 4, True, False) & "|" & String(3, " ") & Dinhdangso(Hbs(i, 1), 12, 2, True, False) & "|" & String(2, " ") & Dinhdangso(mx(i, 1), 6, 2, True, False) & "|" & String(15, " ") & "|" Next For i = nqt + 1 To nqt + ncs Print #FileNum, String(7, " ") & "|" & " " & Dinhdangso(i, 3, 0, True, False) & "|" & String(5, " ") & Dinhdangso(Tendiem(i), 9, 4, True, False) & "|" & String(3, " ") & Dinhdangso(hgd(i), 12, 2, True, False) & "|" & String(2, " ") & Dinhdangso(mx(i, 1), 6, 2, True, False) & "|" & String(15, " ") & "|" Next Print #FileNum, String(7, " ") & "==============================================================" Print #FileNum, Print #FileNum, String(7, " ") & "Sai so trung phuong trong so don vi :" & " " & Dinhdangso(mo(1, 1), 4, 2, False, False) & "mm" Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 49 - Líp Tr¾c ®Þa B-K48 Print #FileNum, Print #FileNum, Print #FileNum, String(7, " ") & "Chuong trinh duoc viet bang ngon ngu lap trinh VISUALBASIC6.0" Close FileNum RichTextBox1.LoadFile CommonDialog1.FileName End Sub 3.5 sö dông ch­¬ng tr×nh 1. Giao diÖn chÝnh cña ch­¬ng tr×nh: H×nh 3.3. H×nh 3.3. Giao diÖn chÝnh cña ch­¬ng tr×nh 2. Tool bar C¸c c«ng cô cña Toolbar lµ New, Open, Save, Copy, Cut, Past cã thÓ kÕt nèi víi c¸c file d¹ng text hay Excel. H×nh 3.4. C¸c c«ng cô cña Toolbar Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 50 - Líp Tr¾c ®Þa B-K48 3. Menu File cã c¸c tÝnh n¨ng nh­ sau : Chóng ta cã thÓ t¹o mét file sè liÖu míi, save mét file sè liÖu, hoÆc më mét file sè liÖu ®· cã s½n. H×nh 3.5. C¸c tÝnh n¨ng cña Menu File Khi më mét file sè liÖu th× cöa sæ Open xuÊt hiÖn (H×nh 3.6), chóng ta theo ®­êng dÉn ®Ó më file. H×nh 3.6. Cöa sæ Open grid Khi thùc hiÖn Save mét file sè liÖu th× cöa sæ Save grid xuÊt hiÖn (H×nh 3.7), chóng ta thùc hiÖn qu¸ tr×nh Save b×nh th­êng. H×nh 3.7. Cöa sæ Save grid Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 51 - Líp Tr¾c ®Þa B-K48 4. Menu Edit cã c¸c tÝnh n¨ng nh­ sau : Chóng ta cã thÓ thùc hiÖn Cut, Copy , Delete vµ Paste trong b¶ng sè liÖu, H×nh 3.8. C¸c tÝnh n¨ng cña Menu Edit C¸c tÝnh n¨ng cña Menu Edit cã thÓ kÕt nèi víi Ecxel, Word. 5. Menu Service cã c¸c tÝnh n¨ng nh­ sau : H×nh 3.9. C¸c tÝnh n¨ng cña Menu Service 3.6 tÝnh to¸n thùc nghiÖm §Ó minh chøng cho tÝnh ®óng ®¾n khi øng dông ph­¬ng ph¸p b×nh sai l­íi tr¾c ®Þa tù do ®Ó xö lý sè liÖu quan tr¾c lón c«ng tr×nh, vµ tÝnh ®óng ®¾n cña ch­¬ng tr×nh mµ chóng t«i ®· thµnh lËp, chóng t«i sö dông ch­¬ng tr×nh ®Ó tÝnh to¸n b×nh sai mét chu kú cña l­íi quan tr¾c lón nhµ §¬n nguyªn I - thuéc L« 5 khu chung c­ §Þnh c«ng - Hµ Néi. L­íi quan tr¾c gåm 4 ®iÓm c¬ së ( M1, TC , M2, M3), vµ 25 ®iÓm quan tr¾c Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 52 - Líp Tr¾c ®Þa B-K48 (thø tù tªn ®iÓm tõ N1 ®Õn N25). §o ®¹c l­íi b»ng m¸y Ni 004 vµ mia Invar. §é cao b×nh sai cña c¸c mèc c¬ së trong chu kú 10 thu ®­îc ë b¶ng 3.1 S¬ ®å l­íi khèng chÕ c¬ së: H×nh 3.10 MÆt b»ng bè trÝ mèc ®o lón: H×nh 3.11 B¶ng 3.1 §é cao sau b×nh sai c¸c ®iÓm l­íi c¬ së ë chu kú 10. H×nh 3.10 S¬ ®å l­íi khèng chÕ c¬ së TT Tªn ®iÓm §é cao sau b×nh sai (mm) 1 M1 8618.26 2 TC 8641.70 3 M2 8429.74 4 M3 8389.90 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 53 - Líp Tr¾c ®Þa B-K48 H×nh 3.11 S¬ ®å bè trÝ mèc ®o lón Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 54 - Líp Tr¾c ®Þa B-K48 File sè liÖu : 4 25 38 0.95 1 2 29.07 1 2 3 -70.18 1 3 4 17.62 1 4 5 63.59 1 5 6 -87.44 1 6 7 95.13 1 7 8 0.06 1 8 9 -100.65 1 9 10 -28.69 1 10 11 29.08 1 11 12 45.63 1 12 13 -44.86 1 13 14 -31.49 1 14 15 43.9 1 15 16 -78.92 1 16 17 -2.3 1 17 1 120.32 1 20 1 91.49 3 20 18 -21.77 1 18 19 193.36 2 19 13 -131.81 2 20 10 9.87 1 10 23 -56.43 1 23 9 84.45 1 20 22 63.75 1 22 21 61.74 1 21 24 -114.7 2 24 25 178.66 1 25 6 -146.43 3 28 5 537.41 2 29 1 537.08 1 27 10 206.51 3 26 28 -186.26 4 28 29 -39.46 2 26 29 -226.03 5 26 27 22.47 3 29 27 248.81 3 28 27 209.58 4 1 N1 2 N2 3 N3 4 N4 5 N5 6 N6 7 N7 8 N8 9 N9 10 N10 11 N11 12 N12 13 N13 14 N14 15 N15 16 N16 17 N17 18 N18 19 N19 20 N20 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 55 - Líp Tr¾c ®Þa B-K48 21 N21 22 N22 23 N23 24 N24 25 N25 26 M1 8618.26 27 TC 8641.7 28 M2 8429.74 29 M3 8389.9 Sö dông ch­¬ng tr×nh xö lý sè liÖu chóng t«i viÕt ®Ó tÝnh to¸n b×nh sai, chóng t«i ®­îc kÕt qu¶ nh­ sau: KET QUA TINH TOAN BINH SAI LUOI DO CAO DO LUN ************************************************* CONG TRINH: DON NGUYEN II - THUOC LO 5 - KDTM DINH CONG CHU KY: 11 - NGAY DO: 07/09/2006 CAC CHI TIEU KY THUAT CUA LUOI : ================================ 1.So luong diem quan trac: 25 2.So luong diem DINH VI : 4 3.So luong chenh cao do: 38 PHAN TICH DO ON DINH CUA CAC DIEM CO SO: ======================================== Phan tich lan : 1 =================== BANG: 1 ============================================================== | TEN | DO CAO | DO LECH | KET QUA | | DIEM | (mm) | CAO DO | PHAN TICH | ============================================================== | M1 | 8618.26 | -.88 | On dinh | | TC | 8641.70 | -1.56 | Khong on dinh | | M2 | 8429.74 | +1.05 | Khong on dinh | | M3 | 8389.90 | +1.39 | Khong on dinh | ============================================================== Loai diem TC ra khoi danh sach cac diem DINH VI Phan tich lan : 2 =================== BANG: 2 ============================================================== | TEN | DO CAO | DO LECH | KET QUA | | DIEM | (mm) | CAO DO | PHAN TICH | ============================================================== | M1 | 8618.26 | -1.40 | Khong on dinh | | M2 | 8429.74 | +.53 | On dinh | | M3 | 8389.90 | +.87 | On dinh | ============================================================== Loai diem M1 ra khoi danh sach cac diem DINH VI Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 56 - Líp Tr¾c ®Þa B-K48 Phan tich lan : 3 =================== BANG: 3 ============================================================== | TEN | DO CAO | DO LECH | KET QUA | | DIEM | (mm) | CAO DO | PHAN TICH | ============================================================== | M2 | 8429.74 | -.17 | On dinh | | M3 | 8389.90 | +.17 | On dinh | ============================================================== Lay 2 diem o bang 3 de dinh vi luoi! ==================================== BINH SAI VA VA DANH GIA DO CHINH XAC: ===================================== Tri do sau binh sai va do chinh xac cua tri do : ================================================ BANG: 4 ============================================================== | T. | TUYEN DO |CH. CAO DO| Vh | SO |CH. CAO SAU| Mh | | T. | DAU-CUOI | (mm) | (mm) |TRAM | BINH SAI | (mm) | ============================================================== | 1 | N1 N2 | +29.07 | -.06 | 1 | +29.01 | .16 | | 2 | N2 N3 | -70.18 | -.06 | 1 | -70.24 | .16 | | 3 | N3 N4 | +17.62 | -.06 | 1 | +17.56 | .16 | | 4 | N4 N5 | +63.59 | -.06 | 1 | +63.53 | .16 | | 5 | N5 N6 | -87.44 | -.07 | 1 | -87.51 | .16 | | 6 | N6 N7 | +95.13 | +.01 | 1 | +95.14 | .16 | | 7 | N7 N8 | +.06 | +.01 | 1 | +.07 | .16 | | 8 | N8 N9 | -100.65 | +.01 | 1 | -100.64 | .16 | | 9 | N9 N10 | -28.69 | +.23 | 1 | -28.46 | .14 | | 10 | N10 N11 | +29.08 | +.02 | 1 | +29.10 | .16 | | 11 | N11 N12 | +45.63 | +.02 | 1 | +45.65 | .16 | | 12 | N12 N13 | -44.86 | +.02 | 1 | -44.84 | .16 | | 13 | N13 N14 | -31.49 | +.02 | 1 | -31.47 | .17 | | 14 | N14 N15 | +43.90 | +.02 | 1 | +43.92 | .17 | | 15 | N15 N16 | -78.92 | +.02 | 1 | -78.90 | .17 | | 16 | N16 N17 | -2.30 | +.02 | 1 | -2.28 | .17 | | 17 | N17 N1 | +120.32 | +.02 | 1 | +120.34 | .17 | | 18 | N20 N1 | +91.49 | -.12 | 3 | +91.37 | .21 | | 19 | N20 N18 | -21.77 | .00 | 1 | -21.77 | .17 | | 20 | N18 N19 | +193.36 | -.01 | 2 | +193.35 | .22 | | 21 | N19 N13 | -131.81 | -.01 | 2 | -131.82 | .22 | | 22 | N20 N10 | +9.87 | -.04 | 1 | +9.83 | .15 | | 23 | N10 N23 | -56.43 | +.22 | 1 | -56.21 | .14 | | 24 | N23 N9 | +84.45 | +.22 | 1 | +84.67 | .14 | | 25 | N20 N22 | +63.75 | +.09 | 1 | +63.84 | .17 | | 26 | N22 N21 | +61.74 | +.09 | 1 | +61.83 | .17 | | 27 | N21 N24 | -114.70 | +.17 | 2 | -114.53 | .23 | | 28 | N24 N25 | +178.66 | +.09 | 1 | +178.75 | .17 | | 29 | N25 N6 | -146.43 | +.26 | 3 | -146.17 | .26 | | 30 | M2 N5 | +537.41 | -.02 | 2 | +537.39 | .20 | | 31 | M3 N1 | +537.08 | -.04 | 1 | +537.04 | .16 | | 32 | TC N10 | +206.51 | +.15 | 3 | +206.66 | .23 | ============================================================== Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 57 - Líp Tr¾c ®Þa B-K48 BANG: 4 (Tiep) ============================================================== | T. | TUYEN DO |CH. CAO DO| Vh | SO |CH. CAO SAU| Mh | | T. | DAU-CUOI | (mm) | (mm) |TRAM | BINH SAI | (mm) | ============================================================== | 33 | M1 M2 | -186.26 | -.33 | 4 | -186.59 | .23 | | 34 | M2 M3 | -39.46 | -.04 | 2 | -39.50 | .18 | | 35 | M1 M3 | -226.03 | -.06 | 5 | -226.09 | .24 | | 36 | M1 TC | +22.47 | +.29 | 3 | +22.76 | .23 | | 37 | M3 TC | +248.81 | +.04 | 3 | +248.85 | .20 | | 38 | M2 TC | +209.58 | -.23 | 4 | +209.35 | .21 | ============================================================== Do cao binh sai va sai so cua chung : ===================================== BANG: 5 ============================================================== | T. | TEN DIEM | DO CAO | Mx | GHI | | T. | | (mm) | (mm) | CHU | ============================================================== | 1 | N1 | 8927.11 | .17 | | | 2 | N2 | 8956.12 | .21 | | | 3 | N3 | 8885.88 | .23 | | | 4 | N4 | 8903.43 | .22 | | | 5 | N5 | 8966.96 | .20 | | | 6 | N6 | 8879.45 | .23 | | | 7 | N7 | 8974.59 | .25 | | | 8 | N8 | 8974.67 | .26 | | | 9 | N9 | 8874.03 | .24 | | | 10 | N10 | 8845.58 | .22 | | | 11 | N11 | 8874.68 | .26 | | | 12 | N12 | 8920.34 | .28 | | | 13 | N13 | 8875.50 | .27 | | | 14 | N14 | 8844.03 | .28 | | | 15 | N15 | 8887.95 | .28 | | | 16 | N16 | 8809.05 | .26 | | | 17 | N17 | 8806.77 | .23 | | | 18 | N18 | 8813.97 | .27 | | | 19 | N19 | 9007.32 | .30 | | | 20 | N20 | 8835.75 | .23 | | | 21 | N21 | 8961.41 | .30 | | | 22 | N22 | 8899.58 | .27 | | | 23 | N23 | 8789.36 | .25 | | | 24 | N24 | 8846.88 | .32 | | | 25 | N25 | 9025.62 | .31 | | | 26 | TC | 8638.92 | .18 | | | 27 | M1 | 8616.16 | .22 | | | 28 | M2 | 8429.74 | .09 | | | 29 | M3 | 8389.90 | .09 | | ============================================================== Sai so trung phuong trong so don vi : .18mm Chuong trinh duoc viet bang ngon ngu lap trinh VISUALBASIC6.0 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 58 - Líp Tr¾c ®Þa B-K48 KÕt luËn Tõ kÕt qu¶ nghiªn cøu lý thuyÕt vµ c¸c tÝnh to¸n øng dông thùc tiÔn ®· ®­îc tr×nh bµy trong néi dung cña ®å ¸n ®· cho phÐp chóng t«i rót ra nh÷ng kÕt luËn sau ®©y: 1. Khi xö lý sè liÖu ®o lón c«ng tr×nh cÇn kÕt hîp ph©n tÝch ®é æn ®Þnh cña c¸c mèc ®é cao c¬ së. ChØ nh÷ng ®iÓm ®é cao c¬ së æn ®Þnh míi ®­îc tham gia vµo viÖc ®Þnh vÞ m¹ng l­íi quan tr¾c. Nh­ vËy ®é lón thùc tÕ cña c¸c c«ng tr×nh míi ®­îc ph¶n ¸nh mét c¸ch kh¸ch quan. 2.ViÖc sö dông ph­¬ng ph¸p b×nh sai tù do ®Ó xö lý hÖ thèng l­íi ®é cao ®o lón lµ mét gi¶i ph¸p toµn diÖn, phï hîp víi ®Æc ®iÓm vµ b¶n chÊt cña d¹ng l­íi nµy. 3. C¸c kÕt qu¶ nghiªn cøu lý thuyÕt ®· ®­îc chóng t«i cô thÓ ho¸ b»ng viÖc x©y dùng mét ch­¬ng tr×nh phÇn mÒm ®Ó gi¶i quyÕt c¸c nhiÖm vô ®Æt ra trªn m¸y tÝnh. So s¸nh kÕt qu¶ b×nh sai cña ch­¬ng tr×nh chóng t«i viÕt vµ c¸c kÕt qu¶ cña c¸c thÇy c« trong bé m«n lµ hoµn toµn gièng nhau. §iÒu nµy cã thÓ kh¼ng ®Þnh thuËt to¸n, quy tr×nh xö lý vµ ch­¬ng tr×nh b×nh sai mµ chóng t«i ®· tr×nh bµy trong ®å ¸n lµ ®óng ®¾n, ®¸p øng ®­îc ®Çy ®ñ nh÷ng yªu cÇu kü thuËt vÒ xö lý sè liÖu quan tr¾c lón c«ng tr×nh. Xö lý sè liÖu ®o biÕn d¹ng c«ng tr×nh nãi chung vµ ®o lón c«ng tr×nh nãi riªng cã vai trß rÊt quan träng, nªn vÊn ®Ò nµy ®ang ®­îc nhiÒu ng­êi quan t©m. NhËn thøc ®­îc tÇm quan träng cña vÊn ®Ò nµy, trong thêi gian lµm ®å ¸n tèt nghiÖp kho¸ häc em ®· cè g¾ng rÊt nhiÒu ®Ó hoµn thiÖn quy tr×nh c«ng nghÖ c«ng t¸c tr¾c ®Þa trong xö lý sè liÖu quan tr¾c lón c«ng tr×nh. Hy väng nhËn ®­îc sù gãp ý cña c¸c thÇy c« vµ c¸c b¹n ®ång nghiªp ®Ó ch­¬ng tr×nh cña chóng t«i hoµn thiÖn h¬n. Mét lÇn n÷a em xin ch©n thµnh c¶m ¬n thÇy gi¸o TS. NguyÔn Quang Phóc, c¸c thÇy c« gi¸o trong khoa tr¾c ®Þa vµ c¸c b¹n ®ång nghiÖp ! Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 59 - Líp Tr¾c ®Þa B-K48 Tµi liÖu tham kh¶o [1]. Phan V¨n HiÕn (1997), Quan tr¾c chuyÓn dÞch vµ biÕn d¹ng c«ng tr×nh,. Tr­êng ®¹i häc Má-§Þa chÊt, Hµ Néi. [2]. NguyÔn Quang Phóc (2001), Nghiªn cøu ph­¬ng ph¸p ph©n tÝch ®é æn ®Þnh cña c¸c mèc chuÈn vµ xö lý sè liÖu ®o lón c«ng tr×nh, LuËn v¨n Th¹c sü kü thuËt, Tr­êng §¹i häc Má-§Þa chÊt, Hµ Néi. [3]. NguyÔn Quang Phóc (2001), Nghiªn cøu ph­¬ng ph¸p xö lý sè liÖu ®o lón c«ng tr×nh x©y dùng , TuyÓn tËp c¸c c«ng tr×nh khoa häc- TËp 32, Tr­êng ®¹i häc Má-§Þa chÊt, Hµ Néi. [4]. NguyÔn Quang Phóc (2001), Tiªu chuÈn æn ®Þnh cña c¸c ®iÓm ®é cao c¬ së trong ®o lón c«ng tr×nh , TuyÓn tËp c¸c c«ng tr×nh khoa häc- TËp 33, Tr­êng ®¹i häc Má-§Þa chÊt, Hµ Néi. [5]. NguyÔn Quang Phóc (2007), Quan tr¾c vµ ph©n tÝch biÕn d¹ng c«ng tr×nh, Bµi gi¶ng dïng cho häc viªn cao häc, Tr­êng ®¹i häc Má-§Þa chÊt, Hµ Néi. [6]. TrÇn Kh¸nh (1996), Nghiªn cøu øng dông b×nh sai tù do trong xö lý sè liÖu tr¾c ®Þa c«ng tr×nh, LuËn ¸n PTS khoa häc kü thuËt, Tr­êng ®¹i häc Má-§Þa chÊt, Hµ Néi. [7]. Tiªu chuÈn x©y dùng ViÖt Nam 203 - 1997 – Bé X©y Dùng. Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 60 - Líp Tr¾c ®Þa B-K48 Phô lôc 1 M· nguån modul Dù b¸o lón Private Sub mnuPhuthuoc_Click() Dim tencongtrinh As String tencongtrinh = txtTencongtrinh Dim Chukydo As String Chukydo = txtChukydo Dim Ngaydo As String Ngaydo = txtngaydo Dim Dg As Long, nDqt As Long, nh As Long Dg = VSGrid.TextMatrix(0, 0) nDqt = VSGrid.TextMatrix(0, 1) nh = VSGrid.TextMatrix(0, 2) Dim id() As Long Dim ic() As Long Dim h() As Double Dim nt() As Long Dim i, j As Long ReDim id(1 To nh), ic(1 To nh), h(1 To nh), nt(1 To nh) For i = 1 To nh id(i) = VSGrid.TextMatrix(i, 0) ic(i) = VSGrid.TextMatrix(i, 1) h(i) = VSGrid.TextMatrix(i, 2) nt(i) = VSGrid.TextMatrix(i, 3) Next Dim Tendiem() As String ReDim Tendiem(1 To Dg + nDqt) For i = 1 To (Dg + nDqt) Tendiem(i) = VSGrid.TextMatrix(nh + i, 1) Next Dim hgd() As Double ReDim hgd(1 To Dg + nDqt) Dim A As Double A = -99999999 For i = 1 To nDqt hgd(i) = A Next j = 0 For i = (nDqt + 1) To (nDqt + Dg) j = j + 1 hgd(i) = VSGrid.TextMatrix(nh + nDqt + j, 2) Next Dim dem As Long dem = 0 1: For i = 1 To nh If (hgd(id(i)) A) And (hgd(ic(i)) = A) Then hgd(ic(i)) = hgd(id(i)) + h(i) dem = dem + 1 End If If (hgd(id(i)) = A) And (hgd(ic(i)) A) Then hgd(id(i)) = hgd(ic(i)) - h(i) dem = dem + 1 End If Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 61 - Líp Tr¾c ®Þa B-K48 Next If dem < nDqt Then GoTo 1 End If Dim mtA() As Double ReDim mtA(1 To nh, 1 To nDqt + Dg) For i = 1 To nh For j = 1 To nDqt + Dg mtA(i, j) = 0 Next mtA(i, id(i)) = -1 mtA(i, ic(i)) = 1 Next Dim mtB() As Double ReDim mtB(1 To nh, 1 To nDqt) For i = 1 To nh For j = 1 To nDqt mtB(i, j) = mtA(i, j) Next Next Dim l() As Double ReDim l(1 To nh, 1 To 1) For i = 1 To nh l(i, 1) = (hgd(ic(i)) - hgd(id(i))) - h(i) Next Dim p() As Double ReDim p(1 To nh, 1 To nh) For i = 1 To nh p(i, i) = 1 / nt(i) Next Dim N() As Double ReDim N(1 To nDqt, 1 To nDqt) N = Nhan_MaTran(chuyenvi_matran(mtB(), nh, nDqt), _ Nhan_MaTran(p(), mtB(), nh, nh, nDqt), nDqt, nh, nDqt) Dim M() As Double ReDim M(1 To nh, 1 To 1) M = Nhan_MaTran(chuyenvi_matran(mtB(), nh, nDqt), _ Nhan_MaTran(p(), l(), nh, nh, 1), nDqt, nh, 1) Dim Nnd() As Double ReDim Nnd(1 To nDqt, 1 To nDqt) Nnd = nghichdao_matran(N(), nDqt) Dim x() As Double ReDim x(1 To nDqt, 1 To 1) x = DoiDau_MaTran(Nhan_MaTran(Nnd(), M(), nDqt, nDqt, 1), nDqt, 1) Dim Hbs() As Double ReDim Hbs(1 To nDqt, 1 To 1) Hbs = cong_matran(ho(), x(), nDqt, 1) Dim V() As Double ReDim V(1 To nh, 1 To 1) V = cong_matran(Nhan_MaTran(mtA(), x(), nh, nDqt, 1), l(), nh, 1) Dim hhbs() As Double ReDim hhbs(1 To nh, 1 To 1) hhbs = cong_matran(hh(), V(), nh, 1) Dim vpv() As Double Dim mo() As Double ReDim vpv(1 To 1, 1 To 1) Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 62 - Líp Tr¾c ®Þa B-K48 ReDim mo(1 To 1, 1 To 1) vpv = Nhan_MaTran(Nhan_MaTran(chuyenvi_matran(V(), nh, 1), p(), 1, nh, nh), V(), 1, nh, 1) For i = 1 To 1 For j = 1 To 1 mo(i, j) = Sqr(vpv(i, j) / (nh - nDqt)) Next Next Dim Rxx() As Double ReDim Rxx(1 To nDqt, 1 To 1) For i = 1 To nDqt Rxx(i, 1) = Sqr(Nnd(i, i)) Next Dim mx() As Double ReDim mx(1 To nDqt, 1 To 1) mx = Nhan_MaTran(Rxx(), mo(), nDqt, 1, 1) Dim fRf() As Double ReDim fRf(1 To nh, 1 To nh) fRf = Nhan_MaTran(Nhan_MaTran(mtA(), Nnd(), nh, nDqt, nDqt), _ chuyenvi_matran(mtA(), nh, nDqt), nh, nDqt, nh) Dim cfRf() As Double ReDim cfRf(1 To nh, 1 To 1) For i = 1 To nh cfRf(i, 1) = Sqr(fRf(i, i)) Next Dim mf() As Double ReDim mf(1 To nh, 1 To 1) mf = Nhan_MaTran(cfRf(), mo(), nh, 1, 1) With CommonDialog1 .DialogTitle = "In file nao" .CancelError = True .Filter = "Text file *.txt|*.txt|All Files | *.*" .ShowSave End With str = CommonDialog1.FileName Dim FileNum2 As Long FileNum2 = FreeFile Open str For Output As FileNum2 .............................. Print #FileNum2, String(5, " ") & "Chuong trinh duoc viet bang ngon ngu lap trinh VISUALBASIC 6.0" Close FileNum2 RichTextBox1.LoadFile CommonDialog1.FileName End Sub Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 63 - Líp Tr¾c ®Þa B-K48 Phô lôc 2 M· nguån modul TÝnh c¸c th«ng sè chuyÓn dÞch Private Sub mnuCacthongsochuyendich_Click() Dim tencongtrinh As String tencongtrinh = txtTencongtrinh Dim Chukydo As String Chukydo = txtChukydo Dim Ngaydo As String Ngaydo = txtngaydo Dim sodiemquantrac As Long Dim sochukyquantrac As Long Dim i, j As Long i = 1 Do While VSGrid.TextMatrix(1, i) "" i = i + 1 Loop sochukyquantrac = i - 1 i = 1 Do While VSGrid.TextMatrix(i, 1) "" i = i + 1 Loop sodiemquantrac = i - 1 Dim khoangthoigiantinhtuCkdau() As Double ReDim khoangthoigiantinhtuCkdau(1 To 1, 1 To sochukyquantrac) Dim khoangthoigiantinh2Ck() As Double ReDim khoangthoigiantinh2Ck(1 To 1, 1 To sochukyquantrac) For i = 1 To sochukyquantrac khoangthoigiantinhtuCkdau(1, i) = VSGrid.TextMatrix(0, i) Next For i = 1 To sochukyquantrac - 1 khoangthoigiantinh2Ck(1, i) = khoangthoigiantinhtuCkdau(1, i + 1) - khoangthoigiantinhtuCkdau(1, i) Next Dim Tendiem() As String ReDim Tendiem(1 To sodiemquantrac) For i = 1 To sodiemquantrac Tendiem(i) = VSGrid.TextMatrix(i, 0) Next Dim docao() As Double ReDim docao(1 To sodiemquantrac, 1 To sochukyquantrac) For i = 1 To sodiemquantrac For j = 1 To sochukyquantrac docao(i, j) = VSGrid.TextMatrix(i, j) Next Next Dim DoLun() As Double ReDim DoLun(1 To sodiemquantrac, 1 To 1) Dim Dlgiua2chukylientiep() As Double ReDim Dlgiua2chukylientiep(1 To sodiemquantrac, 1 To sochukyquantrac - 1) Dim Dlsovoickdautien() As Double ReDim Dlsovoickdautien(1 To sodiemquantrac, 1 To sochukyquantrac - 1) For i = 1 To sodiemquantrac For j = 1 To sochukyquantrac - 1 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 64 - Líp Tr¾c ®Þa B-K48 DoLun(i, 1) = docao(i, j + 1) - docao(i, j) Dlgiua2chukylientiep(i, j) = DoLun(i, 1) Dlsovoickdautien(i, j) = docao(i, j + 1) - docao(i, 1) Next Next Dim TongDocaotungck() As Double ReDim TongDocaotungck(1 To 1, 1 To sochukyquantrac) Dim TG() As Double ReDim TG(1 To sochukyquantrac, 1 To sodiemquantrac) TG() = chuyenvi_matran(docao(), sodiemquantrac, sochukyquantrac) Dim A, b, c As Double For i = 1 To sochukyquantrac b = 0 For j = 1 To sodiemquantrac c = b A = c + TG(i, j) b = A Next TongDocaotungck(1, i) = b NextDim TongDltheo2cklientiep() As Double ReDim TongDltheo2cklientiep(1 To 1, 1 To sochukyquantrac - 1) Dim TongDltheoSovoick0() As Double ReDim TongDltheoSovoick0(1 To 1, 1 To sochukyquantrac - 1) For i = 1 To sochukyquantrac - 1 TongDltheo2cklientiep(1, i) = TongDocaotungck(1, i + 1) - TongDocaotungck(1, i) TongDltheoSovoick0(1, i) = TongDocaotungck(1, i + 1) - TongDocaotungck(1, 1) Next Dim DlTrungbinhtheo2cklientiep() As Double ReDim DlTrungbinhtheo2cklientiep(1 To 1, 1 To sochukyquantrac - 1) Dim DlTrungbinhtheoSovoick0() As Double ReDim DlTrungbinhtheoSovoick0(1 To 1, 1 To sochukyquantrac - 1) For i = 1 To sochukyquantrac - 1 DlTrungbinhtheo2cklientiep(1, i) = TongDltheo2cklientiep(1, i) / sodiemquantrac DlTrungbinhtheoSovoick0(1, i) = TongDltheoSovoick0(1, i) / sodiemquantrac Next Dim Tocdoluntheo2Cklientiep() As Double ReDim Tocdoluntheo2Cklientiep(1 To 1, 1 To sochukyquantrac - 1) Dim TocdoluntheosovoiCk0() As Double ReDim TocdoluntheosovoiCk0(1 To 1, 1 To sochukyquantrac - 1) For i = 1 To sochukyquantrac - 1 Tocdoluntheo2Cklientiep(1, i) = DlTrungbinhtheo2cklientiep(1, i) / khoangthoigiantinh2Ck(1, i) TocdoluntheosovoiCk0(1, i) = DlTrungbinhtheoSovoick0(1, i) / khoangthoigiantinhtuCkdau(1, i + 1) Next Dim Dllechtheo2Cklientiep() As Double ReDim Dllechtheo2Cklientiep(1 To 1, 1 To sochukyquantrac - 1) Dim DllechtheoSovoiCk0() As Double ReDim DllechtheoSovoiCk0(1 To 1, 1 To sochukyquantrac - 1) Dim Smax1() As Double Dim Smax2() As Double Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 65 - Líp Tr¾c ®Þa B-K48 Dim Smin1() As Double Dim Smin2() As Double ReDim Smax1(1 To 1, 1 To sochukyquantrac - 1) ReDim Smax2(1 To 1, 1 To sochukyquantrac - 1) ReDim Smin1(1 To 1, 1 To sochukyquantrac - 1) ReDim Smin2(1 To 1, 1 To sochukyquantrac - 1) For i = 1 To sochukyquantrac - 1 Smax1(1, i) = CviDl2Cklt(i, 1) Smax2(1, i) = CviDlSvCk0(i, 1) For j = 1 To sodiemquantrac If CviDl2Cklt(i, j) >= Smax1(1, i) Then Smax1(1, i) = CviDl2Cklt(i, j) Else Smax1(1, i) = Smax1(1, i) End If If CviDlSvCk0(i, j) >= Smax2(1, i) Then Smax2(1, i) = CviDlSvCk0(i, j) Else Smax2(1, i) = Smax2(1, i) End If Next Next For i = 1 To sochukyquantrac - 1 Smin1(1, i) = CviDl2Cklt(i, 1) Smin2(1, i) = CviDlSvCk0(i, 1) For j = 1 To sodiemquantrac If CviDl2Cklt(i, j) <= Smin1(1, i) Then Smin1(1, i) = CviDl2Cklt(i, j) Else Smin1(1, i) = Smin1(1, i) End If If CviDlSvCk0(i, j) <= Smin2(1, i) Then Smin2(1, i) = CviDlSvCk0(i, j) Else Smin2(1, i) = Smin2(1, i) End If Next Next For i = 1 To sochukyquantrac - 1 Dllechtheo2Cklientiep(1, i) = Smax1(1, i) - Smin1(1, i) DllechtheoSovoiCk0(1, i) = Smax2(1, i) - Smin2(1, i) Next With CommonDialog1 .DialogTitle = "In file nao" .CancelError = True .Filter = "Text file *.txt|*.txt|All Files | *.*" .ShowSave End With str = CommonDialog1.FileName Dim FileNum3 As Long FileNum3 = FreeFile Open str For Output As FileNum3 ................................. Close FileNum3 RichTextBox1.LoadFile CommonDialog1.FileName End Sub Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 66 - Líp Tr¾c ®Þa B-K48 Phô lôc 3 M· nguån modul TÝnh Dù b¸o lón Private Sub mnuDubaoLun_Click() Dim Tendiem As String Tendiem = VSGrid.TextMatrix(0, 0) Dim nDl As Long nDl = VSGrid.TextMatrix(0, 1) Dim Si() As Double ReDim Si(1 To nDl) Dim i, j, k As Long Dim thutuDl() As Double ReDim thutuDl(1 To nDl) Dim t() As Double ReDim t(1 To nDl) Dim Ngay1 As Double Dim Thang1 As Double Dim Nam1 As Double Ngay1 = VSGrid.TextMatrix(0, 2) Thang1 = VSGrid.TextMatrix(0, 3) Nam1 = VSGrid.TextMatrix(0, 4) Dim Ngay() As Double ReDim Ngay(1 To nDl) Dim Thang() As Double ReDim Thang(1 To nDl) Dim Nam() As Double ReDim Nam(1 To nDl) For i = 1 To nDl Ngay(i) = VSGrid.TextMatrix(i, 2) Thang(i) = VSGrid.TextMatrix(i, 3) Nam(i) = VSGrid.TextMatrix(i, 4) Next For i = 1 To nDl thutuDl(i) = VSGrid.TextMatrix(i, 0) Si(i) = VSGrid.TextMatrix(i, 1) Next For i = 1 To nDl t(i) = (((Nam(i) - Nam1) * 12 + (Thang(i) - Thang1)) * 30 + (Ngay(i) - Ngay1)) / 30 Next With CommonDialog1 .DialogTitle = "In file nao" .CancelError = True .Filter = "Text file *.txt|*.txt|All Files | *.*" .ShowSave End With str = CommonDialog1.FileName Dim FileNum5 As Long FileNum5 = FreeFile Open str For Output As FileNum5 Print #FileNum5, String(17, " ") & "KET QUA XAY DUNG MO HINH DU BAO LUN THEO HAM DA THUC" Print #FileNum5, String(15, " ") & "*********************************************************" Print #FileNum5, Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 67 - Líp Tr¾c ®Þa B-K48 Print #FileNum5, Print #FileNum5, String(5, " ") & "1. Ket qua quan trac trong " & nDl & " chu ky :" Print #FileNum5, String(5, " ") & "-------------------------" For i = 1 To nDl Print #FileNum5, String(5, " ") & "- Chu ky " & Dinhdangso(i, 2, 0, True, False) & " : " _ & "Ngay " & Dinhdangso(Ngay(i), 3, 0, False, False) _ & " Thang " & Dinhdangso(Thang(i), 3, 0, False, False) & " Nam " & Dinhdangso(Nam(i), 5, 0, False, False) & " ; Gia tri lun : " & Dinhdangso(Si(i), 7, 2, False, True) & " mm" Next Dim thoigiancandubao As Double thoigiancandubao = 2 ................................. Dim bacdathuc As Double bacdathuc = 0 1: Dim mtA() As Double ReDim mtA(1 To nDl, 1 To bacdathuc + 1) For i = 1 To nDl For j = 1 To bacdathuc + 1 mtA(i, j) = t(i) ^ (j - 1) Next Next Dim mtB() As Double ReDim mtB(1 To 1, 1 To bacdathuc + 1) For i = 1 To 1 For j = 1 To bacdathuc + 1 mtB(1, j) = TgDb ^ (j - 1) Next Next Dim mtL() As Double ReDim mtL(1 To nDl, 1 To 1) For i = 1 To nDl mtL(i, 1) = -Si(i) Next Dim ATA() As Double ReDim ATA(1 To bacdathuc + 1, 1 To bacdathuc + 1) ATA() = Nhan_MaTran(chuyenvi_matran(mtA(), nDl, bacdathuc + 1), mtA(), bacdathuc + 1, nDl, bacdathuc + 1) Dim ATL() As Double ReDim ATL(1 To bacdathuc + 1, 1 To 1) ATL() = Nhan_MaTran(chuyenvi_matran(mtA(), nDl, bacdathuc + 1), mtL(), bacdathuc + 1, nDl, 1) Dim ATAnd() As Double ReDim ATAnd(1 To bacdathuc + 1, 1 To bacdathuc + 1) ATAnd() = nghichdao_matran(ATA(), bacdathuc + 1) Dim Xtg() As Double ReDim Xtg(1 To bacdathuc + 1, 1 To 1) Xtg() = Nhan_MaTran(ATAnd(), ATL(), bacdathuc + 1, bacdathuc + 1, 1) Dim x() As Double ReDim x(1 To bacdathuc + 1, 1 To 1) x() = DoiDau_MaTran(Xtg(), bacdathuc + 1, 1) Dim DlDb() As Double ReDim DlDb(1 To 1, 1 To 1) Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp L­¬ng Anh TuÊn - 68 - Líp Tr¾c ®Þa B-K48 DlDb() = Nhan_MaTran(mtB(), x(), 1, bacdathuc + 1, 1) Dim AX() As Double ReDim AX(1 To nDl, 1 To 1) AX() = Nhan_MaTran(mtA(), x(), nDl, bacdathuc + 1, 1) Dim V() As Double ReDim V(1 To nDl, 1 To 1) V() = cong_matran(AX(), mtL(), nDl, 1) Dim VTV() As Double ReDim VTV(1 To 1, 1 To 1) VTV() = Nhan_MaTran(chuyenvi_matran(V(), nDl, 1), V(), 1, nDl, 1) Dim Saisomohinh As Double Saisomohinh = Sqr(VTV(1, 1) / (nDl - bacdathuc - 1)) If Saisomohinh <= min Then min = Abs(Saisomohinh) bacdathucphuhop = bacdathuc c(bacdathucphuhop) = DlDb(1, 1) Else min = min End If Dim Stinh() As Double ReDim Stinh(1 To nDl, 1 To 1) Stinh() = Nhan_MaTran(mtA(), x(), nDl, bacdathuc + 1, 1) Dim SiTruStinh() As Double ReDim SiTruStinh(1 To nDl, 1 To 1) Dim Sido() As Double ReDim Sido(1 To nDl, 1 To 1) For i = 1 To nDl Sido(i, 1) = Si(i) Next For i = 1 To nDl SiTruStinh(i, 1) = Sido(i, 1) - Stinh(i, 1) Next Dim Ssmh() As Double ReDim Ssmh(1 To 1, 1 To 1) Ssmh() = Nhan_MaTran(chuyenvi_matran(SiTruStinh(), nDl, 1), SiTruStinh(), 1, nDl, 1) ........................ Close FileNum5 End Sub