›

Bài giảng Tin học ứng dụng trong Hóa học - Chương 6: Mô hình hóa và tối ưu hóa

pdf 54 trang huongle 3090
Download
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Tin học ứng dụng trong Hóa học - Chương 6: Mô hình hóa và tối ưu hóa", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfbai_giang_tin_hoc_ung_dung_trong_hoa_hoc_chuong_6_mo_hinh_ho.pdf

Nội dung text: Bài giảng Tin học ứng dụng trong Hóa học - Chương 6: Mô hình hóa và tối ưu hóa

  1. TRƢỜNG ĐẠI HỌC CÔNG NGHIỆP TP.HCM KHOA CÔNG NGHỆ HOÁ HỌC TIN HỌC ỨNG DỤNG TRONG HÓA HỌC CHƢƠNG 6: MÔ HÌNH HÓA & TỐI ƢU HÓA
  2. NỘI DUNG 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH 2. BÀI TOÁN TRUYỀN NHIỆT 3. BÀI TOÁN CHƢNG CẤT 4. TỐI ƢU HÓA THỰC NGHIỆM
  3. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH HPT tuyến tính n phương trình, n ẩn số: a11.x1 + a12.x2 + + a1n.xn = b1 a21.x1 + a22.x2 + + a2n.xn = b2 . an1.x1 + an2.x2 + + ann.xn = bn Hay viết dưới dạng ma trận: A. X = B HPT có nghiệm khi det A 0, khi đó nghiệm của hệ xác định theo phương pháp ma trận : X = A-1.B
  4. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Một số hàm trong Excel: Tính định thức ma trận A: MDETERM(A) Tìm ma trận nghịch đảo A-1: MINVERSE(A) Nhân 2 ma trận A-1 và B: MMULT(A-1,B)
  5. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Giải hệ phương trình: 2,75X1 + 1,78X2 + 1,11X3 = 13,62 3,28X1 + 0,71X2 + 1,15X3 = 17,98 1,15X1 + 2,70X2 + 3,58X3 = 39,72 giải bằng phương pháp ma trận
  6. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 1: lập bảng số liệu
  7. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 2: tính det (A)
  8. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 3: Tính ma trận A-1 Ấn ba phím đồng thời Shift + Ctrl + Enter
  9. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 4: Tính nghiệm X Ấn ba phím đồng thời Shift + Ctrl + Enter
  10. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 5: Nhập giá trị XT và tính BT
  11. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 6: tính BT
  12. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Phƣơng pháp ma trận Bước 7: tính X nhận kết quả nghiệm.
  13. 1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Kết quả giải bằng Solver:
  14. BÀI TẬP 1. Giải các phương trình sau với độ chính xác 10-5: a. ln(8x) – x – 0,5 =0 b. ln (7x) – 3x + 1 = 0 c. ln(6x) – x – 0,4 = 0 2. Giải phương trình sau: a. x6 + 4x4 - 3x - 5 = 0 với x thuộc đoạn [1;2] b. x5 + 5x – 2 = 0 với x thuộc đoạn [0;1]
  15. 2. BÀI TOÁN TRUYỀN NHIỆT Thành lò đốt có 3 lớp: . Gạch chịu nhiệt dày 120 mm, hệ số dẫn nhiệt 0,81 W/ m.K . Gạch cách nhiệt dày 65 mm, hệ số dẫn nhiệt 0,23 W/ m.K . Thép chịu lực dày 10 mm, hệ số dẫn nhiệt 45 W/ m.K Nhiệt độ lò: 8000C, hệ số cấp nhiệt trong lò: 69,6 W/ m2.K Nhiệt độ không khí: 350C, hệ số cấp nhiệt không khí: 13,9 W/ m2.K Yêu cầu: Xác định các nhiệt độ bề mặt các lớp.
  16. 2. BÀI TOÁN TRUYỀN NHIỆT Thép chịu lực Gạch chịu 3 nhiệt 1 2 Gạch cách nhiệt 1= 0,12m, 1= 0,81 W/ m.K 0 2 t1= 800 C, 1=69,6 W/ m .K 2= 0,065 m, 1= 0,23 W/ m.K 0 2 t2= 35 C, 2=13,9 W/ m .K 3= 0,01 m, 3= 45 W/ m.K
  17. 2. BÀI TOÁN TRUYỀN NHIỆT 1 1 q1= 1(t1-tT1)  q2= (tT1-tT2) 1 2   q4= (tT3-tT4) q3= (tT2-tT3)  3 2 3 2 q5= 2(tT4-t2) Lƣu ý: q1 = q2 = q3 = q4 = q5 = q
  18. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Bước 1: Lập bảng tính như sau:
  19. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Bước 2: Nhập giá trị tT1 và tính q1
  20. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Bước 3: tính giá trị tT2, tT3, tT4 thông qua q1
  21. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Bước 4: tính giá trị q5 thông qua tT4 và t2
  22. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Bước 5: Lập biểu thức so sánh q5 với q1
  23. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Bước 6: Thay đổi tT1 để biểu thức so sánh 2%
  24. 2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Kết quả:
  25. 2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 1: chuyển điều kiện về dạng hệ phương trình tuyến tính: q1= 1(t1-tT1) q + 1.tT1 = 1.t1  q2= (tT1-tT2) 1  .q-  .t + .t =0 1 T1 T2  q3= (tT2-tT3) 2 .q -  .tT2 + .tT3 =0 2  .q -  .t + .t =0  3 T3 T4 q4= (tT3-tT4) 3 q - 2.tT4 = 2.t2 q5= 2(tT4-t2)
  26. 2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 2: lập bảng số liệu:
  27. 2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 3: nhập giá trị biến và tính VT
  28. 2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 4: gọi hàm Solver
  29. 2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 5: cài đặt VT = VP
  30. 2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Kết quả nghiệm:
  31. 3. BÀI TOÁN CHƢNG CẤT Chưng cất hỗn hợp nước – Acetic với nồng độ nhập liệu 20% (KL), sản phẩm đỉnh 95%(kl), sản phẩm đáy 0,5% (kl), chỉ số hồi lưu bằng 4. Yêu cầu: Xác định số đĩa lý thuyết.
  32. 3. BÀI TOÁN CHƢNG CẤT Các công thức sử dụng: • Chuyển nồng độ phân mol: 푿 푴 xA = 푿 −푿 + 푴 푴
  33. 3. BÀI TOÁN CHƢNG CẤT Các công thức sử dụng: • Phương trình đoạn cất: 푅 y = + = 0,8x + 19,698 L 푅+1 푅+1 • Phương trình đoạn chưng: 푅+ −1 y = − x = 1,2419x – 0,3986 c 푅+1 푅+1 w • Phương trình đường cân bằng (xây dưng từ số liệu thực nghiệm): y*=4.10-5x3-0,0106x2+1,6853x-0,5523 f = F/ D
  34. 4. TỐI ƢU HÓA THỰC NGHIỆM NỘI DUNG 4.1 PP QUI HOẠCH THỰC NGHIỆM 4.2 KIỂM TRA SỰ TƢƠNG HỢP CỦA PTHQ 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL
  35. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM Mô hình thực nghiệm của đối tượng nghiên cứu: Biến đầu vào: nồng độ, áp suất, nhiệt độ, độ pH, Hàm mục tiêu: hiệu suất, chất lượng sản phẩm, CPSX,
  36. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM Ta cần thiết lập quan hệ: y = f(x1, x2, ,xk) +  hay y = f(X) +  Triển khai hàm dưới dạng chuỗi Taylor: k k k 2 f(x1, x2, . . . , xk) = β0 + βj xj + βjj xj + βij xixj+. . . j=1 j=1 j≠i Với 1 ≤ i ≤ j ≤ k PTHQ thực nghiệm: y = f x1, x2, . . . , xk
  37. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM Bố trí thực nghiệm theo ma trận biến đầu vào - ra: 1 x11 x12 x1k y1 1 x x x2k y X = 21 22 Y = 2 ⋮ 1 xn1 xn2 xnk yn (n dòng, k+1 cột) b0 b Ma trận các hệ số hồi qui tuyến tính có: B = 1 ⋮ bk Theo phương pháp bình phương cực tiểu, ta có: −1 XTX . B = XTY B = XTX . (XTY) Với: XT là ma trận chuyển vị của ma trận X (XTX)-1 là ma trận nghịch của ma trận XTX
  38. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM Ma trận trực giao X có những tính chất sau: + Tính trực giao: tích vô hướng của hai vectơ cột bất kỳ của X bằng 0. n ximxij=0 với j, m= 0, k i=1 + Tính chất đối xứng: tổng các phần tử trong một cột bất kỳ đều bằng 0. n xij=0 với j ≠ 0. i=1
  39. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM Một vài phương pháp qui hoạch thực nghiệm: + Qui hoạch thực nghiệm yếu tố toàn phần 2k ( k là các yếu tố, n là số mức thì số thí nghiệm N = nk) Ma trận X có thêm tính chất chuẩn hóa: n 2 xij=N với j= 0, k i=1 + Qui hoạch thực nghiệm yếu tố từng phần 2k – p (p là giá trị đặc trưng cho độ từng phần)
  40. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM + Qui hoạch trực giao cấp 2: xây dựng ma trận trực giao X bao gồm ba loại thí nghiệm: - Phần cơ sở gồm n = 2k thí nghiệm theo qui hoạch thực nghiệm yếu tố toàn phần. - Phần điểm “*” gồm nk = 2k điểm nằm trên các trục tọa độ của không gian k yếu tố và cách tâm phương án khoảng cách α > 0. - Phần tâm gồm n0 (n0 ≥ 1) thí nghiệm ở tâm phương án dùng để xác định phương sai tái hiện trong công thức kiểm tra ý nghĩa của các hệ số hồi qui. k Tổng số thí nghiệm trong phương án là N = 2 + 2k + n0
  41. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.1 PP QUI HOẠCH THỰC NGHIỆM Một vài dạng phương trình hồi qui: k + hồi qui bậc 1: y = b0 + j=1 bj xj k k + hồi qui bậc 1 đầy đủ: y = b0 + j=1 bj xj + j≠i bij xixj 2 ( với Ck hệ số bij) + hồi qui bậc 2 đầy đủ: k k k 2 y = b0 + j=1 bj xj + j≠i bij xixj + j=1 bjj xj
  42. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.2 KIỂM TRA SỰ TƢƠNG HỢP CỦA PTHQ Bước 1: kiểm tra ý nghĩa của các hệ số phương trình hồi qui bằng tiêu chuẩn Student tα (với α: mức ý nghĩa, α = 0,05) bj Chọn thống kê: tbj = sbj Với sbj: độ lệch quân phương của hệ số thứ i. Nếu tbi > tα(fth) thì hệ số bi được giữ lại trong phương trình hồi qui.(fth = n0 – 1 : bậc tự do tái hiện) Nếu tbi < tα(fth) thì hệ số bi bị loại khỏi phương trình hồi qui.
  43. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.2 KIỂM TRA SỰ TƢƠNG HỢP CỦA PTHQ 2 Phương sai sbj được xác định theo công thức: 2 2 2 sth sth sbj = n 2 = i=1 xi 2 Phương sai tái hiện sth: n0 2 1 2 sth = (yi0 − y o) n0 − 1 i=1
  44. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.2 KIỂM TRA SỰ TƢƠNG HỢP CỦA PTHQ Bước 2: kiểm tra sự tương thích của phương trình hồi qui theo tiêu chuẩn Fisher: 2 sdư F = 2 sth 2 Phương sai sdư được xác định theo công thức: n 1 s2 = (y − y )2 dư N − L i i=1 Với L: số hệ số có ý nghĩa trong phương trình hồi qui. Nếu F < Fα(α, ftt, fth) thì mô hình thống kê phù hợp với số liệu thực nghiệm.(ftt = N – L , fth = n0 – 1)
  45. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Để giải bài toán tối ưu hóa thực nghiệm chúng ta cần tiến hành các bước sau: + Bước 1: chọn phương án tiến hành thí nghiệm. + Bước 2: lập ma trận thực nghiệm X. + Bước 3: tiến hành thí nghiệm để xác định giá trị biến đầu ra Y. + Bước 4: xác định các hệ số trong phương trình hồi qui. + Bước 5: đánh giá phương trình hồi qui thu được. + Bước 6: xác định chế độ thực nghiệm tối ưu.
  46. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Một số hàm thường dùng trong tối ưu hóa: + Nhân hai ma trận: MMULT(array1, array2). + Tính định thức của ma trận: MDETERM(array). + Tính ma trận nghịch đảo: MINVERSE(array). + Tính ma trận chuyển vị: TRANPOSE(array). + Tính giá trị trung bình các số hạng: AVERAGE(number1, number2, ). + Tính tổng bình phương x2: SUMSQ(number1, number2, )
  47. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Một số hàm thường dùng trong tối ưu hóa: 2 + Tính tổng bình phương độ lệch (y푖 − ) : SUMXMY2(array_x, array_y). n (x −x )2 + Tính độ lệch chuẩn của mẫu i=1 i : n−1 STDEV(number1, number2, ). + Tra chuẩn số Student: TINV(p1, p2). + Tra chuẩn số Fisher: FINV(α, p1, p2).
  48. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Hãy tìm mối quan hệ giữa hàm mục tiêu y và các biến Z1, Z2, Z3 theo mô hình trực giao cấp 1 với số liệu thu được như sau: n = 3 n = 2k 0 Biến thực 8 9 1 2 3 4 5 6 7 10 11 225 Z1 150 300 150 300 150 300 150 300 225 225 60 60 60 Z2 30 30 90 90 30 30 90 90 45 30 30 30 Z3 15 15 15 15 45 45 45 12,0 13,8 13,2 y 3,0 6,0 10,0 12,0 15,0 23,0 12,0 18,0
  49. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Bƣớc 1: lập 1 bảng tính excel với các thông tin sau: Zj0 : mức cơ sở. Zj0 : D2 = AVERAGE(B2:C2). Khoảng biến thiên Zj. Zj : E2 = (C2-B2)/2
  50. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Bƣớc 2: chuyển đổi các biến thực sang biến mã hóa: G8=(C8-$D$2)/$E$2 풁 −풁 x1= 풁 x0= 1
  51. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL −1 Bƣớc 3: xác định các hệ số B = XTX . (XTY)
  52. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Bƣớc 4: kiểm tra ý nghĩa của các hệ số hồi qui: STDEV(J16:J18)^2 TINV(0.05,2) 2 2 n0 bj 2 sth sth 2 1 2 tbj = sbj = n 2 = ; sth = (yi0 − y o) s x n0 − 1 bj i=1 i i=1
  53. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL PTHQ có dạng: y = 12,375 + 2,375 x1 + 4,625 x3
  54. 4. TỐI ƢU HÓA THỰC NGHIỆM 4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL Bƣớc 5: kiểm tra sự tương hợp của phương trình hồi qui: 2 sdư F = 2 sth n 1 s2 = (y − y )2 dư N − L i i=1 Fα(α, ftt, fth) = Fα(0.05,5,2): FINV(0.05,5,2) PTHQ có dạng: y = 12,375 + 2,375 x1 + 4,625 x3