Thiết kế tối ưu nam châm điện dạng solenoid

Thiết kế tối ưu nam châm ñiện dạng solenoid ðại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ðẠI HỌC BÁCH KHOA KIM ðIỀN QUANG THIẾT KẾ TỐI ƯU NAM CHÂM ðIỆN DẠNG SOLENOID Chuyên ngành: Thiết bị mạng nhà máy ñiện Mã số ngành: 2.02.01 LUẬN VĂN THẠC SĨ Tp HỒ CHÍ MINH, tháng 10 năm 2006 Thiết kế tối ưu nam châm điện dạng solenoid CƠNG TRÌNH ðƯỢC HỒN THÀNH TẠI TRƯỜNG ðẠI HỌC BÁCH KHOA ðẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH Cán hướng dẫn khoa học: PGS.TS Nguyễn Chu Hùng Cán chấm nhận xét 1: Cán chấm nhận xét 2: Luận văn thạc sĩ ñược bảo vệ HỘI ðỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ðẠI HỌC BÁCH KHOA, ngày …… tháng …… năm 2006 Thiết kế tối ưu nam châm điện dạng solenoid LỜI CẢM ƠN ðể hồn thành luận văn tốt nghiệp này, xin chân thành cảm ơn Quý thầy cô trường ðại học Bách Khoa Thành phố Hồ Chí Minh tận lực đáp ứng nhu cầu sinh hoạt học tập, trang bị kiến thức khoa học kỹ thuật quý báo cho trình học Cao học ngành “Thiết bị, mạng nhà máy điện” ðặc biệt, tơi xin chân thành cảm ơn Thầy Nguyễn Chu Hùng người tận tình hướng dẫn, truyền ñạt kinh nghiệm giúp ñỡ ñiều kiện để hồn thành tập luận văn ðể có thành ngày hơm nay, tơi chân thành cảm ơn người thân, anh em, bạn bè đặc biệt cha, mẹ tơi, người bạn thân tốt bụng ln ln giúp đỡ động viên tơi suốt q trình thực luận văn Xin chân thành cảm ơn Học viên thực Kim ðiền Quang Thiết kế tối ưu nam châm ñiện dạng solenoid LỜI NÓI ðẦU Cùng với phát triển vượt bậc thành tựu ñạt ñược lĩnh vực thiết kế, chế tạo ñối tượng kỹ thuật, vấn đề tối ưu hóa đối tượng trở nên quan trọng ðó phù hợp theo tiêu chuẩn đặt như: Tiêu chuẩn công nghệ chế tạo, Tiêu chuẩn chất lượng, Tiêu chuẩn kinh tế …Nó cho phép tạo đối tượng với tính chất ưu việt hơn, gần đạt tới thơng số lý tưởng cơng nghệ chế tạo q trình làm việc, có giá thành rẻ Luận văn “Thiết kế tối ưu nam châm ñiện dạng solenoid” ñược thực nhằm mục đích trình bày phương pháp thiết kế tối ưu thiết bị điện tính tốn tối ưu hóa phần mềm matlab Như biết nam châm ñiện dạng solenoid ñược ứng dụng nhiều cơng nghiệp Trong nhiều chế tạo van ñiện từ dùng để đóng ngắt thiết bị lị hơi, khí nén, chất lỏng, gas, v.v… Việc lập chương trình tính tốn tối ưu hóa nam châm điện dạng solenoid matlab giúp có nhiều chọn lựa nam châm điện dạng solenoid theo gam cơng suất mà cần Nó cho thơng số kích thước hình học, lực hút điện từ, từ cảm tương ứng với gam cơng suất Nội dung tập luận văn trình bày chương sau: Chương 1: Sơ lược thiết kế tối ưu nam châm ñiện dạng solenoid Chương 2: Bài toán thiết kế tối ưu Chương 3: Phương pháp quy hoạch thực nghiệm Chương 4: Thiết kế tối ưu nam châm điện dạng solenoid Chương 5: Chương trình thiết kế tối ưu nam châm ñiện dạng hút xoay Chương 6: Kết luận kiến nghị Do thời gian có hạn, ñiều kiện tiếp cận thực tế chưa nhiều nên luận văn khơng tránh nhiều sai sót, kính mong ñược giúp ñỡ thông cảm quý thầy cô Học viên thực Kim ðiền Quang Thiết kế tối ưu nam châm ñiện dạng solenoid MỤC LỤC Chương SƠ LƯỢC VỀ THIẾT KẾ TỐI ƯU NAM CHÂM ðIỆN DẠNG SOLENOID 1.1 Mạch từ thiết bị kỹ thuật điện 01 1.2 ðặc tính vật liệu từ 02 1.3 ðặt vấn ñề tối ưu thiết bị ñiện từ 05 1.4 Sơ lược nam châm ñiện dạng solenoid 07 Chương BÀI TOÁN THIẾT KẾ TỐI ƯU 2.1 Các vấn đề chung tốn thiết kế tối ưu 17 2.1.1 Khái niệm mô hình tốn tối ưu 17 2.1.2 Miền giá trị cho phép, giá trị tối ưu tham số thiết kế 21 2.1.3 Lập hàm mục tiêu, tốn đa mục tiêu 21 2.1.4 Phân loại hàm mục tiêu 23 2.1.5 Tính tương đối giải pháp lựa chọn hàm mục tiêu 26 2.1.6 Cơ sở ñể thành lập sơ đồ phân tích để lựa chọn hàm mục tiêu 27 2.1.7 Các mức ñộ thiết kế tối ưu thiết bị ñiện 28 2.2 Phân loại phương pháp giải toán thiết kế tối ưu 28 2.3 Phép quy hoạch phi tuyến 36 2.3.1 Bài toán phép quy hoạch phi tuyến 36 2.3.2 Một vài khái niệm 36 2.3.3 Các tốn biến 38 2.2.4 Các phương pháp tìm kiếm thụ động 39 2.3.5 Các phương pháp tìm kiếm thứ tự 37 2.3.6 Bài tốn đa biến khơng ràng buộc 40 2.3.7 Các phương pháp Gradient 43 Thiết kế tối ưu nam châm điện dạng solenoid 2.3.8 Bài tốn đa biến có mặt ràng buộc 46 Chương QUY HOẠCH THỰC NGHIỆM 3.1 Khái niệm 52 3.2 Quy hoạch thực nghiệm bậc 54 3.3 Quy hoạch thực nghiệm bậc 58 3.3.1 Các phương án cấu trúc có tâm 58 3.3.2 Những phương án trực giao cấp 61 3.3.3 Phương án quay bậc Box Hunter 63 Chương THIẾT KẾ TỐI ƯU NAM CHÂM ðIỆN DẠNG SOLENOID 4.1 Những vấn ñề chung thiết kế tối ưu 69 4.2 ðặc tính nam châm điện dạng Solenoid thơng qua kích thước hình học 71 4.3 Bài toán thiết kế tối ưu nam châm điện dạng Solenoid 78 4.3.1 Chọn thơng số đầu vào 78 4.3.2 Hình thành tốn thiết kế tối ưu 79 4.3.3 Dùng phương pháp QHTN xấp xỉ hóa hàm đầu 80 Chương ỨNG DỤNG PHẦN MỀM MATLAB ðỂ LẬP TRÌNH TÍNH TỐN THIẾT KẾ TỐI ƯU NAM CHÂM ðIỆN DẠNG SOLENOID 5.1 Sơ lược phần mềm Matlab 84 5.1.1 Khái quát 84 5.1.2 Các hộp công cụ Matlab 86 5.1.3 Môi trường làm việc Matlab 88 5.2 Quy hoạch phi tuyến bị ràng buộc 5.2.1 Bài toán 01 mục tiêu 89 89 Thiết kế tối ưu nam châm ñiện dạng solenoid 5.2.2 Cực tiểu hàm vơ hướng nhiều biến có ràng buộc hàm fmincon 90 5.3 Giải thuật phương án tính tốn Matlab 91 5.31 Cơ sở lập trình 91 5.3.2 Giải thuật điều khiển Matlab 93 5.3.3 Kết tính tốn tốn quy hoạch phi tuyến bị ràng buộc 96 Phụ lục chương trình tính tốn thiết kế tối ưu nam châm điện dạng solenoid 100 Chương KẾT LUẬN VÀ KIẾN NGHỊ 6.1 Tóm tắt cơng trình nghiên cứu 108 6.2 Những dự kiến phát triển ñề tài 108 6.2.1 Tương lai việc thiết kế tối ưu thiết bị ñiện 109 6.2.2 Về phía người nghiên cứu 109 Thiết kế tối ưu nam châm ñiện dạng solenoid TÀI LIỆU THAM KHẢO 1-Các phương pháp tối tưu hóa – tác giả Bùi Thế Tâm Trần Vũ Thiệu, NXB GTVT - 1998 2- Thiết kế tối ưu hoá ứng dụng – tác giả Nguyễn Viết Trung, NXB xây dựng Hà Nội- 2003 3- Matlab & Simulink – tác giả Nguyễn Phùng Quang, NXB khoa học kỹ thuật Hà Nội – 2005 4- Giải toán tối ưu hoá ứng dụng Matlab – maple tối ưu hố tĩnh điều khiển tối ưu – tác giả Nguyễn Nhật Lệ Phan Mạnh Dần, NXB Khoa học Kỹ thuật Hà Nội 5-Kỹ thuật ñiện I - tác giả Nguyễn Chu Hùng Tôn Thất Cảnh Hưng, NXB ðHQG TP.HCM - 2001 6- Kỹ thuật ñiện 2-Máy ñiện quay - tác giả Nguyễn Hữu Phúc, NXB đại học Quốc gia TP Hồ Chí Minh - 2003 7- Thiết kế máy ñiện – Tác giả Trần Khánh Hà, NXB khoa học kỹ thuật Hà Nội – 1997 8- Xác suất thống kê Quy hoạch thực nghiệm - tác giả Bùi Minh Trí, NXB Khoa học Kỹ thuật 9- Quy hoạch thực nghiệm - tác giả Nguyễn Minh Tuyển, NXB Khoa học Kỹ thuật 10- Quy hoạch thực nghiệm - tác giả Nguyễn Cảnh, NXB ðại học Quốc gia Tp.HCM Thiết kế tối ưu nam châm ñiện dạng solenoid LÝ LỊCH TRÍCH NGANG: Họ tên: Kim ðiền Quang Ngày, tháng, năm sinh: 26/03/1982 Nơi sinh: Tỉnh Bình Dương ðịa liên lạc: ðiện lực Bình Dương, Số 03 Cách mạng tháng Tám, Phường Phú Thọ, Thị xã Thủ Dầu Một, Tỉnh Bình Dương Q TRÌNH ðÀO TẠO Năm 2004 tốt nghiệp ðại học Bách Khoa Tp.HCM (ngành ðiện - ðiện tử), sau thi đậu vào khố cao học 15 trường ðại học Bách Khoa Tp HCM tiếp tục học Q TRÌNH CƠNG TÁC Sau tốt nghiệp ðại học, công tác ðiện lực Bình Dương Chương 1: Sơ lược thiết kế tối ưu nam châm ñiện dạng solenoid Chương SƠ LƯỢC VỀ THIẾT KẾ TỐI ƯU NAM CHÂM ðIỆN DẠNG SOLENOID 1.1 MẠCH TỪ TRONG CÁC THIẾT BỊ KỸ THUẬT ðIỆN: Mạch từ thiết bị kỹ thuật ñiện tập hợp vật thể mơi trường nhằm mục đích tạo thành đường khép kín cho từ thơng Mạch từ thường cấu tạo vật liệu từ có độ từ thẩm cao Do độ từ thẩm cao vật liệu từ, ñường sức từ chủ yếu chạy phần dẫn từ vật thể, điều giống dịng điện chạy ñoạn dẫn ñiện mạch ñiện Mạch từ bao gồm lõi sắt từ có hay khơng có khe hở khơng khí từ thơng đóng kín qua chúng Về phương diện kết cấu, mạch từ thiết bị kỹ thuật ñiện ñược phân biệt theo ba loại sau : a) Mạch từ tĩnh, mạch từ thường có máy biến áp, trường hợp lý tưởng xem khơng có khe hở khơng khí Mặc đây, chuyển đổi lượng khơng phải điện – cơ, trường hợp đóng vai trị quan trọng chuyển đổi lượng nói chung b) Mạch từ có phần ứng chuyển ñộng tịnh tiến chuyển ñộng xoay ðó loại mạch từ thường có thiết bị đóng - cắt mạch điện contactor, áptomát, relay, máy ngắt cao áp, Ở khe hở khơng khí đóng vai trị việc chuyển đổi lượng ñiện – chuyển ñổi lượng ñiện ln kèm với thay đổi lớn khe hở khơng khí c) Mạch từ có phần ứng phần cảm quay ðó loại mạch từ thường gặp máy ñiện quay Trong mạch từ loại này, biến ñổi lượng Chương 5: Ứng dụng phần mềm Matlab ñể lập trình tính tốn thiết kế tối ưu nam châm điện dạng Solenoid No Sdu F = < Fα S th (α = 0.05) Yes Ước luợng ñược phương trình hồi quy thực nghiệm (bậc 2) Tính tốn tối ưu với hàm mục tiêu C(x, y, z, v) ràng buộc B, F Tìm giá trị tối ưu tương ứng x*, y*, z*, v* - Tìm giá trị: c, χ, η, λ - Tìm giá trị C, B, F, P Kết thúc 95 Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm ñiện dạng Solenoid 5.2.3 Kết tính tốn tốn quy hoạch hàm phi tuyến ràng buộc Kết sau lập trình tính tốn tối ưu ta xác định thơng số x, y z, v, từ tính giá trị F*, C*, B*, P* sau: Biến x y z v Giá trị 1.2 0.6648 1.5 3.1 Từ giá trị hệ ñơn vị tương ñối (x, y, z, v) ta tính thơng số nam châm dạng solenoid bao gồm: - Kích thước: c, δ, η, λ, l; - Giá thành C, lực ñiện từ F, từ cảm B, công suất P (thay biến x, y, z, v vừa tìm vào biểu thức C* = C ⇒ C = C* × ( γ Fe δ G Fe ) γ Fe δ G Fe F* = F.ρ τ F* × k τ τk 3µ o δ ⇒ F = k τ τk 3µ o δ ρτ B B = µo * P* = ρτ ⇒B= 2k τ τk δ B* ì o 2k k P ⇒ P = P * × k τ τ.δ 2 k τ τ.δ Trong đó: thơng số kτ, k3, τ, ρτ, µ0, γFe, GFe xác định mục 4.3.1 (Chương 4) Ví dụ: Cho rc = 0.012m, áp dụng công thức - δ η λ c = x , = y , = z , = v ta ñược rc rc rc rc c = x.rc = 1.2*0.012 = 0.0144m 96 Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm ñiện dạng Solenoid - δ = y*rC = 0.6648*0.012 = 0.00798m - η = z*rC = 1.5*0.012 = 0.018m - λ = v*rC = 3.1*0.012 = 0.0372m - l = δ + η + λ = 0.0632m * Giá thành tương ứng: C* = C γ Fe δ G Fe π π πx ⇒ C =  x − ( z + y + v )k 3σ + (2 z + y + 2v + 1) + y y y ( )   xγ Fe δ G Fe   3.14 π π * 1.2 2 ( ) ( ) 1 6648 * 10 57 * * * 1 =  − + + + + + + + 0.66483 0.6648  0.6648 = 17617(VND) ( ) * Lực từ tương ứng δ−1: (x + β )(x − 1) ×   π   (x + 1)   y sin α    sin α sin α  sin α    F ρτ  *   * k τk µ δ  z + vchy + F = ⇒ F =   shy  ln x ln x  ln x    τ o kτ τk µ oδ   x  sin α    shy   ln x      (1.2 + 1.7 )(1.2 − 1) × π     (1.2 + 1)   0.6648 sin 45    sin 45    sin 45 sin 45   1.5 + 3.1 * ch(0.6648 * +  * =   sh(0.6648 * ln 1.2    ln 1.2 ln 1.2    x  sin 45    sh(0.6648 *   ln 1.2      −7 *11 * 85 * 0.5 * 4π *10 * 0.00798 = 5.36( N ) * Tương tự ta có lực từ tương ứng δ−2: F = 59191 (N) 97 Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm ñiện dạng Solenoid * Từ cảm tương ứng δ−1: B* = B µo ρτ 2kτ τk 3δ  (x + β )(x − 1) ×  (x + 1)  y sin α    sin α sin α  sin α  z + vchy ⇒ B =   shy + ln x ln x  ln x  x   sin α shy     ln x     (1.2 + 1.7 )(1.2 − 1) ×   (1.2 + 1)   0.6648 sin 45   sin 45    sin 45  sin 45   1.5 + 3.1ch(0.6648 )    * )+ =   sh(0.6648  ln ln ln    x    sin 45   ) sh(0.6648 y    ln      * * 3.14 *10 −7 211 * 85 * 0.5 * 0.00798 2.65 * 10 −8τ = 0.88(T ) * Tương tự từ cảm tương ứng δ−2: B =20.64(T) * Công suất tương ứng : P* =  2π  P ⇒ P =  ( x + β )( z + y + v ) * kτ τ δ 2 kτ τ δ y   2π (1.2 + 1.7)(1.5 + 0.6648 + 3.1) *11 * 85 * 0.00798 = 12.9(W) =   0.6648 98      2kτ τk 3δ    * µo  ρτ       4.4 5.7 7.3 9.0 10.8 12.9 15.2 17.6 20.2 23.0 25.9 29.0 32.4 35.9 39.5 43.4 47.4 51.6 Công suất P (W) 0.0070 0.0080 0.0090 0.0100 0.0110 0.0120 0.0130 0.0140 0.0150 0.0160 0.0170 0.0180 0.0190 0.0200 0.0210 0.0220 0.0230 0.0240 BK cuộn dây rc (m) 0.0084 0.0096 0.0108 0.0120 0.0132 0.0144 0.0156 0.0168 0.0180 0.0192 0.0204 0.0216 0.0228 0.0240 0.0252 0.0264 0.0276 0.0288 BK cuộn dây c (m) 0.00465 0.00532 0.00598 0.00665 0.00731 0.00798 0.00864 0.00931 0.00997 0.01064 0.01130 0.01197 0.01263 0.01330 0.01396 0.01463 0.01529 0.01596 Kích thước δ−1 (m) 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 0.00005 Kích thước δ-2 (m) 0.01050 0.01200 0.01350 0.01500 0.01650 0.01800 0.01950 0.02100 0.02250 0.02400 0.02550 0.02700 0.02850 0.03000 0.03150 0.03300 0.03450 0.03600 Kích thước η (m) 0.02630 0.03007 0.03383 0.03760 0.04136 0.04513 0.04889 0.05266 0.05642 0.06019 0.06395 0.06772 0.07148 0.07525 0.07901 0.08278 0.08654 0.09031 Kích thước λ (m)-2 99 0.02170 0.02480 0.02790 0.03100 0.03410 0.03720 0.04030 0.04340 0.04650 0.04960 0.05270 0.05580 0.05890 0.06200 0.06510 0.06820 0.07130 0.07440 Kích thước λ (m) Các kích thước nam châm điện dạng Solenoid 0.0369 0.0421 0.0474 0.0526 0.0579 0.0632 0.0684 0.0737 0.0790 0.0842 0.0895 0.0948 0.1000 0.1053 0.1106 0.1158 0.1211 0.1264 Kích thước l (m) (khi hàm mục tiêu giá thành cực tiểu) 3497 5220 7432 10195 13569 17617 22398 27974 34407 41758 50087 59456 69926 81558 94413 108554 124039 140932 Giá thành (VNð) Hàm mục tiêu NAM CHÂM ðIỆN DẠNG SOLENOI 1.06 1.59 2.26 3.10 4.13 5.36 6.82 8.51 10.47 12.71 15.24 18.10 21.28 24.82 28.74 33.04 37.75 42.89 Lực từ F (N) ứng giá trị δ-1 BẢNG TĨM TẮT KẾ QUẢ TÍNH TỐN THIẾT KẾ TỐI ƯU Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm ñiện dạng Solenoid 3998 7795 14047 23788 38311 59191 88321 127933 180633 249425 337740 449468 588983 761174 971470 1225873 1530984 1894031 Lực từ F (N) ứng giá trị δ-2 0.67 0.72 0.76 0.81 0.85 0.88 0.92 0.95 0.99 1.02 1.05 1.08 1.11 1.14 1.17 1.20 1.22 1.25 Từ cảm B (T) ứng giá trị δ-1 9.52 11.51 13.62 15.85 18.19 20.64 23.18 25.82 28.55 31.38 34.29 37.28 40.36 43.52 46.76 50.07 53.46 56.92 Từ cảm B (T) ứng giá trị δ-2 Chương 5: Ứng dụng phần mềm Matlab ñể lập trình tính tốn thiết kế tối ưu nam châm điện dạng Solenoid Phụ lục CHƯƠNG TRÌNH TÍNH TỐN THIẾT KẾ TỐI ƯU NAM CHÂM ðIỆN DẠNG SOLENOID (Cực tiểu hàm giá thành vật liệu tích cực) Hàm Est_pthq: Xây dựng phương trình hồi quy thực nghiệm %% %clear; load temp; load temp1; load table_dis; syms x1 x2 x3 x4 ; %% %Thong so Bai toan QHTT k=4; n0=1;%So nghiem trung tam %% %Uoc luong cac he so tinh toan(thoa man dieu kien truc giao) %Cac he so anfa va lamda dc tinh toan cho ma tran mat_X mang tinh chat %cua ma tran truc giao N=2^k+n0+2*k; anfa=sqrt(sqrt(N*2^(k-2))-2^(k-1)); lamda=1/N*(2^k+2*anfa^2); %mat_X(2^k+1,2)=anfa; %mat_X(2^k+2,2)=-anfa; %mat_X(2^k+2,3)=anfa; %mat_X(2^k+3,3)=-anfa; %% %matrix P1 for i=1:k for j=1:2*k mat_X(2^k+j+1,i+11)=-lamda; end; end; %matrix P2 for i=1:k mat_X(2^k+2*i,i+1)=anfa; 100 Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm ñiện dạng Solenoid mat_X(2^k+2*i+1,i+1)=-anfa; mat_X(2^k+2*i,i+11)=anfa^2-lamda; mat_X(2^k+2*i+1,i+11)=anfa^2-lamda; end; %matrix P3 for i=1:k for j=1:n0 mat_X(2^k+j,i+11)=-lamda; end; end; %matrix P4 for i=1:2^k for j=1:k mat_X(i,j+11)=1-lamda; end; end; %% %Column tri thuc nghiem %var=[var1;var0;var2]; var=Cost; %matrix X X=[mat_X var]; %% %Uoc luong matrix thong so xac suat B(bij) b0=sum(var)/N; b=zeros(k,k); %Part1 for i=1:k bi(:,i)=mat_X(:,i+1)'*var; end; bi=bi/(2^k+2*anfa^2); %Part2 for i=1:6 bij(:,i)=mat_X(:,i+k+1)'*var; end; bij=bij/(2^k); %Part3 for i=1:k bii(:,i)=mat_X(:,i+11)'*var; end; bii=bii/(2*anfa^4); 101 Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm điện dạng Solenoid % for i=1:k bi(2,i)=0; end; %Part2 for i=1:6 bij(2,i)=0; end; %Part3 for i=1:k bii(2,i)=0; end; %% %Phan: bo sung them cac thi nghiem o tam(co the lay theo lay theo % thi nghiem o cac tam n0) % y_0=sum(y0)/max(size(y0)); %% %Uoc luong tham so phuong sai tai sinh %s2ts=0; %[r0 c0]=size(y0); %y_0=sum(y0)/max(size(y0)); %for i=1:r0 %s2ts=s2ts+(y0(i,1)-y_0)^2; %%end; %s2ts=s2ts/(r0-1); % -s2ts=0; [r0 c0]=size(y0); for i=1:r0 y0(i,1)=var(17,1)*(1+0.2*i) end; y_0=sum(y0)/max(size(y0)); for i=1:r0 s2ts=s2ts+(y0(i,1)-y_0)^2; end; s2ts=s2ts/(r0-1); %% %Cac he so phuong sai khac 102 Chương 5: Ứng dụng phần mềm Matlab để lập trình tính tốn thiết kế tối ưu nam châm ñiện dạng Solenoid sb0=sqrt(s2ts/N); sbi=sqrt(s2ts/(2^k+2*anfa^2)); sbij=sqrt(s2ts/(2^k)); sbii=sqrt(s2ts/(2*anfa^4)); %% %Uoc luong cac he so thong ke tij t0=b0/sb0; ti=bi/sbi; tij=bij/sbij; tii=bii/sbii; %% %Tra bang phan phoi Student [r1 c1]=size(table_stu); %Chon he so anfa_index=0.05 anfa_index=0.05; for i=1:r1 for j=1:c1 if (((r0-1)==table_stu(i,1)) && ((1anfa_index)==table_stu(1,j))) t_anfa=table_stu(i,j); end; end; end; %% %Xay dung phuong trinh hoi quy thuc nghiem if abs(t0)