M Ở ĐẦU
3.1. Xây dựng thuật toán mô phỏng phân bố ngẫu nhiên
thống của kỹ thuật quét gamma phân đoạn. Tìm cách biểu diễn kết quả đo được sao cho hợp lí và chính xác. Mục tiêu là đưa ra phương pháp đánh giá và câu trả lời cho vấn đề đã đặt ra: Khi xử lý kết quả của hệ SGS trong một trường hợp đo đạc thì lấy sai số bao nhiêu là hợp lý?
3.1. Xây dựng thuật toán mô phỏng phân bố ngẫu nhiên của các nguồn phóng xạ. xạ.
Việc bỏ rác ngẫu nhiên vào thùng trên một phân đoạn dạng đĩa được mô phỏng
bằng thuật toán sau (lưu đồ thuật toán ở hình 3.1):
Giai đoạn 1: Lập trình và xuất số liệu trên Fortran (xem phụ lục 2).
B1.1. Sử dụng hàm nội tại (Intrinsic function) của Fortran để gieo một dãy số ngẫu
nhiên phân bố đều từ 0 đến 1.
B1.2. Tạo các cặp số ngẫu nhiên ( , )x y mô tả tọa độ của nguồn bằng cách lấy số ngẫu
nhiên vừa tạo thành trừ đi 0.5 rồi nhân với 58 . Như vậy ta có một bộ các cặp số
( )x y, với − ≤ ≤29 x 29, − ≤ ≤29 y 29 xác định một tập hợp các điểm phân bố đều trên hình vuông giới hạn bởi các đường thẳng x= −29, x=29, y= −29, y=29.
B1.3. Kiểm tra xem điểm vừa gieo có nằm trong phân đoạn hay không ( 2 2 2)
29
y
x + ≤
. Thực hiện tính toán giá trị Is Id cho các điểm nằm trong phân đoạn, tính sự phân bố
xác suất của các khoảng giá trị I Is d , phân bố xác suất của các khoảng sai số. Đối với
các điểm không nằm trong phân đoạn thì gán Is Id =0.
B1.4. Xuất các tập tin số liệu (xem hình 3.2).
21
B2.1. Nhập (Import) các tập tin số liệu vào Origin.
B2.2. Sao chép số liệu qua Excel, sử dụng chức năng sắp xếp dữ liệu (sort data) của Excel để sắp xếp số liệu theo các giá trị Is/Id.
B2.3. Loại bỏ các giá trị ứng với Is/Id = 0.
B2.4. Chuyển số liệu qua Origin và vẽ đồ thị.
22
Hình 3.2.Các tập tin số liệu được xuất.
Hình 3.3.Các tọa độ ( )x y và giá tr, ịIs Id tương ứng.
Tính đúng đắn của thuật toán được minh họa bằng việc gieo một số lượng
23
s d
I I tương ứng, đồ thị nhận được là những hình tròn đồng tâm, nghiệm đúng biểu
thức I Is d = f r( ). Hình 3.3 ứng với trường hợp 1
0.03 cm
µ = − , K =87 cm.
3.2. Đánh giá sai số hệ thống của hệ SGS bằng mô phỏng phân bố ngẫu nhiên của nguồn phóng xạ.