Phát triển phương pháp giải phương trình đạo hàm riêng bằng công nghệ Mạng nơ ron tế bào
3.3.5. Mô phỏng thuật toán trên Matlab
Để mô phỏng thuật toán, ta phải chọn các giá trị tham số phù hợp với thực tế đảm bảo kết quả đúng với các hiện tượng vật lý xảy ra, đây là một việc phải đo đạc, khảo sát rất khó khăn. Giả sử các tham số có giá trị như trong bảng 3.3
Bảng 3.3 Các tham số vật lý cho bài toán thủy lực hai chiều
Tham số Giá trị
g 9,8 (m/s2)
Kx=Ky 20 (s-1)
x=y 100 (m)
Mô phỏng việc tính toán theo thuật toán CNN, ta nhập các giá trị tham số, giá trị ban đầu và mẫu sau đó thực hiện tính toán. Trên cửa sổ ta có thể chọn các giá trị nhập và kiểu điều kiện biên (ta có thể chọn một trong ba kiểu biên). Tuy nhiên ở đây mô phỏng để kiểm tra thuật toán nên chỉ tính toán cho lưới 20x20 điểm, nếu triển khai trên chip CNN thì số điểm sẽ nhiều hơn. Thuật toán thực hiện như sau:
Bước 1: Đưa các giá trị khởi tạo (điều kiện ban đầu);
Bước 2: Tính toán các mẫu;
Bước 3: Thực hiện các phép tính toán;
Bước 4: Đưa ra kết quả tính toán;
Bước 5: Kết thúc thuật toán.
Bước tính toán phải dùng vòng lặp để tính lần lượt cho mọi điểm trong lưới sai phân (MxN) điểm. Giao diện để nhập các tham số và mẫu như trong Hình 3.26
Giả sử, giá trị ban đầu của các hàm như minh hoạ trong Hình 3.27.
Việc tính toán trên máy PC dù mô phỏng theo thuật toán CNN nhưng không thể hiện tốc độ tính toán mà chỉ cho kết quả tính toán của mô hình. Sau thời gian xử lý chúng ta nhận được kết quả tính toán sau thời gian t (Hình 3.28).