Toán học - Bài tập quy hoạch thực nghiệm
Bạn đang xem 20 trang mẫu của tài liệu "Toán học - Bài tập quy hoạch thực nghiệm", để 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:
- toan_hoc_bai_tap_quy_hoach_thuc_nghiem.ppt
Nội dung text: Toán học - Bài tập quy hoạch thực nghiệm
- BỘ CÔNG THƯƠNG TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP THỰC PHẨM TPHCM KHOA CÔNG NGHỆ SINH HỌC VÀ KỸ THUẬT MÔI TRƯỜNG TIỂU LUẬN: GVHD: DƯƠNG HOÀNG KIỆT NHÓM 9: 1.HOÀNG THỊ HOA 2208110022 2.NGUYỄN THỊ HỒNG 2208110023 3.NGUYỄN NGỌC HÙNG 2208110025 4.TRẦN THỊ HƯƠNG 2208110027 5.NGUYỄN THANH THIÊN HƯƠNG 2208110028 6.NGUYỄN THỊ LIÊN 2208110033 7.LÊ THỊ KIM LIÊN 2208110034 09/12/2012
- KẾ HOẠCH LÀM VIỆC NHÓM stt Nội dung công việc Thời gian Người thực hiện hoàn thành 1 Xác định hồi quy thực nghiệm (bài 2.6). 18/04/2012 Nguyễn Thị Liên 2 Kiểm định sự đồng nhất của pp thí 23/04/2012 Nguyễn Thị Hồng nghiệm và tính phương sai tái hiện(bài 2.6), Sildes trang bìa. 3 Kiểm định hệ số hồi quy(bài 2.6) . 09/05/2012 Trần Thị Hương 4 Kiểm định sự phù hợp của mô hình vừa 15/04/2012 Ngọc Hùng tìm được(bài 2.6) , tài liệu tham khảo. 5 Tìm mô hình phù hợp (bài 3.6), Đặt vấn 21/04/2012 Nguyễn Thanh Thiên Hương đề. 6 Kiểm định giả thuyết Ho (bài 3.6), kết 04/05/2012 Lê Thị Kim Liên luận và kiến nghị. 7 Kiểm định sự phù hợp (bài 3.6) , 09/05/2012 Hoàng Thị Hoa Soạn slides và tổng hợp hoàn chỉnh .
- ĐẶT VẤN ĐỀ Phương trình hồi qui gần đúng phụ thuộc vào phương pháp tính dùng để tính các hệ số hồi qui. Phương pháp bình phương nhỏ nhất xác định hệ số phương trình hồi qui sao cho gần đúng với kỳ vọng toán học của thực nghiệm.
- BÀI TẬP 1: PHƯƠNG PHÁP BÌNH PHƯƠNG CỰC TIỂU GIẢI BÀI TẬP 2.6
- CHO BẢNG SỐ LIỆU THỰC NGHIỆM DƯỚI ĐÂY, THỰC HIỆN N=7 (THÍ NGHIỆM LẶP M=4 N X0 X1 X2 X3 X4 Y1 Y2 Y3 Y4 1 1 1.2 1.3 1.4 1.1 -0.0159 -0.49 -0.53 -0.54 2 1 0.7 1.1 0.5 0.4 1.98 2.04 2.05 2.01 3 1 0.5 0.6 0.8 0.4 -0.2767 -1.35 -1.31 -1.34 4 1 0.8 0.5 1.1 0.6 -3.63 -3.72 -3.68 -3.71 5 1 0.9 0.7 0.7 0.9 0.27 0.25 0.262 0.29 6 1 1.1 0.9 0.8 0.7 2.62 2.57 2.65 2.67 7 1 0.6 0.8 1.2 0.5 -3.05 -3.62 -3.04 -3.11
- 1.Tìm PTHQ dạng:ŷ=b0+b1x1+b2x2+b3x3+b4x4 1 1.2 1.3 1.4 1.1 -0.394 1 0.7 1.1 0.5 0.4 2.020 bo -1.319 b 1 0.5 0.6 0.8 0.4 Y= 1 X= B= b2 1 0.8 0.5 1.1 0.6 -3.685 b3 0.268 b 1 0.9 0.7 0.7 0.9 4 2.628 1 1.1 0.9 0.8 0.7 -3.205 1 0.6 0.8 1.2 0.5
- 7 5.8 5.9 6.5 4.6 XtX = 5.8 5.2 5.13 5.54 4.16 5.9 5.13 5.45 5.57 4.07 6.5 5.54 5.57 6.63 4.51 4.6 4.16 4.07 4.51 3.44 3.246 -2.661 -0.912 -1.057 1.932 -2.661 12.092 -2.449 1.089 -9.595 (XtX)-1= -0.912 -2.449 3.079 -0.114 0.688 -1.507 1.089 -0.114 2.287 -2.165 1.932 -9.595 0.688 -2.165 11.336
- -3.688 -4.320 bo= -4.320 -1.458 5.813 b1 = 5.813 t X Y= -0.936 3.462 B= b2 = 3.462 -6.207 -6.035 b3 = -6.035 -1.886 -1.723 b = -1.723 VẬY PTHQ THỰC NGHIỆM TÌM ĐƯỢC 4 LÀ ŷ = -4.320 + 5.813x1 + 3.462x2 -6.035 x3 – 1.723x4
- 2.KIỂM ĐỊNH GIẢ THUYẾT HO: Thí nghiệm 1: y11= -0.0159 ; y12= -0.49 ; y13=-0.53 ; y14=-0.54 1 y1 = (-0.0159-0.49-0.53-0.54)= -0.394 4 2 Sth1 = 1 = [(-0.0159+0.394)2 + (-0.49+0.394)2 + 3 (0.53+0.394)2 +(-0.54+0.394)2] = 0.064 2 2 2 Tương tự: Sth2 = 0.0010 ; Sth3 = 0.0011;Sth4 = 0.0016 2 2 2 Sth5 = 0.00028 ; Sth6 = 0.0019 ; Sth7 = 0.0775
- Giả thuyết cần kiểm định Ho: “phương sai tái hiện của từng thí nghiệm bằng nhau” 3,7 Vì G < G 5% Chấp nhận Ho tức là phương sai tái hiện của từng thí nghiệm như nhau. Khi đó phương sai tái hiện của từng cuộc thí nghiệm là
- 3. KIỂM ĐỊNH GIẢ THUYẾT H0 : ΒJ = 0 Theo câu 2, ta có 5% Tra bảng Student t 21 = 2.080, tbj> tα j bj kết luận 0 -4.32 3.246 0.068 16.524 b0 # 0 0.261 1 5.813 12.092 0.255 11.520 b1# 0 0.505 2 3.462 3.079 0.065 13.597 b2 # 0 0.255 3 -6.035 2.287 0.048 27.501 b3 # 0 0.219 4 -1.723 11.336 0.239 3.527 b4 # 0 0.489 Vậy ŷ = -4.320 + 5.813x + 3.462x Với 1 2 -6.035 x3 – 1.723x4
- 4.KIỂM ĐỊNH SỰ PHÙ HỢP Tra bảng Fisher ta có: = 5.780 Từ PT hồi quy ta có: Y Ŷi -0.394 -3.188 2.020 -0.149 -1.319 -4.854 -3.685 -5.611 0.268 -2.440 2.628 -0.844 -3.205 -6.166 Vậy mô hình ŷ = -4.320 + 5.813x1 + 3.462x2-6.035 x3 – 1.723x4 không phù hợp với số liệu thực nghiệm
- BÀI TẬP 2: QUY HOẠCH TRỰC GIAO GIẢI BÀI TẬP 3.6
- 3.6 Lập mô hình bậc hai nghiên cứu ảnh hưởng của ba yếu tố vào thông số tối ưu hoá y được cho bằng số liệu thực nghiệm sau đây (n0 =1): N x0 x1 x2 x3 y 1 1 1 1 1 4.89 2 1 -1 1 1 4.20 3 1 1 -1 1 2.48 4 1 -1 -1 1 2.22 5 1 1 1 -1 4.70 6 1 -1 1 -1 3.80 7 1 1 -1 -1 2.65 8 1 -1 -1 -1 2.16 9 1 0 0 0 2.30 10 1 α 0 0 3.55 11 1 -α 0 0 4.50 12 1 0 α 0 1.80 13 1 0 -α 0 5.15 14 1 0 0 α 2.32 15 1 0 0 -α 2.56 Ba thí nghiệm ở tâm phương án:
- 1.Tìm mô hình phù hợp Ta có: Đặt Ta có mô hình Khi đó ta có ma trận sau:
- Khi đó ta có ma trận bố trí thí nghiệm sau : N x0 x1 x2 x3 x1x x1x3 x2x3 x’1 x’2 x’3 y 2 1 1 1 1 1 1 1 1 0.27 0.27 0.27 4.89 2 1 -1 1 1 -1 -1 1 0.27 0.27 0.27 4.20 3 1 1 -1 1 -1 1 -1 0.27 0.27 0.27 2.48 4 1 -1 -1 1 1 -1 -1 0.27 0.27 0.27 2.22 5 1 1 1 -1 1 -1 -1 0.27 0.27 0.27 4.70 6 1 -1 1 -1 -1 1 -1 0.27 0.27 0.27 3.80 7 1 1 -1 -1 -1 -1 1 0.27 0.27 0.27 2.65 8 1 -1 -1 -1 1 1 1 0.27 0.27 0.27 2.16 9 1 0 0 0 0 0 0 -0.73 -0.73 -0.73 2.30 10 1 1.215 0 0 0 0 0 0.747 -0.73 -0.73 3.55 11 1 -1.215 0 0 0 0 0 0.747 -0.73 -0.73 4.50 12 1 0 1.215 0 0 0 0 -0.73 0.747 -0.73 1.80 13 1 0 -1.215 0 0 0 0 -0.73 0.747 -0.73 5.15 14 1 0 0 -1.215 0 0 0 -0.73 -0.73 0.747 2.32 15 1 0 0 -1.215 0 0 0 -0.73 -0.73 0.747 2.56
- 1.Tìm mô hình phù hợp Tính bj
- 1.Tìm mô hình phù hợp Vậy phương trình hồi quy là Hay
- 2.Kiểm định H0 : Ba thí nghiệm ở tâm phương án: Với m là số thí nghiệm ở tâm phương án
- 3.Kiểm định giả thuyết H0 : Tra bảng student với mức ý nghĩa α = 5%, m = 3 (số thí nghiệm ở tâm phương án) Tính tương tự ta có
- 3.Kiểm định giả thuyết H0 : Tính tương tự ta có Tính tương tự ta có
- j bj tj Kết luận 0 2.827 86.896 b0≠0 1 0.108 2.837 b1=0 2 0.365 9.588 b2≠0 3 0.017 0.447 b3=0 12 0.105 2.357 b12=0 13 0.088 1.975 b13=0 23 -0.055 1.235 b23=0 11 0.692 11.466 b11≠0 22 0.319 5.286 b22≠0 33 -0.382 6.329 b33≠0 Các hệ số không thoả là b1 , b 2, b3 , b13 , b23 . Nên ta loại nhân tố x1 , x3 , x12 , x23 , x13 và nhận các nhân tố x0 , x2 , x11 , x22 , x33 Vậy phương trình quy hoạch thực nghiệm là: (1)
- Tra bảng Fisher với mức ý nghĩa α = 5%, N=15, l=5, m= 3 (Với m là số thí nghiệm ở tâm phương án) Thế các nhân tố vào (1) ta có: Tương tự ta tính các giá trị:
- Tính phương sai dư: Ta có: Vậy mô hình trên không phù hợp
- KIẾN NGHỊ • Giải bài toán bằng hương pháp bình phương cực tiểu thì các đối tượng thực được biểu diễn một cách gần đúng. như vậy cần có một mô hình đưa ra được một cách chính xác về các đối tượng thực này • Cần đánh giá về sai số của phương pháp, nghĩa là cần có những kết luận thống kê về lời giải. • Các công thức tham khảo của các sách có sự ký hiệu khác nhau đòi hỏi người đọc phải nắm chắc được kiến thức.
- TÀI LIỆU THAM KHẢO 1. PGS.TS. Bùi Minh Trí, Xác xuất thống kê và quy hoạch thực nghiệm, nhà xuất bản khoa học và kỹ thuật, 2006. 2. Nguyễn Cảnh, Quy hoạch thực nghiệm, NXB Đại học quốc gia tp HCM, 2004. 3. Dương Hoàng Kiệt, bài tập quy hoạch thực nghiệm, lưu hành nội bộ trường ĐH CNTP tp HCM, 2012.