Giáo trình Lý thuyết sai số Trường ĐH Công nghiệp Quảng Ninh

BỘ CÔNG THƯƠNG TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP QUẢNG NINH - -Chủ biên: TS Bùi Ngọc Hùng Tham gia: ThS Nguyễn Thị Mai Anh GIÁO TRÌNH LÝ THUYẾT SAI SỐ (LƯU HÀNH NỘI BỘ) Quảng Ninh – 2019 Chương Lý thuyết sai số bình sai trực tiếp dãy trị đo đại lượng 1.1 Khái niệm đo đạc Phép thử nhằm xác định đại lượng trắc địa gọi phép đo Phép đo so sánh đại lượng cần xác định với đại lượng loại chọn làm đơn vị Các điều kiện tạo nên phép đo gọi điều kiện đo Kết thu phép đo theo quy trình định, thời điểm cụ thể gọi trị đo Các điều kiện phép đo bao gồm điều kiện đặc trưng cho môi trường ngoại cảnh, điều kiện đặc trưng cho độ tin cậy máy móc, dụng cụ đo điều kiện đặc trưng cho người đo Ví dụ đo dài trực tiếp thước thép Quy trình đo lúc tiến hành từ thao tác định tuyến, phân đoạn, kéo thước, đọc số đo đại lượng cần thiết để hiệu chỉnh vào đoạn nhiệt độ, hiệu chênh cao hai đầu thước,…Các điều kiện đặc trưng cho môi trường đo kể đến ảnh hưởng nhiệt độ, độ lồi lõm địa hình, tốc độ gió; Đặc trưng cho dụng cụ máy móc độ xác đọc số thước thép, máy móc dùng để đo nhiệt, đo hiệu chênh cao; Đặc trưng cho người đo cách tổ chức nhóm đo, lực kéo hai đầu thước khả đọc số Dễ nhận thấy, thao thác quy trình đo ảnh hưởng đến kết phép đo Trị đo thu cuối tổng hợp thông tin thu từ điều kiện đặc trưng điều kiện đo khoảng thời gian định Trên thực tế phép đo tiến hành phương pháp trực tiếp hay gián tiếp 1.1.1 Đo trực tiếp đo gián tiếp Đo trực tiếp: Là so sánh trực tiếp đại lượng cần đo với đơn vị đo tương ứng Ví dụ: Đo độ dài đoạn thẳng thước thép Đo gián tiếp: Là đại lượng cần xác định tính tốn thơng qua đại lượng đo trực tiếp Ví dụ: Muốn xác định diện tích hình chữ nhật ta cần phải đo tực tiếp chiều dài a chiều rộng b hay trị đo gián tiếp xác định thông qua hai đại lượng đo trực tiếp a b 1.1.2 Đo độ xác khơng độ xác Đo độ xác: Là kết đo nhận điều kiện đo Đo khơng độ xác: Là kết đo nhận điều kiện đo khác Ví dụ, kết đo góc β1, β2, β3 tam giác người đo dùng máy đo theo phương pháp với số lần đo điều kiện thời tiết ổn định kết đo góc nhận độ xác, ngược lại, điều kiện khác đi, chẳng hạn dùng máy khác, người đo khác áp dụng cách đo khác kết đo góc khơng độ xác 1.1.3 Trị đo cần thiết trị đo thừa Trị đo cần thiết (t): Là số lượng đo tối thiểu để xác định đồ hình lưới trắc địa giải toán trắc địa Trị đo thừa (r): Là đại lượng đo dư để kiểm tra nâng cao độ xác Trong bình sai điều kiện trị đo thừa số phương trình điều kiện ràng buộc trị đo với nhau, số trị đo thừa xác định là: r=n-t 1.2 SAI SỐ VÀ PHÂN LOẠI SAI SỐ 1.2.1 Sai số đo Bất kỳ phép đo dù hoàn chỉnh đến đâu tồn sai số, tức xác định trị thực đại lượng cần đo Nếu gọi Li trị đo, L0 trị thực chênh lệch trị thực trị đo gọi sai số thực trị đo, kí hiệu i i = L0 - Li (2.2) Sai số thực ln tồn ta phải tìm quy luật để xác định nó, thực tế ta khơng thể tìm trị thực mà xác định trị xác suất thông qua số hiệu chỉnh Vi: L’ i = L i + V i (2.3) Trong đó: L’i gọi trị xác suất hay trị sau bình sai 1.2.2 Ngun nhân gây nên sai số Có nhiều nguyên nhân gây sai số đo ta phân thành ngun nhân sau: + Sai số dụng cụ đo máy móc: Là sai số sử dụng dụng cụ máy móc khơng xác + Sai số người đo: Do giác quan người không chuẩn xác gây nên Ví dụ, mắt người phân biệt hai điểm cách 0.2mm nên đo chiều dài đoạn thẳng ta ước đọc không phần lẻ mm có sai số đọc số + Sai số ngoại cảnh: Do ảnh hưởng điều kiện thời tiết, điều kiện địa hình, địa vật gây nên Ví dụ, đo chiều dài đoạn thẳng thước thép đo nhiệt độ môi trường thay đổi làm chiều dài thước bị co dãn địa hình gập ghềnh thước bị cong, vênh,… kết đo xác 1.2.3 Phân loại sai số Dựa vào tính chất quy luật xuất xử lý số liệu đo đạc người ta phân làm loại sai số: 1.2.3.1 Sai số thô (sai lầm) Là sai số nhầm lẫn người đo đạc tính tốn, sai số thơ thường có giá trị lớn khơng xuất theo quy luật phát loại bỏ nhờ phép đo thừa 2.1.3.2 Sai số hệ thống Là sai số thường có trị số dấu khơng đổi biến đổi theo quy luật Các nguyên nhân sinh sai số hệ thống dụng cụ máy móc khơng hồn chỉnh, thói quen người đo điều kiện ngoại cảnh,… Ví dụ: Giả sử dùng thước thép 20m để đo đường thẳng đó, chiều dài thật thước 20,001m Như kết lần kéo thước có chứa +1mm sai số gọi sai số hệ thống Sai số hệ thống loại bỏ giảm bớt biết nguyên nhân quy luật xuất dùng phương pháp kiểm định tìm trị số để cải vào kết đo 2.1.3.3 Sai số ngẫu nhiên Sai số ngẫu nhiên sinh tổng hợp nhiều nguồn sai số, chúng tồn kết đo, xuất biến thiên phức tạp dấu trị số Giả sử đo đoạn thẳng thước thép sai số khắc vạch sai số giãn nở nhiệt mang tính chất hệ thống, cịn sai số lực kéo khơng đều, gió thổi,…tất sai số đồng thời tác động khoảnh khắc theo độ lớn (trị số) chiều (dấu) khác lên chúng sai số ngẫu nhiên Sai số ngẫu nhiên xuất không theo quy luật ta khơng khắc phục mà tìm cách hạn chế ảnh hưởng tới kết đo Vì sai số thơ sai số hệ thống tìm biện pháp loại trừ giảm bớt nên coi sai số ngẫu nhiên thành phần chủ yếu sai số đối tượng nghiên cứu lý thuyết sai số 1.3 Đặc tính sai số ngẫu nhiên 1.3.1 Ví dụ khảo nghiệm Sai số ngẫu nhiên có trị số dấu xuất khơng theo quy luật, điều kiện đo định, sai số ngẫu nhiên xuất theo quy luật Để nghiên cứu đặc tính sai số ngẫu nhiên, người ta tiến hành thực nghiệm đo góc 162 tam giác, xác định sai số thực sau: bảng 1-1 khoảng sai số -∆ +∆ tổng số lượng % số lượng % số lượng % 0'' 00 -:-  0''20 21 13,0 21 13,0 42 26,0 0'' 20 -:-  0''40 19 11,7 19 11,7 38 23,4 0'' 40 -:-  0''60 11 6,8 16 9,9 27 16,7 0'' 60 -:-  0''80 5,0 13 8,0 21 13,0 0'' 80 -:-  1''00 4,3 5,6 16 9,9 1'' 00 -:-  1''20 3,7 3,1 11 6,8 1'' 20 -:-  1''40 1,8 0,6 2,4 1'' 40 -:-  1''60 1,2 0,6 1,8  1''60 0 0 0 tổng 77 47,5 85 52,5 162 100,0 Để dễ nhận biết quy luật phân bố sai số, ta dùng số liệu bảng1-1 để vẽ đồ thị: Trên trục hoành biểu diễn khoảng sai số tương ứng với bảng diện tích hình chữ nhật biểu diễn số phần trăm xuất khoảng sai số ghi bảng 1-1 ta xác định chiều cao để vẽ hình chữ nhật sau: ví dụ sai số xuất khoảng 0'' 00 -:-  0''20 chiếm 13% chiều cao hình chữ nhật tương ứng là: h= 13 = 0.65 100x 0.2 Nếu nối điểm cạnh đáy phía hình chữ nhật ta đường gãy khúc đối xứng qua trục tung tiệm cận với trục hồnh 06 05 04 03 hình 1-1 đồ thị tính chất sai số ngẫu nhiên 0+02+04+06+08+1+12+14+16 01 -16-14-12-10 -08 -06-04 -02 1.3.2 Các đặc tính sai số ngẫu nhiên Qua thực nghiệm ta thấy sai số ngẫu nhiên có đặc tính sau đây: a Đặc tính giới hạn: Trong điều kiện đo đạc cụ thể, trị tuyệt đối sai số ngẫu nhiên không vượt giới hạn định b Đặc tính tập trung: Sai số ngẫu nhiên có trị tuyệt đối nhỏ, có khả xuất nhiều c.Đặc tính đối xứng: Sai số ngẫu nhiên dương âm với trị số tuyệt đối bé có số lần xuất gần d Đặc tính bù trừ: Khi số lần đo tiến tới vô cùng, số trung bình cộng sai số đo đạc ngẫu nhiên đại lượng tiến tới không lim n→ 1 +  + +  n  = = lim n→ n n 1.4 Tiêu chuẩn đánh giá độ xác Một nhiệm vụ chủ yếu công tác xử lý số liệu trắc địa đánh giá độ tin cậy đại lượng đo Đánh giá độ tin cậy xây dựng đại lượng đặc trưng cho kết đo kết xử lý, làm tiêu chuẩn để so sánh gọi tiêu chuẩn đánh giá độ xác Trong thực tế, sau đo nhiều lần ta thu trị trung bình cộng trường hợp với số lần đo khác khác Nhưng trị trung bình cộng khác trị thực Nảy sinh vấn đề trị trung bình cộng trị trung bình coi đáng tin cậy Giải mâu thuẫn người ta dựa vào tiêu chuẩn sau để đánh giá 1.4.1 Sai số trung bình Khi đo đại lượng n lần, giả sử sai số thực n lần đo : 1 , 2 3 , 4 n Vì tổng đại số sai số đo theo tính chất thứ tự sai số ngẫu nhiên triệt tiêu nên khơng dùng để đánh giá độ xác kết đo được.Vì ta dùng số trung bình để đánh giá độ xác kết đo Sai số trung bình trị trung bình cộng trị tuyệt đối sai số thực thành phần, xác định công thức: =   (2.5) n Trong đó:  sai số thực n số lần đo Tuy nhiên dùng sai số trung bình θ nhiều trường hợp chưa đánh giá xác chưa phản ánh độ biến động sai số ngẫu nhiên Ví dụ: Đo đại lượng điều kiện khác (khơng độ xác) sau: Điều kiện 1: : -5”, -3”, +7”, +1” Điều kiện 2: : +5”, -4”, -3”, -4” Sai số trung bình: Điều kiện I: Điều kiện II: + + +1 = 4" 5+ 4+3+ = 4"  II = I = Ta thấy I = II có nghĩa hai điều kiện đo khác cho độ xác Với kết ta lầm tưởng hai kết qủa đo hai điều kiện đo đạt độ xác nhau, Nhưng điều kiện II tốt biến động sai số điều kiện II < điều kiện I (hình 2.1) Khoảng biến động điều kiện II -5 -4 Khoảng biến động điều kiện I +5 1.4.2 Sai số trung phương (m) +7 Để phản ánh mức độ dao động sai số người ta dùng sai số trung phương để đánh giá độ xác đại lượng đo Khi có dãy sai số thực, n vơ lớn sai số trung phương tính : Cơng thức gần Gauss: m=  =  n 21 + 22 + + 21n n (2-6) Trong đó:  sai số thực: i = L0 - Li (2.7) n số lần đo Vậy sai số trung phương (m): Là bậc hai trị trung bình cộng bình phương sai số thực thành phần Ví dụ: Với số liệu mục (2.2.1), thay vào (2-6): mI = mII = 25 + + 49 + = 4."58 25 + 16 + + 16 = 4."06 Từ kết ta thấy rõ ràng điều kiện II tốt điều kiện I 1.4.3 Sai số xác suất (r) Sai số xác suất trị số sai số ngẫu nhiên mà sai số ngẫu nhiên khác có trị tuyệt đối lớn nhỏ có khả xuất Nếu xếp sai số ngẫu nhiên theo thứ tự tăng dần coi sai số nằm sai số xác suất Bằng lý thuyết xác suất người ta chứng minh mối quan hệ sai số xác suất sai số trung phương m: r = 0.6745m  m  = 0.7979m  m (2-8) (2-9) 1.4.4 Sai số trung phương tương đối Sai số trung bình, sai số trung phương, sai số giới hạn sai số tuyệt đối Trong đo dài dùng sai số tương đối phản ánh rõ mức độ xác kết đo Sai số tương đối tỷ số sai số đo giá trị đại lượng đo Trong tử số ln nhận mẫu số làm tròn đến bội số 10 Mẫu số sai số tương đối biểu thị cho chất lượng đo đạc, mẫu số lớn độ xác đo cao ngược lại mL = T L (2-10) Ví dụ: Đo chiều dài đoạn thẳng thước thép kết sau: m S = 5mm S1 = 50m; mS = 5mm S2 = 250m; Nừu nhìn vào trị số tuyệt đối mS ta tưởng hai cạnh đo với độ xác nhau, dùng sai số tương đối ta thấy cạnh S2 có độ xác cao Sai số trung phương tương đối: mS1 S1 mS2 S 21 = = ; 50 000 10 000 = = 250 000 50 000 1.4.5 Sai số giới hạn Theo đặc tính thứ sai số ngẫu nhiên nêu rõ: điều kiện đo định, trị số tuyết đối sai số ngẫu nhiên không vượt giới hạn định Như vậy, dãy kết đo, trị đo có sai số vượt q giới hạn trị đo không phù hợp với quy luật xem khơng đảm bảo độ xác khơng dùng để cải vào kết đo Giá trị sai số giới hạn phụ thuộc chủ yếu vào điều kiện đo.Trong lý thuyết xác suất người ta chứng minh sai số đo phù hợp với quy luật phân phối chuẩn có 3% sai số ngẫu nhiên có giá trị lớn lần sai số trung phương, 5% sai số ngẫu nhiên có giá trị vượt lần sai số trung phương Trong thực tế số lần đo không nhiều, sai số ngẫu nhiên lớn lần sai số trung phương có khả xuất Vì ta thường lấy lần sai số trung phương làm giới hạn sai số ngẫu nhiên Gọi max sai số giới hạn, ta có: max = 3m (2.11) Trong trắc địa cơng trình với u cầu độ xác cao, thường quy định: (2.12) max = 2m 1.5 Sai số trung phương hàm trị đo Trong toán trắc địa đại lượng cần xác định thường hàm số mà đối số đại lượng đo trực tiếp Vì ta cần phải xác định sai số trung phương hàm đối số trị đo trực tiếp độc lập Ví dụ: Để xác định gia số toạ độ X, Y ta phải xác định thơng qua trị đo cạnh S góc định hướng , góc  lại xác định thơng qua góc  1.5.1 Sai số trung phương hàm số dạng tổng quát Giả sử ta có hàm số dạng tổng quát sau: F = f(x,y,…,u) (2.13) Trong đó: x, y, …,u đại lượng đo độc lập, có sai số trung phương tương ứng m x, my, , mu Ta cần xác định sai số trung phương hàm trị đo mf Gọi X, Y, ,U trị thực đại lượng đo độc lập, sai số thực tương ứng là: X = xi + xi →xi = X - xi Y = yi + yi →yi = Y - yi ………………………… U = ui + ui →ui = U - ui Sai số thực hàm: Fi = F - fi Hay: Fi = F ( X, Y,…, U) – fi (x, y,…,u) Hoặc: Fi = F (xi + xi, yi + yi,…, ui + ui) – f( xi, yi,…, ui) (2.15) Khai triển Taylor bỏ qua số hạng bậc cao cách lấy đạo hàm riêng hàm số F theo đối số x, y,…, u, ta có:  F   F   F  y i + +  Fi =  xi +  u i  x   u   y  (2.16) (với i = 1,2, ,n) Bình phương vế (2-16) sau lấy tổng vế chia cho số lần đo n: F F  =  F  x.x +  F  y.y  + + 2 F . F  x.y  + +  F  u.u  n    x   y    n     x   y  n n    u  n Theo giả thiết đại lượng x, y,…, u đại lượng đo độc lập sai số x, y,…, u sai số ngẫu nhiên độc lập Theo tính chất sai số ngẫu nhiên định nghĩa sai số trung phương ta có:  F   F   F   m y + +  =  m u  m x +   x   u   y  m F (2.17) Hay:  F   F   F   m y + +  mF =   m x +   m u  x   u   y  2 (2.18) Ví dụ: Trong tam giác đo góc: 1 = 500 00’ 00”; 2 = 600 00’ 00” Với sai số m1 = 3”; m2 = 4” Hãy tính: m3 3? Giải: Ta có: f = 3 =1800 - (1 + 2) = 700 00’ 00” 𝜕𝑓 2 Vì ( ) = 𝑛ê𝑛 𝑚𝛽 = √𝑚𝛽1 + 𝑚𝛽2 = √16 + = ±5′′ 𝜕𝛽 1.5.2 Sai số trung phương số hàm đơn giản * Sai số trung phương hàm dạng: Giả sử có hàm số dạng: F = a1L1 + a2L2 + +anLn (2.19) Trong đó: số Li đại lượng đo với sai số trung phương mi Vi phân hàm F theo đối số đại lượng đo chuyển vi phân thành sai số thực ta được: (2.20) F = a1L1 + a2L2 + +anLn Và chuyển sang sai số trung phương: m F = a 21m L + a 2 m L + + a n m L n Hay: m F = a 21 m L1 + a 2 m L2 + + a n m Ln (2.21) Nếu a1 = a2 = = an =1: → m F = m L + m L + + m L n Nếu đo độ xác: m L = m L = = m L = m L ta có: n mF = mL n (2.22) * Sai số trung phương hàm số dạng: Giả sử hàm số có dạng: F= L1 L2 L3 Áp dụng cơng thức tính sai số trung phương hàm tổng quát: m F L =   L3 2 LL L    m L1 +   m L2 +   L3  L3     m L3  (2.23) Khi toán vừa hàm có trị đo góc vừa có trị đo chiều dài chênh cao sử dụng công thức hàm tổng quát ta phải biểu thị sai số trung phương đo góc dạng mβ/ρ’’ Ví dụ: Xác định chênh cao h phương pháp đo cao lượng giác: h = S tgV + i - t Với trị đo sai số trung phương là: S = 143.5m; V =2030 ; i = 1.5m; l = 2.2m, mS = ±1’; mi =ml= (coi sai số), ta có hàm: h = 143.5tg2030 +1.5 - 2.2 = 5.76m 𝑚ℎ = √𝑡𝑔2 𝑉 𝑚𝑆2 + 𝑠 𝑚𝑣2 = ±0.048𝑚 𝑐𝑜𝑠 𝑉 𝜌2 Ví dụ: Một cạch đường chuyền S đo thành đoạn liền kết sau: L1= 91,24m  0,20m ; L2= 79,15m  0,15m L3= 102,05m  0,35m ; L4= 117,35m  0,18m Hãy tính sai số trung phương sai số trung phương tương đối cạnh đường chuyền - Chiều dài cạnh đường chuyền : S = L1+L2+L3+L4 = 390.39m - Sai số trung phương : m S =  m L  + m L  + m L  + m L  = ,m Vậy  mS ,  = =  T S ,  1.5 Công thức Fererô 10 + Tính β0 - Phương trình số hiệu chỉnh cạnh đo Trong hình vẽ Si chiều dài cạnh đo hai điểm cần xác định j K Phương trình trị bình sai cạnh đo là: S’ = S + vS = S0 + dS k Vậy vS = dS + S0 – S = dS +lS Phương trình số hiệu chỉnh cạnh đo là: v S = a jk dx j + b jk dy j − a jk dxk − b jk dy k + l S Trong đó: a jk = − x jk S jk ; b jk = − αjk y jk J S jk 3.2.2.Các ví dụ lập hệ phương trình số hiệu chỉnh + Ví dụ 1: Cho lưới độ cao hình vẽ Kết đo hiệu chênh cao ghi bảng sau Các điểm A B có độ cao cho trước HA= 0.0 m; HB = 8,000m Hãy lập phương trình số hiệu chỉnh cho mạng lưới TT hi(m) 3.001 4.002 1.002 5.004 h2 P Q h4 h1 h3 B Giải: A - Chọn ẩn số: t= 4-2 = Ẩn số: HP’; HQ’ - Tính trị gần đúng: HP0= 3,0m; HQ0= 7,0m Phương trình số hiệu chỉnh: v1 = dHP + l1 = dHP -1 v2 = - dHP + dHQ +l2 = - dHP + dHQ -2 v3 = - dHQ +l3 = - dHQ -2 v4 = - dHP +l4 = - dHP -4 Trong đó: l1= HP0 – HA- h1= 3-0-3.001 = -1(mm) l2= HQ0 – HP0- h2= 7-3-4.002 = -2(mm) l3= HB – HQ0- h3= 8-7-1.002 = -2(mm) l4= HB – HP0- h4= 8-3-5.004 = -4(mm) 74 +Ví dụ 2: Viết phương trình số hiệu chỉnh cho mạng lưới sau: VI I II Cạnh đo I II = 663,389m, góc phương vị đo αI II = 1910 23’28’’ TT Góc đo TT Góc đo 117 50 57.8 38 57 44.58 25 13 6.03 22 35 00 17 03 53.67 19 41 23.25 98 45 42.83 19 52 19.5 Tên điểm X(mm) Y(mm) I 1000000 1000000 Giải: +t=6 + Ẩn số: XII’, YII’, XIII’, YIII’, XIV’, YIV’ + Tính trị gần đúng: - Tên cạnh S(mm) I II 663389 II III 288130,8357 I IV 226153,2359 I III 477867,0694 III IV 333766,241 Phương vị : αI III 50 57 10,75 αI V 129 50 43,42 75 III - Tọa độ gần tính theo tốn trắc địa thuận: tên điểm X Y II 349678,551 868977,1842 III 531188,837 1092748,037 IV 855099,507 1173635,057 - Tính hệ số hướng phương trình số hiệu chỉnh + Tọa độ gần Tên điểm X I 1000000 II III IV Y 1000000 349678,6 868977,1842 531188,8 1092748,037 855099,5 1173635,057 Tên hướng dx dy I II -650321,45 -131022,816 I III -468811,16 92748,03711 I IV -144900,49 173635,0574 II I 650321,449 131022,8158 II IV 505420,956 304657,8732 II III 181510,286 223770,853 III II -181510,29 -223770,853 III I 468811,163 -92748,0371 III IV 323910,67 80887,02025 IV III -323910,67 -80887,0203 IV II -505420,96 -304657,873 IV I 144900,493 -173635,057 76 + Hệ số phương trình số hiệu chỉnh Tên góc Tên hướng dx dy s a b III I 468811,163 -92748,037 477897,589 -0,084 -0,423 III II -181510,286 -223770,853 288130,836 -0,556 0,451 III IV 323910,670 80887,020 333857,503 0,150 -0,599 III I 468811,163 -92748,037 477897,589 -0,084 -0,423 IV II -505420,956 -304657,873 590141,307 -0,180 0,299 IV III -323910,670 -80887,020 333857,503 -0,150 0,599 IV I 144900,493 -173635,057 226153,236 -0,700 -0,584 IV II -505420,956 -304657,873 590141,307 -0,180 0,299 I III -468811,163 92748,037 477897,589 0,084 0,423 I IV -144900,493 173635,057 226153,236 0,700 0,584 I II -650321,449 -131022,816 663389,000 -0,061 0,305 I III -468811,163 92748,037 477897,589 0,084 0,423 II IV 505420,956 304657,873 590141,307 0,180 -0,299 II I 650321,449 131022,816 663389,000 0,061 -0,305 II III 181510,286 223770,853 288130,836 0,556 -0,451 II IV 505420,956 304657,873 590141,307 0,180 -0,299 cạnh gốc I II -650321,449 -131022,816 663389,000 0,980 0,198 Pvij gốc -650321,449 -131022,816 663389,000 -0,061 0,305 I II + Giá trị góc gần đúng: Tính theo phương pháp tính góc phương vị sau tính góc Ví dụ góc 10= α31- α32 Giá trị góc gần TT 117 51 22,73 25 12 42,67 17 34,56 98 45 52,72 38 57 50,06 22 34 54,52 19 41 22,70 19 52 20,05 77 + Phương trình số hiệu chỉnh: TT dx2 dy2 dx3 dy3 dx4 dy4 l -0,556 0,451 0,472 -0,874 0,000 0,000 24,93 0,000 0,000 0,233 -0,176 -0,150 0,599 -23,36 0,180 -0,299 -0,150 0,599 -0,031 -0,300 -19,11 -0,180 0,299 0,000 0,000 -0,520 -0,884 9,89 0,000 0,000 -0,084 -0,423 0,700 0,584 5,48 0,061 -0,305 0,084 0,423 0,000 0,000 -5,48 0,119 0,005 3,000 3,000 -0,180 0,299 -0,55 0,376 -0,152 -0,556 0,451 0,180 -0,299 0,55 pvi gốc 0,061 0,305 0,000 0,000 0,000 0,000 0,00 cạnh gốc -0,980 -0,198 0 0 0,00 3.3 Thành lập hệ phương trình chuẩn giải phương trình chuẩn 3.3.1 Thành lập hệ phương trình chuẩn Số lượng phương trình chuẩn bình sai gián tiếp số ẩn số t Muốn thành lập hệ phương trình chuẩn ta làm sau: TT dx1 dx2 … dxt li Si vi a1 b1 … t1 l1 S1 v1 a2 b2 … t2 l2 S2 v2 a3 b3 … t3 l3 S3 v3 … … … … … … … … n an bn … tn ln Sn [a] [b] … [t] [l] [S] dx1 dx2 … dxt [paa] [paa] … [pat] [pal] [pbb] … [pbt] … S’ f Sj [pas] S1’ f1 S1 [pbl] [pbs] S2’ f2 S2 … … … [ptt] [ptl] [pts] [pl] [p.s] … … St’ … ft St Cách xác định tương tự bình sai điều kiện ta lập phương trình chuẩn sau: 78  paadx1 +  pabdx2 + +  patdxt +  pal =  pabdx1 +  pbbdx2 + +  pbtdxt +  pbl =  patdx1 +  pbtdx2 + +  pttdxt +  ptl = 3.3.2 Giải phương trình chuẩn A 34 B 79 I II Kí hiệu dx1 a dx2 b dx3 c [.l] [.s] e Q1j e a [aa] [ab] [ac] [al] [as] E1 -1  aa  b E1bxa -  ab  =E1b  aa  - [bb] -  ac = E1c  aa  - [bc]  ab  [ab]  aa  - c E1cxa E2cx(b.1) c.2 E3 [bb.1] -1 −  ab  [ac]  aa  bc.1 = E bb.1 2c [cc] E1c.[ac] E2cx[bc.1] - [bl] - [bc.1] b.1 E2  al = E1l  aa   ab  [al]  aa  [bl.1] -  bl.1 =E  bb.1 2l [cl] E1c.[al] E2cx[bl.1]  as   aa  - [bs] -  ab  [as]  aa  - [bs.1] -  bs.1  bb.1 [cs] E1c.[as] E2cx[bs.1]  ab   aa  E1b - E1b  bb.1 E1c E2cxE1b e Q3j fj Sj 0 f1 S1 0 0 1  bb.1 0 0 E2c 0 Q2j - f1 = E1f  aa  − f2 -  ab  f2  aa  [f2.1] - S1 aa  S2 -  ab  S  aa  [ S 1] f 1 =E  bb.1 2f S2 1 -   bb.1 f3 E1cf1 E2cx[f2.1] S3 E1c S1 E2cx[ S 1] [cc.2] -1 [cl.2] - cl.2 cc.2 [cs.2] cs.2 cc.2 80 E2c E E1c + E2cxE1b − 2c − cc.2  cc.2 E1c+ E2cxE1b [f3.2] - cc.2 - f3 2 cc.2 [ S3 2] - [S3 2] cc.2 dx1 dx2 dx3 [ll] [sl] 81 3.4 Đánh giá độ xác bình sai gián tiếp Như ta biết để đánh giá độ xác ẩn số, phải xác định phần tử ma trận trọng số đảo ẩn số Qx Sử dụng thuật tốn Gauss, thơng thường áp dụng phương pháp cột phụ Gauss, phương pháp Ganzen hay phương pháp Enke để xác định trọng số đảo Nội dung phương pháp sau: Phương pháp cột phụ Gauss Cơ sở toán học phương pháp cột phụ Gauss xuất phát từ đẳng thức: N.N-1 = N.Qx = Qx.N = E * Hay N Qx - E = Trong phương trình N ma trận hệ số hệ phương trình chuẩn, E ma trận đơn vị kích thước txt với cột tương ứng là: Cột Cột Cột t 1    0 0   . 0   0   1 0   . 0   0   0 0   . 1   … Ta thấy hệ phương trình chuẩn phương trình * có ma trận hệ số N, khác ma trận ẩn số ma trận số hạng tự Vì cho phép kết hợp việc tìm phần tử ma trận trọng số đảo ẩn số sơ đồ giải hệ phương trình chuẩn Gauss Việc tiến hành gắn thêm vào phía phải sơ đồ giải hệ phương trình chuẩn mảng cột phụ ứng với ma trận E Ứng với cột ta ký hiệu cột e/Q ij Cụ thể cột có dạng: e e Q2 i Qti Q1i ; ;…; 0 0 1       1 0 0             0 t 0 e Dựa vào cột phụ này, sơ đồ Gauss kết hợp với việc giải hệ phương trình chuẩn , tìm trọng số đảo Qik ẩn số Phương pháp Ganzen Cơ sở phương pháp Ganzen quan hệ: ** N.Qx = E dựa sở hệ số có sẵn sơ đồ Gauss giải hệ phương trình chuẩn t ( t+1) t ( t+1) trọng số đảo cần tìm Trên sở hệ phương trình hệ số trọng số, ta chọn Thực phép nhân ** ta có hệ t2 phương trình hệ số trọng số ẩn số chứa phương trình giải hệ phương trình tìm trọng số đảo Để hiểu rõ nội dung phương pháp Ganzen , ta lấy ví dụ với hệ ẩn số Các hệ số trọng số càn tìm lúc Q11, Q12, Q13, Q22, Q23, Q33 Từ phương trình hệ số trọng số dạng: 82  paa  Q13 +  pab Q23 +  pac Q33 = 0   pbb.1 Q23 +  pbc.1 Q33 =0    pcc.2 Q33 =   paa  Q12 +  pab Q22 +  pac Q23 = 0   pbb.1 Q22 +  pbc.1 Q23 =1    pcc.2 Q23 =   paa  Q11 +  pab Q12 +  pac Q13 = 1   pbb.1 Q12 +  pbc.1 Q13 =0    pcc.2 Q13 =  Ta chọn phương trình để tìm trọng số đảo:  paa  Q13 +  pab Q23 +  pac Q33 = 0   pbb.1 Q23 +  pbc.1 Q33 =0  (a)   pcc.2 Q33 =   paa  Q12 +  pab Q22 +  pac Q23 = 0  (b)  pbb.1 Q22 +  pbc.1 Q23 =1    paa  Q11 +  pab Q12 +  pac Q13 = 1 (c) Từ nhóm phương trình a ta tìm được:  pbc.1 Q =E Q Q33 = ; Q23 = pcc.2  pbb.1 33 2c 33 Q13 =-  pab Q -  pac Q =E Q +E Q  paa  23  paa  33 1b 23 1c 33 Từ nhóm phương trình (b) ta có Q22 =  pbc.1 Q = +E Q  pbb.1  pbb.1 23  pbb.1 2c 23 Q12 =-  pab Q -  pac Q =E Q +E Q  paa  22  paa  23 1b 22 1c 23 Từ phương trình c ta có: Q11 =  pac Q = +E Q +E Q  pab Q12  paa  13  paa  1b 12 1c 13  paa   paa  Phương pháp Enke Phương pháp Enke dùng để tính trọng số đảo ẩn số cuối gần cuối hệ phương trình chuẩn Bởi ẩn số cần đánh giá độ xác phải xếp theo thứ tự ẩn số cuối gần cuối Chính đặc điểm phương pháp Enke, nên phương pháp tối ưu, lưới cần đánh giá độ xác ẩn số trở lại Chúng ta xây dựng nội dung phương pháp Enke cụ thể với trường hợp hệ phương trình chuẩn có ẩn số, ẩn số cuối dx3, cịn ẩn số gần cuối dx2 Để tìm trọng số đảo Q33 Q22 ta tiến hành sau: a Từ nội dung phương pháp Ganzen ta có: 83 Q33 = hay P33=[pcc.2]  pcc.2 b Để tìm Q22 ta sử dụng phương trình hệ số trọng số dạng:  paa  Q12 +  pab Q22 +  pac Q23 = 0   pab Q12 +  pbb Q22 +  pbc Q23 = 1  pac Q12 +  pbc Q22 +  pcc Q23 =  Sau khử lần theo Gauss, ta thu phương trình:  pbb.1 Q22 +  pbc.1 Q23 =1 (a)    pbc.1 Q22 +  pcc.1 Q23 =0 (b)  pbc.1 cộng vào đẳng thức b ta có: Nhân vế biểu thức a với  pbb.1  pbc.1 hay Q =-  pbc.1 Q  pcc.2 Q23 =23  pbb.1  pbb.1 33 Thay Q23 vào phương trình b, ta có: Q22 =-  pcc.1 Q =  pcc.1 Q  pbc.1 23  pbb.1 33 Vì [pcc.1] khơng có sơ đồ giải hệ phương trình chuẩn phương pháp Gauss, nên tính từ:  pbc.1  pcc.1 =  pcc.2 +  pbb.1 3.5.Bài tập ứng dụng Bình sai gián tiếp Cho lưới độ cao cho trước hình vẽ Hãy lập giải hệ phương trình chuẩn lưới theo phương pháp bình sai gián tiếp Cho số liệu đo sau: K TT Chênh cao h(m) Chiều dài S(m) 0.93 2.51 4.24 1.5 6.78 h1 S1 A S4 Tên điểm Độ cao H(m) A 0.000 B 7.720 84 S2 h4 J B Số liệu gốc: Với mạng lưới ta có ẩn số HK’ HJ’ h2 h3 S3 Trị gần ẩn số: HK0 = HA+h1=0.930m HJ0 = HB -h3=3.480m Ta lập phương trình số hiệu chỉnh sau: v1 = dH K v = - dH K + dH J +4cm v3 = -dH J v = -dH K +1cm Từ ta có ma trận hệ số sau: 1  -1 A=  0   -1 0 0    1 ; L=   0 -1     4x2   x1 Nếu coi trị đo độc lập chọn trọng số Phi theo công thức: Pi = viết ma trận trọng số trị đo: 1  0 P=  0  0  0 0 ta Si        2  -1   −2.5  T  ; M = A PL =    -1   2.0  Ta tìm ma trận N = A T PA =  Và N −1 = 12    ; X= 25  2   dH K  -1   = -N M =  dH J   0.92cm     -1.32cm  Từ tìm độ cao điểm sau bình sai: HK’ = HK0 + dHK = 0.930+0.0092 = 0.9392m HJ’ = HJ0 + dHJ = 3.480+0.0132 = 3.4668m + Giải theo sơ đồ Gauss: Từ hệ phương trình số hiệu chỉnh ta lập hệ phương trình chuẩn sau: 2dH K − 0.5dH J − 2.5 =   −0.5dH K + dH J + = Giải hệ phương trình chuẩn theo sơ đồ Gauss ta được: 85 Kí hiệu dx1 a a E11 dx2 [.l] b 2.000 -0.500 -2.500 -1.000 0.250 1.250 2.000 -0.625 1.375 -1.320 0.92 1.167 -0.125 1.042 -1.000 -1.320 b E1bxa b.1 E2 Vậy độ cao điểm sau bình sai là: HK’ = HK0 + dHK = 0.930+0.0092 = 0.9392m HJ’ = HJ0 + dHJ = 3.480+0.0132 = 3.4668m Bài tập: Cho mạng lưới tứ giác hình vẽ Biết điểm gốc A, B, C có tọa độ sau A(1000,0); B(0,0); C(0, 1000) Các góc đo với số liệu sau: STT Trị đo STT 45 00 00 90 00 01 Trị đo 45 00 00 90 00 00 45 00 00 45 00 01 A B Hãy lập phương trình số hiệu chỉnh mạng lưới Bài giải Bước + Chọn ẩn số: Là tọa độ điểm cần xác định X’I, Y’I + Tính tọa độ gần điểm I0 (1000; 1000) + Lập phương trình số hiệu chỉnh v1 = −a BI dxI − bBI dyI + l1 v2 = a AI dxI + bAI dyI + l v3 = (a IA − a IB )dxI + (bIA − bIB )dy I + l3 v4 = a BI dxI + bBI dyI + l v5 = −aCI dxI − bCI dy I + l5 v6 = (a IB − a IC )dxI + (bIB − bIC )dy I + l6 86 I C Tên điểm A B C I0 x 1000 0 1000 y 0 1000 1000 Bảng tính hệ số hướng: Góc Hướng BI BA AB AI IA IB BC BI CI CB IB IC ∆X 1000 1000 -1000 0 -1000 1000 1000 -1000 -1000 ∆Y 1000 0 1000 -1000 -1000 1000 1000 -1000 -1000 S 1414.214 1000 1000 1000 1000 1414.214 1000 1414.214 1000 1000 1414.214 1000 Trị gần góc: = αBI – αBA = 45-0 = 450 = αAB – αAI = 180 – 90 = 900 = αIA – αIB = 270 – 225 = 450 = αBC – αBI = 45-0 = 450 = αCI – αCB = 360 – 270 = 900 = αIB – αIC = 225 – 180 = 450 Hệ phương trình số hiệu chỉnh: v1 = −103.13dxI + 103.1dyI + v2 = 206.27dxI + 0dyI −1 v3 = −103 2dxI − 103 2dyI + v4 = 103.13dxI −103.13dyI + v5 = 0dxI + 206 27 dyI + v6 = −103 1dxI − 103 1dyI − Bước + Thành lập hệ phương trình chuẩn 87 α 45 180 90 270 225 90 45 270 225 180 a 103.133 0.000 0.000 206.265 -206.265 -103.133 206.265 103.133 0.000 -206.265 -103.133 0.000 b -103.133 -206.265 206.265 0.000 0.000 103.133 0.000 -103.133 -206.265 0.000 103.133 206.265 - Lập bảng hệ số phương trình số hiệu chỉnh TT dxI dyI li Si vi -103.1 103.1 0 v1 206.2 -1 205.2 v2 -103.1 -103.1 -206.2 v3 103.1 -103.1 0 … 206.2 206.2 -103.1 -103.1 -1 0 -2 -207.2 -2 + Hệ số phương trình chuẩn: dx1 dx2 [pl] [p.s] [paa] [paa] [pal] [pas] [pbb] [pbl] [pbs] dx1 dx2 [pl] [p.s] 85036.88 -103.1 84933.78 85036.88 103.1 85139.98 Vậy phương trình chuẩn là: 85036.88dxI + 0dyI − 103.1 = 0dxI + 85036.88dyI + 103.1 = 88 ... Sai số xác suất (r) Sai số xác suất trị số sai số ngẫu nhiên mà sai số ngẫu nhiên khác có trị tuyệt đối lớn nhỏ có khả xuất Nếu xếp sai số ngẫu nhiên theo thứ tự tăng dần coi sai số nằm sai số. .. tới kết đo Vì sai số thơ sai số hệ thống tìm biện pháp loại trừ giảm bớt nên coi sai số ngẫu nhiên thành phần chủ yếu sai số đối tượng nghiên cứu lý thuyết sai số 1.3 Đặc tính sai số ngẫu nhiên... 1.2.3 Phân loại sai số Dựa vào tính chất quy luật xuất xử lý số liệu đo đạc người ta phân làm loại sai số: 1.2.3.1 Sai số thô (sai lầm) Là sai số nhầm lẫn người đo đạc tính tốn, sai số thơ thường