Header Page of 161 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC NÔNG LÂM THÀNH PHỐ HỒ CHÍ MINH  TIỂU LUẬN TỐT NGHIỆP NGHIÊN CỨU VỀ XỬ LÝ SONG SONG TRONG GIS VÀ XÂY DỰNG ỨNG DỤNG SONG SONG HÓA THUẬT TOÁN ĐỊNH DÕNG CHẢY TRÊN BỀ MẶT Sinh viên thực hiện: TRẦN CÔNG HUẤN Ngành: HỆ THỐNG THÔNG TIN MÔI TRƢỜNG Niên khóa: 2010 – 2014 Thành phố Hồ Chí Minh, Tháng năm 2014 Footer Page of 161 Header Page of 161 NGHIÊN CỨU VỀ XỬ LÝ SONG SONG TRONG GIS VÀ XÂY DỰNG ỨNG DỤNG SONG SONG HÓA THUẬT TOÁN ĐỊNH DÕNG CHẢY TRÊN BỀ MẶT Tác giả TRẦN CÔNG HUẤN Giáo viên hƣớng dẫn: ThS KHƢU MINH CẢNH Tháng năm 2014 [I] Footer Page of 161 Header Page of 161 LỜI CẢM ƠN Trong trình thực đề tài này, nhận đƣợc hƣớng dẫn giúp đỡ thầy hƣớng dẫn, quý thầy cô thuộc khoa môi trƣờng tài nguyên nói chung môn thông tin địa lý ứng dụng nói riêng thuộc trƣờng đại học Nông Lâm Thành phố Hồ Chí Minh Qua xin gửi lời cảm ơn chân thành tới: - ThS.Khƣu Minh Cảnh, công tác Trung tâm Ứng dụng Hệ thống Thông tin Địa lý – Sở Khoa học Công nghệ TP.HCM, ngƣời hƣớng dẫn trực tiếp trình làm đề tài - Thầy Phó giáo sƣ, tiến sĩ Nguyễn Kim Lợi, chủ nhiệm Bộ môn Thông tin địa lý ứng dụng - Trƣờng Đại học Nông lâm Tp HCM - KS Lê Hoàng Tú Kĩ sƣ GIS - Nghiên cứu viên, Trung tâm Nghiên cứu Biến đổi Khí hậu, Trƣờng Đại học Nông Lâm TP.HCM - KS Nguyễn Duy Liêm Kĩ sƣ GIS - Giảng viên Khoa Môi trƣờng – Tài nguyên, Trƣờng Đại học Nông Lâm TP.HCM - Đặc biệt xin cảm ơn đến ba mẹ sinh thành, nuôi nấng giáo dục đến ngày hôm Gia đình bạn bè động viên giúp đỡ, tạo điều kiện thuận lợi cho trình học tập nhƣ thời gian thực đề tài, đồng thời cho gửi lời cảm ơn sâu sắc tới Sở Khoa Học Công Nghệ Thành Phố Hồ Chí tạo điều kiện cho thực tập, giúp đỡ, cung cấp số liệu giúp thực tốt đề tài TP.HCM, ngày tháng năm 2014 Trần Công Huấn Bộ môn Tài nguyên GIS Khoa Môi trƣờng & Tài nguyên Trƣờng Đại học Nông Lâm Tp.Hồ Chí Minh Footer Page of 161 [II] Header Page of 161 TÓM TẮT Đề tài nghiên cứu “Nghiên cứu xử lý song song GIS xây dựng ứng dụng song song hóa thuật toán định dòng chảy bề mặt” đƣợc thực hoàn thành Phòng kỹ thuật – Trung tâm ứng dụng Hệ Thống Thông Tin Địa Lý – Sở Khoa học Công nghệ Thành phố Hồ Chí Minh Thời gian thực từ 2/3/2014 đến 31/5/2014 Gồm nội dung nhƣ sau: - Tìm hiểu mô hình, công nghệ, họ máy tính song song - Tìm hiểu thuật toán song song, cấu trúc thuật toán song song - Tìm hiểu thuật toán Floyd – Warshall thuật toán tìm tích lũy dòng chảy - Tìm hiểu công cụ tìm dòng ArcGis - Tìm hiểu thuật toán xác định hƣớng dòng chảy theo D8 cài đặt ngôn ngữ C sharp - Tìm hiểu dạng mô hình liệu DEM - Tìm hiểu gói phát triển ArcEngine ArcGis Visual studio Trên xây dựng công cụ hỗ trợ hiển thị, cập nhật, phân tích liệu chuyển dạng liệu Các kết thu đƣợc: - Xây dựng đƣợc công cụ xác định hƣớng dòng chảy (theo D8) tính tích lũy dòng chảy bề mặt - Ứng dụng khai thác tốc độ tối đa máy tính cho kết nhanh nhằm tiết kiệm thời gian cho ngƣời sử dụng - Xây dựng công cụ chuyển liệu dạng file text sang dạng raster hiển thị liệu raster lên form ứng dụng Footer Page of 161 [III] Header Page of 161 MỤC LỤC LỜI CẢM ƠN II TÓM TẮT III MỤC LỤC IV DANH MỤC TỪ VIẾT TẮT VII DANH MỤC BẢNG BIỂU VIII DANH MỤC HÌNH ẢNH IX PHẦN MỞ ĐẦU 1.1 Tính cấp thiết đề tài 1.2 Mục tiêu nghiên cứu đề tài 1.3 Đối tƣợng phạm vi nghiên cứu 1.3.1 Đối tƣợng nghiên cứu 1.3.2 Phạm vi nghiên cứu PHẦN TỔNG QUAN VỀ THUẬT TOÁN VÀ TÍNH TOÁN SONG SONG 2.1 Đại cƣơng tính toán song song 2.1.1 Một số khái niệm thuật ngữ 2.1.2 Các mức độ song song (Level of parallelism) .4 2.1.3 Phân loại kiến trúc song song .6 2.1.4 Mô hình SIMD (PRAM) .8 2.1.5 Dùng công nghệ EREW mô kiến trúc CRCW, CREW 2.1.6 Họ máy MIND .10 2.1.6.1 Hệ đa xử lý với nhớ phân tán (Multi processor system with Distributed Memory) 11 Footer Page of 161 [IV] Header Page of 161 2.1.6.2 Hệ đa xử lý dùng chung nhớ (Multi processor system with Shared Memory) 12 2.1.6.3 Hệ đa xử lý với nhớ dùng chung phân tán (Multi processor system with distributed shared memory) 13 2.1.7 2.2 Ngôn ngữ mô tả thuật toán song song 13 Các mô hình tính toán song song minh họa .15 2.2.1 Mô hình nhị phân (Bimary Tree Model) 15 2.2.2 Mô hình mạng .19 2.2.3 Thuật toán k-cube-Min 23 2.2.4 Thuật toán song song tính tích ma trận 24 2.2.5 Đánh giá hiệu thuật toán song song 26 2.3 Tính toán song song NET minh họa 34 2.3.1 Task 36 2.3.2 Vòng lặp song song (Parallel Loops) 37 2.3.3 Parallel LINQ 38 2.4 Thuật toán Floyd – Warshall toán tìm đƣờng ngắn cặp đỉnh đồ thị .38 PHẦN DỮ LIỆU, NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU .41 3.1 Dữ liệu .41 3.1.1 Mô hình liệu DEM 41 3.1.2 File text độ cao 43 3.2 Thuật toán định dòng chảy bề mặt địa hình 44 3.2.1 Giới thiệu thuật toán phân tích dòng chảy D8 44 3.2.1.1 Xác định hƣớng dòng chảy 44 3.2.1.2 Tính toán tích lũy dòng chảy 45 3.2.2 Giới thiệu công cụ tìm dòng ArcGIS .46 Footer Page of 161 [V] Header Page of 161 3.3 3.2.2.1 ArcSWAT .46 3.2.2.2 Bộ công cụ tìm dòng chảy tích lũy ArcGIS .49 Cài đặt thuật toán D8 (tuần tự) 52 3.3.1 Đọc liệu (đọc file text độ cao) .52 3.3.2 Xác định hƣớng dòng chảy theo D8 53 3.3.3 Tính toán tích lũy dòng chảy (D8) 56 3.4 Tại phải cài đặt thuật toán song song .59 3.5 Cài đặt thuật toán song song D8 63 3.5.1 Đọc liệu 64 3.5.2 Xác định song song hƣớng dòng chảy theo D8 65 3.5.3 Tính toán song song tích lũy dòng chảy theo D8 .67 PHẦN CÁC KẾT QUẢ NGHIÊN CỨU 74 4.1 Giới thiệu liệu thử nghiệm 74 4.2 Nhóm công cụ xây dựng chƣơng trình 75 4.3 Các kết thực đƣợc ứng dụng phân tích hƣớng dòng chảy76 PHẦN KẾT LUẬN VÀ KIẾN NGHỊ .87 5.1 Kết luận 87 5.2 Kiến nghị 87 PHẦN TÀI LIỆU THAM KHẢO .88 PHỤ LỤC .90 Footer Page of 161 [VI] Header Page of 161 DANH MỤC TỪ VIẾT TẮT GIS: Geographic Information System (Hệ thống thông tin địa lý) ESRI: Economic and Social Research Institute (Viện nghiên cứu hệ thống môi trƣờng) TP.HCM: Thành phố Hồ Chí Minh DEM: Digital Elevation Models (Mô hình độ cao số) CPU: Central Processing Uint (Bộ xử lý trung tâm) Thuật toán D8 (Thuật toán xác định hƣớng dòng chảy đơn) Footer Page of 161 [VII] Header Page of 161 DANH MỤC BẢNG BIỂU Bảng 2.2.2: so sánh số đặc trƣng k- cute với đồ thị đầy đủ .21 Bảng 2.2.5: So sánh thời gian thuật toán Boolean-AND, Boolean-AND – 1, Boolean-AND - 26 Bảng 3.2.1.1: Hƣớng dòng chảy tính lƣu vực 45 Bảng 3.2.1.2: Sự tích lũy dòng chảy lƣu vực 46 Bảng 3.4: Thống kê thời gian phép toán song song 62 Footer Page of 161 [VIII] Header Page 10 of 161 DANH MỤC HÌNH ẢNH Hình 1.1: DEM raster Hình 1.2: Cấp độ xám DEM .1 Hình 2.1.2 Các mức độ song song Hình 2.1.3 Phân loại kiến trúc song song .7 Hình 2.1.6.1: Hệ đa xử lý với nhớ phân tán .11 Hình 2.1.6.2: Hệ đa xử lý dùng chung nhớ 12 Hình 2.1.6.3: Hệ đa xử lý với nhớ dùng chung phân tích 13 Hình 2.2.1.1: Mô hình nhị phân cộng số 16 Hình 2.2.1.2: Cây nhị phân thực tính toán 18 Hình 2.2.3: Tính tổng với xử lý .18 Hình 2.2.2.1: Mạng 3-cube 20 Hình 2.2.2.2: Phân bổ đầu vào 22 Hình 2.3: Mô tả thuật toán song song NET Framework 4.0 .35 Hình 3.1.1.1 Quy trình chuyển đổi DEM .41 Hình 3.1.1.2: Những điểm lỗi có DEM 42 Hình 3.1.1.3: Hàm fill ArcToolbox 42 Hình 3.1.2.1: Công cụ Raster to Ascii ArcToolbox .43 Hình 3.1.2.2: Text file dùng để xử lý 44 Hình 3.2.2.2.1: Sơ đồ tạo DEM từ đồ địa hình ArcGis 49 Hình 3.2.2.2.2: Sơ đồ tính hƣớng, tích lũy tìm dòng ArcGis 50 Hình 3.2.2.2.3 Sơ đồ tìm liên kết dòng cửa xả ArcGis 51 Hình 3.4.1: Hiệu suất CPU thực phép toán 62 Hình 3.4.2: Hiệu suất CPU thực phép toán song 63 Hình 3.5.1: Chia liệu thành mảng 64 Hình 3.5.1: Form chuyển liệu sang dạng raster 72 Hình 4.1: Bộ liệu thử nghiệm phần mềm phân tích dòng chảy 74 Hình 4.2 :Form phần mềm phân tích song song dòng chảy theo D8 75 Hình 4.3.1 Mở file text độ cao .76 Footer Page 10 of 161 [IX] Header Page 92 of 161 Kết tích lũy hiển thị Textbox Hình 4.3.8: Kết tích lũy in Textbox Lƣu kết tích lũy, nhƣ phần tính thuật toán D8 ta chọn nút lƣu file text kết khung thao tác xuất hộp thoại cho ta chọn ổ đĩa lƣu lƣu tên file tex kết Xong bƣớc ta đƣợc file text kết với tên đƣờng dẫn nhƣ ta lƣu Để chuyển file text kết sang dạng raster ta làm nhƣ sau Trong khung thao tác ta chọn chuyển dạng raster để chuyển đổi định dạng file text sang dạng raster Chọn chuyển dạng raster Hình 4.3.9: Chọn nút chuyển file text sang dạng raster Footer Page 92 of 161 [81] Header Page 93 of 161 Lúc phần mềm chuyển sang form ConvertToRaster để ta dễ dàng thực Trong Form ta làm theo bƣớc đƣợc dẫn Các bƣớc nhƣ sau: Nhấn vào chọn file text cần chuyển Nhấn vào chọn nơi lƣu raster sau chuyển Nhấn vào để chƣơng trình bắt đầu chuyển Nhấn vào để trở fom Hình 4.3.10: Sử dụng chức chuyển sang dạng raster Thực theo bƣớc, chƣơng trình tự chuyển liệu sang dạng raster, chuyển xong phần mềm báo “OK” ta cần vào đƣờng dẫn mà ta chọn trƣớc để xem kết Sau convert xong, để xem kết dạng raster form ta làm nhƣ sau Trên form khung công cụ ta chọn nút hiển thị kết Để hiển thị kết ta chọn nút hiển thị kết Hình 4.3.11: Chọn nút hiển thị raster lên Form Footer Page 93 of 161 [82] Header Page 94 of 161 Sau ta chọn nút hiển thị kết quả, phần mềm tự động chuyển sang Form Trong form ta làm nhƣ sau Trong form ta chọn nút để add DEM vào form Xuất hộp thoại ta chọn DEM hộp thoại Chọn nút để add DEM Chọn DEM convert trƣớc để hiển thị Chọn nơi lƣu DEM mục Nhấn nút Open để mở DEM Chọn hiển thị file raster Hình 4.3.12: Chọn file raster hiển thị lên form Nhấn nút Open, phần mềm thông báo nhấn OK Kết thuật toán tính hƣớng D8 liệu độ cao trƣờng đại học Nông Lâm Tp.HCM dƣới dạng raster hiển thị Form nhƣ sau: Footer Page 94 of 161 [83] Header Page 95 of 161 Hình 4.3.13: Kết thuật toán D8 dạng raster Kết hiển thị hƣớng theo D8, nhìn rõ form này, ô có hƣớng hiển thị cấp độ xám, theo hƣớng đƣợc xếp theo thứ tự lớn dần từ đến 128, hƣớng lớn có cấp độ xám nhỏ Ngƣợc lại ô nhỏ có cấp độ xám lớn (màu đen) Còn kết thuật toán tìm tích lũy dòng chảy liệu DEM trƣờng đại học Nông Lâm Tp.HCM Footer Page 95 of 161 [84] Header Page 96 of 161 Hình 4.3.14: Kết tích lũy dạng raster hiển thị form Kết cho ta thấy rõ dòng tích lũy Màu tối thể cho nơi ô chảy Màu sáng thể điểm tập trung nhiều ô chảy góc dƣới bên trái có tọa độ thể rõ vị trí điểm ta muốn xét cách rê chuột lên vị trí muốn xét, điều đặc biệt Dựa ứng dụng ta nhận thấy rõ dòng chảy lƣu vực lớn mà ta muốn xét, từ ta áp dụng cho công trình nhƣ xây dựng đƣờng xá, cống rãnh hay vùng thƣờng xuyên có lũ lụt xảy ra, từ ta ý đặc biệt cảnh báo cho vùng vị trí Bảng so sánh kết thởi gian thực thuật toán đƣợc thực core i5 – 3230M (4 Cpu – 2.6 Ghz) Footer Page 96 of 161 [85] Header Page 97 of 161 Thuật toán Thuật toán Thuật toán song song Tính hƣớng D8 giây 626 giây 260 Tính tích lũy phút 53 giây 243 phút 20 giây 854 Tổng thời gian phút 55 giây 869 phút 22 giây 124 Phép toán Bảng 4.3: So sánh thời gian thực thuật toán song song Bảng 4.3 so sánh thời gian thực thuật toán song song kết cho ta thấy khác biệt, thời gian tiết kiệm đƣợc thuật toán song song lớn so với thuật toán Đặc biệt, liệu lớn máy đƣợc trang bị lõi cpu đời cho ta kết tốt hơn, khả quan hơn, thời gian tiết kiệm đƣợc lên đến hàng đồng hồ, chí với liệu lớn siêu lớn kết tăng lên gấp bội, lên đến hàng ngày, hàng tuần Footer Page 97 of 161 [86] Header Page 98 of 161 PHẦN KẾT LUẬN VÀ KIẾN NGHỊ 5.1 Kết luận Kết quả, nghiên cứu xây dựng đƣợc ứng dụng phân tích song song thuật toán định hƣớng dòng chảy bề mặt với tính ƣu điểm sau: - Ứng dụng có khả đọc mở file text ArcGis, lƣu liệu kết tính đƣợc vào file text riêng để ta dễ dàng theo dõi tiếp tục tính toán sau - Ứng dụng có khả sử dụng file text ArcGis nhập vào để xác định hƣớng dòng chảy (theo D8) - Ứng dụng có khả dùng kết hƣớng dòng chảy xác định trƣớc để tính toán tích lũy dòng chảy ô - Ứng dụng có khả chuyển liệu dạng file text ArcGis sang liệu dạng raster hiển thị khác kiểu liệu nhƣ liệu đồ (shapefile), geodatabases, đặc biệt raster lên form ứng dụng - Ứng dụng đƣợc sử dụng điều chỉnh thuật toán song song cho ứng dụng khai thác tối đa khả máy tính, luôn đạt tốc độ lớn nhất, tiết kiệm thời gian ngắn cho ngƣời sử dụng 5.2 Kiến nghị Kết tính tích lũy ứng dụng cho ta tranh tổng quát tích lũy dòng chảy lƣu vực lớn, kết hình thành đƣờng tích lũy mà dòng chảy qua Dựa vào ta áp dụng xây dựng cho công trình giao thông, cầu cống cho hiệu cao cảnh báo vùng có nguy lũ lụt tránh rủi ro ngƣời Ứng dụng đƣợc phát triển lên thành phiên với công cụ lấy kết từ tính hƣớng D8 tích lũy để xây dựng thêm nhƣ công cụ tìm dòng (Stream line), công cụ tìm liên kết dòng (Stream link), tìm cửa xả (watershed) để ứng dụng thêm đa dạng phục vụ đƣợc nhiều ngành, nhiều lĩnh vực từ đời sống nguòi ngày đƣợc nâng cao Footer Page 98 of 161 [87] Header Page 99 of 161 PHẦN TÀI LIỆU THAM KHẢO Tiếng việt: - [1] Pgs Ts Nguyễn Kim Lợi, ThS Lê Cảnh Định, Ths Trần Thống Nhất, 2009, Hệ thống thông tin địa lý nâng cao, NXB Nông Nhiệp - [2] Ks Nguyễn Duy Liêm, 2013,Chuyên đề Swat, biên soạn - [3] Ks Lê Minh Hoàng, 2012, Giải thuật lập trình, giảng chuyên đề, Đại học Sƣ phạm Hà Nội - [4] Pgs Ts Nguyễn Đức Nghĩa, 2008, Tính toán song song, giảng môn học, NXB Bác Khoa Hà Nội Tiếng Anh - [5] Ole Mark,Terry van Kalken,K Rabbi, Jesper Kjelds, A mouse GIS study of the drainage in Dhake city - [6] Peter Van Capelleveen, Urban drainage network modeling better analyzed using ArcView 3D analyst - [7] Jurgen Garbrecht and Lawrence W Martz, June 1997 The assignment of drainage direction over flat surfaces in raster digital elevation models Journal of Hydrology, 193:204–213 - [8] D.M Mark, 1988, Modelling in Geomorphological Systems, chapter Network models in geomorphology - [9] John O’Callaghan and David Mark, December 1984, The extraction of drainage networks from digital elevation data Computer Vision, Graphics, and Image Processing - [10] Chase Wallis, Dan Watson, David Tarboton, and Robert Wallace, 2009, Parallel Flow-Direction and Contributing Area Calculation for Hydrology Analysis in Digital Elevation Models In International Conference on Parallel and Distributed Processing Techniques and Applications Website - [11] www.esri.com Footer Page 99 of 161 [88] Header Page 100 of 161 - [12] www.msdn.microsoft.com - [13] www.hcmuaf.edu.vn - [14] www.hochiminhcity.gov.vn - [15] www.stackoverflow.com - [16] www.swat.tamu.edu - [17] www.proceedings.esri.com - [18] www.help.arcgis.com Footer Page 100 of 161 [89] Header Page 101 of 161 PHỤ LỤC Phụ lục 1: Lý thuyết tìm đƣờng ngắn Hiện tại, nhiều thuật toán song song tìm đƣờng ngắn đồ thị Chúng ta nêu thuật toán, cụ thể thuật toán tiến sĩ Lyudmil Aleksandrov đề xuất Giả sử, có đồ thị G, n đỉnh Khi đó, thuật toán là:  Giai đoạn tiền xử lý: - Bƣớc 0: chia đồ thị n đỉnh thành r đồ thị Mỗi đồ thị có n/r đỉnh đồ thị có - n đỉnh biên đồ thị r Bƣớc 1: Tải đồ thị vào p xử lý Hình PL1: Hình minh họa việc phân bố r=5 đồ thị vào p=3 xử lý - Bƣớc 1: (tiếp theo) Sau đó, xử lý chuyển toàn đỉnh biên đến tất xử lý khác để hình thành đồ thị G* có cạnh nhƣ G (nhƣng có đỉnh biên ngoài) - Bƣớc 2: Tính toán xử lý Mỗi xử lý tính toán đƣờng ngắn Footer Page 101 of 161 [90] Header Page 102 of 161 - Bƣớc 3: Mỗi xử lý trau đổi thông tin biên Khi giá trị nhỏ đƣợc cập nhật đồ thị G* - Bƣớc 4: Mỗi xử lý tính toán giá trị ngắn từ đỉnh lại đến đỉnh đồ thị G*  Giai đoạn xử lý: Khi có yêu cầu tìm đƣờng ngắn nhất: - Bƣớc 1: Bộ xử lý trung tâm xác định xác vị trí đỉnh xử lý chứa đỉnh nguồn đích Hình PL2: Nhập điểm nguồn điểm đích - Bƣớc 2: Tính giá trị trọng số ngắn đỉnh đồ thị chứa điểm nguồn đến đồ thị chứa điểm đích Hình PL3 : Tìm khoảng cách ngắn tập điểm biên đồ thị chứa điểm nguồn đến đồ thị chứa điểm đích Footer Page 102 of 161 [91] Header Page 103 of 161 - Bƣớc 3: Đƣờng đƣợc hình thành cách kết hợp đƣờng ngắn từ điểm nguồn đến đỉnh đồ thị chứa điểm nguồn với đƣờng tìm đƣợc bƣớc đƣờng ngắn từ điểm đích đến đỉnh đồ thị chứa điểm đích Kết thúc thuật toán Hình PL4 : Xác định đƣợc đƣờng ngắn kết thúc thuật toán Để cài đặt đƣợc thuật toán trên, phải xây dựng đối tƣợng phƣơng thức nhƣ sau: - Đối tƣợng đồ thị G, - Xác định đối tƣợng đỉnh trong, đỉnh đồ thị - Đồ thị G* (bao gồm tập đỉnh đồ thị cạnh đồ thị G), - Các node liên kết Thuật toán đƣợc chứng minh vấn đề sau: - Thuật toán giải đƣợc vấn đề tìm đƣờng ngắn - Việc truyền thông node (bộ xử lý) cực tiểu - Từ đây, với đồ, xây dựng thành đồ thị với đỉnh giao đƣờng cạnh đƣờng Đối với liệu raster, phân chia theo không gian Footer Page 103 of 161 [92] Header Page 104 of 161 Phụ lục 2: Thuật toán Fox nhân ma trận Phép toán nhân ma trận phép toán có số lƣợng tính toán lớn Với hai ma trận A=(aij) B = (bij) với kích thƣớc giả định nxn Tích ma trận A.B = C = (cij) ma trận kích thƣớc nxn với phần tử ma trận C đƣợc xác định nhƣ sau: cij = ∑aikbkj, với k=0,n-1 i=0,1,…,n-1 j = 0,1,…,n-1 Với thuật toán nhân ma trận đƣợc cài đặt theo định nghĩa, ta có độ phức tap thuật toán O(n3) Để thực thuật toán Floyd, ta thực công việc sau: - Ma trận a b ma trận kề đồ thị - Thay Phép nhân Phép cộng - Thay Phép cộng phép lấy giá trị cực tiểu (min) Vì phép nhân ma trận cần tính toán cho n2 phần tử, đó, ta cố gắng thực thuật toán giảm mức độ phức tạp thành O(n2) Thuật toán Fox nhân ma trận với ý tƣởng thực n bƣớc (từ bƣớc đến bƣớc n-1) Trong bƣớc, giá trị cij đƣợc cập nhật Cụ thể nhƣ sau: - Bƣớc 0: cij  aii x bij, - Bƣớc k: (1≤k