Đề 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

®Þnh vÞ lµ qu¸ tr×nh tÝnh lÆp, ®­îc thùc hiÖn ®ång thêi víi qu¸ tr×nh b×nh sai. L­íi ®é cao ®o lón c«ng tr×nh thùc chÊt lµ mét m¹ng l­íi do lÆp ë nhiÒu thêi ®iÓm (mçi thêi ®iÓm lµ mét chu kú). ViÖc xö lý l­íi quan tr¾c t¹i mét thêi ®iÓm nµo ®ã thùc chÊt lµ x¸c ®Þnh ®é cao c¸c ®iÓm cña l­íi dùa vµo nh÷ng ®iÓm ®é cao gèc æn ®Þnh t¹i thêi ®iÓm ®ã. Nh­ vËy, khi xö lý hÖ thèng l­íi ®é cao ®o lón cÇn ph¶i gi¶i quyÕt 2 nhiÖm vô ®ã lµ: - Thø nhÊt, 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. - Thø hai, b×nh sai l­íi quan tr¾c, x¸c ®Þnh ®é cao cña c¸c ®iÓm ®o lón g¾n trªn c«ng tr×nh trong hÖ thèng c¸c ®iÓm gèc æn ®Þnh. Nghiªn cøu m« h×nh cña bµi to¸n b×nh sai tù do víi ­u ®iÓm tËp hîp sè liÖu gèc chØ tham gia vµo qu¸ tr×nh ®Þnh vÞ, mµ kh«ng tham gia vµo qu¸ tr×nh b×nh sai, v× vËy vector trÞ b×nh sai nhËn ®­îc kh«ng chÞu ¶nh h­ëng cña sai sè sè liÖu gèc vµ nh÷ng chuyÓn dÞch (nÕu cã) cña c¸c sè liÖu gèc. Vµ theo nh÷ng yªu cÇu nãi trªn cña viÖc ®Þnh vÞ hÖ thèng l­íi ®é cao ®o lón, chóng t«i nhËn thÊy r»ng viÖc sö dông thuËt to¸n b×nh sai tù do ®Ó xö lý hÖ thèng løíi ®é cao lµ mét gi¶i ph¸p toµn dÞªn, phï hîp víi ®Æc ®iÓm vµ b¶n chÊt cña d¹ng l­íi nµy. Sù phï hîp ®ã kh«ng chØ dõng l¹i ë ý nghÜa ®iÒu kiÖn ®Þnh vÞ . + = 0, mµ cßn cho phÐp xö lý ®ång thêi 2 cÊp l­íi, to¹ ra kh¶ n¨ng tù ®éng ho¸ cao khi sö dông m¸y tÝnh ®iÖn tö. Do ®ã, khi nghiªn cøu ph­¬ng ph¸p xö lý sè liÖu ®o lón c«ng tr×nh chóng t«i nhËn thÊy cÇn gi¶i quyÕt hai vÊn ®Ò sau: X©y dùng mét quy tr×nh hîp lý ®Ó cã thÓ gi¶i quyÕt ®­îc ®ång thêi 2 nhiÖm vô kÓ trªn khi xö lý c¸c m¹ng l­íi ®é cao ®o lón c«ng tr×nh. X©y dùng mét ch­¬ng tr×nh phÇn mÒm theo thuËt to¸n b×nh sai tù do phï hîp víi quy tr×nh ®· nªu ®Ó tù ®éng ho¸ xö lý sè liÖu ®o lón. VÊn ®Ò nµy sÏ ®­îc chóng t«i gi¶i quyÕt trong ch­¬ng sau. Ch­¬ng 3 xö lý sè liÖu quan tr¾c lón c«ng tr×nh theo ph­¬ng ph¸p b×nh sai l­íi tr¾c ®Þa tù do Trong ch­¬ng 2 khi xem xÐt m« h×nh cña bµi to¸n b×nh sai tù do vµ nh÷ng yªu cÇu cña viÖc ®Þnh vÞ hÖ thèng l­íi ®é cao ®o lón, cã thÓ kh¼ng ®Þnh r»ng, 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. Trong ch­¬ng nµy, chóng t«i x©y dùng mét quy tr×nh hîp lý tõ ®ã x©y dùng mét ch­¬ng tr×nh phÇn mÒm m¸y tÝnh theo m« h×nh bµi to¸n b×nh sai l­íi tù do ®Ó cã thÓ tù ®éng ho¸ qu¸ tr×nh ph©n tÝch x¸c ®Þnh hÖ thèng c¸c ®iÓm l­íi c¬ së æn ®Þnh ®ång thêi b×nh sai hÖ thèng l­íi ®é cao ®o lón c«ng tr×nh. C¸c thuËt to¸n ®­îc x©y dùng theo néi dung cña ph­¬ng ph¸p b×nh sai gi¸n tiÕp kÌm ®iÒu kiÖn. Chóng t«i sÏ sö dông thuËt to¸n b×nh sai tù do ®Ó b×nh sai hçn hîp hai bËc l­íi. 3.1 thuËt to¸n 1. Chän Èn sè Èn sè ®­îc chän lµ ®é cao b×nh sai cña tÊt c¶ c¸c ®iÓm trong l­íi, bao gåm c¸c ®iÓm cña l­íi c¬ së vµ c¸c ®iÓm cña l­íi quan tr¾c. 2. LËp ph­¬ng tr×nh sè hiÖu chØnh cña c¸c trÞ ®o Ph­¬ng tr×nh sè hiÖu chØnh cña trÞ ®o thø m trong l­íi gi÷a hai ®iÓm i vµ k (h×nh 3.1) ®­îc biÔu diÔn d­íi d¹ng tuyÕn tÝnh: , (3.1) k i hm H×nh 3.1 Chªnh cao gi÷a hai ®iÓm Trong ®ã, nh­ ®· biÕt: lµ sè h¹ng tù do. vµ lµ sè hiÖu chØnh vµo ®é cao gÇn ®óng cña ®iÓm k vµ i. lµ träng sè cña chªnh cao t­¬ng øng. §iÒu ®¸ng l­u ý ë ®©y lµ ®é cao gÇn ®óng vµ ø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è 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; 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: 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: (3.2) Trong ®ã: 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: (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 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) Trong ®ã, B lµ ma trËn hÖ sè cña phÐp chuyÓn ®æi to¹ ®é Helmert, trong ®ã: §èi víi l­íi ®é cao tù do, 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): TÝnh ma trËn gi¶ nghÞch ®¶o: (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: (3.7) Dùa vµo vector nghiÖm 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. 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: (3.8) - Vector ®é cao b×nh sai: (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Þ: (3.10) - Sai sè trung ph­¬ng cña c¸c Èn sè: (3.11) - Sai sè trung ph­¬ng cña hµm c¸c Èn sè: (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Þ 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 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 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): , HÖ ph­¬ng tr×nh sè hiÖu chØnh viÕt d­íi d¹ng ma trËn: 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. 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 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: . + = 0 Víi 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 . 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: = 1 øng víi c¸c ®iÓm ®é cao c¬ së. = 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: 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. 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) 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 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: 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. 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 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 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) 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 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 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 |" 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" 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 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 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 (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 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 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ë H×nh 3.11 S¬ ®å bè trÝ mèc ®o lón 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 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 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 | ============================================================== 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 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 ! 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. 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 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) 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 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 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 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 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, 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) 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