BÀI LUẬN BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN NĂNG LƯỢNG NGUYÊN TỬ VIỆT NAM NGUYỄN THỊ THU HÀ NGHIÊN CỨU VÀ PHÁT TRIỂN VẬT LIỆU LITHIUM ALUMINATE (LiAlO2) ĐỂ ĐO LIỀU PHOTON LUẬN ÁN TIẾ[.]
BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN NĂNG LƯỢNG NGUYÊN TỬ VIỆT NAM - NGUYỄN THỊ THU HÀ NGHIÊN CỨU VÀ PHÁT TRIỂN VẬT LIỆU LITHIUM ALUMINATE (LiAlO2) ĐỂ ĐO LIỀU PHOTON LUẬN ÁN TIẾN SĨ VẬT LÝ Hà Nội – 2023 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN NĂNG LƯỢNG NGUYÊN TỬ VIỆT NAM - NGUYỄN THỊ THU HÀ NGHIÊN CỨU VÀ PHÁT TRIỂN VẬT LIỆU LITHIUM ALUMINATE (LiAlO2) ĐỂ ĐO LIỀU PHOTON LUẬN ÁN TIẾN SĨ VẬT LÝ Chuyên ngành: Vật lý nguyên tử hạt nhân Mã số: 9.44.01.06 Người hướng dẫn khoa học: TS Trịnh Văn Giáp TS Nguyễn Trọng Thành Hà Nội - 2023 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu tơi hướng dẫn TS Trịnh Văn Giáp TS Nguyễn Trọng Thành Các kết nghiên cứu trình bày luận án trung thực, khách quan chưa bảo vệ học vị Tôi xin cam đoan giúp đỡ cho việc thực luận án cám ơn, thơng tin trích dẫn luận án rõ nguồn gốc Tác giả luận án Nguyễn Thị Thu Hà i LỜI CẢM ƠN Để hoàn thành chương trình tiến sĩ viết luận án, tác giả nhận quan tâm, giúp đỡ tận tình tổ chức, cá nhân Lời đầu tiên, tác giả xin bày tỏ lòng biết ơn sâu sắc tới hai thầy hướng dẫn TS Trịnh Văn Giáp TS Nguyễn Trọng Thành tận tình hướng dẫn truyền đạt kiến thức, kinh nghiệm quý báu chuyên môn, động viên tạo điều kiện thuận lợi cho tác giả suốt trình học tập, nghiên cứu hoàn thiện luận án Tác giả xin gửi lời cảm ơn thầy cô, cán tham gia giảng dạy công tác Trung tâm Đào tạo hạt nhân - Viện Năng lượng nguyên tử Việt Nam tận tình giảng dạy, giúp đỡ hỗ trợ thủ tục cần thiết cho tác giả trình học tập, nghiên cứu thực luận án Trong trình thực luận án “Nghiên cứu phát triển vật liệu lithium aluminate (LiAlO2) để đo liều photon”, tác giả nhận nhiều giúp đỡ, tạo điều kiện tập thể ban lãnh đạo, nhà khoa học, cán Viện Khoa học Kỹ thuật Hạt nhân, Viện Khoa học vật liệu, Trung tâm Đào tạo Hạt nhân Tôi xin bày tỏ lòng cảm ơn chân thành giúp đỡ Cuối tác giả muốn giành tình cảm chân thành gửi lời cảm ơn bố mẹ, người thân đặc biệt người chồng yêu quý bên cạnh động viên, hỗ trợ giúp đỡ tơi học tập, làm việc hồn thành luận án Luận án khơng tránh khỏi cịn nhiều khiếm khuyết, thiếu sót, tác giả mong muốn nhận ý kiến đóng góp q thầy cơ, nhà khoa học, đồng nghiệp người quan tâm, để tác giả tiếp tục hồn thiện luận án Hà Nội, ngày tháng năm 2023 Nghiên cứu sinh Nguyễn Thị Thu Hà ii MỤC LỤC DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG BIỂU vii DANH MỤC CÁC HÌNH VẼ ix MỞ ĐẦU Chương NGHIÊN CỨU TỔNG QUAN 1.1 Tương tác xạ với vật chất 1.1.1 Sự ion hóa trực tiếp gián tiếp 1.1.2 Tương tác xạ ion hóa với vật chất 1.2 Các đại lượng, đơn vị phương pháp đo liều xạ 1.2.1 Các đại lượng đơn vị đo liều xạ 1.2.2 Các phương pháp đo liều xạ 14 1.2.3 Phương pháp đo liều nhiệt phát quang 15 1.2.4 Các tham số động học đo liều nhiệt phát quang 23 1.2.5 Vật liệu dùng đo liều nhiệt phát quang 28 1.3 Tổng quan vật liệu LiAlO2 29 1.3.1 Vật liệu LiAlO2 ứng dụng đo liều xạ 29 1.3.2 Tổng quan nghiên cứu vật liệu LiAlO2 giới nước 35 1.4 Tiểu kết chương 40 Chương PHƯƠNG PHÁP NGHIÊN CỨU, CHẾ TẠO VẬT LIỆU LiAlO2 41 2.1 Phương pháp chế tạo vật liệu LiAlO2 41 2.1.1 Dụng cụ thiết bị dùng cho chế tạo vật liệu 41 2.1.2 Phương pháp chế tạo vật liệu LiAlO2 42 2.2 Khảo sát đặc trưng cấu trúc hình thái học vật liệu LiAlO2 45 2.2.1 Phương pháp nhiễu xạ tia X (XRD) 45 2.2.2 Phương pháp kính hiển vi điện tử quét (SEM) 46 2.3 Nghiên cứu khảo sát đặc tính đo liều vật liệu LiAlO2 48 2.3.1 Nguồn chiếu xạ 48 2.3.2 Thiết bị đo nhiệt phát quang 48 2.4 Phân tích đường cong TL LiAlO2 phương pháp giải chập 49 2.5 Tiểu kết chương 51 iii Chương KẾT QUẢ NGHIÊN CỨU VÀ PHÁT TRIỂN VẬT LIỆU LiAlO2 53 3.1 Kết chế tạo khảo sát đặc trưng cấu trúc vật liệu LiAlO2 53 3.1.1 Ảnh hưởng nhiệt độ thiêu kết giá trị pH đến cấu trúc pha mẫu chế tạo 53 3.1.2 Kết đo phổ nhiễu xạ tia X 55 3.1.3 Kết đo kính hiển vi điện tử quét 59 3.2 Kết nghiên cứu khảo sát đặc tính đo liều vật liệu γ-LiAlO2 60 3.2.1 Phông giới hạn phát LiAlO2 60 3.2.2 Độ đồng vật liệu LiAlO2 sau chế tạo 61 3.2.3 Độ tuyến tính 62 3.2.4 Độ nhạy tín hiệu nhiệt phát quang LiAlO2 63 3.2.5 Nghiên cứu, khảo sát khả tái sử dụng vật liệu LiAlO2 65 3.2.6 Đáp ứng liều vật liệu LiAlO2 65 3.2.7 Độ suy giảm tín hiệu TL sau chiếu xạ 68 3.3 Nghiên cứu, khảo sát đường cong TL tham số bẫy vật liệu LiAlO2 72 3.3.1 Dạng đường cong ứng với loại xạ khác 72 3.3.2 Nghiên cứu khảo sát tín hiệu TL vật liệu LiAlO2 thay đổi tốc độ gia nhiệt 74 3.4 Phân tích đường cong TL LiAlO2 phương pháp giải chập 79 3.4.1 Đường cong TL mô so sánh với đường cong thực nghiệm vật liệu LiAlO2 chế tạo phương pháp sol-gel 80 3.4.2 Đường cong TL mô so sánh với đường cong thực nghiệm vật liệu LiAlO2 chế tạo phương pháp sol-gel kết hợp với EDTA 84 3.5 Tiểu kết chương 87 Chương ÁP DỤNG MẠNG NƠRON NHÂN TẠO ĐỂ NHẬN DẠNG, ĐÁNH GIÁ LIỀU CỦA VẬT LIỆU LiAlO2 90 4.1 Giới thiệu mạng nơron nhân tạo 90 4.1.1 Khái niệm mạng nơron nhân tạo 90 4.1.2 Phân loại mạng nơron nhân tạo 90 4.1.3 Mơ hình mạng nhiều lớp 93 4.1.4 Quá trình học mạng nơron 93 iv 4.1.5 Thuật toán lan truyền ngược 95 4.2 Xây dựng cấu trúc mơ hình mạng 96 4.3 Xây dựng sở liệu dùng cho việc huấn luyện mạng 97 4.4 Quá trình huấn luyện 98 4.5 Kết áp dụng mạng nơron nhân tạo để nhận dạng, đánh giá liều vật liệu LiAlO2 100 4.5.1 Kết áp dụng ANN để nhận dạng, đánh giá liều vật liệu LiAlO2 sau 52 lần học 100 4.5.2 Kết áp dụng ANN để nhận dạng, đánh giá liều vật liệu LiAlO2 sau 109 lần học 102 4.5.3 Viết chương trình thu thập, tính tốn liệu, ứng dụng mạng nơron 105 4.6 Tiểu kết chương 109 KẾT LUẬN VÀ KIẾN NGHỊ 110 TÀI LIỆU THAM KHẢO 115 PHỤ LỤC 123 PHỤ LỤC 128 v DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT Kí hiệu Tên tiếng Việt Tên tiếng Anh ANN Artificial Neural Network Mạng nơron nhân tạo E Activation Energy Năng lượng kích hoạt EDTA Etylen Diamin Tetra Acetic Hợp chất hóa học EDTA EPR Electron Paramagnetic Resonance Cộng hưởng thuận từ điện tử FOM Figure Of Merit Hệ số làm khớp GOT General Of Trap Một bẫy tổng quát ICRP International Commission on Radiation Protection ICRU International Commission on Radiation Units and Measurements Ủy ban quốc tế đơn vị đo lường xạ LiAlO2 Lithium aluminate Vật liệu LiAlO2 LiF Lithi fluoride Hợp chất vô LiF MLP Multilayer perceptron Mạng truyền thẳng nhiều lớp OSL Optical Stimulated Luminesence Phát quang cưỡng quang s Frequency Factor Tần số thoát SEM Scanning Electron Microscope Kính hiển vi điện tử quét TL Thermo-Luminescence Nhiệt phát quang TLD Thermo-Luminescence Detector Đầu dò nhiệt phát quang TSC Thermal Stimulated Conductivity Độ dẫn điện cưỡng nhiệt TSL Thermally Stimulated Luminescence XRD X-Ray Diffraction Ủy ban quốc tế An toàn xạ Phát quang cưỡng nhiệt Nhiễu xạ tia X vi DANH MỤC CÁC BẢNG BIỂU STT Tiêu đề bảng Trang 1.1 Giá trị trọng số xạ (wR) số loại xạ ion hóa có lượng khác 12 1.2 Giá trị trọng số mô wT ICRP khuyến cáo 13 1.3 Một vài tính chất lithium aluminate 30 1.4 Đặc tính đo liều số vật liệu nhiệt phát quang 39 3.1 Cấu trúc pha vật liệu bị ảnh hưởng việc điều chỉnh giá trị pH thiêu kết nhiệt độ 600 0C 53 3.2 Cấu trúc pha vật liệu bị ảnh hưởng việc điều chỉnh giá trị pH thiêu kết nhiệt độ 800 0C 54 3.3 Cấu trúc pha vật liệu bị ảnh hưởng việc điều chỉnh giá trị pH thiêu kết nhiệt độ 900 0C 54 3.4 Dữ liệu cấu trúc vật liệu γ- LiAlO2 sau chế tạo 59 3.5 Giá trị phông giới hạn phát LiAlO2 60 3.6 Kết đánh giá độ đồng LiAlO2 62 3.7 Đánh giá độ tuyến tính đo liều LiAlO2 63 3.8 Kết kiểm tra khả tái sử dụng vật liệu LiAlO2 chế tạo phương pháp sol-gel kết hợp với EDTA sau chiếu xạ gamma 4mGy 65 3.9 Suy giảm tín hiệu TL số vật liệu khác sử dụng đo liều xạ vật liệu chế tạo 71 3.10 Dịch chuyển đỉnh phổ TL LiAlO2 theo tốc độ gia nhiệt 74 3.11 Các tham số bẫy xác định phương pháp dạng đỉnh LiAlO2 chiếu xạ gamma 0,3 Gy với tốc độ gia nhiệt 0C/s 77 vii 3.12 Các tham số bẫy xác định phương pháp dạng đỉnh LiAlO2 chiếu xạ beta 30 phút với tốc độ gia nhiệt 10 C/s 77 3.13 Các tham số động học tính phương pháp vị trí đỉnh LiAlO2 chiếu xạ chiếu xạ gamma 0,3 Gy với tốc độ gia nhiệt 0C/s 78 3.14 Các tham số động học tính phương pháp vị trí đỉnh LiAlO2 chiếu xạ beta với tốc độ gia nhiệt 10 0C/s 78 3.15 Các tham số bẫy LiAlO2 chế tạo sol-gel chiếu xạ gamma 0,3 Gy 81 3.16 Các tham số bẫy LiAlO2 chế tạo sol-gel với EDTA chiếu xạ gamma 0,3 Gy 85 4.1 Bảng tóm tắt kết tính tốn mạng nơron sau 52 lần học 101 4.2 Bảng tóm tắt kết tính tốn mạng nơron sau 109 lần học 104 P2.1 Bảng kết tính tốn mạng nơron sau 52 lần học 128 P2.2 Bảng kết tính tốn mạng nơron sau 109 lần học 131 viii 65 S.W Kwon, S.B Park (1997), Effect of precursors on the preparation of lithium aluminate Mater Sci 246, 131-138 66 J.C Mittani, M Prokic, E.G Yukihara (2008), Opticaly stimulated luminescence and thermoluminescence of terbium-activated silicates and aluminates Radiat Meas., 43, 323-326 67 J.I Lee, A.S Pradhan, J.L Kim, I Chang, B.H Kim, K.S Chung (2012), Preliminary study on development and characterization of high sensitivity LiAlO2 optically stimulated luminescence material Radiat Meas 47, 837-840 68 B Dhabekar, E Alagu Raja, S Menon, T.K Gundu Rao, R.K Kher, B.C Bhatt (2008), ESR, PL and TL studies of LiAlO2: Mn/Ce phosphor Radiat Meas 43, 291-294 69 H Teng, S Zhou, H Lin, T Jia, X Hou, J Wang (2010), Growth and characterization of high-quality Mn-doped LiAlO2 single crystal Chin Opt Lett 8, 414-417 70 M Sadel, P Bilski, J Swakon (2013), Relative TL and OSL ef ficiency to protons of various dosimetric materials Radiat Prot Dosim 71 J.I Lee, A.S Pradhan, J.L Kim, I Chang, B.H Kim, K.S Chung (2013), Characteristics of LiAlO2- radioluminescence and optically stimulated luminescence Radiat Meas 56, 217-222 72 P Twardak, Bilski, B Marczewska, W Gieszczyk (2014), Analysis of TL and OSL kinetics of lithium aluminate Radiation Measurement, 1-5 73.Bế Kim Giáp (2005), Nghiên cứu huỳnh quang cưỡng nhiệt, cưỡng quang vật liệu Al2O3 ứng dụng đo liều xạ, Luận án Tiến sĩ Vật lý, Viện Vật lý Điện tử, Hà Nội 74.Ngô Quang Thành, Nguyễn Trọng Thành, Vũ Xuân Quang, Nguyễn Thế Khôi (2006), “Các đặc trưng quang học đo liều nhiệt phát quang vật liệu Li2B4O7:Cu”, Tuyển tập báo cáo Hội nghị “Advances in Optics Photonics Spectroscopy and Application”, TP Cần Thơ, Tr.394-399 75 Lê Văn Tuất, Ngô Khoa Quang, Ngô Quang Thành (2007), “Chế tạo liều kế nhiệt phát quang CaSO4: Dy dạng viên nén dùng đo liều xạ”, Tạp chí Khoa học Đại học Huế: Khoa học Tự nhiên, 35, pp 61-68 76 T Ngọc, V.X Quang, N.M Sơn (2004), “Phát quang cưỡng quang thạch anh tự nhiên Việt Nam”, Tạp chí Khoa học Đại học Huế, số 25, tr 61-67 120 77 Ngô Quang Thành (2009), Nghiên cứu tượng phát quang cưỡng nhiệt số vật liệu rắn nhằm ứng dụng đo liều xạ, Luận án tiến sĩ Vật lý, Viện Vật lý, VAST 78 Trần thị Thanh Lam (2020), Chế tạo vật liệu K2GdF5: Tb phương pháp phản ứng pha rắn ứng dụng đo liều, Luận văn thạc sĩ Vật lý, Học viện Khoa học Công nghệ, VAST 79 J.V Dam, G Marinello (1994), Methods for in vio dosimetry in external radiotherapy, Physics for clinical radiotherapy, Booklet No.1 1994 80 D.J Huntley, D.L Godfrey-Smith, M.L.W Thewalt (1985), Optical dating of sediments Nature, 313, 105–107 81 S.W.S McKeever, M Moscovitch (2003), On the advantages and disadvantages of optically stimulated luminescence dosimetry and thermoluminescence dosimetry Radiat Prot Dosim., 104, 263–270 82 M.S Akselrod, S.W.S McKeever (1999), A radiation dosimetry method using pulsed optically stimulated luminescence Radiat Prot Dosim., 81, 167–176 83 L Brian et.al (2000), Dosimetry measurement, Copyright CRC Press LLC 84 A Halperin, and A A Braner (1960), Evaluation of thermal activation energies from glow curves Phys Rev., 117: 408-424 85 I Broser, and R Broser-Warminsky (1954), Luminescence and electrical conductivity of crystal phosphors Br J Appl Phys., Suppl., 4, 90-94 86 G.A Dussel, and R.H Bube (1967), Theory of thermally stimulated conductivity in a previously photoexcited crystal Phys Rev., 155, 764779 87 R Visocekas, (1985), Tunneling radiative recombination in labradorite: its association with anomalous fading of thermoluminescence, Nucl Tracks, 10, 521–529 88 R.A Ribeiro, G.G Silva, N.D.S Mohallem (2001), The influences of heat treatment on the structural properties of lithium aluminates, Journal of Physics and Chemistry of Solids, 62, 857–864 89 IEC 62387 Ed.1, 2012 90 K.S.V Nambi (1977), Thermoluminescence, Its understanding and applications, Instituti De Energia Atomica, Sao Paulo-Brasil 91 Karan Kumar Gupta, R.M Kadam, N.S Dhoble, S.P Lochab, S.J Dhoble (2007), A comparative investigation of Ce3ỵ/Dy3ỵ and Eu2ỵ doped 121 LiAlO2 phosphors for high dose radiation dosimetry: Explanation of defect recombination mechanism using PL, TL and EPR study, Journal of Luminescence, 188, 81–95 92 Nguyễn Dỗn Phước (2016), Tối ưu hóa điều khiển điều khiển tối ưu, Nhà xuất bách khoa Hà Nội 122 PHỤ LỤC CHƯƠNG TRÌNH LÀM KHỚP THEO MƠ HÌNH FOK, SOK, GOK BẰNG MATLAB % MATLAB SCRIPT % XOA TAT CA clear all; % XAY DUNG DU LIEU DUONG CONG THEO MO HINH GOK disp ('CHUONG TRINH XAY LAM KHOP THEO MO HINH FOK, SOK, GOK'); % DOC FILE DU LIEU DAU VAO filename = 'DatavatlieuLiAlO2.xls'; DATA = xlsread(filename); % LAY DU NHIEU NHIET DO ,DON VI K TK= DATA(:,2); % LAY DU LIEU CUONG DO IS1= DATA(:,3); IS2= DATA(:,4); IS3= DATA(:,5); IS4= DATA(:,6); IS5= DATA(:,7); IS6= DATA(:,8); IS7= DATA(:,9); IS8= DATA(:,10); % VE DU LIEU PHO DO DC plot(TK, IS1,'*','MarkerIndices',1:5:length(IS1)); plot(TK, IS2,'*','MarkerIndices',1:5:length(IS2)); plot(TK, IS3,'*','MarkerIndices',1:5:length(IS3)); plot(TK, IS4,'*','MarkerIndices',1:5:length(IS4)); plot(TK, IS5,'*','MarkerIndices',1:5:length(IS5)); plot(TK, IS6,'*','MarkerIndices',1:5:length(IS6)); plot(TK, IS7,'*','MarkerIndices',1:5:length(IS7)); plot(TK, IS8,'*','MarkerIndices',1:5:length(IS8)); %======================================================================= % XAY DUNG HAM FTTING THEO MO MINH FOK T= TK; I1=Im1*exp(1+(E1./(8.617385*10^-5.*T)).*((T-Tm1)/Tm1)-((T.*T)/(Tm1.*Tm1)).* exp((E1./(8.617385*10^-5*T)).*((T-Tm1)/Tm1).*(1 - (2.*8.617385*10^-5.*Tm1)/E1))- (2.*8.617385*10^- 123 5.*Tm1)/E1); I2=Im2*exp(1+(E2./(8.617385*10^-5.*T)).*((T-Tm2)/Tm2)-((T.*T)/(Tm2.*Tm2)).* exp((E2./(8.617385*10^-5*T)).*((T-Tm2)/Tm2).*(1 - (2.*8.617385*10^-5.*Tm2)/E2))- (2.*8.617385*10^5.*Tm2)/E2); I3=Im3*exp(1+(E3./(8.617385*10^-5.*T)).*((T-Tm3)/Tm3)-((T.*T)/(Tm3.*Tm3)).* exp((E3./(8.617385*10^-5*T)).*((T-Tm3)/Tm3).*(1 - (2.*8.617385*10^-5.*Tm3)/E3))- (2.*8.617385*10^5.*Tm3)/E3); I4=Im4*exp(1+(E4./(8.617385*10^-5.*T)).*((T-Tm4)/Tm4)-((T.*T)/(Tm4.*Tm4)).* exp((E4./(8.617385*10^-5*T)).*((T-Tm4)/Tm4).*(1 - (2.*8.617385*10^-5.*Tm4)/E4))- (2.*8.617385*10^5.*Tm4)/E4); I5=Im5*exp(1+(E5./(8.617385*10^-5.*T)).*((T-Tm5)/Tm5)-((T.*T)/(Tm5.*Tm5)).* exp((E5./(8.617385*10^-5*T)).*((T-Tm5)/Tm5).*(1 - (2.*8.617385*10^-5.*Tm5)/E5))- (2.*8.617385*10^5.*Tm5)/E5); I6=Im6*exp(1+(E6./(8.617385*10^-5.*T)).*((T-Tm6)/Tm6)-((T.*T)/(Tm6.*Tm6)).* exp((E6./(8.617385*10^-5*T)).*((T-Tm6)/Tm6).*(1 - (2.*8.617385*10^-5.*Tm6)/E6))- (2.*8.617385*10^5.*Tm6)/E6); I7=Im7*exp(1+(E7./(8.617385*10^-5.*T)).*((T-Tm7)/Tm7)-((T.*T)/(Tm7.*Tm7)).* exp((E7./(8.617385*10^-5*T)).*((T-Tm7)/Tm7).*(1 - (2.*8.617385*10^-5.*Tm7)/E7))- (2.*8.617385*10^5.*Tm7)/E7); I8=Im8*exp(1+(E8./(8.617385*10^-5.*T)).*((T-Tm8)/Tm8)-((T.*T)/(Tm8.*Tm8)).* exp((E8./(8.617385*10^-5*T)).*((T-Tm8)/Tm8).*(1 - (2.*8.617385*10^-5.*Tm8)/E8))- (2.*8.617385*10^5.*Tm8)/E8); IT =I1+I2+I3+I4+I5+I6; % XAY DUNG HAM FTTING THEO MO MINH SOK %======================================================================= T= TK; I1=4.*Im1.*(exp((E1./(8.617385*10^-5.*T)).*((T-Tm1)/Tm1))).*(((T.*T)/(Tm1.*Tm1)).*(1 (2.*8.617385*10^-5.*Tm1)/E1).*( exp((E1./(8.617385*10^-5.*T)).*((T-Tm1)/Tm1)))+(1+ (2.*8.617385*10^-5.*Tm1)/E1)).^-2; I2=4.*Im2.*(exp((E2./(8.617385*10^-5.*T)).*((T-Tm2)/Tm2))).*(((T.*T)/(Tm2.*Tm2)).*(1 (2.*8.617385*10^-5.*Tm2)/E2).*( exp((E2./(8.617385*10^-5.*T)).*((T-Tm2)/Tm2)))+(1+ (2.*8.617385*10^-5.*Tm2)/E2)).^-2; I3=4.*Im3.*(exp((E3./(8.617385*10^-5.*T)).*((T-Tm3)/Tm3))).*(((T.*T)/(Tm3.*Tm3)).*(1 (2.*8.617385*10^-5.*Tm3)/E3).*( exp((E3./(8.617385*10^-5.*T)).*((T-Tm3)/Tm3)))+(1+ (2.*8.617385*10^-5.*Tm3)/E3)).^-2; I4=4.*Im4.*(exp((E4./(8.617385*10^-5.*T)).*((T-Tm4)/Tm4))).*(((T.*T)/(Tm4.*Tm4)).*(1 (2.*8.617385*10^-5.*Tm4)/E4).*( exp((E4./(8.617385*10^-5.*T)).*((T-Tm4)/Tm4)))+(1+ (2.*8.617385*10^-5.*Tm4)/E4)).^-2; I5=4.*Im5.*(exp((E5./(8.617385*10^-5.*T)).*((T-Tm5)/Tm5))).*(((T.*T)/(Tm5.*Tm5)).*(1 (2.*8.617385*10^-5.*Tm5)/E5).*( exp((E5./(8.617385*10^-5.*T)).*((T-Tm5)/Tm5)))+(1+ (2.*8.617385*10^-5.*Tm5)/E5)).^-2; I6=4.*Im6.*(exp((E6./(8.617385*10^-5.*T)).*((T-Tm6)/Tm6))).*(((T.*T)/(Tm6.*Tm6)).*(1 (2.*8.617385*10^-5.*Tm6)/E6).*( exp((E6./(8.617385*10^-5.*T)).*((T-Tm6)/Tm6)))+(1+ (2.*8.617385*10^-5.*Tm6)/E6)).^-2; 124 IT =I1+I2+I3+I4+I5+I6; % XAY DUNG HAM FTTING THEO MO MINH GOK %%=============================================================== I1=Im1.*(b1.^( b1./(b1-1))).*exp((E1./(8.617385*10^-5.*T)).*((T-Tm1)/Tm1)).*( (b1-1).*(1 (2.*8.617385*10^-5.*Tm1)/E1).*((T.*T)/(Tm1.*Tm1)).*(exp((E1./(8.617385*10^-5.*T)).*((TTm1)/Tm1)))+1+((2.*8.617385*10^-5.*Tm1).*(b1-1))./E1).^ -( b1./(b1-1)); I2=Im2.*(b2.^( b2./(b2-1))).*exp((E2./(8.617385*10^-5.*T)).*((T-Tm2)/Tm2)).*( (b2-1).*(1 (2.*8.617385*10^-5.*Tm2)/E2).*((T.*T)/(Tm2.*Tm2)).*(exp((E2./(8.617385*10^-5.*T)).*((TTm2)/Tm2)))+1+((2.*8.617385*10^-5.*Tm2).*(b2-1))./E2).^ -( b2./(b2-1)); I3=Im3.*(b3.^( b3./(b3-1))).*exp((E3./(8.617385*10^-5.*T)).*((T-Tm3)/Tm3)).*( (b3-1).*(1 (2.*8.617385*10^-5.*Tm3)/E3).*((T.*T)/(Tm3.*Tm3)).*(exp((E3./(8.617385*10^-5.*T)).*((TTm3)/Tm3)))+1+((2.*8.617385*10^-5.*Tm3).*(b3-1))./E3).^ -( b3./(b3-1)); I4=Im4.*(b4.^( b4./(b4-1))).*exp((E4./(8.617385*10^-5.*T)).*((T-Tm4)/Tm4)).*( (b4-1).*(1 (2.*8.617385*10^-5.*Tm4)/E4).*((T.*T)/(Tm4.*Tm4)).*(exp((E4./(8.617385*10^-5.*T)).*((TTm4)/Tm4)))+1+((2.*8.617385*10^-5.*Tm4).*(b4-1))./E4).^ -( b4./(b4-1)); I5=Im5.*(b5.^( b5./(b5-1))).*exp((E5./(8.617385*10^-5.*T)).*((T-Tm5)/Tm5)).*( (b5-1).*(1 (2.*8.617385*10^-5.*Tm5)/E5).*((T.*T)/(Tm5.*Tm5)).*(exp((E5./(8.617385*10^-5.*T)).*((TTm5)/Tm5)))+1+((2.*8.617385*10^-5.*Tm5).*(b5-1))./E5).^ -( b5./(b5-1)); I6=Im6.*(b6.^( b6./(b6-1))).*exp((E6./(8.617385*10^-5.*T)).*((T-Tm6)/Tm6)).*( (b6-1).*(1 (2.*8.617385*10^-5.*Tm6)/E6).*((T.*T)/(Tm6.*Tm6)).*(exp((E6./(8.617385*10^-5.*T)).*((TTm6)/Tm6)))+1+((2.*8.617385*10^-5.*Tm6).*(b6-1))./E6).^ -( b6./(b6-1)); %======================================================================= % LUA CHON HAM FITTING, s = fitoptions('Method','NonlinearLeastSquares', 'Lower',[0,0], 'Upper',[Inf,max(cdate)], 'Startpoint',[1 1]); % FITTING HAM f = fittype(IT,'problem'); [kqparameter ,parameter] = fit(cdate,pop,f,'problem',2) %=========CHAY KET QUA FITTING==================================== E1 = kqparameter.E1; E2 = kqparameter.E2; E3 = kqparameter.E3; E4 = kqparameter.E4; E5 = kqparameter.E5; E6 = kqparameter.E6; Tm1 = kqparameter.Tm1; Tm2 = kqparameter.Tm2; 125 Tm3 = kqparameter.Tm3; Tm4 = kqparameter.Tm4; Tm5 = kqparameter.Tm5; Tm6 = kqparameter.Tm6; Im1 = kqparameter.Im1; Im2 = kqparameter.Im2; Im3 = kqparameter.Im3; Im4 = kqparameter.Im4; Im5 = kqparameter.Im5; Im6 = kqparameter.Im6; b1 = kqparameter.b1; b2 = kqparameter.b2; b3 = kqparameter.b3; b4 = kqparameter.b4; b5 = kqparameter b5; b6 = kqparameter.b6; % VE DUONG CONG FITING VA THUC NGHIEM p=plot(T,IT ,T,I1,T,I2,T,I3,T,I4,T,I5,T,I6,T,I7,T,I8); p(1).LineWidth = 2; hold; xlabel('Temperature(K)'); ylabel('Intensity(a.u)'); grid on xlim([321 650]) ylim([-inf inf]) legend('10','8','7','5','3','2'); legend boxoff p(1).LineWidth = 1; p(1).LineWidth = 1; p(2).LineWidth = 1; p(3).LineWidth = 1; p(4).LineWidth = 1; p(5).LineWidth = 1; p(6).LineWidth = 1; %TINH HE SO THOAT S Beta = 10; %TOC DO GIA NHIET CUA MAY 126 delta1=(2*8.617385*10^-5.*Tm1)/E1; delta2=(2*8.617385*10^-5.*Tm2)/E2; delta3=(2*8.617385*10^-5.*Tm3)/E3; delta4=(2*8.617385*10^-5.*Tm4)/E4; delta5=(2*8.617385*10^-5.*Tm5)/E5; S1 = ((Beta.*E1)./(8.617385*10^-5.*Tm1.*Tm1)).*(exp(E1./(8.617385*10^-5.*Tm1)))*(1./(1+delta1)); S2 = ((Beta.*E2)./(8.617385*10^-5.*Tm2.*Tm2)).*(exp(E2./(8.617385*10^-5.*Tm2)))*(1./(1+delta2)); S3 = ((Beta.*E3)./(8.617385*10^-5.*Tm3.*Tm3)).*(exp(E3./(8.617385*10^-5.*Tm3)))*(1./(1+delta3)); S4 = ((Beta.*E4)./(8.617385*10^-5.*Tm4.*Tm4)).*(exp(E4./(8.617385*10^-5.*Tm4)))*(1./(1+delta4)); S5 = ((Beta.*E4)./(8.617385*10^-5.*Tm5.*Tm3)).*(exp(E5./(8.617385*10^-5.*Tm5)))*(1./(1+delta5)); %======================================================================= % TINH TOAN SAI SO FOM TongFIT =sum(IT); HieuFIT =sum( abs(IS1-IT)); FOM = (HieuFIT /TongFIT)*100 ; %======================================================================= 127 PHỤ LỤC Bảng P2.1: Bảng kết tính tốn mạng nơron sau 52 lần học STT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 ĐẦU VÀO HỌC Tốc độ Liều Đỉnh gia nhiệt chiếu (0C/s) (Gy) (K) 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 368,00 369,00 372,00 373,00 374,00 376,00 378,00 382,00 385,00 388,00 390,00 392,00 394,00 395,00 397,00 398,00 400,00 401,00 402,00 403,00 404,00 405,00 371,00 372,00 374,00 375,00 376,00 378,00 381,00 385,00 388,00 390,00 392,00 394,00 Đỉnh (K) 477,00 479,00 482,00 483,00 484,00 486,00 490,00 494,00 498,00 501,00 504,00 506,00 508,00 510,00 512,00 514,00 515,00 517,00 518,00 520,00 521,00 522,00 478,00 479,00 482,00 483,00 485,00 487,00 490,00 494,00 497,00 500,00 503,00 505,00 128 KẾT QUẢ HỌC AI (52 LẦN HỌC) Tốc độ Liều Đỉnh Đỉnh gia nhiệt chiếu (0C/s) (Gy) (K) (K) 1,48 1,57 1,77 1,87 1,97 2,18 2,50 3,00 3,48 3,96 4,45 4,95 5,47 6,00 6,53 7,06 7,58 8,08 8,56 9,03 9,47 9,89 1,50 1,59 1,79 1,89 1,99 2,18 2,47 2,97 3,48 3,99 4,49 4,99 0,27 0,28 0,29 0,30 0,30 0,31 0,31 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,29 0,29 0,28 0,28 0,88 0,89 0,90 0,91 0,91 0,91 0,91 0,90 0,89 0,88 0,87 0,88 373,00 373,70 375,03 375,66 376,27 377,42 379,02 381,39 383,53 385,54 387,46 389,34 391,15 392,91 394,59 396,19 397,71 399,13 400,47 401,72 402,89 403,97 369,77 371,12 373,42 374,43 375,38 377,12 379,48 382,94 386,03 388,91 391,63 394,22 483,58 484,42 486,04 486,81 487,56 488,99 490,98 493,96 496,62 499,05 501,30 503,39 505,34 507,14 508,80 510,33 511,73 513,01 514,18 515,24 516,21 517,10 477,29 478,09 480,03 481,09 482,18 484,40 487,70 492,86 497,39 501,27 504,56 507,35 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 396,00 398,00 400,00 401,00 402,00 404,00 405,00 406,00 407,00 408,00 379,00 380,00 382,00 383,00 384,00 386,00 388,00 391,00 393,00 395,00 397,00 399,00 400,00 402,00 403,00 404,00 405,00 406,00 407,00 408,00 409,00 410,00 374,00 376,00 378,00 379,00 380,00 382,00 384,00 388,00 507,00 509,00 511,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 478,00 480,00 483,00 484,00 485,00 487,00 490,00 494,00 498,00 500,00 503,00 505,00 507,00 509,00 511,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 479,00 480,00 483,00 484,00 485,00 488,00 491,00 495,00 129 5,49 5,99 6,49 6,99 7,49 7,99 8,49 8,99 9,49 9,98 1,52 1,58 1,82 1,94 2,06 2,26 2,53 2,96 3,45 3,98 4,52 5,04 5,54 6,03 6,50 6,98 7,46 7,95 8,45 8,96 9,48 10,01 1,60 1,60 1,75 1,87 1,99 2,21 2,52 3,00 0,88 0,89 0,90 0,90 0,91 0,91 0,91 0,90 0,90 0,89 1,21 1,19 1,18 1,18 1,18 1,19 1,19 1,20 1,21 1,20 1,20 1,19 1,19 1,19 1,19 1,20 1,20 1,20 1,21 1,20 1,20 1,19 1,50 1,49 1,49 1,50 1,51 1,53 1,53 1,52 396,67 398,95 401,04 402,94 404,64 406,15 407,49 408,67 409,70 410,60 379,24 381,15 383,89 384,89 385,72 387,08 388,70 390,99 393,08 395,05 396,89 398,60 400,18 401,61 402,92 404,10 405,16 406,11 406,97 407,75 408,44 409,06 372,46 375,01 378,93 380,33 381,39 382,80 383,98 385,92 509,73 511,75 513,48 514,96 516,21 517,28 518,19 518,95 519,58 520,10 479,75 480,60 482,95 484,18 485,35 487,50 490,25 494,00 497,15 499,91 502,41 504,69 506,79 508,73 510,52 512,17 513,68 515,05 516,29 517,38 518,32 519,12 482,82 482,49 483,19 483,89 484,68 486,36 488,80 492,49 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,40 2,40 2,40 2,40 391,00 393,00 395,00 397,00 399,00 401,00 402,00 403,00 405,00 406,00 407,00 408,00 409,00 410,00 367,00 368,00 371,00 372,00 373,00 376,00 378,00 382,00 386,00 389,00 391,00 394,00 396,00 398,00 399,00 401,00 402,00 404,00 405,00 406,00 408,00 409,00 367,00 368,00 371,00 372,00 498,00 501,00 504,00 506,00 508,00 510,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 522,00 475,00 476,00 479,00 481,00 482,00 484,00 487,00 492,00 496,00 499,00 502,00 504,00 506,00 508,00 510,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 476,00 477,00 480,00 481,00 130 3,49 3,98 4,48 4,98 5,49 6,00 6,50 7,00 7,50 8,00 8,49 8,98 9,48 9,97 1,49 1,50 1,75 1,89 2,00 2,21 2,50 3,00 3,49 3,99 4,48 4,99 5,49 6,00 6,50 6,99 7,48 7,99 8,49 9,00 9,49 9,98 1,49 1,59 1,79 1,89 1,50 1,49 1,48 1,48 1,49 1,49 1,50 1,51 1,51 1,51 1,51 1,50 1,49 1,48 2,09 2,12 2,12 2,11 2,10 2,09 2,09 2,10 2,09 2,11 2,13 2,13 2,12 2,11 2,10 2,10 2,09 2,10 2,10 2,10 2,10 2,10 2,41 2,39 2,39 2,40 388,70 391,88 394,97 397,67 399,90 401,65 403,00 404,05 404,88 405,55 406,12 406,62 407,06 407,47 368,06 368,40 370,17 371,24 372,34 374,61 378,40 384,28 387,10 388,09 389,45 391,77 394,46 396,96 399,09 400,91 402,53 404,00 405,32 406,43 407,29 407,88 367,35 368,74 370,50 371,27 495,89 499,13 502,22 505,11 507,73 510,06 512,14 514,02 515,75 517,35 518,85 520,23 521,50 522,64 474,82 476,34 478,90 480,12 481,32 483,61 486,77 491,77 496,25 500,36 504,00 506,95 509,01 510,24 511,01 511,74 512,70 513,98 515,49 517,08 518,56 519,79 476,06 477,25 479,83 481,21 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 373,00 375,00 377,00 381,00 384,00 387,00 389,00 391,00 393,00 395,00 397,00 398,00 399,00 401,00 402,00 403,00 404,00 405,00 482,00 485,00 488,00 492,00 495,00 499,00 501,00 504,00 506,00 508,00 509,00 511,00 513,00 514,00 515,00 517,00 518,00 519,00 2,00 2,20 2,49 3,00 3,50 3,98 4,49 5,01 5,50 5,98 6,46 6,96 7,47 7,98 8,50 9,00 9,49 9,96 2,40 2,41 2,40 2,39 2,40 2,42 2,41 2,40 2,39 2,38 2,38 2,39 2,39 2,40 2,40 2,40 2,40 2,40 372,10 374,02 377,31 382,06 384,62 386,30 388,81 391,74 394,32 396,34 397,93 399,24 400,38 401,40 402,29 403,02 403,56 403,90 482,59 485,15 488,04 491,03 494,51 499,06 502,99 505,54 507,00 507,94 508,82 509,90 511,28 512,93 514,78 516,69 518,53 520,21 Bảng P2.2: Bảng kết tính tốn mạng nơron sau 109 lần học ĐẦU VÀO HỌC STT 10 11 12 13 14 15 Tốc độ gia nhiệt (0C/s) Liều chiếu (Gy) 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 Đỉnh Đỉnh (K) (K) 368,00 369,00 372,00 373,00 374,00 376,00 378,00 382,00 385,00 388,00 390,00 392,00 394,00 395,00 397,00 477,00 479,00 482,00 483,00 484,00 486,00 490,00 494,00 498,00 501,00 504,00 506,00 508,00 510,00 512,00 131 KẾT QUẢ HỌC AI (109 LẦN HỌC) Tốc độ Liều Đỉnh Đỉnh gia nhiệt chiếu (0C/s) (Gy) (K) (K) 1,50 1,60 1,80 1,90 2,00 2,20 2,49 2,99 3,51 4,01 4,51 5,00 5,50 5,99 6,49 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,30 368,17 369,43 371,80 372,90 373,96 375,91 378,51 382,14 385,12 387,67 389,93 391,95 393,79 395,47 397,00 477,27 478,75 481,53 482,83 484,07 486,38 489,50 493,96 497,69 500,87 503,62 506,05 508,23 510,20 512,03 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 398,00 400,00 401,00 402,00 403,00 404,00 405,00 371,00 372,00 374,00 375,00 376,00 378,00 381,00 385,00 388,00 390,00 392,00 394,00 396,00 398,00 400,00 401,00 402,00 404,00 405,00 406,00 407,00 408,00 379,00 380,00 382,00 383,00 384,00 386,00 388,00 391,00 393,00 395,00 397,00 514,00 515,00 517,00 518,00 520,00 521,00 522,00 478,00 479,00 482,00 483,00 485,00 487,00 490,00 494,00 497,00 500,00 503,00 505,00 507,00 509,00 511,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 478,00 480,00 483,00 484,00 485,00 487,00 490,00 494,00 498,00 500,00 503,00 132 6,99 7,50 8,01 8,51 9,01 9,50 9,99 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,62 1,81 1,91 2,00 2,20 2,50 3,00 3,50 4,00 4,50 0,30 0,30 0,30 0,30 0,30 0,30 0,30 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 0,90 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 398,39 399,67 400,82 401,88 402,83 403,70 404,49 370,43 371,87 374,38 375,49 376,51 378,37 380,83 384,35 387,39 390,05 392,40 394,50 396,37 398,06 399,59 401,01 402,35 403,62 404,81 405,91 406,91 407,78 379,16 380,53 382,65 383,52 384,32 385,78 387,73 390,61 393,19 395,50 397,48 513,72 515,31 516,80 518,21 519,53 520,78 521,96 477,86 479,48 482,32 483,58 484,75 486,87 489,65 493,62 497,07 500,12 502,81 505,17 507,27 509,13 510,82 512,39 513,87 515,29 516,63 517,87 519,00 519,98 478,56 480,25 482,98 484,15 485,24 487,26 490,03 494,12 497,66 500,67 503,16 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 2,10 2,10 2,10 2,10 2,10 2,10 2,10 399,00 400,00 402,00 403,00 404,00 405,00 406,00 407,00 408,00 409,00 410,00 374,00 376,00 378,00 379,00 380,00 382,00 384,00 388,00 391,00 393,00 395,00 397,00 399,00 401,00 402,00 403,00 405,00 406,00 407,00 408,00 409,00 410,00 367,00 368,00 371,00 372,00 373,00 376,00 378,00 505,00 507,00 509,00 511,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 479,00 480,00 483,00 484,00 485,00 488,00 491,00 495,00 498,00 501,00 504,00 506,00 508,00 510,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 522,00 475,00 476,00 479,00 481,00 482,00 484,00 487,00 133 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,19 2,47 3,00 3,58 4,10 4,54 5,00 5,50 6,00 6,51 7,00 7,50 8,00 8,50 9,01 9,50 9,97 1,49 1,60 1,80 1,90 2,00 2,20 2,50 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,20 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,49 1,49 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 1,50 2,11 2,10 2,10 2,10 2,10 2,10 2,10 399,13 400,51 401,71 402,80 403,87 404,95 406,03 407,06 408,00 408,80 409,43 374,21 375,73 378,20 379,22 380,13 381,76 383,93 387,34 390,58 393,37 395,65 397,57 399,17 400,51 401,73 402,97 404,32 405,75 407,13 408,34 409,28 409,93 366,76 368,37 371,03 372,17 373,24 375,30 378,26 505,25 507,10 508,83 510,54 512,23 513,88 515,45 516,87 518,09 519,05 519,75 478,35 480,38 483,39 484,57 485,64 487,54 490,17 494,56 498,72 502,00 504,47 506,61 508,55 510,30 511,89 513,41 514,93 516,48 518,01 519,46 520,79 521,99 475,25 476,43 479,13 480,52 481,88 484,38 487,51 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 1,50 1,60 1,80 1,90 2,00 2,20 2,50 3,00 3,50 4,00 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 10,00 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 382,00 386,00 389,00 391,00 394,00 396,00 398,00 399,00 401,00 402,00 404,00 405,00 406,00 408,00 409,00 367,00 368,00 371,00 372,00 373,00 375,00 377,00 381,00 384,00 387,00 389,00 391,00 393,00 395,00 397,00 398,00 399,00 401,00 402,00 403,00 404,00 405,00 492,00 496,00 499,00 502,00 504,00 506,00 508,00 510,00 512,00 514,00 515,00 517,00 518,00 519,00 520,00 476,00 477,00 480,00 481,00 482,00 485,00 488,00 492,00 495,00 499,00 501,00 504,00 506,00 508,00 509,00 511,00 513,00 514,00 515,00 517,00 518,00 519,00 134 3,00 3,50 3,97 4,50 5,00 5,50 6,00 6,50 7,00 7,50 8,00 8,50 9,00 9,50 9,97 1,50 1,60 1,80 1,90 2,00 2,21 2,50 3,00 3,50 4,00 4,55 5,00 5,50 6,01 6,50 6,99 7,49 8,00 8,50 9,00 9,50 10,00 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,10 2,09 2,09 2,10 2,10 2,11 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 2,40 382,58 386,03 388,71 391,15 393,46 395,65 397,71 399,50 401,07 402,66 404,32 405,77 406,68 406,98 406,83 367,02 368,16 370,62 371,79 372,90 374,93 377,57 381,27 384,18 386,74 389,39 391,67 393,56 395,03 396,20 397,33 398,63 400,15 401,74 403,19 404,38 405,29 491,79 495,65 498,70 501,61 504,33 506,49 508,27 509,96 511,79 513,72 515,56 517,10 518,26 519,20 520,19 475,48 477,17 480,03 481,27 482,43 484,59 487,60 492,25 495,56 498,09 501,11 503,57 505,75 507,82 509,72 511,43 512,97 514,32 515,48 516,50 517,50 518,67