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
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
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
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
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)
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
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
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)
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
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
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
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
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.
1. GIẢI HỆ PHƢƠNG TRÌNH TUYẾN TÍNH Kết quả giải bằng Solver:
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]
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.
2. BÀI TOÁN TRUYỀN NHIỆT Thép chịu lực Gạch chịu nhiệt
3
1 2
Gạch cách nhiệt 1= 0,12m, 1= 0,81 W/ m.K 2= 0,065 m, 1= 0,23 W/ m.K = 0,01 m, = 45 W/ m.K
t1= 8000C, 1=69,6 W/ m2.K t2= 350C, 2=13,9 W/ m2.K
2. BÀI TOÁN TRUYỀN NHIỆT 1
1
q1= 1(t1-tT1)
q2= (tT1-tT2) 1
2
q3= (tT2-tT3) 2
q4= (tT3-tT4) 3
3
2
q5= 2(tT4-t2)
Lƣu ý: q1 = q2 = q3 = q4 = q5 = q
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:
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
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 q 1
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
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
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%
2. BÀI TOÁN TRUYỀN NHIỆT a. Sử dụng hàm Goalseek Kết quả:
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)
q2= (tT1-tT2) 1
q3= (tT2-tT3) 2
q4= (tT3-tT4) 3
q5= 2(tT4-t2)
q
+ 1.tT1
= 1.t1
1 .q-
=0 .tT1 + .tT2 - .tT2 + .tT3 =0 2 .q - .tT3 + .tT4=0 3 .q
q
- 2.tT4 = 2.t2
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:
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
2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 4: gọi hàm Solver
2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Bước 5: cài đặt VT = VP
2. BÀI TOÁN TRUYỀN NHIỆT b. Sử dụng hàm Solver Kết quả nghiệm:
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.
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 = +
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:
= 0,8x + 19,698 yL = + + • Phương trình đoạn chưng:
+ − yc = + +xw = 1,2419x – 0,3986 • 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
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
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,…
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
j=
j=
j≠i
f(x,x,...,xk) β0 βj xj βjj xj βij xixj... Với 1 ≤ i ≤ j ≤ k
PTHQ thực nghiệm:
y f x,x,...,xk
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:
y 1 x x … xk … x y x 1 x k X … Y ⋮ y 1 x x … xk (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 ⋮ bk Theo phương pháp bình phương cực tiểu, ta có: − T T T X X .B X Y B X X .(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
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
i=1
với j ≠ 0.
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
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)
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.
Tổng số thí nghiệm trong phương án là N = 2k + 2k + n0
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: + hồi qui bậc 1:
y b0 kj= bj xj
+ hồi qui bậc 1 đầy đủ:
y b0 kj= bj xj kj≠i bij xixj ( với
Ck hệ số b ) ij
+ hồi qui bậc 2 đầy đủ:
y b0 kj= bj xj kj≠i bij xixj kj= bjj xj
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) Chọn thống kê:
tj
Với sbj: độ lệch quân phương của hệ số thứ i. Nếu tbi > tα(f th) thì hệ số bi được giữ lại trong phương trình hồi qui.(f th = n0 – 1 : bậc tự do tái hiện) Nếu tbi < tα(f th) thì hệ số bi bị loại khỏi phương trình hồi qui.
4. TỐI ƢU HÓA THỰC NGHIỆM 4.2 KIỂM TRA SỰ TƢƠNG HỢP CỦA PTHQ
sj được xác định theo công thức: h h s s sj x i= i : Phương sai tái hiện sh 1 sh n0 1 (yi0 yo) i= Phương sai
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:
Phương sai
ư F
sdư được xác định theo công thức: 1 sdư N L (yi y) i=
Với L: số hệ số có ý nghĩa trong phương trình hồi qui. Nếu F < Fα(α, f tt, f th) thì mô hình thống kê phù hợp với số liệu thực nghiệm.(f tt = N – L , f th = n0 – 1)
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.
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
x: SUMSQ(number1, number2,…)
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: + Tính tổng bình phương độ lệch
(y ):
SUMXMY2(array_x, array_y). + Tính độ lệch chuẩn của mẫu
(−) : −
STDEV(number1, number2,…). + Tra chuẩn số Student: TINV(p1, p2). + Tra chuẩn số Fisher: FINV(α, p1, p2).
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: Biến thực
n0 = 3
n = 2k 8
1
2
3
4
5
6
7
Z1
150
300
150
300
150
300
150
300
Z2
30
30
90
90
30
30
90
90
Z3
15
15
15
15
45
45
45
y
3,0
6,0
10,0
12,0
15,0
23,0
12,0
45
18,0
9
10
11
225
225
60
60
60
30
30
30
12,0
13,8
13,2
225
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:
Z j0 : mức cơ sở.
Z j0 : D2 = AVERAGE(B2:C2).
Khoảng biến thiên Z j.
Z j : E2 = (C2-B2)/2
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
4. TỐI ƢU HÓA THỰC NGHIỆM
4.3 TỐI ƢU HÓA THỰC NGHIỆM BẰNG EXCEL − T Bƣớc 3: xác định các hệ số B X X .(X T Y)
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)
b tj sjj
s s 1 h h sj x ;sh n0 1 (yi0 yo) i= i i=
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 x 4,625 x