Sự ổn định của một vài lớp phương trình sai phân dạng phân thức Sự ổn định của một vài lớp phương trình sai phân dạng phân thức Sự ổn định của một vài lớp phương trình sai phân dạng phân thức luận văn tốt nghiệp,luận văn thạc sĩ, luận văn cao học, luận văn đại học, luận án tiến sĩ, đồ án tốt nghiệp luận văn tốt nghiệp,luận văn thạc sĩ, luận văn cao học, luận văn đại học, luận án tiến sĩ, đồ án tốt nghiệp
ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Phạm Đình Tùng THIẾT KẾ MẶT ĐÁP ỨNG 3-MỨC VỚI TÍNH CHẤT HIỆU ỨNG BẬC HAI TRỰC GIAO LUẬN ÁN TIẾN SĨ TOÁN HỌC Hà Nội - 2016 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Phạm Đình Tùng THIẾT KẾ MẶT ĐÁP ỨNG 3-MỨC VỚI TÍNH CHẤT HIỆU ỨNG BẬC HAI TRỰC GIAO Chuyên ngành: Lý thuyết xác suất Thống kê Toán học Mã số: 62 46 01 06 LUẬN ÁN TIẾN SĨ TOÁN HỌC Người hướng dẫn khoa học: TS NGUYỄN KỲ NAM GS.TS NGUYỄN HỮU DƯ Hà Nội - 2016 Lời cam đoan Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết trình bày luận án trung thực chưa công bố cơng trình khác Nghiên cứu sinh Phạm Đình Tùng Lời cảm ơn Luận án thực hoàn thành hướng dẫn khoa học TS Nguyễn Kỳ Nam GS.TS Nguyễn Hữu Dư, hai người thầy vơ mẫu mực, tận tình giúp đỡ đường khoa học Hai thầy dìu dắt tơi đường tốn học, đưa tơi bước vào lĩnh vực tốn học đầy thú vị, tạo thử thách giúp tự học hỏi, tìm tịi sáng tạo, tơi may mắn tiếp nhận từ hai người thầy đáng kính Tơi xin bày tỏ lịng biết ơn sâu sắc đến hai thầy Trong trình học tập nghiên cứu để hồn thành luận án, tơi nhận nhiều giúp đỡ quý báu thầy cô Bộ môn Xác suất Thống kê Khoa Toán-Cơ-Tin học, Trường Đại học Khoa học Tự Nhiên-ĐHQG Hà Nội Tôi xin trân trọng giúp đỡ thầy cô Tôi muốn bày tỏ cảm ơn chân thành đến Ban Giám hiệu, Ban Chủ nhiệm Khoa Tốn-Cơ-Tin học, Phịng sau đại học phòng ban chức Trường Đại học Khoa học Tự Nhiên-ĐHQG Hà Nội tạo điều kiện thuận lợi cho học tập nghiên cứu Tôi xin bày tỏ lịng biết ơn đến gia đình tồn thể bạn bè ln khuyến khích, động viên để tơi vững bước đường tốn học chọn Mục lục Trang Lời cam đoan Lời cảm ơn Thuật ngữ Anh Việt Danh mục ký hiệu chữ viết tắt Danh mục bảng 10 Danh mục hình vẽ 11 Mở đầu 12 Chương Giới thiệu phương pháp bề mặt đáp ứng 17 1.1 Thiết kế thí nghiệm 17 1.2 Phương pháp bề mặt đáp ứng 19 1.2.1 Mơ hình bề mặt đáp ứng 20 1.2.2 Tối ưu bề mặt đáp ứng 22 1.2.3 Thiết kế thí nghiệm phù hợp mơ hình bề mặt đáp ứng 23 1.3 Thí nghiệm sinh tổng hợp enzym Lipase 26 Chương Thiết kế Box-Behnken nhỏ với nhóm trực giao 32 2.1 Thiết kế Box-Behnken 32 2.2 Thiết kế nhóm khơng đầy đủ với cỡ nhóm khác 34 2.3 Phương pháp xây dựng thiết kế Box-Behnken nhỏ 36 2.4 Thiết kế Box-Behnken nhỏ với nhóm trực giao 38 2.5 Chia nhóm trực giao thiết kế bậc hai dựa IBD 40 2.6 Kết luận chương 41 Chương Thiết kế ước lượng tương đương SPLIT-PLOT 42 3.1 Thiết kế SPLIT-PLOT 43 3.2 Mơ hình bề mặt đáp ứng cho thiết kế SPLIT-PLOT tổng quát 45 3.3 Cấu trúc mong muốn ma trận thông tin SPD 46 3.4 Thuật toán SPLIT 47 3.5 Kết đánh giá 50 3.5.1 Các ví dụ minh chứng 51 3.5.2 Kết tính tốn 54 3.6 Kết luận chương 57 Kết luận 60 Danh mục cơng trình khoa học tác giả liên quan đến luận án 62 Tài liệu tham khảo 63 Phụ lục 73 Thuật ngữ anh việt Augmented-Pair Design Thiết kế tăng cặp Balanced Incomplete Block Design Thiết kế nhóm khơng đầy đủ cân Block Nhóm Block Effect Hiệu ứng nhóm Blocking Factor Nhân tố chia nhóm Box-Behnken Design Thiết kế Box-Behnken Box-Plot Đồ thị hộp ria mèo Balanced Lattice Lưới cân Center Point/Run Điểm tâm Central Composite Design Thiết kế tổng hợp trung tâm Concurrence Matrix Ma trận lặp lại Coordinate-exchange Algorithm Thuật tốn đổi tọa độ Contour Plot Đồ thị đường vịng Cyclic Tuần hoàn Design Matrix Ma trận thiết kế D-optimal D-tối ưu D-efficient D-hiệu Easy-to-Change Factor Nhân tố dễ thay đổi mức giá trị Effect Hiệu ứng Equivalent-Estimation Split- Plot Design Thiết kế ước lượng tương đương SPLIT-PLOT Expanded Design Matrix Ma trận thiết kế mở rộng Experimental Region Miền thí nghiệm Factor Nhân tố Factorial Design Thiết kế nhân tố Fractional Factorial Design Thiết kế phần nhân tố First-order Model Mơ hình đa thức bậc First-order Orthogonal Design Thiết kế bậc trực giao Generalized Least Square Method Phương pháp bình phương tối thiểu tổng quát Hard-to-Change Factor Nhân tố khó thay đổi mức giá trị Incomplete Block Design Thiết kế nhóm khơng đầy đủ Interaction Effect Hiệu ứng tương tác Interchange Algorithm Thuật toán hốn đổi Intrablock Trong nhóm Main/First-order Effect Hiệu ứng hiệu ứng bậc Least Square Method Phương pháp bình phương tối thiểu Level Mức giá trị Orthogonal Quadratic Effect Property Tính chất hiệu ứng bậc hai trực giao Orthogonal Block Nhóm trực giao Partial Balanced Incomplete Block Design Thiết kế nhóm khơng đầy đủ bán cân Placket-Burman Design Thiết kế Placket-Burman Point-exchange Algorithm Thuật toán hoán đổi điểm Quadratic Effect Hiệu ứng bậc hai Qualitative Định tính Quantitative Định lượng Randomization Ngẫu nhiên hóa Regression Technique Kỹ thuật hồi quy Regular Graph Design Thiết kế đồ thị quy Replicate Set Nhóm lặp lại Replication Lặp lại phép thử Resolvable Tách Response Biến đáp ứng Response Function Hàm đáp ứng Resovable Incomplete Block Design Thiết kế nhóm khơng đầy đủ tách Response Surface Design Thiết kế bề mặt đáp ứng Response Surface Model Mơ hình bề mặt đáp ứng Response Surface Methodology Phương pháp bề mặt đáp ứng Rotatability Tính xoay quanh Robust Tính vững Run Phép thử thí nghiệm Runs Số phép thử thí nghiệm Screening Design Thiết kế sàng lọc Second-order Design Thiết kế bậc hai Second-order Split-Plot Design Thiết kế SPLIT-PLOT bậc hai Second-order Model Mơ hình đa thức bậc hai Small Composite Design Thiết kế tổng hợp nhỏ Stationary Point Điểm dừng Surface Plot Đồ thị bề mặt Subset Design Thiết kế tập Supplementary Set Design Thiết kế tập phần phụ Split-Plot Design Thiết kế SPLIT-PLOT Sub-Plot Ô nhỏ Try Lần thử Trial and Error Thử sai Whole-Plot Ô lớn White Noise Nhiễu trắng Danh mục ký hiệu chữ viết tắt APD Thiết kế tăng cặp BBD Thiết kế Box-Behnken BIBD Thiết kế nhóm khơng đầy đủ cân CCD Thiết kế tổng hợp trung tâm EE-SPD Thiết kế ước lượng tương đương SPLIT-PLOT ETC Dễ thay đổi mức giá trị FD Thiết kế nhân tố đầy đủ FFD Thiết kế phần nhân tố đầy đủ GLS Bình phương tối thiểu tổng quát HTC Khó thay đổi mức giá trị IBD Thiết kế nhóm khơng đầy đủ OLS Bình phương tối thiểu thông thường PBD Thiết kế Placket-Burman PBIBD Thiết kế nhóm khơng đầy đủ bán cân RGD Thiết kế đồ thị quy RSD Thiết kế bề mặt đáp ứng RSM Phương pháp bề mặt đáp ứng SCD Thiết kế trung tâm nhỏ SBBD Thiết kế Box-Behnken nhỏ SPD Thiết kế SPLIT-PLOT [49] Vining, G G., Kowalski, S M & Montgomery, D C (2005), "Response surface designs within a split-plot structure", Journal of Quality Technology, 37, pp.115-129 [50] Wong, W.K (1994), "Comparing robust properties of A, D, E and Goptimal designs", Computational Statistics & Data Analysis, 18, pp.441-448 [51] Xiao, L L., Lin, D K J & Bai, F S (2012), "Constructing Definitive Screening Designs Using Conference Matrices" Journal of Quality Technology, 44, pp.2-8 [52] Zang, T.F, Yang,J.F & Lin, D.K.J (2011), "Small Box-Behnken design", Statistics and Probability Letters, 81, pp.1027-1033 68 Phụ lục Ma trận thiết kế cho SBBDs với m = 6, · · · , 16 Ma trận cho SBBDs với m = 6, 8, 12 15 nhận cách xóa bỏ cột cuối ma trận với m + Các dòng đánh dấu † sử dụng nửa thiết kế nhân tố 24 Các hàng đánh dấu ‡ sử dụng hai lần m=7 ±1 ±1 ±1 0 ±1 ±1 ±1 0 ±1 ±1 ±1 0 ±1 ±1 ±1 0 ±1 ±1 ±1 0 ±1 ±1 ±1 0 0 ±1 ±1 ±1 m=9 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 ±1 ±1 0 0 0 0 ±1 ±1 ±1 0 0 0 0 ±1 ±1 ±1 0 ±1 0 0 ±1 ±1 0 ±1 ±1 0 ±1 ±1 ±1 0 0 ±1 ±1 ±1 ±1 0 0 0 ±1 ±1 ±1 69 ±1 ±1 ±1 ±1 m = 10 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 ±1 ±1 0 0 0 0 ±1 ±1 ±1 0 0 0 0 ±1 ±1 ±1 0 ±1 0 0 ±1 ±1 0 ±1 ±1 0 ±1 ±1 ±1 0 0 ±1 ±1 ±1 0 0 ±1 ±1 ±1 ±1 0 0 0 ±1 ±1 ±1 0 ±1 0 0 ±1 ±1 0 0 0 ±1† ±1 ±1 ±1 0 m = 11 ±1 0 ±1 0 ±1 0 0 ±1 0 ±1 0 ±1 0 0 ±1 0 ±1 0 ±1 0 ±1 ±1 ±1 0 0 0 0 0 ±1 ±1 ±1 0 0 0 0 0 ±1 ±1 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 ±1 0† 0 ±1 ±1 0 ±1 ±1 0 ±1 0 ±1 ±1 0 ±1 ±1 0 0 ±1 ±1 0 ±1 0 0 ±1 ±1 0 ±1† ±1 ±1 0 0 ±1 ±1 0 ±1 ±1 ±1 0 ±1 0 0 0 0 ±1 ±1‡ 70 m = 13 ±1 0 ±1 ±1 0 0 0 ±1 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 0 ±1 ±1 0 0 0 ±1 ±1 ±1 ±1 m = 16 ±1 ±1 ±1 ±1 0 0 0 0 0 0 0 ±1 ±1 ±1 ±1 0 0 0 0 0 0 0 ±1 ±1 ±1 ±1 0 0 0 0 0 0 0 ±1 ±1 ±1 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 0 ±1 0 0 ±1 0 0 ±1 0 ±1 0 0 0 ±1 0 ±1 0 ±1 0 0 ±1 ±1 0 0 ±1 0 0 ±1 0 ±1 0 ±1 0 ±1 0 ±1 0 0 ±1 0 0 ±1 ±1 0 ±1 0 0 ±1 0 ±1 ±1 0 0 ±1 ±1 0 0 ±1 0 0 0 0 ±1 ±1 0 ±1 0 0 ±1 ±1 0 0 0 ±1 ±1 0 0 ±1 ±1 0 0 ±1 ±1 0 0 0 0 ±1 0 ±1 0 0 ±1 ±1 0 0 ±1 ±1 0 0 ±1 0 ±1 71 ±1 ±1 ±1 ±1 ±1 Điều kiện tương đương ước lượng bình phương tối thiểu tổng qt thơng thường Như trình bày mục 3.2, EE-SPD SPD có ước lượng bình phương tổi thiểu thơng thường tổng quát giống nhau, nghĩa βˆ OLS = βˆ GLS βˆ OLS = (X X)−1 X Y βˆ GLS = (X V−1 X)−1 X V−1 Y Khi βˆ OLS = βˆ GLS , ước lượng βˆ OLS βˆ GLS tương đương Định lý 3.1 Với σγ2 = σε2 = 1, ước lượng βˆ OLS βˆ GLS tương đương trace(C C) = với C C = kB − B(X X)−1 B, B = X Jb X Chứng minh: Với giả thiết σγ2 = σε2 = 1, ma trận hiệp phương sai mơ hình V = In + Jb có nghịch đảo V−1 = In − k+1 Jb với n số phép thử thí nghiệm SPD, b số lớn SPD k = n/b cỡ ô lớn Đặt P = X(X X)−1 X , ước lượng bình phương tối thiểu tổng qt khai triển dạng sau: βˆ GLS = (X V−1 X)−1 X V−1 (PY + (In − P)Y) = (X X)−1 X Y + (X V−1 X)−1 X V−1 (In − P)Y = βˆ OLS + (X V−1 X)−1 X V−1 (In − P)Y (3.12) Khi đó, tính chất ước lượng tương đương xảy (X V−1 X)−1 X V−1 (In − P)Y = O, ∀ Y ⇔ X V−1 (In − P) = O (3.13) Thay ma trận V−1 , P (3.13) làm đơn giản, ta (In − X(X X)−1 X )Jb X = O (3.14) Việc kiểm tra phần tử ma trận C = (In − X(X X)−1 X )Jb X tương đương với tổng bình phương phần tử C hay trace(C C) = Mặt khác, sử dụng tính chất lũy linh ma trận A = In − X(X X)−1 X (nghĩa 72 A2 = A ) dùng ma trận B thay cho X Jb X, ta có C C = X Jb AAJb X = X Jb AJb X = X Jb (In − X(X X)−1 X )Jb X = kX Jb X − X Jb X(X X)−1 X Jb X = kB − B(X X)−1 B 73 (3.15) ... thuộc vào giá trị nhân tố dạng y = f (x1 , x2 , · · · , xm ) + ε, (1.1) hàm đáp ứng f chưa biết sai số ε sai số thí nghiệm giả sử có phân bố chuẩn với trung bình phương sai σ Khi đó, E(y) = f... · · · , βm−1m ] ước lượng phương pháp bình phương tối thiểu βˆ = (X X)−1 X Y, có ma trận hiệp phương sai cov(βˆ ) = σ (X X)−1 Đại lượng ˆ x)) Nvar(Y( σ2 gọi phương sai dự báo mơ hình điểm ˆ... đưa công trình Vining et al [49] năm 2005 Loại thiết kế bắt nguồn từ tên ước lượng tham số mơ hình phương pháp bình phương tối thiểu giống với ước lượng phương pháp bình phương tối thiểu tổng quát