Từ khoá: Số Hiệu, Tiêu đề hoặc Nội dung ngắn gọn của Văn Bản...

Đăng nhập

Đang tải văn bản...

Tiêu chuẩn quốc gia TCVN 12741:2019 về Xác định và sử dụng hàm hiệu chuẩn đường thẳng

Số hiệu: TCVN12741:2019 Loại văn bản: Tiêu chuẩn Việt Nam
Nơi ban hành: *** Người ký: ***
Ngày ban hành: Năm 2019 Ngày hiệu lực:
ICS:03.120.30 Tình trạng: Đã biết

1

Thu thập thông tin về độ không đảm bảo và hiệp phương sai gắn với dữ liệu đo - mặc dù phụ thuộc vào lĩnh vực đo cụ thể, các ví dụ được cung cấp trong tiêu chuẩn này;

2

Cung cấp ước lượng tốt nhất của các tham số đường thẳng;

3

Xác định giá trị sử dụng của mô hình, cả về dạng chức năng (dữ liệu có phản ánh mối quan hệ đường thẳng không"?) và về mặt thống kê (độ trải của dữ liệu có nhất quán với độ không đảm bảo kèm theo của chúng không?) bằng cách sử dụng kiểm nghiệm Khi-bình phương;

4

Thu được độ không đảm bo chuẩn và hiệp phương sai gắn với ước lượng của các tham số đường thẳng;

5

S dụng hàm hiệu chuẩn để dự đoán, tức là, xác định ước lượng x ca biến X và độ không đảm bảo kèm theo tương ứng với giá trị đo được y của biến Y và độ không đảm bảo kèm theo.

Các bước trên đây được trình bày bằng sơ đồ trên Hình 1.

Mục đích chính của tiêu chuẩn này là để xem xét các bước từ 2 đến 5. Do đó, như một phần bước 1, trước khi s dụng tiêu chuẩn này, người s dụng sẽ cần cung cấp độ không đảm bảo chuẩn, và hiệp phương sai nếu liên quan, gắn với giá trị Y đo được và, khi thích hợp, độ không đảm bảo chuẩn, và hiệp phương sai nếu liên quan, gắn với giá trị X đo được, Cần tính đến các nguyên tắc của GUM khi đánh giá độ không đảm bo này dựa trên mô hình đo cụ thể cho lĩnh vực liên quan.

TCVN 9598:2013 (ISO 11095:1996) [14] liên quan đến hiệu chuẩn tuyến tính s dụng mẫu chuẩn. Những khác biệt so với tiêu chuẩn này được trình bày trong Bảng 1.

Phương pháp số được cho dựa trên Tài liệu tham khảo [6].

Hình 1 - Tóm tắt các bước trong việc xác định và sử dụng hàm hiệu chuẩn đường thẳng

Bảng 1 - Khác biệt giữa TCVN 9598:2013 (ISO 11095:1996) và tiêu chuẩn này

Đặc điểm

TCVN 9598:2013 (ISO 11095:1996)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Đề cập cụ thể đến mẫu chuẩn

Tổng quát hơn

X-giá trị giả định là biết chính xác

Thông tin độ không đảm bảo chung hơn

Tất cả các giá trị đo thu được một cách độc lập

Thông tin độ không đảm bảo chung hơn

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Không

Số kiểu cấu trúc độ không đảm bảo được xử lý

Hai

Năm, bao gồm trường hợp tổng quát nhất

Chỉ độ không đảm bo gắn với sai số ngẫu nhiên

Thông tin độ không đảm bảo tổng quát hơn

Kiểm nghiệm tính nhất quán

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Khi-bình phương

Độ không đảm bảo gắn với các dự đoán

Đặc biệt

Tương thích với GUM

 

XÁC ĐỊNH VÀ SỬ DỤNG HÀM HIỆU CHUẨN ĐƯỜNG THẲNG

Determination and use of straight-line calibration functions-

1  Phạm vi áp dụng

Tiêu chuẩn này đề cập đến hàm hiệu chuẩn tuyến tính, tức hàm hiệu chuẩn đường thẳng, mô tả mối quan hệ giữa hai biến XY, gọi là hàm có dạng Y = A + BX. Mặc dù có nhiều nguyên tắc áp dụng cho các loại hàm hiệu chuẩn tổng quát hơn, nhưng bất kỳ khi nào có thể, các cách tiếp cận được mô tả đều khai thác dạng đơn giản của hàm hiệu chuẩn đường thẳng.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Ước lượng của các tham số A B được xác định bằng cách sử dụng phương pháp bình phương tối thiểu. Điểm nhấn mạnh ca tiêu chuẩn này là việc chọn phương pháp bình phương tối thiểu thích hợp nhất cho loại dữ liệu đo, cụ thể là các phương pháp phản ánh độ không đảm bảo kèm theo. Loại ma trận hiệp phương sai tổng quát nhất gắn với dữ liệu đo được xử lý, nhưng những trường hợp đặc biệt quan trọng dẫn đến các tính toán đơn giản hơn sẽ được mô tả chi tiết.

Đối với tất c các trường hợp xem xét, phương pháp xác nhận việc sử dụng hàm hiệu chuẩn đường thẳng và đánh giá độ không đảm bảo, hiệp phương sai kèm theo các ước lượng tham số sẽ được đưa ra.

Tiêu chuẩn này cũng mô tả việc sử dụng các ước lượng tham số của hàm hiệu chuẩn và độ không đảm bảo, hiệp phương sai kèm theo của chúng đ dự đoán giá trị của X và độ không đảm bảo chun kèm theo dựa vào giá trị đo được của Y và độ không đảm bảo chuẩn kèm theo.

CHÚ THÍCH 1: Tiêu chun này không đưa ra cách xử lý chung các giá trị bất thường trong dữ liệu đo, mặc dù kiểm nghiệm hiệu lực đã cho có th sử dụng làm cơ sở để xác định dữ liệu không thống nhất.

CHÚ THÍCH 2: Tiêu chuẩn này mô tả phương pháp đánh giá độ không đảm bảo gắn với dữ liệu đo trong trường hợp độ không đảm bo ch được biết dựa vào hệ số thang đo (Phụ lục E).

2  Tài liệu viện dẫn

Các tài liệu viện dẫn trong tiêu chuẩn này rất cần thiết cho việc áp dụng tiêu chuẩn. Đối với các tài liệu có ghi năm công bố thì áp dụng bản được nêu. Đối với các tài liệu không ghi năm công bố thì áp dụng phiên bản mới nhất, bao gồm cả các sửa đổi.

TCVN 6165:2009 (ISO/IEC Guide 99:2007), Từ vựng quốc tế về đo lường học - Khái niệm, thuật ngữ chung và cơ bản (VIM)

TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008), Độ không đảm bảo đo - Phần 3: Hướng dẫn trình bày độ không đảm bảo đo (GUM:1995)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Tiêu chuẩn này áp dụng các thuật ngữ và định nghĩa trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008), TCVN 6165:2009 (ISO/IEC Guide 99:2007) và các thuật ngữ và định nghĩa dưới đây.

Giải thích các ký hiệu chính được nêu trong Phụ lục G.

3.1

Giá trị đại lượng đo được (measured quantity value)

Giá trị đại lượng đại diện cho kết quả đo.

[TCVN 6165:2009 (ISO/IEC Guide 99:2007), 2.10]

3.2

Độ không đảm bảo đo (measurement uncertainty)

Tham số không âm đặc trưng cho sự phân tán của các giá trị đại lượng được quy cho đại lượng đo, trên cơ sở thông tin đã sử dụng.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

3.3

Độ không đảm bảo đo chuẩn (standard measurement uncertainty)

Độ không đảm bảo đo được thể hiện là độ lệch chuẩn.

[TCVN 6165:2009 (ISO/IEC Guide 99:2007), 2.30]

3.4

Hiệp phương sai gắn với hai giá trị đại lượng (covariance associated with two quantity values)

Tham số đặc trưng cho sự phụ thuộc lẫn nhau của các giá trị đại lượng được quy cho hai đại lượng đo, trên cơ sở thông tin đã sử dụng.

3.5

Ma trận hiệp phương sai đo (measurement covariance matrix)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Ma trận kích thước N x N gắn với ước lượng véc tơ của đại lượng véc tơ kích thước N x 1, chứa trên đường chéo của nó các bình phương của độ không đảm bảo chuẩn gắn với các thành phần tương ứng của ước lượng véc tơ của đại lượng véc tơ, và, trên các vị trí ngoài đường chéo của nó, là các hiệp phương sai gắn với các cặp thành phần của ước lượng véc tơ của đại lượng véc tơ.

CHÚ THÍCH 1: Ma trận hiệp phương sai Ux kích thước N x N gắn với ước lượng véc tơ x của đại lưng véc tơ X có dạng

trong đó cov(xi, xj) = u2(xi) là phương sai (độ không đm bảo chuẩn bình phương) gắn với xi và cov(xi, xj) là hiệp phương sai gắn với xixj. cov(xi, xj) = 0 nếu các thành phần XiXj, của X là không tương quan.

CHÚ THÍCH 2: Các hiệp phương sai còn gọi là độ không đảm bảo tương hỗ.

CHÚ THÍCH 3: Ma trận hiệp phương sai còn gọi là ma trận phương sai-hiệp phương sai.

CHÚ THÍCH 4: Định nghĩa lấy từ [ISO/IEC Guide 98-3:2008/Suppl. 1:2008, định nghĩa 3.11 [13].

3.6

Mô hình đo (measurment model)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

[TCVN 6165:2009 (ISO/IEC Guide 99:2007), 2.48]

3.7

Mô hình hàm (functional model)

Mô hình thng kê liên quan đến các sai số gắn với biến phụ thuộc.

3.8

Mô hình cấu trúc (structural model)

Mô hình thống kê liên quan đến các sai số gắn với biến độc lập và biến phụ thuộc.

3.9

Hiệu chuẩn (calibration)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 1: Hiệu chuẩn có thể diễn tả bằng một tuyên bố, hàm hiệu chuẩn, biểu đồ hiệu chun, đường cong hiệu chuẩn, hoặc bảng hiệu chuẩn. Trong một số trường hợp nó có thể bao gồm sự hiệu chính cộng hoặc nhân của số chỉ với độ không đảm bảo đo kèm theo.

CHÚ THÍCH 2: Không được nhầm lẫn hiệu chuẩn với hiệu chỉnh hệ thống đo, thường gọi sai là “tự hiệu chuẩn”, cũng không được nhầm lẫn với kiểm định của hiệu chuẩn.

CHÚ THÍCH 3: Thông thường chỉ riêng bước đầu tiên trong định nghĩa trên được hiểu là hiệu chuẩn.

[TCVN 6165:2009 (ISO/IEC Guide 99:2007), 2.39]

3.10

Phân bố xác suất (probability distribution)

<biến ngẫu nhiên> hàm đưa ra xác suất mà biến ngẫu nhiên lấy giá trị bất kỳ đã cho hoặc thuộc tập hợp các giá trị đã cho.

CHÚ THÍCH 1: Xác suất trên toàn bộ tập hợp giá trị của biến ngẫu nhiên bằng 1.

CHÚ THÍCH 2: Phân bố xác suất được gọi là đơn biến khi nó liên quan đến một biến ngẫu nhiên (vô hướng), và đa biến khi nó liên quan đến véc tơ của các biến ngẫu nhiên. Phân bố xác suất đa biến còn được mô tả là phân bố đồng thời.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 4: Định nghĩa và Chú thích 1 lấy từ TCVN 8244-1:2010 (ISO 3534-1:2006), định nghĩa 2.11 và TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008), định nghĩa C.2.3; Chú thích 2 và 3 lấy từ ISO/IEC Guide 98-3:2008/Suppl. 1:2008, định nghĩa 3.1 [13].

3.11

Phân b chuẩn (normal distribution)

Phân bố xác suất của biến ngẫu nhiên liên tục X có hàm mật độ xác suất

với -∞ < x < +∞

CHÚ THÍCH 1: µ là kỳ vọng và σ là độ lệch chuẩn của X.

CHÚ THÍCH 2: Phân bố chuẩn còn được gọi là phân bố Gauss.

CHÚ THÍCH 3: Định nghĩa và Chú thích 1 lấy từ TCVN 8244-1:2010 (ISO 3534-1:2006), định nghĩa 2.50; Chú thích 2 lấy từ TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008), định nghĩa C.2.14.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Phân b t (t-distribution)

Phân bố xác suất của biến ngẫu nhiên liên tục X có hàm mật độ xác suất

với -∞ < x < +∞, với tham số v, số nguyên dương, là bậc tự do của phân bố, trong đó

là hàm gamma.

[ISO/IEC Guide 98-3:2008/Suppl. 1:2008, định nghĩa 3.5].

3.13

Phân bố Khi-bình phương (chi-squared distribution)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Phân bố xác suất của biến ngẫu nhiên liên tục X có hàm mật độ xác suất

với 0 ≤ x < +∞, với tham số v, số nguyên dương, trong đó  là hàm gamma.

CHÚ THÍCH: Tổng bình phương của n biến chun chuẩn hóa độc lập là biến ngẫu nhiên c2 với tham số n; khi đó n được gọi là bậc tự do.

3.14

Ma trận xác định dương (positive definite matrix)

Ma trận M kích thước n x n có tính chất zTMz > 0 đối với mọi véc tơ z khác 0 kích thước n x 1.

3.15

Ma trận nửa xác định dương (positive semi-definite matrix)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4  Quy ước và ký hiệu

Với mục đích của tiêu chuẩn này, các quy ước và ký hiệu sau được sử dụng.

4.1  X được gọi là biến độc lập và Y là biến phụ thuộc ngay cả khi hiểu biết về XY là ‘có thể hoán đổi’ ví dụ như trong Điều 7.

4.2  Đại lượng AB được gọi là các tham số của hàm hiệu chuẩn đường thẳng Y = A + BX. A và B còn được dùng để ký hiệu các biến (giả) trong biểu thức liên quan đến các tham số hàm hiệu chuẩn.

4.3  Đại lượng Xi và Yi được dùng làm biến (giả) để ký hiệu cho tọa độ của điểm dữ liệu thứ i.

4.4  Các hằng số A*B* là giá trị (chưa biết) của AB xác định hàm hiệu chuẩn đường thẳng Y = A* + B*X đối với hệ thống đo cụ thể đang xét.

4.5  Các hằng số  và  là tọa độ (chưa biết) của điểm dữ liệu thứ i do hệ thống đo cung cấp thỏa mãn .

4.6  xiyi là giá trị đo được của tọa độ của điểm dữ liệu thứ i.

4.7  ab là các ước lượng của tham số hàm hiệu chuẩn đối với hệ thống đo.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.9  Véc tơ kích thước m x 1 được ký hiệu là:

và ma trận kích thước m x n được ký hiệu là:

Kích thước của véc tơ hoặc ma trận luôn được quy định đ tránh nhầm lẫn có thể có.

4.10  T ký hiệu cho chuyển vị.

4.11  Ma trận 0 được ký hiệu bằng 0 và véc tơ đơn vị được ký hiệu là 1.

4.12  Một số ký hiệu có nhiều hơn một nghĩa. Bối cảnh sẽ làm rõ việc sử dụng.

4.13  Các số trình bày trong các bảng tới số chữ số thập phân xác định là các đại diện được làm tròn chính xác của các số lưu giữ ở độ chụm cao hơn, ví dụ như trong trường hợp bảng tính. Do đó, có th nhận thấy sự không thống nhất nhỏ giữa tổng cột hiển thị và tổng cột các số hiển thị.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.15  Trong các ví dụ, trong khi giá trị dữ liệu được cung cấp với độ chụm đã cho thì các kết quả tính toán được cung cấp với độ chụm cao hơn để cho phép người sử dụng so sánh kết quả khi thực hiện tính toán.

5  Các nguyên tắc hiệu chuẩn đường thẳng

5.1  Khái quát

5.1.1  Điều này xem xét hệ thức Y = A + BX mô tả biến phụ thuộc Y (còn gọi là ‘đáp ứng’) như hàm của biến độc lập X (còn gọi là mô phỏng’) có thể được xác định từ dữ liệu đo như thế nào. Trong bối cảnh hiệu chuẩn, dữ liệu đo phát sinh khi phương tiện đo được quy định bi các giá trị A*B* (chưa biết) của tham số hàm hiệu chuẩn được mô phng’ bằng vật mẫu với giá trị hiệu chuẩn của Xi được cho theo đơn vị của chuẩn, của tính chất của vật mẫu, và đáp ứng’ tương ứng hoặc số chỉ Yi của phương tiện được ghi lại. Hệ thức này cung cấp đáp ứng Y của hệ thống dựa vào vật mẫu có đại lượng được hiệu chuẩn X. Quá trình này được gọi là ‘đánh giá tiếp’. Hữu ích hơn trong thực tiễn là hệ thức này cho phép chuyển đổi đáp ứng đo được y của Y thành ước lượng x, theo các đơn vị của chuẩn, của tính chất X của vật mẫu. Quá trình này được gọi là ‘đánh giá ngược’ hoặc ‘dự đoán’.

5.1.2  Việc hiệu chuẩn hệ thống đo cần tính đến độ không đảm bảo đo, và, nếu có, hiệp phương sai gắn với dữ liệu đo. Đầu ra của quy trình hiệu chuẩn là hàm hiệu chuẩn được sử dụng cho dự đoán (và, nếu cần, đánh giá tiếp). Đầu ra cũng bao gồm độ không đảm bo chuẩn và hiệp phương sai gắn với ước lượng a b của tham số mô tả hàm hiệu chuẩn, được sử dụng để đánh giá độ không đảm bảo chuẩn gắn với dự đoán (và đánh giá tiếp).

5.2  Đầu vào để xác định hàm hiệu chuẩn

5.2.1  Dữ liệu đo

Thông tin yêu cầu để xác định hàm hiệu chuẩn đường thẳng là dữ liệu đo và độ không đảm bảo chuẩn kèm theo cùng với hiệp phương sai. Trong tiêu chuẩn này, dữ liệu đo được ký hiệu bằng (xi, yi), i = 1,…, m, tức là m cặp giá trị đại lượng đo được của XY. Giả định rằng m ít nhất bằng hai và các giá trị của xi không bằng nhau.

CHÚ THÍCH: Độ không đảm bảo kèm theo các ước lượng ab thường giảm khi m tăng. Do đó, hiệu chuẩn cần hướng tới sử dụng càng nhiều đim dữ liệu đo được nếu khả thi về kinh tế.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Độ không đảm bảo chuẩn gắn với xi yi được ký hiệu tương ứng là u(xi) và u(yi). Hiệp phương sai gắn với xixj được ký hiệu là cov(xi, xj). Tương tự, giá trị gắn với yiyj và với xiyj được ký hiệu tương ứng là cov(yi, yj) và cov(xi, yi). Phụ lục D chỉ ra cách độ không đảm bảo và hiệp phương sai gắn với biến đáp ứng và biến mô phỏng đo được có thể được đánh giá và đưa ra giải thích về thông tin độ không đảm bảo đó. Thông tin hoàn chỉnh về độ không đảm bảo được trình bày bằng dãy các phần tử (ma trận) U kích thước 2m x 2m chứa các phương sai (độ không đảm bảo chuẩn bình phương) u2(xi) và u2(yi) và các hiệp phương sai:

Đối với nhiều ứng dụng, một số hoặc tất cả các hiệp phương sai được lấy bằng 0 (xem 5.3).

CHÚ THÍCH: Tiêu chuẩn này liên quan đến các vấn đề trong đó u(xi) hoặc u(yi) thường khác nhau.

5.3  Xác định hàm hiệu chuẩn

5.3.1  Đầu vào để xác định hàm hiệu chuẩn là dữ liệu đo, độ không đảm bảo kèm theo và các hiệp phương sai có thể có. Cho trước các tham số AB, đầu vào có thể sử dụng để cung cấp thước đo sai lệch của điểm dữ liệu thứ i (xi, yi) so với đường thẳng Y= A + BX. Các ước lượng ab được xác định bằng cách tối thiểu hóa tng bình phương các sai lệch này, hoặc thước đo chung hơn khi mọi hiệp phương sai khác 0. Điều này đạt được như thế nào tùy thuộc vào ‘cấu trúc độ không đảm bảo’ gắn với dữ liệu đo. Cấu trúc này liên quan đến câu trả lời cho các câu hỏi sau đây:

i) Độ không đảm bảo gắn với giá trị đo được xi là không đáng kể?

ii) Hiệp phương sai gắn với cặp giá trị đo được là không đáng kể?

5.3.2  Trong tiêu chuẩn này xem xét các trường hợp sau đây, được cho theo thứ tự phức tạp tăng dần và tùy thuộc vào câu trả lời cho các câu hỏi ở 5.3.1:

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

b) Độ không đảm bảo gắn với giá trị đo được xiyi và tất cả các hiệp phương sai gắn với dữ liệu được coi là không đáng kể (Điều 7);

c) Độ không đảm bảo gắn với giá tr đo được xiyi và chỉ các hiệp phương sai gắn với cặp (xi, yi) (Điều 8);

d) Ch có độ không đảm bảo gắn với giá trị đo được yi và chỉ các hiệp phương sai gắn với yiyj (i ≠ j) (Điều 9);

e) Trường hợp chung nhất trong đó có độ không đảm bảo gắn với giá trị đo được xi và yi và hiệp phương sai gắn với tất cả các cặp giá trị của xi, xj, và yi (Điều 10).

5.3.3  Đối với từng trường hợp trong 5.3.2 được cho

a) dữ liệu đo và cấu trúc độ không đảm bảo quy định,

b) mô hình thống kê tương ứng,

c) vấn đề bình phương tối thiểu nhắm đến,

d) các bước tính toán,

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

f) giá trị sử dụng của mô hình,

g) tổ chức tính toán cho máy tính, khi thích hợp,

h) thuật toán số, khi thích hợp, và

i) một hoặc nhiều ví dụ thực tế.

5.4  Xử lý số

Phụ lục C đưa ra cách tiếp cận chung cho trường hợp tổng quát nhất là e) trong 5.3.2. Nó có thể sử dụng để xử lý tất cả các trường hợp khác và sử dụng các phương pháp tinh vi, ổn định về số học. Tuy nhiên, các trường hợp từ a) đến c) trong 5.3.2 có thể được xử lý bằng cách sử dụng các thao tác cơ bản, ví dụ, có thể thực hiện trong bảng tính. Các trường hợp d) và e) trong 5.3.2 đòi hỏi một số thao tác ma trận, dễ áp dụng theo ngôn ngữ máy tính hỗ trợ số học ma trận, nhưng không phù hợp lắm cho tính bằng bảng tính.

5.5  Độ không đảm bảo và hiệp phương sai gắn với tham số hàm hiệu chuẩn

5.5.1  Đối với tất cả các trường hợp xem xét, các ước lượng của tham số hàm hiệu chuẩn có thể được thể hiện (rõ ràng hoặc hàm ý) như hàm của dữ liệu đo. Các nguyên tắc của GUM [TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008)] có thể được áp dụng để lan truyền độ không đảm bảo và hiệp phương sai gắn với dữ liệu đo thông qua các hàm này để thu được giá trị gắn với các ước lượng tham số này. Theo cách này, dữ liệu đo được sử dụng đ cung cấp ước lượng ab của tham số hàm hiệu chuẩn, và để đánh giá độ không đảm bảo chuẩn u(a) và u(b) và hiệp phương sai cov(a,b) gắn với các ước lượng này. Đối với trường hợp a) và d) trong 5.3.2, sự lan truyền là chính xác vì ước lượng tham số có thể được biểu thị như tổ hợp tuyến tính của các đầu vào yi. Đối với các trường hợp khác, trong đó ước lượng tham số không thể biểu thị như vậy thì lan truyền là gần đúng, dựa trên tuyến tính hóa về ước lượng tham số. Vì nhiều mục đích, giá trị gần đúng thu được bằng phép tuyến tính hóa sẽ đủ chính xác.

CHÚ THÍCH: Khi lan truyền độ không đảm bảo là gần đúng, và đặc biệt nếu độ không đảm bảo liên quan lớn (ví dụ, trong một số lĩnh vực đo sinh học), thì có thể sử dụng cách tiếp cận dựa trên lan truyền phân bố. Cách tiếp cận này [ISO/IEC Guide 98-3:2008/Suppl. 1:2008] s dụng phương pháp Monte Carlo (không đề cập trong tiêu chuẩn này).

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó u(a) và u(b) là độ không đảm bảo chuẩn gắn với ab, tương ứng, và cov(a,b) = cov(b,a) là hiệp phương sai gắn với ab.

5.6  Xác nhận giá trị sử dụng của mô hình

5.6.1  Trong việc xác định ước lượng ab của tham số hàm hiệu chuẩn đường thẳng, giả định rằng mô hình Y = A + BX là hợp lý và độ không đảm bảo gắn với dữ liệu đo đưa ra thước đo tin cậy về độ lệch của dữ liệu đo so với đường thẳng. Khi ab đã được xác định, độ lệch thực tế của điểm dữ liệu so với hàm hiệu chuẩn khớp tốt nhất có thể được đánh giá dựa trên độ lệch dự đoán. So sánh này liên quan đến thước đo tổng hợp độ lệch biểu thị bằng tng các bình phương  của m phần dư có trọng số, phần dư có trọng số th i là thước đo độ lệch của điểm dữ liệu thứ i so với đường thẳng, hoặc, khi hiệp phương sai gắn với điểm dữ liệu thứ i (xi, yi) là khác 0 so với một dạng tổng quát hơn. Nếu  lớn hơn nhiều so với dự kiến thì trên cơ sở thống kê, có lý do để đặt ra câu hỏi về hiệu lực của các giả định mô hình.

5.6.2  Từ quan điểm thống kê, dữ liệu đo có thể được coi như các thể hiện của biến ngẫu nhiên. Nếu đã biết phân bố xác suất đặc trưng cho các biến ngẫu nhiên này thì về nguyên tắc có thể xác định phân bố xác suất cho thước đo tổng hợp độ lệch trong 5.6.1. Khi đó, xác suất có th được tính là , coi như được lấy từ phân bố tổng hợp này, vượt quá phân vị cụ thể bất kỳ của phân bố. Tuy nhiên, vì thông tin về các đại lượng này thường giới hạn ở bản thân các giá trị đo được và phương sai kèm theo (lấy là kỳ vọng và phương sai, tương ứng, của biến ngẫu nhiên đặc trưng bi các phân bố này), nên không có đủ thông tin để xác định phân bố cho thước đo này. Thay vào đó, việc đánh giá hiệu lực được thực hiện với giả định rằng phân bố của các đại lượng là phân bố chuẩn. Với giả định này, duy trì từ giờ về sau, ít nhất là cho mục đích xác nhận giá trị sử dụng, phân bố cho thước đo này là  với n = m - 2 bậc tự do. Theo đó, xác suất mà  vượt quá phân vị cụ thể bất kỳ của  có thể được xác định (xem 6.3, 7.3, 9.3 và 10.3). Sử dụng phân vị 95 %.

CHÚ THÍCH 1: Nếu  vượt quá phân vị 95 % của  thì hàm hiệu chuẩn đường thẳng có thể được coi như không giải thích cho dữ liệu đủ tt vì mục đích thực tế. Trong trường hợp này, dữ liệu và độ không đảm bảo kèm theo cần được kiểm tra về khả năng sai lỗi. Hàm hiệu chuẩn bao gồm đa thức của X có bậc 2 hoặc cao hơn hoặc một dạng toán học khác có thể được bàn đến; xem xét như vậy nằm ngoài phạm vi tiêu chuẩn này.

CHÚ THÍCH 2: Có khả năng mô hình 'quá tốt' chỗ giá trị quan trắc được  nhỏ hơn nhiều so với giá trị kỳ vọng. Khả năng này thường ứng với độ không đảm bảo kèm theo dữ liệu đo được trích ra quá lớn và sẽ không được xem xét thêm trong tiêu chuẩn này.

5.6.3  Để thu được càng nhiều giá trị càng tốt từ hàm hiệu chuẩn, mong muốn là độ không đảm bảo đầu vào cần được rút ra trước khi xác định tham số hàm hiệu chuẩn chứ không phải được đánh giá khi khớp với dữ liệu đã được xác định, với độ không đảm bảo kèm theo ước lượng từ dữ liệu hoặc biết theo hệ số thang đo. Trường hợp đề cập sau được xem xét trong Phụ lục E.

5.6.4  Nếu trong trường hợp cụ thể bất kỳ, việc xác nhận giá trị của mô hình không đạt, tức là  vượt quá phân vị 95 % của  (xem 5.6.2), thì độ không đảm bảo chuẩn u(a) và u(b), hiệp phương sai cov(a,b) tính được (xem 5.5.2) cần được coi là không đáng tin cậy, khi độ không đảm bảo gắn với giá trị dự đoán (xem 5.7).

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

5.7.1  Hàm hiệu chuẩn thường được sử dụng cho việc dự đoán (đánh giá ngược) trong đó, dựa vào ước lượng của Y và độ không đảm bảo chuẩn kèm theo, giá trị tương ứng của X được ước lượng và độ không đảm bảo chuẩn kèm theo được đánh giá. Việc đánh giá độ không đảm bảo đề cập sau s dụng độ không đảm bảo chuẩn gắn với các ước lượng ab cũng như hiệp phương sai kèm theo của chúng. Xem 11.1.

5.7.2  Đôi khi cần có đánh giá tiếp trong đó giá trị tương ứng của Y thu được cùng với độ không đảm bảo chuẩn kèm theo, dựa vào ước lượng của X và độ không đảm bảo kèm theo, ví dụ, khi so sánh các hiệu chuẩn của một tập hợp các phương tiện đo tương tự nhau. Xem 11.2.

CHÚ THÍCH: Giả định rằng điều kiện đo duy trì trong quá trình hiệu chuẩn giữ nguyên ở thời điểm hàm hiệu chuẩn được sử dụng về sau. Nếu không, sẽ cần một hiệu chuẩn mới hoặc thực hiện việc điều chnh thích hợp để tính đến thay đổi bất kỳ như độ trôi có th xảy ra (và độ không đảm bo kèm theo bất kỳ cũng được x lý). Biểu đồ kim soát có thể hữu ích cho mục đích này.

5.8  Xác định đường thẳng bình phương tối thiểu thông thường khớp nhất với dữ liệu

5.8.1  Đường thẳng bình phương tối thiểu thông thường (không trọng số) khớp nhất với dữ liệu được xác định bằng các giá trị ab của tham số AB tối thiểu hóa

(2)

Các giá trị này thỏa mãn công thức đã cho bằng cách cho đạo hàm riêng bậc một của biểu thức (2) đối với AB bằng 0.

5.8.2  Có thể tính ab theo các bước sau đây:

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

5.8.3  Giá trị x0y0 sao cho đường khớp nhất với các điểm dữ liệu chuyển đổi  đi qua gốc và có cùng độ dốc như với đường khớp nhất với các điểm dữ liệu gốc (xi, yi).

CHÚ THÍCH: Về mặt toán học, các tham số khớp nhất được xác định bằng cách giải cặp phương trình tuyến tính liên quan đến ma trận kích thước 2 x 2. Đối với các điểm dữ liệu chuyn đổi, ma trận này là đường chéo, cho phép dễ dàng xác định các tham s nghiêm. Chuyn đổi dữ liệu cũng có tác dụng hữu ích về độ chính xác số học của giải pháp tính toán [4, trang 33].

5.8.4  Các phương pháp mô tả trong Điều 6 đến Điều 10 dưới đây tạo thành phần mở rộng của các các tính toán trong 5.8.2 có tính đến thông tin về độ không đảm bảo quy định.

6  Mô hình dùng cho độ không đảm bảo gắn với yi

6.1  Khái quát

6.1.1  Điều này xem xét trường hợp 5.3.2 a), tức là khi thông tin dưới đây được cung cấp cho i = 1,..., m:

a) dữ liệu đo (xi, yi), và

b) độ không đảm bảo chuẩn u(yi) gắn với yi.

Phụ lục D cung cấp hướng dẫn về việc thu được các độ không đảm bảo này. Tất cả các độ không đảm bảo và hiệp phương sai khác gắn với dữ liệu được coi là không đáng kể.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

yi = A* + B*xi + ei1

i = 1,..., m,

(3)

trong đó ei là thể hiện của các biến ngẫu nhiên độc lập với kỳ vọng 0 và phương sai u2(yi) [9, trang 1]. A*B* là các giá trị (chưa biết) của tham số hàm hiệu chuẩn cho hệ thống đo yêu cầu hiệu chuẩn và cung cấp dữ liệu đo. Mô hình này được gọi là mô hình hàm, không có độ không đảm bảo gắn với xi.

6.1.3  Lấy wi = 1/u(yi), i = 1,..., m. Các ước lượng a và b là giá tr tối thiểu hóa tổng các bình phương có trọng số

(4)

đối với AB. Vấn đề tối thiểu hóa này được gọi là vấn đề bình phương tối thiểu có trọng số (WLS). Các ước lượng này thỏa mãn công thức đã cho bằng cách gn các đạo hàm riêng bậc một của biểu thức (4) đối với AB bằng 0.

6.2  Các ước lượng tham số hiệu chuẩn, độ không đảm bảo chuẩn và hiệp phương sai kèm theo

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 1: Các bước 1 đến 5 tương ứng với các bước:

CHÚ THÍCH 2: Các bước từ 1 đến 5 xác định phép tính bình phương tối thiểu cho hệ thống các công thức

wia + wixib = wiyi, i = 1,..., m.

CHÚ THÍCH 3: Nếu u(yi) giống nhau, sao cho wi giống nhau, thì ab giống với giá trị cho bi đường bình phương tối thiểu thông thường khớp nhất trong 5.8.2.

CHÚ THÍCH 4: u2(a), u2(b) và cov(a,b) trong bước 6 thu được bằng cách áp dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008) cho ab như được cung cấp bởi bước 1 đến 5.

6.2.2  Các ước lượng ab xác định trong 6.1.3 có các tính chất sau đây [15] đối với dữ liệu yi theo mô hình (3):

i) Các ước lượng ab được cho bởi tổ hợp tuyến tính các dữ liệu yi.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

iii) Ma trận hiệp phương sai đối với biến ngẫu nhiên trong ii) được quy định là u2(a), u2(b) và cov(a,b) tính được trong 6.2.1.

Tính chất i) nêu rõ ab được rút ra bằng cách sử dụng phương pháp ước lượng tuyến tính. Tính chất ii) nêu rõ rằng phương pháp ước lượng tuyến tính là không chệch. Tính chất ii) và iii) cùng cho thấy rằng phương pháp ước lượng là nhất quán theo nghĩa khi số lượng m điểm dữ liệu tăng lên, thì ước lượng ab hội tụ A*B*, tương ứng.

Phương pháp ước lượng của 6.1.3 có tính chất tối ưu sau đây đối với dữ liệu yi theo mô hình (3):

iv) Các ước lượng  và  cung cấp bởi phương pháp ước lượng tuyến tính không chệch bất kỳ có thể được coi là các thể hiện của biến ngẫu nhiên có phương sai ít nhất là lớn bằng phương sai gắn với phương pháp ước lượng WLS.

Tính chất iv) có thể được giải thích như dưới đây. Đối với hằng số cd, độ không đảm bảo chuẩn  gắn với tổ hợp tuyến tính các ước lượng  và  cung cấp bởi phương pháp ước lượng tuyến tính không chệch bất kỳ ít nhất là lớn bằng u(ca + db). Các tính chất từ i) đến iv) chứng tỏ việc sử dụng phương pháp bình phương tối thiểu cho dữ liệu tương thích với mô hình (3). Lưu ý là khi sử dụng mô hình này, các tuyên bố chỉ được đưa ra về kỳ vọng và phương sai gắn với ei; các phân bố kèm theo không được quy định thêm. Nếu đưa ra giả định bổ sung rằng ei là thể hiện của biến ngẫu nhiên phân bố chuẩn thì khi đó các tính chất khác gắn với phương pháp ước lượng WLS có thể đưa ra:

v) Các biến ngẫu nhiên trong ii) được đặc trưng bởi phân bố chuẩn hai chiều tập trung ở A*B* với ma trận hiệp phương sai quy định bởi u2(a), u2(b) và cov(a,b).

vi) Các ước lượng ab là các ước lượng hợp lý lớn nhất, tương ứng với giá trị hợp lý nhất của AB có thể làm tăng dữ liệu đo quan trắc được yi.

vii) Trong bi cảnh suy luận Bayes, phân bố trong phạm vi kiến thức cho AB, đã cho dữ liệu đo quan trắc yi, là phân bố chuẩn hai chiều tập trung ab với ma trận hiệp phương sai xác định bởi u2(a)u2(b) và cov(a,b).

6.3  Xác nhận giá trị sử dụng của mô hình

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

7  Thiết lập ri = wi(yia bxi), i = 1,..., m;

8  Thiết lập giá trị Khi-bình phương quan trắc  và bậc tự do n = m - 2;

9  Kiểm tra xem  có vượt quá phân vị 95 % của  không, và xem liệu có bác b mô hình đường thẳng hay không.

CHÚ THÍCH: Kim nghiệm Khi-bình phương dựa trên giả định là ei trong mô hình (3) là các th hiện của biến ngẫu nhiên chuẩn độc lập.

6.4  Tổ chức các tính toán

Các tính toán trong 6.2.1 và 6.3 có thể được tổ chức thành một hoặc hai bảng để thực hiện trong bảng tính như trong Bảng 2 và Bảng 3, có thể hợp nhất thành một bảng hoặc bảng tính.

Bảng 2 - Dữ liệu cho hàm hiệu chuẩn đường thẳng bằng bình phương ti thiểu có trọng s

Bảng 3 - Tổ chức tính toán để xác định hàm hiệu chuẩn đường thẳng bằng bình phương tối thiểu có trọng số

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

VÍ DỤ: (TRỌNG SỐ BẰNG NHAU) Bng 4 đưa ra sáu điểm dữ liệu và độ không đảm bảo chuẩn kèm theo. Giá trị đo được xi được lấy là chính xác và độ không đảm bảo chuẩn gắn với giá trị đo được yi u(yi) = 0,5. Do đó, trọng số được lấy là wi = 1/u(yi) = 2,0, i = 1...., 6.

Bảng 4- Dữ liệu thể hiện sáu điểm đo, trọng số bằng nhau

xi

yi

u(yi)

1,0

3,3

0,5

2,0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0,5

3,0

7,1

0,5

4,0

9,3

0,5

5,0

10,7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

6,0

12,1

0,5

Tham số đường thẳng khớp nhất được tính như trong Bảng 5. Từ bảng này, go = 84,000/24,000 = 3,500, h0 = 192,400/24,000 = 8,017, b = 123,000/70,000 = 1,757 và a = 8,017 - (1,757)(3,500) = 1,867.

Bảng 5 - Bảng tính toán gắn với dữ liệu trong Bảng 4

wi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

hi

gihi

ri

 

 

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8.017

 

 

a = 1.867

 

2.000

4.000

4.000

13.200

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-9.433

25.000

47.167

-0.648

0.419

2.000

4.000

8.000

22.400

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-4.833

9.000

14.500

0.438

0.192

2.000

4.000

12.000

28.400

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-1.833

1.000

1.833

-0.076

0.006

2.000

4.000

10.000

37.200

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.567

1.000

2.567

0.810

0.655

2.000

4.000

20.000

42.800

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

5.367

9.000

16.100

0.095

0.009

2.000

4.000

24.000

48.400

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8.167

25.000

40.833

-0.619

0.383

 

24.000

84.000

192.400

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

70.000

123.000

b = 1,757

1.665

Độ không đảm bảo chuẩn và hiệp phương sai gắn với tham số được làm khớp cũng có thể được đánh giá, sử dụng công thức trong 6.2.1, từ thông tin trong Bảng 5:

u2(a) = 1/24,000 +(3,500)2/70,000, như vậy u(a) = 0,465;

u2(b) = 1/70,000, như vậy u(b) = 0,120;

u(a, b) = -3,500/70,000 = -0,050.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Dữ liệu và hàm hiệu chuẩn đường thẳng khớp được trình bày trên Hình 2. Độ không đảm bảo chuẩn gắn với yi được minh họa (và trong các hình tiếp sau) bằng đường thẳng đứng, tập trung vào yi và có các cực trị yi - u(yi)yi + u(yi). Phần dư có trọng s được trình bày trên Hình 3.

Hình 2 - Dữ liệu trong Bảng 4 và hàm hiệu chuẩn đường thẳng khớp thu được trong Bảng 5

Hình 3 - Phần dư có trọng số ri thu được trong Bảng 5

VÍ DỤ: (TRỌNG SỐ KHÔNG BẰNG NHAU) Bảng 6 đưa ra sáu đim dữ liệu và độ không đảm bảo chuẩn kèm theo. Giá tr đo được xi được lấy là chính xác. yi thu được bằng cách sử dụng hai thiết lập phương tiện, sao cho giá trị X càng lớn thì yi càng kém chính xác.

Bảng 6 - Dữ liệu th hiện sáu điểm đo, trọng s không bằng nhau

xi

yi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1,0

3,2

0,5

2,0

4,3

0,5

3,0

7,6

0,5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8,6

1,0

5,0

11,7

1,0

6,0

12,8

1,0

Tham số đường thẳng khớp nhất được tính như trong Bảng 7. Từ bảng này, g0 = 39,000/15,000 = 2,600, h0 = 93,500/15,000 = 6,233, b = 65,000/31,600 = 2,057 và a = 6,233 - (2,057)(2,600) = 0,885.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

w2

gi

hi

gihi

ri

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

 

 

 

2.600

6.233

 

 

a = 0.885

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.000

4.000

4.000

12.800

-3.200

-0,007

10.240

10.413

0.510

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.000

4.000

8.000

17.200

-1.200

-3.867

1.440

4.040

-1.398

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.000

4.000

12.000

30.400

0.800

0.733

0.640

2.187

1.088

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.000

1.000

4.000

8.600

1.400

2.367

1.960

3.313

-0.513

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.000

1.000

5.000

11.700

2.400

5.467

5.760

13.120

0.830

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.000

1.000

6.000

12.800

3.400

6.567

11.560

22.327

-0.427

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

15.000

39.000

93.500

 

 

31.600

63.000

b = 2.057

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Độ không đm bảo chuẩn và hiệp phương sai gắn với tham số được làm khớp cũng có thể được đánh giá, sử dụng công thức trong 6.2.1, từ thông tin trong Bảng 7:

u2(a) = 1/15,000 +(2,600)2/31,600, như vậy u(a) = 0,530;

u2(b) = 1/31,600, như vậy u(b) = 0,178;

u(a, b) = -2,600/31,600 = -0,082.

Giá trị Khi-bình phương quan trắc được là  = 4,131 với v = 4 bậc tự do, như tính được ở Bảng 7 sử dụng 6.3. Vì  không vượt quá phân vị 95 % của  tức là 9,488, nên không có lý do để nghi ngờ sự nhất quán của mô hình đường thẳng và dữ liệu.

Dữ liệu và hàm hiệu chuẩn đường thẳng khớp được trình bày trên Hình 4. Phần dư có trọng số được trình bày trên Hình 5.

Hình 4 - Dữ liệu trong Bảng 6 và hàm hiệu chuẩn đường thẳng khớp thu được trong Bảng 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

7  Mô hình độ không đảm bảo gắn với xi và yi

7.1  Khái quát

7.1.1  Điều này xem xét trường hợp 5.3.2 b), tc là khi thông tin dưới đây được cung cấp cho i = 1,…, m:

a) dữ liệu đo (xi, yi),

b) độ không đảm bảo chuẩn u(xi) gắn với xi, và

b) độ không đảm bảo chuẩn u(yi) gắn với yi.

Phụ lục D cung cấp hướng dẫn về việc thu được các độ không đảm bảo này. Tất cả các hiệp phương sai khác gắn với dữ liệu được coi là không đáng kể.

7.1.2  Trường hợp 5.3.2 b) tương ứng với mô tả bởi mô hình thống kê

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

7.1.3  Vì xi (cùng với yi - xem Điều 6) có độ không đảm bảo kèm theo nên cũng phải tính đến chúng khi xác định hàm hiệu chuẩn đường thẳng. Vấn đề xác định ab trong bối cảnh này là một trong hồi quy khoảng cách trực giao có trọng số (ODR có trọng số) hoặc hồi quy khoảng cách tổng quát hóa (GDR) [2]. Trong tài liệu thống kê nó được gọi là mô hình sai số-trong-biến [7], [9, trang 50], [17, trang 189], Các ước lượng ab là ước lượng tối thiểu hóa tng các bình phương

(6)

đối với A, BXi, i = 1,..., m, các trọng số ni = 1/u(xi) wi = 1/u(yi). Mỗi ước lượng giải , cùng với ab, xác định ước lượng , của  trong mô hình (5).

7.1.4  Cho trước AB, các giá trị  tối thiểu hóa tổng các bình phương (6) đối với Xi được cho bởi

(7)

Sử dụng biểu thức (7) vấn đề tối ưu hóa có thể được đặt ra đối với tham số AB độc lập bằng cách thay thế Xi trong biểu thức (6) bằng ,

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

thì tổng các bình phương (8) tương đương với

Số hạng Ri có giải thích hình học dưới đây. Véc tơ vuông góc với đường thẳng Y = A + BX được cho bởi (-B,1)T/(1 + B2)1/2Ri là bội số trọng số của thành phần có dấu của  theo chiều của véc tơ vuông góc.

CHÚ THÍCH 1: Bình phương tối thiểu thông thường (xem 5.8) và bình phương tối thiểu có trọng số (xem Điều 6), khoảng cách tới đường thẳng được đo ‘theo chiều thẳng đứng’, tức là, theo chiu Y, phản ánh thực tế là độ lệch của đim đo được (xi, yi) so với đường thẳng có thể được tính về số hạng sai số ei kèm theo yi, vì xi được giả định là đã biết chính xác. ODR có trọng số gii quyết trường hợp trong đó còn có độ không đảm bảo gắn với xi.

CHÚ THÍCH 2: Biểu thức (7) được cho bằng cách cho đạo hàm riêng bậc một của biểu thức (6) đối với A, BXi bằng 0.

CHÚ THÍCH 3: Nếu ui(xi) = 0 thì khi đó trong biểu thức (7)  (như vậy  và  . Kết quRi trong biu thức (9) được cho bởi Ri = wi(yi A – Bxi). Do đó, nếu u(xi) = 0 thì số hạng Ri được đánh giá theo cách tương tự như trong biểu thức (4) ở 6.1.3.

CHÚ THÍCH 4: Nếu u(xi) = u(yi) = ui thì khi đó  xác định điểm trên đường thẳng Y= A + BX gần nhất với (xi, yi) và

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

7.1.6  Điều 7.1.3 liên quan đến A, BXi, i = 1,..., m, như các biến trong tối thiểu hóa. Các tính toán nêu trong 7.2.1 thực hiện quá trình tối thiểu hóa này trong phép lặp hai giai đoạn [2]:

i) từ các gần đúng với ab, xác định tối ưu  tương ứng, và

ii) đối với các  này, xác định các gần đúng mới với ab sẽ làm giảm tổng các bình phương (6).

CHÚ THÍCH: Không có phân biệt nào về ký hiệu giữa   phép lặp điển hình và giá trị tính được cuối cùng.

7.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn và hiệp phương sai kèm theo

7.2.1  Các ước lượng ab được tính ở các bước từ 1 đến 6 dưới đây sử dụng quy trình lặp chỉ ra ở 7.1.6; độ không đảm bảo chuẩn u(a), u(b) và hiệp phương sai cov(a,b) được đánh giá ở bước 7 (xem Phụ lục B):

1  Thu được các xấp xỉ ban đầu  và  cho ab, ví dụ, bằng cách xác định đường khớp nhất bình phương tối thiểu có trọng số với dữ liệu (xem 6.2.1 bước 1 đến 5), bỏ qua độ không đảm bảo gắn với ;

4  Xác định giải pháp bình phương tối thiểu (không có trọng số) δa và δb cho hệ thức

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

đó là

5  Cập nhật các xấp xỉ hiện tại cho các tham số và phần dư: , i = 1,..., m

6  Lặp lại các bước từ 2 đến 5 cho đến khi đạt được sự hội tụ. Đặt ;

;

trong đó g0, h0, ...là các giá trị tính được bước 4.

CHÚ THÍCH 1: Các tính toán ở bước 4 tương tự như các tính toán bước 1 đến 5 trong 6.2.1.

CHÚ THÍCH 2: Ở bước 2,  xác định điểm  trên xấp x hiện tại với hàm hiệu chuẩn đường thẳng khớp nhất là gần nhất, như khoảng cách có trọng số, với điểm dữ liệu đo được (xi, yi).

CHÚ THÍCH 3: Ở bước 3, hi tính được đại diện cho giá trị của khoảng cách tổng quát hóa Ri trong biểu thức (9) từ điểm dữ liệu thứ i đến ước lượng hiện tại của hàm hiệu chuẩn đường thẳng. Thuật toán được thiết kế đ tối thiểu hóa tổng các bình phương của các khoảng cách này.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 5: Phần dư tính được c 5 được gắn với phép giải hệ thức tính bước 4. Tại điềm hội tụ, ri tính được bước 5 bằng hi tính được bước 3.

CHÚ THÍCH 6: Hoàn toàn ch yêu cầu các phần dư tính được bước 5 trong lần lặp cuối. Tuy nhiên, để dễ trình bày dạng bảng (Bng 9 7.4), phần dư được tính tại mỗi lần lặp.

CHÚ THÍCH 7: u2(a), u2(b) và cov(a,b) trong bước 7 thu được bằng cách áp dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008) cho ab như được cung cấp bởi các bước 1 đến 6.

7.2.2  Trong khi có thể rút ra các tính chất trong 6.2.2 đối với phương pháp ước lượng WLS thì thực tế là các ước lượng ab xác định bằng cách tối thiểu hóa tổng các bình phương (6) phụ thuộc phi tuyến vào dữ liệu xiyi nghĩa là không thể dễ dàng nêu rõ các tính chất tương ứng đối với ODR có trọng số. Các ước lượng ab xác định trong 7.1.3 có các tính chất sau đây đi với dữ liệu xiyi theo mô hình (5):

i) Các ước lượng ab được cho bởi hàm phi tuyến của các dữ liệu xi và yi.

ii) Các ước lượng ab được coi là các thể hiện của biến ngẫu nhiên có kỳ vọng tương ứng là A*B* gần đúng.

iii) Các thành phần của ma trận hiệp phương sai đối với biến ngẫu nhiên trong ii) xấp xỉ bằng u2(a), u2(b) và cov(a,b) tính được trong 7.2.1.

Các gần đúng trong ii) và iii) sẽ chính xác hơn đối với dữ liệu có độ không đảm bảo kèm theo nhỏ hơn. Tuy nhiên, phương pháp ước lượng có tính chất nhất quán sau đây:

iv) Đối với dữ liệu thỏa mãn mô hình (5), khi số lượng m điểm dữ liệu tăng, các ước lượng ab hội tụ đến A*B*, tương ứng [16].

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Nếu giả định bổ sung được đưa ra là diei là các thể hiện của biến ngẫu nhiên phân bố chuẩn thì các tính chất khác gắn với phương pháp ước lượng ODR có trọng số có thể được nêu ra:

v) Các biến ngẫu nhiên trong ii) được đặc trưng gần đúng bởi phân bố chuẩn hai chiều tập trung ở A*B* với ma trận hiệp phương sai quy định bởi u2(a), u2(b) và cov(a,b).

vi) Các ước lượng ab là ước lượng hp lý lớn nhất, tương ứng với giá trị hợp lý nhất của AB có th làm tăng dữ liệu đo quan trắc được xiyi.

vii) Trong bối cảnh suy luận Bayes, phân bố trong phạm vi kiến thức cho AB, đã cho dữ liệu đo quan trắc xiyi, là xấp xỉ phân bố chuẩn hai chiều tập trung ab với ma trận hiệp phương sai xác định bởi u2(a), u2(b) và cov(a,b).

7.3  Xác nhận giá trị sử dụng của mô hình

Nếu m > 2 thì giá trị sử dụng của mô hình có thể kiểm nghiệm một phần bằng cách sử dụng phần dư ri có trọng số tính được trong bước 5 ở 7.2.1 ở lần lặp cuối cùng, tức là sự hội tụ (tiếp tục từ 7.2.1):

8  Thiết lập giá trị Khi-bình phương quan trắc  và bậc tự do n = m - 2;

9  Kiểm tra xem  có vượt quá phân vị 95 % của  không, và xem liệu có bác bỏ mô hình đường thẳng hay không.

CHÚ THÍCH: Kim nghiệm Khi-bình phương dựa trên giả định là diei trong mô hình (5) là các th hiện của biến ngẫu nhiên chuẩn độc lập và phép gần đúng bậc một.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Các tính toán trong 7.2.1 và 7.3 có thể được tổ chức thành hai bảng liên tiếp, phù hợp để thực hiện trong bảng tính. Bảng thứ nhất (Bảng 8), đưa ra các xấp xỉ  và  (xem 7.2.1 bước 1), tính fi, gihi (xem 7.2.1 bước 3). Bảng thứ hai (Bng 9) sử dụng các giá trị fi, gihi này để tính các hiệu chính δa và δb (xem 7.2.1 bước 4).

Bảng 8 - Tính toán để xác định hàm hiệu chuẩn đường thẳng, đã cho các xấp x  và  với các ước lượng tham số đường thẳng a b

Bng 9 - T chức tính toán để xác định các hiệu chính δa và δb cho hàm hiệu chuẩn đường thng GDR

VÍ DỤ: Bảng 10 đưa ra sáu điểm dữ liệu đo được và độ không đảm bảo chuẩn kèm theo.

Bảng 10 - Sáu điểm dữ liệu đo được và độ không đảm bảo chuẩn tương ứng

xi

w(xi)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

u(yi)

1,2

0,2

3,4

0,2

1,9

0,2

4,4

0,2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0,2

7,2

0,4

4,0

0,2

8,5

0,4

4,7

0,2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0,4

5,9

0,2

13,5

0,4

Để xác định các xấp xỉ  và  ban đầu (7.2.1 bước 1), hàm hiệu chuẩn đường thẳng bình phương tối thiểu có trọng số được xác định. Tiếp theo quy trình mô tả trong 6.2, thu được các dữ liệu cho trong Bảng 11 và Bảng 12.

Bảng 11 - Dữ liệu đại diện cho sáu điểm do

xi

yi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1,2

3,4

0,2

1,9

4,4

0,2

2,9

7,2

0,4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8,5

0,4

4,7

10,8

0,4

5,9

13,5

0,4

Bảng 12 - Bảng tính toán gắn với dữ liệu trong Bảng 11 đ xác định các xấp xỉ  và  ban đầu

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

gi

hi

gi hi

ri

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

 

 

2.573 3

6.186 7

 

 

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

25.000 0

30.000 0

65.000 0

-6.866 7

-13.933 3

47.151 1

95.675 6

0.818 6

0.670 1

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

25.000 0

47.500 0

110.000 0

-3.366 7

-8.933 3

11.334 4

30.075 6

-1.700 6

2.892 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

25.000 0

72.500 0

180.000 0

1.633 3

5.066 7

2.667 8

8.275 6

1.557 7

2.426 4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

6.250 0

25.000 0

53.125 0

3.566 7

5.783 3

12.721 1

20.627 2

-1.879 1

3.531 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

6.250 0

29.375 0

67.500 0

5.316 7

11.533 3

28.266 9

61.318 9

0.111 3

0.012 4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

6.250 0

36.875 0

84.375 0

8.316 7

18.283 3

63.166 9

152.056 4

0.416 3

0.173 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

93.750 0

241.250 0

580.000 0

 

 

171.308 3

368.029 2

b = 2.148 3

9.705 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Vi các giá trị  và  cập nhật này, hai dữ liệu bảng mới được hình thành (Bảng 15 và Bảng 16) đề xác định thêm các hiệu chính δa = = -0,001 0 và δb = 0,000 2. Quá trình này được lặp lại lần thứ ba (Bảng 17 và Bảng 18). Trong trường hợp này, độ lớn của các hiệu chính nhỏ hơn 0,000 05, được điều chỉnh để thành không đáng kể cho mục đích này, và các xấp xỉ cuối cùng cho ước lượng tham số là a = 0,578 8 và b = 2,159 7.

Độ không đảm bảo chuẩn và hiệp phương sai (7.2.1 bước 7) gắn với tham số được làm khớp cũng có thể được đánh giá từ thông tin trong bảng cuối cùng (Bảng 18):

Giá trị Khi-bình phương quan trắc được là  = 2,743 với n = 4 bậc tự do, như tính được Bảng 18 sử dụng 7.3. Vì  không vượt quá phân vị 95 % của , tức là 9,488, nên không có lý do để nghi ngờ sự nhất quán của mô hình đường thẳng và dữ liệu.

Các điểm dữ liệu và hàm hiệu chuẩn đường thẳng ODR có trọng số được trình bày bằng đồ thị trên Hình 6. Đối với mỗi giá trị i, đồ thị cũng đưa ra vị trí của , điểm trên đường thẳng gần nhất về mặt xác suất với điểm dữ liệu (xi,yi). Phần dư có trọng số được minh họa trên Hình 7.

Bảng 13 - Lần lặp đầu tiên để xác định fi, gihi, cho trước  và

xi

u(xi)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

u(yi)

ti

zi

fi

gi

hi

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

0.658 3

2.148 3

 

 

 

 

1.200 0

0.200 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.200 0

4.452 2

1.262 6

0,163 7

2.110 0

2.664 2

0.345 5

1.900 0

0.200 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.200 0

4.452 2

1.769 9

-0.340 1

2.110 0

3.734 5

-0.717 6

2.900 0

0.200 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.200 0

4.452 2

3.019 2

0.311 6

2.110 0

6.370 6

0.657 5

4.000 0

0.200 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.400 0

2.901 9

3.812 6

-0.751 5

1.703 5

6.494 7

-1.280 2

4.700 0

0.200 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0,400 0

2.901 9

4.711 1

0.044 7

1.703 5

8.025 3

0.076 1

5.900 0

0.200 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.400 0

2.901 9

5.941 6

0.166 7

1.703 5

10.121 4

0.284 0

Bảng 14 - Lần lặp đầu tiên để xác định các hiệu chính δa và δb, cho trước fi, gihi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

fihi

ri

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

3.123 9

-0.043 7

 

 

δa = -0.078 4

 

4.452 2

5.621 6

0.729 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.437 8

15.423 0

-1.719 3

0.481 4

0.231 8

4.452 2

7.879 9

-1.514 1

-2.857 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8.162 3

1.786 4

-0.593 5

0.352 3

4.452 2

13.442 2

1.387 4

0.220 9

0.749 8

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.165 6

0.752 3

0.565 9

2.901 9

11.063 6

-2.180 7

1.178 2

-1.205 7

1.376 4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-1.218 7

1.485 2

2.901 9

13.671 0

0.129 7

2.703 8

0.150 6

7.310 8

0.407 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.014 5

2.901 9

17.241 6

0.483 5

4.799 9

0.358 5

23.038 7

1.720 6

0.305 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

22.062 2

68.919 9

-0.964 8

 

 

55.360 6

0.615 2

δb = 0.011 1

2.742 9

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

xi

u(xi)

yi

u(yi)

ti

zi

fi

gi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

 

 

 

0.579 9

2.159 4

 

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.200 0

0.200 0

3.400 0

0.200 0

4.414 6

1.287 3

0.228 9

2.101 1

2.704 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.900 0

0.200 0

4.400 0

0.200 0

4.414 6

1.792 2

-0.282 7

2.101 1

3.705 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.900 0

0.200 0

7.200 0

0.200 0

4.414 6

3.030 5

0.357 9

2.101 1

6.379 9

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.000 0

0.200 0

8.500 0

0.400 0

2.885 8

3.821 2

-0.717 5

1.698 8

6.491 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.700 0

0.200 0

10.800 0

0.400 0

2.885 8

4.717 7

0.070 9

1.698 8

8.014 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

5.900 0

0.200 0

13.500 0

0.400 0

2.885 8

5.944 8

0.179 6

1.698 8

10.098 8

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Bảng 16 - Như Bảng 14 nhưng cho lần lặp thứ hai

figi

fihi

ri

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

 

 

3.141 2

-0.000 3

 

 

δa = -0.001 0

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

5.682 7

1.010 3

-3.895 3

0.481 4

15.173 4

-1.875 1

0.482 3

0.232 6

4.414 6

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-1.248 2

-2.834 4

-0.593 5

8.033 9

1.682 2

-0.592 8

0.351 4

4.414 6

13.404 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-0.220 1

0.752 4

0.048 4

-0.165 6

0.752 5

0.566 2

2.885 8

11.027 1

-2.070 6

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-1.218 4

1.334 2

-1.407 4

-1.218 7

1.485 2

2.885 8

13.614 3

0.204 6

2.678 1

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

7.172 0

0.323 8

0.120 3

0.014 5

2.885 8

17.155 5

0.518 3

4.762 6

0.305 6

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.455 3

0.304 4

0.092 7

21.901 2

68.790 1

-0.005 7

 

 

54.444 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

δb = 0.000 2

2.742 7

Bảng 17 - Như Bảng 13 nhưng cho lần lặp thứ ba

xi

u(xi)

yi

u(yi)

ti

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

fi

gi

hi

 

 

 

 

0.578 8

2.159 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

 

 

1.200 0

0.200 0

3.400 0

0.200 0

4.413 8

1.287 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.100 9

2.705 0

0.482 3

1.900 0

0.200 0

4.400 0

0.200 0

4.413 8

1.792 4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.100 9

3.765 7

-0.592 8

2.900 0

0.200 0

7.200 0

0.200 0

4.413 8

3.036 6

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

2.100 9

6.379 5

0.752 5

4.000 0

0.200 0

8.500 0

0.400 0

2.885 5

3.821 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.698 7

6.490 9

-1.218 7

4.700 0

0.200 0

10.800 0

0.400 0

2.885 5

4.717 6

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.698 7

8.013 7

0.120 3

5.900 0

0.200 0

13.500 0

0.400 0

2.885 5

5.944 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.698 7

10.098 0

0.304 4

Bảng 18 - Như Bảng 14 nhưng cho lần lặp thứ ba

figi

fihi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

ri

 

 

 

 

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

δa = 0.000 0

 

4.413 8

5.682 9

1.013 3

-3.894 7

0.452 3

15.168 5

-1.878 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.232 7

4.413 8

7.911 3

1.245 4

-2.834 0

-0.592 8

8.031 5

1.680 0

-0.592 8

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.413 8

13.402 7

1.580 9

-0.220 2

0.752 5

0.048 5

-0.165 7

0.752 5

0.566 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

11.025 8

-2.070 2

1.154 8

-1.218 7

1.333 5

-1.407 3

-1.218 7

1.485 2

2.885 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.204 3

2.677 6

0.120 3

7.169 5

0.322 0

0.120 3

0.014 5

2.885 5

17.153 1

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.761 9

0.304 4

22.675 6

1.449 6

0.304 4

0.092 7

21.897 7

68.788 4

0.000 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

54.427 1

0.000 1

δb = 0.000 0

2.742 7

 

Hình 6 - Dữ liệu trong Bảng 10 và hàm hiệu chuẩn đường thẳng khớp thu được trong Bảng 11 đến 18

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8  Mô hình độ không đảm bảo gắn với xi và yi và hiệp phương sai gắn với cặp (xi, yi)

8.1  Khái quát

8.1.1  Điều này xem xét trường hợp 5.3.2 c), tức là khi thông tin sau đây được cung cấp cho i = 1,…, m:

a) dữ liệu đo (xi, yi),

b) độ không đảm bảo chuẩn u(xi) gắn với xi,

c) độ không đm bo chuẩn u(yi) gắn với yi, và

d) hiệp phương sai cov(xi, yi) gắn với xiyi.

Phụ lục D cung cấp hướng dẫn về việc thu được các độ không đảm bảo và hiệp phương sai này. Tất cả các hiệp phương sai khác gắn với dữ liệu được coi là không đáng kể.

8.1.2  Trường hợp 5.3.2 c) tương ứng với mô tả bởi mô hình thống kê

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó mỗi cặp (di,ei) là thể hiện của các biến ngẫu nhiên hai chiều với kỳ vọng (0,0)T và ma trận hiệp phương sai có thành phần đường chéo u2(xi)u2(yi) và thành phần ngoài đường chéo cov (xi, yi) = cov(yi,xi), tức là

độc lập với các biến ngẫu nhiên khác.

CHÚ THÍCH: Giả định rng (di,ei) là các th hiện của biến ngẫu nhiên chuẩn hai chiều chỉ cần thiết cho việc xác nhận giá trị sử dụng của mô hình (10).

8.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn và hiệp phương sai kèm theo

8.2.1  Về thuật toán, trường hợp này có thể được xử lý bằng cách mở rộng (xem Phụ lục B) xử lý trong Điều 7. Các tính toán giống hệt như trong Điều 7, ngoại trừ bước 2 trong 7.2.1 được thay bằng

8.2.2  Tất cả các tính chất nêu trong 7.2.2 áp dụng cho dữ liệu được tạo ra theo mô hình (10) và tuân thủ theo phần còn lại của Điều 7.

9  Mô hình độ không đảm bảo và hiệp phương sai gắn với yi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

9.1.1  Điều này xem xét trường hợp 5.3.2 d), tức là khi thông tin sau đây được cung cấp cho i = 1,..., m:

a) dữ liệu đo (xi, yi),

b) độ không đảm bảo chuẩn u(yi) gắn với yi,

c) hiệp phương sai cov(yi, yj) gắn với cặp (yi, yj) j = 1,... , m, j ≠ i.

9.1.2  Độ không đảm bảo chuẩn bình phương và hiệp phương sai bao gồm ma trận hiệp phương sai

kích thước m x m gắn với y = (yi,..., ym)T. Phụ lục D cung cấp hướng dẫn về việc thu được các độ không đảm bảo và hiệp phương sai này. Tất cả các độ không đảm bảo và hiệp phương sai khác gắn với dữ liệu được coi là không đáng kể.

9.1.3  Trường hợp 5.3.2 d) tương ứng với mô tả bởi mô hình thống kê

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó e = (e1,..., em)T là thể hiện của biến ngẫu nhiên đa biến với kỳ vọng véc tơ bằng véc tơ 0 có kích thước m x 1 và ma trận hiệp phương sai có kích thước m x m bằng Uy [21].

9.1.4  Các ước lượng ab là giá trị tối thiểu hóa tng các bình phương tổng quát [8]

trong đó e = y - A1 - Bx, đối với AB. Vấn đề xác định ab trong bối cảnh này được gọi là một trong hồi quy Gauss-Markov (GMR) [2].

CHÚ THÍCH: Đối với trường hợp trong đó Uy là đường chéo, tổng bình phương tng quát hóa (12) đơn giản thành biểu thức (4) trong 6.1.3 dẫn đến vấn đề WLS.

9.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn và hiệp phương sai kèm theo

9.2.1  Nếu Uy hoàn toàn dương, thì tồn tại thừa số Cholesky Ly tam giác dưới kích thước m x m của Uy = LyLyT [10] (xem thêm A.4), ước lượng ab của AB có thể được tính trực tiếp bằng cách sử dụng cùng một quy trình chung như trong 6.2.1 sau một số tính toán sơ bộ sử dụng thao tác ma trận-véc tơ. Nếu không có thể cần đến nhiều phương pháp số học hơn. Các thao tác này chuyển đổi tổng bình phương tổng quát hóa (12) thành tổng bình phương thông thường (2) như trong 5.8.1, tức là, vấn đề trở thành bình phương tối thiểu không trọng số và không có hiệp phương sai.

9.2.2  Ước lượng tham số ab được tính ở các bước từ 1 đến 7 dưới đây; độ không đảm bảo chuẩn u(a)u(b) và hiệp phương sai cov(a,b) được đánh giá trong bước 8:

1  Tính thừa số Cholesky Ly kích thước m x m của Uy = LyLyT; xem A.4.1;

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

9.2.3  Ước lượng ab xác định trong 9.1.4 có các tính chất sau đây [15] đối với dữ liệu yi theo mô hình (11):

i) Các ước lượng ab được cho bởi tổ hợp tuyến tính các dữ liệu yi.

ii) Các ước lượng ab được coi là các thể hiện của biến ngẫu nhiên có kỳ vọng tương ứng là A*B*.

iii) Ma trận hiệp phương sai đối với biến ngẫu nhiên trong ii) được quy định là u2(a), u2(b) và cov(a,b) tính được trong 9.2.2.

Tính chất i) nêu rõ ab được rút ra bằng cách sử dụng phương pháp ước lượng tuyến tính. Tính chất ii) nêu rõ rằng phương pháp ước lượng tuyến tính là không chệch. Tính cht ii) và iii) cùng cho thấy rằng phương pháp ước lượng là nhất quán theo nghĩa khi số lượng m điểm dữ liệu tăng lên, thì ước lượng ab hội tụ A*B*, tương ứng.

Phương pháp ước lượng của 9.1.4 có tính chất tối ưu sau đây đối với dữ liệu yi theo mô hình (11):

iv) Các ước lượng  và  cung cấp bởi phương pháp ước lượng tuyến tính không chệch bất kỳ có thể được coi là các thể hiện của biến ngẫu nhiên có phương sai ít nhất là lớn bằng phương sai gắn với phương pháp ước lượng GMR.

Tính chất iv) có thể được giải thích như dưới đây. Đối với hằng số cd, độ không đảm bảo chuẩn  gắn với tổ hp tuyến tính các ước lượng  và  cung cấp bởi phương pháp ước lượng tuyến tính không chệch bất kỳ ít nhất là lớn bằng u(ca + db). Các tính chất từ i) đến iv) chứng tỏ việc s dụng phương pháp bình phương tối thiểu cho dữ liệu tương thích với mô hình (11). Lưu ý là khi s dụng mô hình này, các tuyên bố chỉ được đưa ra về kỳ vọng và phương sai gắn với ei; các phân bố kèm theo không được quy định thêm. Nếu đưa ra giả định bổ sung rằng ei là thể hiện của biến ngẫu nhiên đặc trưng bởi phân bố chuẩn nhiều chiều thì khi đó các tính chất khác gắn với phương pháp ước lượng GMR có thể đưa ra:

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

vi) Các ước lượng ab là ước lượng hợp lý lớn nhất, tương ứng với giá trị hợp lý nhất của AB có thể làm tăng dữ liệu đo quan trắc được yi.

vii) Trong bối cảnh suy luận Bayes, phân bố trong phạm vi kiến thức cho AB, đã cho dữ liệu đo quan trắc yi, là phân bố chuẩn hai chiều tập trung ở a b với ma trận hiệp phương sai xác định bởi u2(a), u2(b) và cov(a,b).

CHÚ THÍCH 1: Các tính chất liệt kê trên giống với trong phương pháp ước lượng WLS của 6.1.3 đối với dữ liệu tạo ra theo mô hình (3).

CHÚ THÍCH 2: u2(a), u2(b) và cov(a,b) trong bước 8 thu được bằng cách áp dụng định luật lan truyền độ không đm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008) cho ab như cung cấp bởi các bước từ 1 đến 7.

9.3  Xác nhận giá trị sử dụng của mô hình

Nếu m > 2 thì giá trị sử dụng của mô hình có thể kiểm nghiệm một phần bằng cách sử dụng phần dư ri có trọng số (tiếp theo 9.2.2):

9  Thiết lập

10  Thiết lập giá trị Khi-bình phương quan trắc được  và bậc tự do n = m - 2;

11  Kiểm tra xem  có vượt quá phân vị 95 % của  không, và xem liệu có bác bỏ mô hình đường thẳng hay không.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

9.4  Tổ chức tính toán

Các tính toán trong 9.2.2 và 9.3 có thể được tổ chức thành các bảng như trong các bảng từ 19 đến 21. Bảng 20 chứa fi, gihi tính được ở bước 1 và 2 trong 9.2.2 đối với thừa số hóa Cholesky Ly của ma trận hiệp phương sai Uy. Bảng 21 sử dụng các fi, gihi này để tính các ước lượng ab của tham số hàm hiệu chuẩn đường thẳng.

Bảng 19 - Dữ liệu cho hàm hiệu chuẩn đường thẳng Gauss-Markov

Bảng 20 - Tính toán ban đầu để xác định hàm hiệu chuẩn đường thẳng Gauss-Markov

Bảng 21 - Tổ chức tính toán để xác định hàm hiệu chuẩn đường thẳng Gauss-Markov

VÍ DỤ: Bảng 22 đưa ra mười điểm dữ liệu đo được (xi, yi) và độ không đm bo chuẩn gắn với yi.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Bảng 22 - Dữ liệu thể hiện mười điểm đo, yi có ma trận hiệp phương sai kèm theo

xi

yi

1,0

1,3

2,0

4,1

3,0

6,9

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

7,5

5,0

10,2

6,0

12,0

7,0

14,5

8,0

17,1

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

19,5

10,0

21,0

Ma trận hiệp phương sai Uy kích thước 10 x 10 gắn với yi

Thừa số Cholesky Ly kích thước 10 x 10 của Uy = LyLyT, tính được sử dụng một trong hai thuật toán mô tả trong A.4.1, là

Các véc f, g và h trong Bảng 23 được tính theo bước 2 trong 9.2.2.

Bảng 23 - Tính toán ban đầu bảng kèm theo dữ liệu trong Bảng 22

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

gi

hi

0,707 1 0,408 2 0,288 7 0,223 6 0,182 6 0,447 2 0,149 1 0,092 5 0,067 3 0,052 9

0,707 1 1,224 7 1,732 1 2,236 1 2,738 6 2,683 3 1,639 8 1,849 0 2,219 8 2,646 3

0,919 2 2,816 9 4,416 7 3,957 8 5,696 3 5,366 6 3,652 2 4,428 4 5,320 8 5,536 0

Tham số đường thẳng khớp nhất trong Bảng 24 được tính theo Bảng 21. Từ Bảng 24, g0 = 4,404 8/1,071 4 = 4,111 1, h0 = 9,0048/1,071 4 = 8,404 4, b = 54,218 5/24,629 6 = 2,201 4 và a = 8,404 4 (2,201 4)(4,111 1) = -0,645 6.

Bảng 24 - Bảng tính toán kèm theo với dữ liệu trong Bảng 22

figi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

ri

 

 

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8.404 4

 

 

a = -0,645 6

 

0.500 0

0.500 0

0.650 0

-2.199 9

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4.839 5

11.051 4

-0.180 9

0.032 7

0.166 7

0.500 0

1.150 0

-0.453 6

-0.614 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.278 6

0.394 4

0.147 7

0.083 3

0.500 0

1.275 0

0.545 3

1.990 6

0.297 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.790 2

0.624 5

0.050 0

0.500 0

0.885 0

1.316 8

2.078 5

1.734 0

2.737 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.672 7

0.033 3

0.500 0

1.040 0

1.988 0

4.161 9

3.952 3

8.273 9

-0.214 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.200 0

1.2000

2.400 0

0.844 7

1.608 0

0.713 6

1.358 3

-0.251 6

0.063 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.244 4

0.544 4

1.020 9

2.399 4

1.054 6

2.464 0

0.138 7

0.019 2

0.008 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.409 4

1.468 9

3.651 4

2.457 8

5.363 6

0.417 7

0.174 5

0,004 5

0.149 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1.943 3

4.755 5

3.776 3

9.241 2

0.477 7

0.228 2

0.002 8

0.140 1

0.293 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

5.091 2

5.898 6

12.365 0

-0.255 2

0.065 1

1.071 4

4.404 8

9.004 8

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

24.629 6

54.218 5

b = 2.201 4

2.074 0

Các độ không đảm bảo chuẩn và hiệp phương sai gắn với ab được đánh giá từ Bảng 24 sử dụng bước 8 trong 9.2.2:

Giá trị Khi-bình phương quan trắc được là  = 2,074 với 8 bậc tự do, như tính được trong Bảng 24 sử dụng 9.3. Vì  không vượt quá phân vị 95 % của , tức là 15,507, nên không có lý do để nghi ngờ sự nhất quán của mô hình đường thẳng và dữ liệu.

Các điểm dữ liệu và hàm hiệu chuẩn đường thẳng khớp như trình bày trên Hình 8. Phần dư có trọng số được minh họa trên Hình 9.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Hình 9 - Phần dư có trọng số ri tính được trong Bảng 24

10  Mô hình độ không đảm bảo và hiệp phương sai gắn với xi và yi

10.1  Khái quát

10.1.1  Điều này xem xét trường hợp 5.3.2 e), tức là trường hợp chung nhất trong đó tất cả các dữ liệu đo đều kèm theo độ không đảm bảo và hiệp phương sai. Phụ lục D cung cấp hướng dẫn để thu được các độ không đảm bảo và hiệp phương sai này.

10.1.2  Độ không đảm bảo chuẩn và hiệp phương sai bao gồm ma trận hiệp phương sai

kích thước 2m x 2m gắn với véc tơ (x1,..., xm, y1,..., ym)T của dữ liệu đo kích thước 2m x 1.

10.1.3  Trường hợp 5.3.2 e) tương ứng với mô tả bởi mô hình thống kê

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

10.1.4  Các ước lượng ab là giá trị tối thiểu hóa tổng các bình phương tổng quát

trong đó d = x - X e = y - A1 - Bx, đối với A, BXi, i = 1,..., m. Vấn đề xác định a và b trong bối cảnh này được gọi là một trong hồi quy Gauss-Markov tổng quát hóa (GGMR) [2].

10.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn và hiệp phương sai kèm theo

10.2.1  Nếu U hoàn toàn dương, thì tồn tại thừa số Cholesky L tam giác dưới kích thước 2m x 2m của U = LLT [10] (xem thêm A.4), ước lượng ab của AB có thể được tính theo quy trình lặp sử dụng thao tác ma trận-véc tơ. Nếu không có thể cần đến nhiều phương pháp số học hơn. Các thao tác này chuyển đổi tổng bình phương tng quát hóa (14) thành tổng bình phương thông thường (2) như trong 5.8.1, tức là, vấn đề trở thành bình phương tối thiểu không trọng số và không có hiệp phương sai. Quy trình lặp cũng liên quan đến các xấp xỉ , xác định các điểm  trên đường gần nhất với điểm dữ liệu đo được (xi, yi), trong đó độ gần được đo theo khoảng cách có trọng số, có tính đến thông tin về độ không đảm bảo xác định bởi U.

10.2.2  Ước lượng ab được tính các bước từ 1 đến 10 dưới đây sử dụng quy trình lặp như nêu trong 6.2.1; độ không đảm bảo chuẩn u(a)u(b) và hiệp phương sai cov(a,b) được đánh giá trong bước 11:

1  Thu được các xấp xỉ ban đầu cho tham số ;

2  Tính véc tơ kích thước 2m x 1,

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó  và  được trích từ ước lượng  hiện tại của véc tơ tham số;

3  Tính thừa số Cholesky L kích thước 2m x 2m của U = LLT [10]; xem A.4.1;

4  Giải hệ tam giác dưới

để xác định véc tơ chuyển đổi  kích thước 2m x 1 và ma trận chuyển đổi  kích thước 2m x (m + 2); xem A.4.3;

5  Từ véc tơ  có kích thước (m + 2) x 1 và ma trận  kích thước (m + 2) x (m + 2);

6  Xác định thừa số Cholesky M, ma trận tam giác dưới kích thước (m + 2) x (m + 2) trong ; xem A.4.1;

7  Giải hệ tam giác dưới Mq = -g để xác định véc tơ q kích thước (m + 2) x 1 ; xem A.4.3;

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

9  Cập nhật các xấp xỉ hiện tại cho các tham số: ;

10  Lặp lại các bước từ 2 đến 9 cho đến khi đạt được sự hội tụ. Đặt  và  (các thành phần m + 1 và m + 2 của );

11  Phân vùng M thu được trong bước 6 là

trong đó

là ma trận con tam giác dưới phía dưới bên phải kích thước 2 x 2 của M. Khi đó

CHÚ THÍCH 1: Trong bước 1, các xấp x ban đầu được cung cấp bởi , trong đó a0b0 là giá trị tham số đường thng xác định bởi bình phương tối thiểu có trọng số khớp với dữ liệu; xem 6.2.1.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 3: Trong bước 8, hiệu chính δt được cho bởi phép giải bình phương tối thiểu của công thức ma trận

cách giải được xác định bởi công thức chuẩn

CHÚ THÍCH 4: Các bước từ 5 đến 8 giải công thức chuẩn sử dụng thừa số Cholesky. Một cách tiếp cận bằng số ổn định hơn là sử dụng thừa số QR [10] của  (xem A.5.1). Quy trình mô tả trong C.2 sử dụng thừa số QR và tránh các tính toán liên quan đến số nghịch đảo L như nêu trong bước 4.

CHÚ THÍCH 5: Về mặt ma trận, ma trận hiệp phương sai gắn với các ước lượng ab

CHÚ THÍCH 6: Một cách tiếp cận tổng quát hơn và ổn định hơn bng số để giải quyết vấn đề hồi quy Gauss-Markov tổng quát được nêu trong C.2. Cách tiếp cận trên giả định rằng ma trận U là hoàn toàn dương và không thể hiện bất kỳ mối tương quan mạnh nào.

CHÚ THÍCH 7: u2(a), u2(b) và cov(a,b) trong bước 11 thu được bằng cách áp dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008) cho ab như cung cấp bởi các bước từ 1 đến 10.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

i) Các ước lượng ab được cho bởi hàm phi tuyến các dữ liệu xiyi.

ii) Các ước lượng ab có thể được coi là các thể hiện của biến ngẫu nhiên có kỳ vọng tương ứng là A*B*.

iii) Các thành phần của ma trận hiệp phương sai đối với biến ngẫu nhiên trong ii) được lấy xấp xỉ bởi u2(a), u2(b) và cov(a,b) tính được trong 10.2.2.

Các xấp xỉ trên sẽ chính xác hơn đối với dữ liệu có độ không đảm bảo kèm theo nhỏ hơn. Tuy nhiên, phương pháp ước lượng có tính chất nhất quán sau đây:

iv) Đối với dữ liệu thỏa mãn mô hình (13), khi số m điểm dữ liệu tăng thì các ước lượng ab hội tụ A*B*, tương ứng [16].

Nếu đưa ra giả định bổ sung rằng diei là thể hiện của biến ngẫu nhiên đặc trưng bởi phân bố chuẩn nhiều chiều thì khi đó các tính chất khác gắn với phương pháp ước lượng GGMR có th đưa ra:

v) Các biến ngẫu nhiên trong ii) đưc đặc trưng gần đúng bởi phân bố chuẩn hai chiều tập trung ở A*B* với ma trận hiệp phương sai quy định bởi u2(a), u2(b) và cov(a,b).

vi) Các ước lượng ab là ước lượng hợp lý lớn nhất, tương ứng với giá trị hợp lý nhất của AB, có thể làm tăng dữ liệu đo quan trắc được xiyi.

vii) Trong bối cảnh suy luận Bayes, phân bố trong phạm vi kiến thức cho AB, đã cho dữ liệu đo quan trắc xiyi, xấp xỉ phân bố chuẩn hai chiều tập trung ab với ma trận hiệp phương sai xác định bởi u2(a), u2(b) và cov(a,b).

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Nếu m > 2 thì giá trị sử dụng của mô hình có thể kiểm nghiệm một phần bằng cách sử dụng phần dư  có trọng số (tiếp theo 10.2.2):

12  Thiết lập giá trị Khi-bình phương quan trắc được  và bậc tự do n = m - 2;

13  Kiểm tra xem  có vượt quá phân vị 95 % của  không, và xem liệu có bác b mô hình đường thẳng hay không.

CHÚ THÍCH: Kiểm nghiệm Khi-bình phương dựa trên giả định là diei trong mô hình (13) là các th hiện của biến ngẫu nhiên đặc trưng bởi phân bố chuẩn nhiều chiều và dựa trên xấp x bậc một.

VÍ DỤ: Bảng 25 đưa ra bảy điểm dữ liệu đo được (xi, yi) thu được bằng cách sử dụng các mô hình đo mô tả trong D.2 và D.4.

Ma trận hiệp phương sai gắn vi yi  được rút ra bằng cách sử dụng mô hình đo (D.1) với uS = 2,0 và uR = 1,0.

Dữ liệu xi và ma trận hiệp phương sai kèm theo được rút ra bằng cách sử dụng mô hình (D.2) với z1 = 50, z2 = 100, z3 = 200, u(z1) = 0,5, u(z2) = u(z3) = 1,0 và uD,I = 0,5.

Bảng 25 - Dữ liệu thể hiện bẩy điểm đo, xiyi có ma trận hiệp phương sai kèm theo

xi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

50,4

52,3

99,0

97,8

149,9

149,7

200,4

200,1

248,5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

299,7

300,9

349,1

349,2

Ma trận hiệp phương sai Ux kích thước 7 x 7 gn với xi

Thừa số Cholesky Lx kích thước 7 x 7 của Ux = LxLxT, tính được sử dụng một trong hai thuật toán mô tả trong A.4.1, là

Ma trận hiệp phương sai Uy kích thước 7 x 7 gắn với yi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Thừa số Cholesky Ly kích thước 7 x 7 của Uy = LyLyT, tính được sử dụng một trong hai thuật toán mô tả trong A.4.1, là

Ma trận hiệp phương sai U kích thước 14 x 14 được cho bởi

CHÚ THÍCH: Đối với ví dụ này, có mối tương quan gắn với từng cặp xixj và từng cặp yiyj nhưng không có mối tương quan gắn với cặp xi yi, tức là, cov(xi, yi) = 0 đối với mọi ij.

Thừa số Cholesky L kích thước 14 x 14 được cho bởi

Bình phương tối thiểu có trọng số khớp với dữ liệu (bước 1 đến 5 trong 6.2.1) cho các xấp xỉ  = 0,270 7 và  = 1,001 1. Quy trình lặp bắt đầu với .

Bảng 26 đưa ra véc tơ ban đầu , các hiệu chính δtk = đối với lần lặp thứ k, k = 1,...,4 và ước lượng cuối .

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

50.400 0

17.253 1

1.258 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.290 4

50.572 7

99.000 0

-43.150 1

-3.214 5

-6.320 1

-0.710 1

98.568 2

149.900 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-3.960 4

-3.888 9

-0.756 4

149.608 0

200.400 0

2.967 7

-10.762 9

-0.602 4

-1.716 5

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

248.500 0

24.039 4

-11.406 4

3.237 8

-1.706 4

248.739 3

299.700 0

-22.251 0

-15.776 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-2.611 0

299.475 9

349.100 0

-20.619 2

-16.621 7

-3.380 5

-2.742 9

348.892 1

0.270 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-33.395 7

0.100 6

-5.301 9

0.342 4

1.001 1

0.011 0

0.211 3

0.007 6

0.033 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Các ước lượng tốt nhất của ABa = 0,342 4 và b = 1,001 2.

Ở lần lặp cuối ma trận M kích thước 9 x 9 là

như vậy M22(10.2.2 bước 11) là

Độ không đảm bảo chuẩn và hiệp phương sai gắn với ab được đánh giá trong bước 11 của 10.2.2 là

Giá trị Khi-bình phương quan trắc được là  = 1,772 với n = 5 bậc tự do, như tính được trong bước 12 10.3. Vì  không vượt quá phân vị 95 % của , tức là 11,070, nên không có lý do để nghi ngờ sự nhất quán của mô hình đường thẳng và dữ liệu.

11  Sử dụng hàm hiệu chuẩn

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

11.1  Dự đoán

11.1.1  Xét nội dung dưới đây được quy định, sau khi áp dụng một trong các điều từ 6 đến 10:

a) ước lượng tham số đường thẳng ab, độ không đảm bảo chuẩn u(a)u(b) và hiệp phương sai cov(a,b) gắn với ab, và

b) giá trị đo được y của Y và độ không đảm bào chuẩn kèm theo u(y).

Coi là y thu được độc lập với dữ liệu đo sử dụng để thiết lập hàm hiệu chuẩn.

11.1.2  Ước lượng x của X tương ứng với y

(15)

11.1.3  Độ không đảm bảo chuẩn u(x) gắn với x được cho bởi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 1: Công thức cho u2(x) được thiết lập bằng cách sử dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008). Xấp xỉ được dựa trên tuyến tính hóa công thức (15). c(a), c(b)c(y) đại diện cho hệ số độ nhạy.

CHÚ THÍCH 2: Công thức ma trận có thể thuận lợi cho mục đích tính toán:

CHÚ THÍCH 3: Trong trường hợp b = 0, tức là, đường thẳng khớp nhất Ià y = a, tương ứng với hàm hiệu chuẩn không chấp nhận được, thì không thể tiến hành dự đoán.

CHÚ THÍCH 4: Giá trị sử dụng của độ không đảm bảo chuẩn u(x) phụ thuộc vào sự thỏa mãn của kiểm nghiệm Khi-bình phương liên quan được nêu trong các điều từ 6 đến 10.

VÍ DỤ 1: Xét ví dụ số về bình phương tối thiểu có trọng số (WLS) có trọng số đã biết bằng nhau mô tả trong Điều 6, tham số đường thẳng khớp nhất và độ không đảm bảo chuẩn, hiệp phương sai kèm theo là

a = 1,867, b = 1,757, u(a) = 0,465, u(b) = 0,120, cov(a,b) = -0,050.

Lấy y = 10,5 là giá trị đo được b sung của Y và độ không đảm bảo chuẩn kèm theo u(y) = 0,5.

Từ 11.1.2, ước lượng giá trị x của X tương ứng với y

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Sử dụng 11.1.3, độ không đảm bảo chuẩn kèm theo u(x) được cho bởi

VÍ DỤ 2: Xét ví dụ số về bình phương tối thiểu có trọng số (WLS) có trọng số không bằng nhau đã biết mô tả Điều 6, tham số đường thẳng khớp nhất và độ không đảm bảo chuẩn, hiệp phương sai kèm theo là

a = 0,885, b = 2,057, u(a) = 0,530, u(b) = 0,178, cov(a,b) = -0,082.

Lấy y = 10,5 là giá trị đo được bổ sung của Y và độ không đảm bảo chun kèm theo u(y) = 1,0.

Từ 11.1.2, ưc lượng giá trị x của X tương ứng với y là

x = (10,5 - 0,885)/2,057 = 4,674.

Sử dụng 11.1.3, độ không đảm bảo chuẩn kèm theo u(x) được cho bởi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

11.2  Đánh giá tiếp

Xét nội dung dưới đây được quy định, sau khi áp dụng một trong các điều từ 6 đến 10:

a) ước lượng tham số đường thẳng ab, độ không đảm bảo chuẩn u(a)u(b) và hiệp phương sai cov(a,b) gắn với ab, và,

b) giá trị đo được x của X và độ không đảm bảo chuẩn kèm theo u(x).

Coi là x thu được độc lập với dữ liệu đo sử dụng để thiết lập hàm hiệu chuẩn.

11.2.1  Ước lượng y của Y tương ứng với x

y = a + bx

(16)

11.2.2  Độ không đảm bảo chuẩn u(y) gắn với y được cho bởi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH 1: Công thức cho u2(y) được thiết lập bằng cách sử dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008). Xấp x được dựa trên tuyến tính hóa công thức (16). c(a), c(b)c(y) đại diện cho hệ số độ nhạy.

CHÚ THÍCH 2: Công thức ma trận có th thuận lợi cho mục đích tính toán:

CHÚ THÍCH 3: Giá trị sử dụng của độ không đảm bảo chuẩn u(y) phụ thuộc vào sự thỏa mãn của kiểm nghiệm Khi-bình phương liên quan được nêu trong các điều từ 6 đến 10.

VÍ DỤ: Xét ví dụ số về bình phương tối thiểu có trọng số (WLS) có trọng số đã biết bằng nhau mô tả trong Điều 6, tham số đường thẳng khớp nhất và độ không đảm bảo chuẩn, hiệp phương sai kèm theo là

a = 1,867, b = 1,757, u(a) = 0,465, u(b) = 0,120, cov(a,b) = -0,050.

Lấy x = 3,5 là giá trị đo được b sung của X và độ không đảm bo chuẩn kèm theo u(x) = 0,2, giả định rằng cov(x, a) = cov(x, b) = 0, tức là, không có mối tương quan gắn với x a, xb.

Từ 11.2.1, ước lượng giá trị y của Y tương ứng với x

y = 1,867 + (1,757)(3,5) = 8,017.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

như vậy u(y) = 0,406.

 

Phụ lục A

(tham khảo)

Phép tính ma trận

A.1  Khái quát

Phụ lục này mô tả các phép tính ma trận được sử dụng trong tiêu chuẩn này.

A.2  Phép tính cơ bản

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

A.2.1  Phép nhân ma trận-véc tơ

Tích ma trận-véc tơ Ad là véc tơ e kích thước m x 1 với phần tử thứ i ei xác định bởi

A.2.2  Phép nhân ma trận-ma trận

Tích ma trận AB là ma trận kích thước m x k có cột thứ j là tích của A và cột thứ j của B.

A.2.3  Chuyển vị ma trận

Chuyển vị AT của ma trận A là ma trận kích thước n x m với phần tử A(j, i) = aij trong hàng thứ j và cột thứ i.

A.2.4  Ma trận đơn vị

Ma trận đơn vị bậc m là ma trận I kích thước m x m sao cho I(j, j) = 1, j = 1,…, m, và tất cả các phần tử khác bằng 0.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Nghịch đảo của C, nếu có, được ký hiệu là C-1 và là ma trận kích thước m x m sao cho

Chuyển vị của C-1 bằng nghịch đo của CT và được ký hiệu là C-T.

A.3  Các định nghĩa cơ bản

Trong các định nghĩa sau đây, C là ma trận (vuông) kích thước m x m với phần tử C(i, j) = cij trong hàng thứ i và cột thứ j.

A.3.1  Ma trận đối xứng

Ma trận C là đối xứng nếu cij = cij, i = 1,..., m, j = 1,..., m, tức là, C = CT.

A.3.2  Ma trận khả nghịch

Ma trận C là khả nghịch nếu tồn tại nghịch đảo C-1 của nó (xem A.2.5).

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Ma trận C là ma trận tam giác dưới nếu cij = 0, i < j, và ma trận tam giác trên nếu cij = 0, i > j.

A.3.4  Ma trận trực giao

Ma trận C là trực giao nếu CTC = I.

A.4  Thừa số Cholesky

Thừa số Cholesky của ma trận hoàn toàn dương đối xứng U kích thước m x m là ma trận tam giác dưới L kích thước m x m sao cho U = LLT [10].

A.4.1  Thuật toán thừa số Cholesky

A.4.1.1  Thuật toán sau đây tính ma trận tam giác dưới L sao cho U = LLT.

CHÚ THÍCH: Đ ghi đè các phần tử tam giác dưới U(i, j), i j của U với thừa số Cholesky của nó, ch thực hiện các bước trong giai đoạn Thừa số hóa của thuật toán trong A.4.1.1, sử dụng U thay cho L.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

CHÚ THÍCH: Để ghi đè các phần tử tam giác dưới U(i, j), i j của U với thừa số Cholesky của nó, chỉ thực hiện các bước trong giai đoạn Thừa số hóa của thuật toán trong A.4.1.2, sử dụng U thay cho L.

A.4.2  Giải thích phép thừa số hóa Cholesky của ma trận hiệp phương sai

A.4.2.1  Giả định Ei, i = 1,…, m là m biến ngẫu nhiên độc lập, mỗi biến có kỳ vọng bằng 0 và phương sai bằng 1, lấy ei là thể hiện của Ei. Lấy

Khi đó, . S phụ thuộc chung của y1y2 vào ei nghĩa là y1y2 có mối tương quan, với hiệp phương sai cov(y1, y2) = l11l21. Tiếp tục, giả định

A.4.2.2  Theo thuật ngữ ma trận, y = Le, với tam giác dưới L. Sự phụ thuộc chung của y1y3 vào e1 nghĩa là có mối tương quan gắn với y1 y3. Tương tự, sự phụ thuộc chung của y2y3 vào e1e2 nghĩa là có mối tương quan gắn với y2y3, và cứ tiếp như vậy.

A.4.2.3  Cho trước ma trận hiệp phương sai U gắn với dữ liệu yi, thừa số Cholesky U = LLT tính các hệ số Iij sao cho ma trận hiệp phương sai có thể giải thích bằng cách giả định rằng yi được định nghĩa trong A.4.2.1 như thể hiện của các tổ hợp tuyến tính xác định bởi Iij của các biến ngẫu nhiên độc lập Ei. Trên thực tế, các ma trận hiệp phương sai thường được xác định về các thừa số U = BBTU đã cho có vô cùng nhiều thừa số B có thể sử dụng để xây dựng U. Thừa s Cholesky, trong đó các tổ hợp tuyến tính được đại diện bởi ma trận tam giác dưới, là duy nhất đến dấu trị số của các cột của L.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

A.4.3.1  Nếu L là ma trận tam giác dưới kích thước m x m sao cho L(j, j) ≠ 0, j = 1,…, m x là véc tơ kích thước m x 1, thuật toán sau đây tính véc tơ y, trong đó y sao cho Ly = x, tức là, y = L-1x.

CHÚ THÍCH: Để ghi đè véc tơ x với lời giải y, ch thực hiện các bước trong giai đoạn Giải của thuật toán trong A.4.3.1, sử dụng x thay cho y.

A.4.3.2  Thuật toán trong A.4.3.1 có thể áp dụng để giải phương trình ma trận LY = X bằng cách áp dụng liên tiếp cho từng cột của X. Về toán học, lời giải được cho bởi Y = L-1X.

A.4.4  Giải hệ tam giác trên

A.4.4.1  Phép giải hệ tam giác trên có thể được xác định theo chuyển vị ma trận tam giác dưới. Nếu L là ma trận tam giác dưới kích thước m x m sao cho L(j, j) ≠ 0, j = 1,..., m x là véc tơ kích thước m x 1, thuật toán sau đây tính véc tơ y, trong đó y sao cho LTj = x, tức Ià, y = L-Tx.

CHÚ THÍCH: Để ghi đè véc tơ x với lời giải y, chỉ thực hiện các bước trong giai đoạn Giải của thuật toán trong A.4.4.1, sử dụng x thay cho y.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

A.5  Thừa số trực giao

Các ma trận trực giao là t hợp các phép quay và phép đối xứng và có tính chất nhân trước véc tơ với ma trận trực giao mà không làm thay đổi độ lớn của véc tơ đó (căn bậc hai tổng các bình phương các phần t của nó). Các cột của ma trận trực giao có thể được coi như xác định hệ trục vuông góc. Điều quan trọng của kỹ thuật thừa số trực giao là nó cho phép giải các phương trình ma trận theo cách số học ổn định. Các thuật toán để tính thừa số trực giao được mô tả trong các tài liệu tham khảo [1,10, 20].

A.5.1  Thừa số QR

Thừa số QR của ma trận A kích thước m x n, với m ≥ n, có thể được viết là

trong đó Q = [Q1 Q2] là ma trận trực giao kích thước m x m, Q1 là ma trận bao gồm n cột đầu tiên của Q, với , và R1 là ma trận tam giác trên kích thước n x n.

CHÙ THÍCH: Cũng có thể thu được thừa số QR của ma trận A kích thước m x n, với m < n. Trong tiêu chuẩn này, vì tất cả các ma trận mà tính toán thừa số QR được yêu cầu có m n, nên thừa số không được cung cấp.

A.5.2  Thừa số RQ

A.5.2.1  Thừa số RQ của ma trận B kích thước m x n, với m n có thể được viết là

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó Z là ma trận trực giao T2 là ma trận tam giác trên.

A.5.2.2  Thừa số RQ của ma trận B kích thước m x n, với m < n có thể được viết là

trong đó Z là ma trận trực giao T2 là ma trận tam giác trên.

 

Phụ lục B

(tham khảo)

Áp dụng thuật toán Gauss-Newton cho hồi quy khoảng cách tổng quát

B.1  Phụ lục này rút ra các thuật toán trong 7.2.1 và 8.2.1 bằng cách sử dụng thuật toán Gauss-Newton.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

B.3  Lấy  là xấp xỉ của các tham số giải a

tương ứng là véc tơ kích thước m x 1 của các giá trị hàm số và ma trận Jacobi kích thước m x n của đạo hàm riêng bậc một đối với các tham số, đánh giá tại  xấp xỉ với các tham số.

B.4  Lấy p giải

Khi đó, ước lượng mới của các tham số giải được cho bởi .

B.5  Đối với các thuật toán trong 7.2.1 và 8.2.1, a = (A, B)T và hàm ƒi(A) là thước đo khoảng cách tổng quát tính từ điểm dữ liệu thứ i(xi, yi) đến đường y - A + Bx.

B.6  Lấy Ui là ma trận hiệp phương sai gắn với điểm dữ liệu thứ i:

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

và lấy , gọi là điểm chân thứ i, giải

(B.2)

hàm của AB.

B.7  Nếu  được xác định bằng

tức là,  đánh giá tại giải , thì khi đó các giá trị của A và B làm tối thiểu

xác định đường hồi quy khoảng cách tổng quát khớp nhất. Áp dụng thuật toán Gauss-Newton đòi hỏi xác định đạo hàm riêng bậc một của ƒi(A, B) đối với A và B đ hình thành ma trận Jacobi J.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

(B.3)

khi đó, biểu thị giá trị hàm và đạo hàm đối với Ui, A, B, xi, yi,

(B.4)

và

(B.5)

B.9  Giải  của vấn đề điểm chân (B.2) được cho bởi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

(B.6)

CHÚ THÍCH: Biểu thức (B.3), (B.4), (B.5) và (B.6) được xác định đối với Ui chứ không phải nghịch đảo của ,. Không có yêu cầu Ui khả nghịch nhưng nTUin phải khác không.

B.10  Các thuật toán trong 7.2.1 và 8.2.1 áp dụng thuật toán Gauss-Newton sử dụng các biểu thức rõ ràng cho ƒi(A,B),  và . Giải cho bước mới p trong biểu thức (B.1) được hình thành như vn đề xác định đường thẳng khớp nhất bình phương tối thiểu có trọng số (xem 6.2.1 các bước từ 1 đến 5) đối với dữ liệu chuyển đổi rút ra từ dữ liệu đo (xi,yi), ma trận hiệp phương sai kèm theo Ui và các gần đúng hiện tại với AB.

 

Phụ lục C

(tham khảo)

Cách tiếp cận thừa số trực giao để giải quyết bài toán Gauss-Markove tổng quát

C.1  Khái quát

Thuật toán lặp mô tả trong 10.2.2 giả định rằng ma trận hiệp phương sai U kích thước 2m x 2m là xác định dương và do đó là khả nghịch. Đặc biệt, tính khả nghịch đòi hỏi mọi u(xi) > 0 và u(yi) > 0. Trong phụ lục này thuật toán tổng quát được mô tả là phù hợp cho mọi ma trận hiệp phương sai U (bán xác định dương đối xứng). Tất cả những gì được yêu cầu là ma trận hiệp phương sai có thể thừa số hóa là U = BBT, trong đó B là ma trận kích thước 2m x p (p m). Thông thường, ma trận hiệp phương sai được bắt nguồn từ thừa số hóa như vậy. Nếu U là khả nghịch thì B có thể là thừa số Cholesky của nó. Thuật toán tiến hành tương tự như mô tả trong 10.2.2 và yêu cầu tính toán số dư ƒ và ma trận Jacobi J, nhưng số gia δt được xác định bằng cách sử dụng hai thừa số trực giao, về toán học, δt tối thiểu hóa

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

C.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn, hiệp phương sai kèm theo

Các ước lượng ab được tính như trong các bước từ 1 đến 9 dưới đây; độ không đảm bảo chuẩn u(a) và u(b) được đánh giá trong bước 10:

1  Thu được các xấp xỉ ban đầu  cho tham số;

2  Tính véc tơ kích thước 2m x 1,

và ma trận (Jacobi) kích thước 2m x (m + 2),

trong đó được trích từ ước lưng  hiện tại của véc tơ tham số;

3  Xác định thừa số QR của J:

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó Q là ma trận trực giao, kích thước 2m x 2 mR1 là ma trận tam giác trên kích thước (m + 2) x (m + 2); xem A.5.1;

4  Hình thành tích ma trận QTB và xác định thừa số RQ của nó

QTB = TZ,

trong đó T ma trận, kích thước 2m x pZ là ma trận trực giao kích thước p x p; xem A.5.2;

5  Đặt và phân vùng  và T.

trong đó  là véc tơ kích thước (m + 2) x 1,  là véc tơ kích thước (m - 2) x 1, T11 là ma trận kích thước (m + 2) x (p - m + 2), T12 là ma trận kích thước (m + 2) x (m - 2) và T22 là ma trận tam giác trên kích thước (m - 2) x m - 2);

6  Giải hệ tam giác trên  để xác định véc tơ  kích thước (m - 2) x 1; xem A.4.4;

7  Giải hệ tam tác trên  đ xác định số gia , xem A.4.4;

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

9  Lặp lại các bước từ 2 đến 8 cho đến khi đạt được sự hội tụ. Đặt  và  (các thành phần m + 1 và m + 2 của );

10  Lấy Ra là ma trận dưới bên phải kích thước 2 x 2 của R1Ta là ma trận dưới bên phải kích thước 2 x 2 của T11. Giải hệ tam giác trên

đối với ma trận tam giác trên Ka kích thước 2 x 2 (xem A.4.4) và đặt Ua = KaKaT. Khi đó

CHÚ THÍCH 1: Cách tiếp cận mô tả trong C.2 đại diện cho cách giải tổng quát nhất đ xác định hàm hiệu chuẩn tuyến tính sử dụng phương pháp bình phương tối thiểu. Tt cả các cách tiếp cận khác mô tả trong tiêu chuẩn này có th được giải như trường hợp đặc biệt.

CHÚ THÍCH 2: Các bước 1, 2, 8 và 9 trong C.2 tương ứng giống với các bước 1, 2, 9 và 10 trong 10.2.2.

C.3  Xác nhận giá trị sử dụng của mô hình

Nếu m > 2 thì giá trị sử dụng của mô hình có thể kiểm nghiệm một phần bằng cách sử dụng thành phần véc tơ  (tiếp theo C.2):

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

12  Kiểm tra xem  vượt quá phân vị 95 % của  không, và xem liệu có bác bỏ mô hình đường thẳng hay không.

CHÚ THÍCH: Kiểm nghiệm Khi-bình phương dựa trên giả định là diei trong mô hình (13) là các thể hiện của biến ngẫu nhiên đặc trưng bởi phân bố chuẩn nhiều chiều và dựa trên xp xỉ bậc một. Với gi định này, véc tơ  kích thước (m - 2) x 1 gắn với phân bố Gauss nhiều chiều có ma trận hiệp phương sai bằng với ma trận đơn vị kích thước (m - 2) x (m - 2) sao cho  gắn với phân bố x2 với m - 2 bậc tự do.

VÍ DỤ 1: Cách tiếp cận thừa số QR có thể áp dụng cho ví dụ số mô tả trong Điều 10.

Ma trận hiệp phương sai Ux sinh ra dưới dạng thừa số (xem D.4) là

Ux = BxBxT,

trong đó

Ma trận hiệp phương sai Uy = ByByT cũng có dạng thừa số với

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Đối với ví dụ này, thuật toán trong C.2 tương đương về toán học với thuật toán trong 10.2.2 và hai cách tiếp cận đều cho kết quả số rt giống nhau.

VÍ DỤ 2: Bảng C.1 cho 7 điểm dữ liệu đo được (xi, yi) thu được bằng cách sử dụng mô hình đo mô tả trong D.2 và D.5.

Ma trận hiệp phương sai gắn với yi được rút ra bằng cách sử dụng mô hình (D.1) với uS = 2,0 và uR = 1,0 và giống như trong Ví dụ 1 Phụ lục C.

Dữ liệu xi và ma trận hiệp phương sai kèm theo được rút ra bằng cách sử dụng mô hình đo (D.3) với z1 = 50, z2 = 100, z3 = 200, u(z1) = 0,5 và u(z2) = u(z3) = 1,0, sao cho

Ux = BxBxT,

trong đó

Ma trận hiệp phương sai hoàn chỉnh U kích thước 14 x 14 được thừa số hóa thành U = BBT, trong đó B là ma trận kích thước 14 x 11

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Đối với ví dụ này, không thể áp dụng thuật toán trong 10.2.2 vì U không xác định dương. Thay vào đó có thể sử dụng thuật toán mô tả trong C.2.

Bảng C.2 đưa ra véc tơ ban đầu , các hiệu chính  = đối với lần lặp th k, k = 1,...,5 và ước lượng cuối .

Bảng C.1 - Dữ liệu thể hiện 7 điểm đo, xiyi có ma trận hiệp phương sai kèm theo

xi

yi

50,5

47,1

99,7

98,4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

153,7

199,5

194,0

249,9

251,9

299,2

297,5

349,7

349,0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

50.500 0

99.700 0

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

199.500 0

249.900 0

299.200 0

349.700 0

30.822 9

55.831 3

86.654 2

-59.071 1

-28.248 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

27.583 1

3.187 4

-13.836 5

-10.649 1

-45,597 6

-45.410 2

-62.434 1

-59.246 7

23.295 7

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

49.409 3

-49.584 9

-26.289 1

-23.471 3

0,175 5

8.112 4

-0.206 3

7.906 1

-44.790 4

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-44.996 7

-36.884 3

16.423 1

15.677 0

32.100 2

-43.047 0

-26.623 7

-27.369 8

-10.946 8

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

100.257 0

151.065 5

198.904 4

249.613 0

299.161 3

349.969 9

-1.852 8

1.004 2

-50.620 3

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-140.085 6

0.857 1

-63.931 6

0.321 7

-100.943 2

0.610 8

-68,134 5

0.372 2

-2.373 1

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

Phụ lục D

(tham khảo)

Quy định về độ không đảm bảo và hiệp phương sai gắn với x-giá trị và y-giá trị đo được

D.1  Khái quát

Phụ lục này chỉ ra cách thức để thu được độ không đảm bảo và hiệp phương sai gắn với giá trị đo đáp ứng và kích thích đo được. Cách tiếp cận dựa trên việc sử dụng mô hình đo gồm các quá trình làm cơ s cho việc xác định dữ liệu đáp ứng và kích thích, và áp dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008). Các ví dụ minh họa được sử dụng cho mục đích này.

D.2  Dữ liệu đáp ứng 1

D.2.1  Khái quát

D.2.1.1  Giả định đại lượng Y đại diện đáp ứng phương tiện đo có thể được biểu thị bằng mô hình đo

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

(D.1)

trong đó Y0 là đại lượng thể hiện bởi đáp ứng chỉ thị và E là đại lượng đại diện cho ảnh hưng hệ thống. Giả định rằng hiểu biết về Y0 được mã hóa bằng phân bố có độ lệch chuần uR. Phân bố này thường dựa trên phân tích số lượng chỉ thị lặp lại của Y. Y0 được ước lượng bằng trung bình của các chỉ thị này và uR là độ không đảm bảo chuẩn gắn với ước lượng này. Giả định rằng hiểu biết về E sao cho E có kỳ vọng 0 (tức là, mọi hiệu chính cần thiết đều được áp dụng) và phương sai  (thu được từ hiểu biết phụ thuộc vào tính chất cụ thể của phương tiện đo).

D.2.1.2  Tiếp sau biểu thức (D.1), bằng việc áp dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008), độ không đảm bảo chuẩn u(yi) gắn với giá trị đo được yi của Y được cho bởi

Ngoài ra, hiệp phương sai gắn với giá trị đo được yi yj của Y là

D.2.1.3  Do đó, ma trận hiệp phương sai trong trường hợp này là

D.2.2  Mô hình đo dùng cho độ không đảm bảo và hiệp phương sai gắn với yj

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó Y0,i là đại lượng thể hiện bởi đáp ứng chỉ thị thứ i và các đại lượng E1 E2 đại diện cho ảnh hưởng hệ thống. Giả định rằng hiểu biết về Y0,i sao cho Y0,i có phương sai  và hiểu biết về Ek sao cho Ek có kỳ vọng bằng 0 và phương sai , đối với k = 1, 2.

D.2.2.2  Độ không đảm bảo chuẩn u(yi) gắn với giá trị đo được yi của Yi được cho bởi

Hiệp phương sai gắn với giá trị đo đưc yiyj

D.2.2.3  Ma trận hiệp phương sai trong trường hợp này là

D.3  Dữ liệu đáp ng 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Y= Y0(1 + D).

D.3.2  Xử lý tương tự như D.2 ngoại trừ là lúc này, sử dụng uD để ký hiệu cho độ không đảm bảo chuẩn tương đối gắn với ước lượng của Y0,

D.3.3  Ma trận hiệp phương sai trong trường hợp này là

D.4  Dữ liệu kích thích 1

D.4.1  Dữ liệu sử dụng trong ví dụ ở Điều 10 được rút ra từ mô hình đo dưới đây, thúc đẩy bởi thực tiễn trong đo lường khối lượng trong đó số các khối lượng được sử dụng để tạo ra nhiều giá trị hiệu chuẩn xi. Dữ liệu kích thích xi là thể hiện của các biến ngẫu nhiên Xj, i = 1,…, 7, được xác định theo các biến ngẫu nhiên Zk, k = 1,2,3 và Di, i = 1,…,7:

Các biến ngẫu nhiên Zk, k = 1,2,3, có các kỳ vọng zk và phương sai u2(zk), trong khi Di có kỳ vọng bằng 0 và phương sai . (Trong hiệu chuẩn khối lượng, giá trị zk là giá trị đã hiệu chuẩn đối với các khối lượng và u(zk) là độ không đm bảo kèm theo.)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

là ma trận độ nhạy kích thước 7 x 10, trong đó CD = I là ma trận đơn vị kích thước 7 x 7 và

D.4.3  Lấy SD là ma trận đường chéo kích thước 7 x 7 với thành phần đường chéo SD(i,i) = uD,i, i = 1,...,7, và Sz là ma trận đường chéo kích thước 3 x 3 với thành phần đường chéo Sz(k,k) = u(zk), k = 1,2,3. Đặt

D.4.4  Khi đó, ước lượng tốt nhất của X được cho bằng x = Czz kích thước 7 x 1 và ma trận hiệp phương sai kèm theo kích thước 7 x 7 được cho bởi

Số hạng  là phương sai đóng góp phát sinh từ Di trong khi số hạng thứ hai là đóng góp từ Zk.

D.5  Dữ liệu kích thích 2

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Các biến ngẫu nhiên Zk, có các kỳ vọng zk và phương sai u2(zk), k = 1, 2, 3.

D.5.2  Độ không đảm bảo u(zk) được lan truyền thông qua mô hình đo trong D.5.1 sang những độ không đảm bảo gắn với các ước lượng xi của Xi sử dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008). Theo sau ký hiệu của D.4, ước lượng tốt nhất gắn với Xi được cho bởi x = Czz kích thước 7 x 1 và ma trận hiệp phương sai kèm theo kích thước 7 x 7 được cho bởi

Trong trường hợp này Ux là không nghịch đảo được.

D.6  Dữ liệu kích thích và đáp ứng

D.6.1  Sự tương quan, tức là, hiệp phương sai khác không, gắn với dữ liệu đo xiyi phát sinh thông qua sự có mặt của các ảnh hưởng chung đến cả hai.

D.6.2  Giả định X Y có thể được biểu thị bằng mô hình đo

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

trong đó X0, Y0T là các biến ngẫu nhiên độc lập với kỳ vọng x0, y0 và 0, còn phương sai tương ứng là u2(x0), u2(y0) và u2(t).

D.6.3  Tiếp sau biểu thức (D.4), bằng cách áp dụng định luật lan truyền độ không đảm bảo trong TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008), độ không đảm bảo chuẩn u(xi) và u(yi) gắn với giá trị đo được xi của X yi của Y được cho bởi

Ngoài ra, hiệp phương sai gắn với xiyi

D.6.4  Nếu thay X Y có thể biểu thị bằng mô hình đo

hiệp phương sai gắn với xiyi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Phụ lục E

(tham khảo)

Độ không đảm bảo đã biết theo hệ số thang đo

E.1  Phụ lục này mô tả phương pháp đ đánh giá độ không đảm bảo gắn với dữ liệu đo trong trường hợp độ không đảm bảo chỉ được biết theo hệ số thang đo.

E.2  Tiêu chuẩn này giả định chung rằng độ không đảm bảo gắn với dữ liệu đo được cung cấp. Đây là tình huống khi các đại lượng (các biến) liên quan đã được mô tả đặc trưng theo các nguyên tắc của TCVN 9595-3:2013 (ISO/IEC Guide 98-3:2008) và ISO/IEC Guide 98-3:2008/Suppl. 1:2008 [13] về phân bố xác suất. Giá trị đo được cho bằng kỳ vọng của biến và phương sai gắn với nó được cho bằng phương sai của biến.

CHÚ THÍCH: Đặc biệt, liên quan đến giá trị đo được y như thể hiện của biến đặc trưng bởi phân bố t có tham số t lệ s và bậc tự do n (n >2), độ không đảm bảo chuẩn gắn với y được cho bởi u(y) = [n/(n - 2)1/2 s, độ lệch chuẩn của phân bố đó.

E.3  Vì hàm hiệu chuẩn thường được sử dụng trong thực tiễn đo nên việc đánh giá độ không đảm bảo gắn với dữ liệu hiệu chuẩn cần hoàn chỉnh và chặt chẽ nhất có th. Các ước lượng của tham số hàm hiệu chuẩn và độ không đảm bảo kèm theo khi đó có thể được sử dụng tin cậy. Sự xem xét này giảm nhẹ đề cập trong tiêu chuẩn này là trường hợp trong đó độ không đảm bảo được biết theo hằng số nhân. Trường hợp phổ biến nhất [9, trang 30] là trường hợp trong đó tin tưng rằng giá trị y đo được có độ không đảm bảo về cơ bản giống nhau nhưng độ không đảm bảo chuẩn chung σ của chúng là chưa biết. (Đây là ví dụ về trường hợp tổng quát hơn trong đó ma trận hiệp phương sai , trong đó  đã cho, nhưng σ chưa biết.) Nếu m > 2 thì có thể cung cấp ước lượng  của σ trên cơ sở sự phân tán của các điểm dữ liệu quanh đường thẳng khớp. Ước lượng này được biết là ước lượng hậu nghiệm của σ, khái niệm hậu nghiệm đề cập đến thực tế là chỉ có thể xác định sau khi đã thu được đường thẳng khớp nhất cho dữ liệu.

E.4  Ước lượng hậu nghiệm được xác định bằng cách sử dụng cùng các khái niệm như đã sử dụng trong việc xác nhận giá trị sử dụng của mô hình. Bằng cách đặt thêm giả định rằng dữ liệu đầu vào là thể hiện của biến phân bố chuẩn nhiều chiều, ước lượng hậu nghiệm σ được chọn sao cho  bằng với m - 2, kỳ vọng của phân bố Khi-bình phương với m - 2 bậc tự do. Trong trường hợp này không thể tiến hành việc xác nhận giá trị sử dụng của mô hình và dữ liệu vì ước lượng hậu nghiêm đã được chọn sao cho chuẩn mực hiệu lực được tự động đáp ứng.

E.5  Do đó, phương pháp này chỉ nên được sử dụng với chú ý đặc biệt. Ví dụ, nếu đồ thị dữ liệu chỉ ra rằng hàm hiệu chuẩn đường thẳng không thích hợp thì không nên sử dụng phương pháp này.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

E.7  Trong trường hợp dữ liệu có thể được coi là thể hiện của phân bố chuẩn nhiều chiều với ma trận hiệp phương sai U đã biết, phân bố hai chiều trong E.6 là phân bố chuẩn với ma trận hiệp phương sai Ua có các thành phần u2(a), u2(b) và hiệp phương sai cov(a,b) như trong biểu thức (1).

E.8  Đối với trường hợp trong E.4 trong đó phân bố chuẩn nhiều chiều có ma trận hiệp phương sai u = σ2 U0, trong đó U0 đã biết và σ chưa biết, thì U0 được sử dụng thay cho U trong tính toán. Ma trận hiệp phương sai

gắn với các ước lượng của tham số hiệu chuẩn đường thẳng có thể được tính. Nếu m > 2 thì có thể sử dụng giá trị Khi-bình phương quan trắc được (xem 6.3) để cung cấp ước lượng hậu nghiệm của hệ số thang đo gắn với độ không đảm bảo đầu vào. Lấy  được tính như trong bước 8 của 6.3, và đặt

E.9  Ma trận hiệp phương sai điều chỉnh tỷ lệ

khi đó được cho bởi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

E.10  Các ước lượng (E.1) được dựa trên việc khớp với một số lượng hữu hạn m điểm dữ liệu và đối với m nhỏ sẽ ước lượng thấp phương sai của phân bố đối với tham số khớp. Với m > 4, ước lượng tốt hơn được xác định [19, chương 8] bằng cách sử dụng

CHÚ THÍCH: Với gi định về tính chuẩn, các ước lượng tham số khi đó được gắn với phân bố t hai chiều với ma trận tỷ lệ  và m - 2 bậc tự do. Với m > 4, ma trận hiệp phương sai của phân bố đó được cho bởi

(E.2)

trong đó hệ số tăng (m - 2)/(m - 4) tính đến thực tế là σ được ước lượng chứ không phải được biết trước.

VÍ DỤ: (CHƯA BIẾT TRỌNG SỐ) Trong ví dụ này, xi được lấy là chính xác và yi có độ không đảm bảo chuẩn bằng nhau nhưng chưa biết, còn ước lượng hậu nghiệm của độ không đảm bảo gắn với tham số khớp được đánh giá từ các phần dư của việc khớp. Việc khớp được xác định bằng cách lấy trọng số bằng 1 (nghĩa là độ không đảm bảo chuẩn u(yi) trên danh nghĩa cũng bằng 1). Dữ liệu được cho trong Bảng E.1.

Bảng E.1 - Dữ liệu đại diện cho 6 điểm đo, với trọng số đặt bằng 1

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

yi

u(yi)

1,000

3,014

1

2,000

5,225

1

3,000

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

1

4,000

9,061

1

5,000

11,201

1

6,000

12,762

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Tham số đường thẳng khớp nhất được tính như trong Bng E.2. Từ bảng này, g0 = 21,000/6,000 = 3,500, h0 = 48,267/6,000 = 8,044, b = 34,363/17,500 = 1,964 và a = 8,044 - (1,964)(3,500) = 1,172.

Bảng E.2 - Bng tính toán đi kèm với dữ liệu trong Bảng E.1

w1

gi

hi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

ri

 

 

 

 

3.500

8.044

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

a = 1.172

 

1.000

1.000

1.000

3.014

-2.000

-5.031

6.250

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-0.122

0.015

1.000

1.000

2.000

5.225

-1.500

-2.819

2.250

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.126

0.016

1.000

1.000

3.000

7.004

-0.500

-1.040

0.250

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-0.059

0.003

1.000

1.000

4.000

9.061

0.500

1.017

0.250

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.035

0.001

1.000

1.000

5.000

11.201

1.500

3.157

2.250

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

0.211

0.045

1.000

1.000

6.000

12.762

2.500

4.718

6.250

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

-0.191

0.037

 

6.000

21.000

48.267

 

 

17.300

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

b = 1.964

0.116

Dữ liệu và hàm hiệu chuẩn đường thẳng khớp được vẽ trên Hình E.1. Số dư có trọng số được minh họa trên Hình E.2. Do u(yi) là tùy ý nên đưa ra giá trị 1, trong trường hp này các vạch độ không đảm bảo vượt quá nhiều các số dư về độ lớn.

Hình E.1 - Dữ liệu trong Bảng E.1 và hàm hiệu chuẩn đường thẳng khớp thu được trong Bảng E.2

Nếu biết trước là u(yi) = 1, i = 1,…, m, thì khi đó độ không đảm bo gắn với các tham số khớp sẽ được tính toán từ thông tin trong Bảng E.2:

Hình E.2 - Số dư có trọng số tính bằng cách sử dụng hàm hiệu chuẩn đường thẳng khớp thu được trong Bảng E.2.

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Giá trị này của  đại diện cho ước lượng của độ không đảm bảo chuẩn u(yi) gắn với yi dựa trên giá trị Khi-bình phương quan trắc được. Với ước lượng hậu nghiệm này, các tính toán có thể lặp lại với u(yi) = 0,171. Các ước lượng cho ab sẽ không thay đổi nhưng giá trị Khi-bình phương quan trắc và độ không đảm bảo sẽ cỏ tỷ lệ như sau:

như vậy , kỳ vọng của phân bố Khi-bình phương với 4 bậc tự do. Sử dụng công thức (E.1),

Các thành phần của  là những thành phần sẽ được đánh giá nếu đã biết, theo tiên nghiệm, là . Tuy nhiên,  là ước lượng của độ không đảm bảo chuẩn gắn với yi. Với m > 4, tỷ lệ tăng của (m - 2)/(m - 4) có thể được kết hợp vào ma trận hiệp phương sai để tính đến độ không đảm bảo b sung phát sinh từ thực tế  là ước lượng rút ra từ m điểm dữ liệu. Sử dụng công thức (E.2),

như vậy

 

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

(tham khảo)

Triển khai phần mềm các thuật toán được mô tả

F.1  Việc triển khai phần mềm các thuật toán mô tả trong tiêu chuẩn này đ xác định và sử dụng hàm hiệu chuẩn đường thẳng đã được Phòng thí nghiệm vật lý quốc gia (NPL - National Physical Laboratory) Anh xây dựng. Phầm mềm có sẵn dưới dạng thư mục nén ZIP trên trang tin điện tử của NPL tại địa chỉ www.npl.co.uk/mathematics-scienctific-computing/software-support-for-metrology/software-downloads-(ssfm) và Tổ chức Tiêu chuẩn hóa Quốc tế tại standards.iso.org/iso/ts/28037/.

F.2  Phần mềm, phát triển trong ngôn ngữ lập trình MATLAB [18], được cung cấp dưới dạng tệp M và tệp html xuất bản bằng MATLAB version 7.10.0 (R2010a). Đối với người sử dụng MATLAB, tệp M có thể chạy trực tiếp và cũng có thể sửa đi để chạy các thuật toán đối với các dữ liệu đo khác nhau. Đối với người sử dụng không tiếp cận được MATLAB, tốt nhất là đọc phần mềm ở dạng tệp html được cung cấp. Có thể sử dụng phần mềm làm cơ sở cho việc chuẩn bị triển khai các thuật toán trong các ngôn ngữ lập trình khác. Trong các tệp, các lệnh được thực hiện cho một số hàm MATLAB cũng bao gồm phần mềm. Ví dụ, hàm algm_gdr1_steps_2_to_5 triển khai các bước 2 đến 5 của quy trình tính toán cho trường hợp 5.3.2 b) (độ không đảm bảo gắn với giá trị đo được xiyi và tất cả các hiệp phương sai gắn với dữ liệu được coi là không đáng kể) quy định trong 7.2.1. Ngoài ra, một số ứng dụng các chức năng tích hợp của MATLAB được đưa ra, ví dụ như để thu được thừa số Cholesky của ma trận. Tập lệnh MATLAB (có đuôi '.m') và tệp html (‘.html’) được cung cấp như sau:

- TS28037_WLS1 chạy ví dụ số về bình phương tối thiểu có trọng số (WLS) với trọng số bằng nhau đã biết mô tả trong Điều 6 và thực hiện dự đoán mô tả trong Ví dụ 1 của 11.1 và đánh giá tiếp mô tả trong 11.2;

- TS28037_WLS2 chạy ví dụ số về bình phương tối thiểu có trọng số (WLS) với trọng số bằng nhau đã biết mô tả trong Điều 6 và thực hiện dự đoán mô tả trong Ví dụ 2 của 11.1;

- TS28037_WLS3 chạy ví dụ s về bình phương tối thiểu có trọng s (WLS) với trọng số không bằng nhau chưa biết mô tả trong Phụ lục E;

- TS28037_GDR1 chạy ví dụ số về hồi quy khoảng cách tổng quát (GDR) mô tả trong Điều 7;

- TS28037_GDR2 chạy ví dụ số để minh họa thuật toán cho hồi quy khoảng cách tổng quát (GDR) mô tả trong Điều 8;

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

- TS28037_GGMR1 chạy ví dụ số về hồi quy Gauss-Markov tổng quát (GGMR) mô tả trong Điều 10;

- TS28037_GGMR2 chạy ví dụ số về GGMR mô tả trong Điều 10 và Ví dụ 1 của Phụ lục C sử dụng cách tiếp cận thừa số trực giao mô tả trong C.2;

- TS28037_GGMR3 chạy ví dụ số về GGMR mô tả trong Ví dụ 2 của Phụ lục C sử dụng cách tiếp cận thừa số trực giao mô tả trong C.2;

Trong khi dự đoán và đánh giá tiếp chỉ thực hiện trong các tập lệnh giải quyết vấn đề WLS, có thể sao chép và dán mã MATLAB tương ứng với các ứng dụng này của hàm hiệu chuẩn vào bất kỳ tập lệnh nào được cung cấp.

F.3  Phần mềm cần được sử dụng cùng với tiêu chuẩn này. Đặc biệt khuyến nghị người sử dụng nghiên cứu tiêu chuẩn này trước khi chạy phần mềm.

F.4  Phần mềm được cung cấp với thỏa thuận cấp phép phần mềm (REF: MSC/L/10/001) và việc sử dụng phần mềm phải tuân thủ các điều khoản được nêu trong thỏa thuận. Bằng việc chạy mã MATLAB, người sử dụng chấp nhận các điều khoản của thỏa thuận. Các thắc mắc về phần mềm cần được gửi đến NPL tại địa chỉ enquiries@npl.co.uk.

 

Phụ lục G

(tham khảo)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

A

phần chắn của hàm hiệu chuẩn đường thẳng

A*

giá trị chưa biết của A đối với một hệ thống đo cụ thể

a

ước lượng của A

a

véc tơ (a, b)T của ước lượng tham số

B

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

B*

giá trị chưa biết của B đi với một hệ thống đo cụ thể

b

ước lượng của B

cov(a,b)

hiệp phương sai gắn với ab

di

, thể hiện của biến ngẫu nhiên có kỳ vọng bằng 0 và phương sai u2(xi)

ei

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

L

ma trận tam giác dưới

m

số điểm đo được

ri

số dư có trọng số hoặc khoảng cách có trọng số đối với đim dữ liệu thứ i liên quan đến ab

Ri

số dư có trọng số hoặc khoảng cách có trọng số đối với điểm dữ liệu thứ i biểu thị về đại số liên quan đến AB

U

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Ua

ma trận hiệp phương sai kích thước 2 x 2 gắn với a

Ux

ma trận hiệp phương sai kích thước m x m gắn với dữ liệu đo xi, i = 1,..., m

Uy

ma trận hiệp phương sai kích thước m x m gắn với dữ liệu đo yi, i = 1,..., m

uR

độ lệch chuẩn của biến ngẫu nhiên có phân bố mã hóa kiến thức về ảnh hưởng ngẫu nhiên

uS

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

u(z)

độ không đảm bảo chuẩn gắn với z, với z ký hiệu cho a, b, xi, yi, ...

vi

nghịch đảo của u(xi)

wi

nghịch đảo của u(yi)

X

biến độc lập (kích thích)

Xi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

giá trị chưa biết của biến độc lập (kích thích) thứ i cung cấp bởi hệ thống đo

x

ước lượng của X (trong trường hợp dự đoán) hoặc giá trị đo được của X (đánh giá tiếp)

xi

giá trị đo được thứ i của X

ước lượng của biến độc lập (kích thích) thứ i

Y

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Yi

biến (đáp ứng) phụ thuộc thứ i

giá trị chưa biết của biến (đáp ứng) phụ thuộc thứ i cung cấp bởi hệ đo

y

ước lượng của Y (trong trường hợp dự đoán) hoặc giá trị đo được của Y (đánh giá tiếp)

yi

giá trị đo được thứ i của Y

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

n

bậc tự do của mô hình, phân bố Khi-bình phương hoặc phân bố t

σ

độ lệch chuẩn của biến ngẫu nhiên đặc trưng bởi phân bố xác suất

ước lượng hậu nghiệm của σ

giá trị Khi-bình phương quan trắc được

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

 

Thư mục tài liệu tham khảo

[1] Anderson, E., Bai, Z., Bischof, C. H., Blackford, S., Demmel, J., Dongarra, J. J., Croz, J. D., I. A. Greenbaum, Hammarling, S., McKenney, A., and Sorensen, D. C. LAPACK Users’ Guide, 3rd ed., SIAM, Philadelphia, PA, 1999, http://www.netlib.org/lapack/lug/

[2] Bartholomew-Biggs, M., Butler, B. P., and Forbes, A. B. Optimisation algorithms for generalised regression on metrology, In Advanced Mathematical and Computational Tools in Metrology IV (Singapore, 2000), P. Ciarlini, A. B. Forbes, F. Pavese, and D. Richter, Eds., World Scienti c, pp. 21-31 (Các thuật toán tối ưu hóa cho hồi quy tổng quát về đo lường, trong các công cụ toán học và tính toán nâng cao trong đo lường IV)

[3] Boggs, P. T., Byrd, R. H., and Schnabel, R. B. A stable and e cient algorithm for nonlinear orthogonal distance regression, SIAM J. Sci. Stat. Comput. 8, 6 (1987), 1052-1078 (Thuật toán ổn định và hiệu quả cho hồi quy khoảng cách trực giao phi tuyến)

[4] Butler, B. P., Cox, M. G., Ellison, S. L. R., and Hardcastle, W. A., Eds., Statistics Software Quali cation: Reference Data Sets, Royal Society of Chemistry, Cambridge, 1996 (Phần mềm thống kê Quali cation)

[5] Carroll, R. J., Ruppert, D., and Stefanski, L. A. Measurement error in nonlinear models, Chapman& Hall/CRC, Boca Raton, 1995 (Lỗi đo lường trong các mô hình phi tuyến tính)

[6] Cox, M. G., Forbes, A. B., Harris, P. M., and Smith, I. M. The classi cation and solution regression problems for calibration, Tech. Rep. CMSC 24/03, National Physical Laboratory, Teddington, UK, 2003 (Các cation phân loại và giải pháp cho các vấn đề hồi quy đối với hiệu chuẩn)

[7] Draper, N. R., and Smith, H. Applied Regression Analysis, Wiley, New York, 1998, Third edition (Phân tích hồi quy ứng dụng)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

[9] Fuller, W. A. Measurement Error Models, Wiley, New York, 1987

[10] Golub, G. H., and Van Loan, C. F. Matrix Computations, North Oxford Academic, Oxford, 1983

[11] TCVN 8244-1:2010 (ISO 3534-1:2006), Thống kê học - Từ vựng và ký hiệu - Phần 1: Thuật ngữ chung về thống kê và thuật ngữ dùng trong xác suất

[12] TCVN 8244-2:2010 (ISO 3534-2:2006), Thống kê học - Từ vựng và ký hiệu - Phần 2: Thống kê ứng dụng

[13] ISO/IEC Guide 98-3/Suppl 1, Uncertainty of measurement - Part 3: Guide to the expression of uncertainty in measurement (GUM:1995) - Supplement 1: Propagation of distributions using Monte Carlo method (Độ không đảm bảo đo - Phần 3: Hướng dẫn trình bày độ không đảm bảo đo (GUM:1995) - Phần bổ sung 1: Lan truyền phân bố sử dụng phương pháp Monte Carlo)

[14] TCVN 9598:2013 (ISO 11095:1996), Hiệu chuẩn tuyến tính sử dụng mẫu chuẩn

[15] Kendall, M. G., and Stuart, A. The Advanced Theory of Statistics, Volume 2: Inference and Relationship. Charles Gri n, London, 1961 (Lý thuyết thống kê nâng cao)

[16] Kukush, A., and Van Huffel, S. Consistency of elementwise-weighted total least squares estimator in a multivariate errors-in-variables model AX = B, Metrika 59, 1 (February 2004), 75-97 (Tính nhất quán của ước lượng tổng bình phương nhỏ nhất theo trọng số trong mô hình sai lỗi đa biến)

[17] Mardia, K. V., Kent, J. T., and Bibby, J. M. Multivariate Analysis, Academic Press, London, 1979 (Phân tích đa biến)

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

[19] Migon, H. S. and Gamerman, D. Statistical Inference: An Integrated Approach, Arnold, London, 1999

[20] Paige, C. C. Fast numerically stable computations for generalized least squares problems, SIAM J. Numer. Anal. 16 (1979), 165-171 (Tính toán ổn định số nhanh cho các bài toán bình phương tối thiểu)

[21] Strang, G., and Borre, K. Linear Algebra, Geodesy and GPS, Wiley, Wellesley-Cambridge Press, 1997

 

MỤC LỤC

Lời nói đầu

Lời giới thiệu

1  Phạm vi áp dụng

2  Tài liệu viện dẫn

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

4  Quy ước và ký hiệu

5  Các nguyên tắc hiệu chun đường thẳng

5.1  Khái quát

5.2  Đầu vào để xác định hàm hiệu chuẩn

5.3  Xác định hàm hiệu chuẩn

5.4  Xử lý số

5.5  Độ không đảm bảo và hiệp phương sai gắn với tham số hàm hiệu chuẩn

5.6  Xác nhận giá trị sử dụng của mô hình

5.7  Sử dụng hàm hiệu chuẩn

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

6  Mô hình dùng cho độ không đảm bảo gắn với yi

6.1  Khái quát

6.2  Các ước lượng tham số hiệu chuẩn, độ không đảm bảo chuẩn và hiệp phương sai kèm theo

6.3  Xác nhận giá trị sử dụng của mô hình

6.4  Tổ chức tính toán

7  Mô hình độ không đảm bảo gắn với xiyi

7.1  Khái quát

7.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chun và hiệp phương sai kèm theo

7.3  Xác nhận giá trị sử dụng của mô hình

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

8  Mô hình độ không đm bảo gắn với xiyi và hiệp phương sai gắn với cặp (xi, yi)

8.1  Khái quát

8.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn và hiệp phương sai kèm theo

9  Mô hình độ không đảm bảo và hiệp phương sai gắn với yi

9.1  Khái quát

9.2  Ước lượng tham số hiệu chuẩn và độ không đảm bảo chuẩn và hiệp phương sai kèm theo

9.3  Xác nhận giá trị sử dụng của mô hình

9.4  Tổ chức tính toán

10  Mô hình độ không đảm bảo và hiệp phương sai gắn với xiyi

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

10.2  Ước lượng tham số hiệu chuẩn và độ không đm bảo chuẩn và hiệp phương sai kèm theo

10.3  Xác nhận giá trị sử dụng của mô hình

11  Sử dụng hàm hiệu chuẩn

11.1  Dự đoán

11.2  Đánh giá tiếp

Phụ lục A (tham khảo) Phép tính ma trận

Phụ lục B (tham khảo) Áp dụng thuật toán Gauss-Newton cho hồi quy khoảng cách tổng quát

Phụ lục C (tham khảo) Cách tiếp cận thừa số trực giao để giải quyết bài toán Gauss-Markove tổng quát

Phụ lục D (tham khảo) Quy định về độ không đảm bảo và hiệp phương sai gắn với x-giá trị và y-giá trị đo được

...

...

...

Bạn phải đăng nhập hoặc đăng ký Thành Viên TVPL Pro để sử dụng được đầy đủ các tiện ích gia tăng liên quan đến nội dung TCVN.

Mọi chi tiết xin liên hệ: ĐT: (028) 3930 3279 DĐ: 0906 22 99 66

Phụ lục F (tham khảo) Triển khai phần mềm các thuật toán được mô t

Phụ lục G (tham khảo) Từ vựng ký hiệu chính

Thư mục tài liệu tham khảo

Văn bản này chưa cập nhật nội dung Tiếng Anh

Bạn Chưa Đăng Nhập Thành Viên!


Vì chưa Đăng Nhập nên Bạn chỉ xem được Thuộc tính của văn bản.
Bạn chưa xem được Hiệu lực của Văn bản, Văn bản liên quan, Văn bản thay thế, Văn bản gốc, Văn bản tiếng Anh,...


Nếu chưa là Thành Viên, mời Bạn Đăng ký Thành viên tại đây


Bạn Chưa Đăng Nhập Thành Viên!


Vì chưa Đăng Nhập nên Bạn chỉ xem được Thuộc tính của văn bản.
Bạn chưa xem được Hiệu lực của Văn bản, Văn bản liên quan, Văn bản thay thế, Văn bản gốc, Văn bản tiếng Anh,...


Nếu chưa là Thành Viên, mời Bạn Đăng ký Thành viên tại đây


Bạn Chưa Đăng Nhập Thành Viên!


Vì chưa Đăng Nhập nên Bạn chỉ xem được Thuộc tính của văn bản.
Bạn chưa xem được Hiệu lực của Văn bản, Văn bản liên quan, Văn bản thay thế, Văn bản gốc, Văn bản tiếng Anh,...


Nếu chưa là Thành Viên, mời Bạn Đăng ký Thành viên tại đây


Bạn Chưa Đăng Nhập Thành Viên!


Vì chưa Đăng Nhập nên Bạn chỉ xem được Thuộc tính của văn bản.
Bạn chưa xem được Hiệu lực của Văn bản, Văn bản liên quan, Văn bản thay thế, Văn bản gốc, Văn bản tiếng Anh,...


Nếu chưa là Thành Viên, mời Bạn Đăng ký Thành viên tại đây


Tiêu chuẩn quốc gia TCVN 12741:2019 (ISO/TS 28037:2010) về Xác định và sử dụng hàm hiệu chuẩn đường thẳng

Bạn Chưa Đăng Nhập Thành Viên!


Vì chưa Đăng Nhập nên Bạn chỉ xem được Thuộc tính của văn bản.
Bạn chưa xem được Hiệu lực của Văn bản, Văn bản liên quan, Văn bản thay thế, Văn bản gốc, Văn bản tiếng Anh,...


Nếu chưa là Thành Viên, mời Bạn Đăng ký Thành viên tại đây


1.793

DMCA.com Protection Status
IP: 52.15.136.223
Hãy để chúng tôi hỗ trợ bạn!