Phương pháp điều khiển trực tiếp tải điều hòa

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH ĐẠI HỌC BÁCH KHOA TP.HCM KHOA ĐIỆN - ĐIỆN TỬ THIẾT BỊ, MẠNG VÀ NHÀ MÁY ĐIỆN -WOX -NIÊN KHÓA 2011 - 2013 LUẬN VĂN CAO HỌC PHƯƠNG PHÁP ĐIỀU KHIỂN TRỰC TIẾP TẢI ĐIỀU HỊA GVHD : PGS-TS PHAN THỊ THANH BÌNH SVTH : NGUYỄN HỒNG KHOA MSHV : 11186121 TP.Hồ Chí Minh, 11 / 2012 ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: Nguyễn Hoàng Khoa MSHV: 11186121 Ngày, tháng, năm sinh: 08/10/1988 Nơi sinh: TPHCM Chuyên ngành: Thiết bị, mạng nhà máy điện Mã số : 605250 I TÊN ĐỀ TÀI: Phương pháp điều khiển trực tiếp tải điều hòa II NHIỆM VỤ VÀ NỘI DUNG: Tìm hiểu phương pháp điều khiển trực tiếp tải điều hịa, thực chương trình trình bày phương pháp điều khiển theo nhóm g-DLC III NGÀY GIAO NHIỆM VỤ : (Ghi theo QĐ giao đề tài) 02/07/2012 IV NGÀY HOÀN THÀNH NHIỆM VỤ: (Ghi theo QĐ giao đề tài) 30/11/2012 V CÁN BỘ HƯỚNG DẪN (Ghi rõ học hàm, học vị, họ, tên): PGS-TS Phan Thị Thanh Bình Tp HCM, ngày 02 tháng 12 năm 2012 CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) CHỦ NHIỆM BỘ MÔN ĐÀO TẠO (Họ tên chữ ký) TRƯỞNG KHOA….……… (Họ tên chữ ký) Học viên thực hiện: Nguyễn Hồng Khoa TĨM TẮT LUẬN VĂN Nhu cầu điện tăng 15% / năm từ đến 2020 nước ta tình trạng thiếu điện, cao điểm Để đáp ứng nhu cầu lượng giảm tổn thất mức 7%, ngành điện Việt Nam cần tiếp cận công nghệ tiên tiến, nâng cấp hệ thống cũ cải thiện tình hình Dựa ý tưởng Chi-Min Chu Tai-Lang Jong, luận văn nêu cách tiếp cận việc điều khiển trực tiếp tải điều hòa g-DLC (group-Direct Load Control), giúp tối ưu hóa việc cắt giảm tải điều hịa mà đảm bảo hai tiêu chí: ngăn ngừa tượng trả lại tải đáp ứng công suất cắt tải điện lực đưa xuống cho nhóm điều khiển FINAL THESIS SUMMARY Electric demand will rise 15% / year from now to 2020 and our country still lacks electric, especially in rush-hours To satisfy the energy demand and reduce electric loss to 7%, EVN needs to reach advanced technologies, upgrade old systems to improve this status Based on idea of Chi-Min Chu and Tai-Lang Jong, this final thesis raise a method for air-conditioning direct load control g-DLC (group-Direct Load Control), to optimize shedding air-conditioning loads but remaining ensure two conditions: prevent transient payback load phenomenon and satisfy load shedding demand for each control group Học viên thực hiện: Nguyễn Hoàng Khoa LUẬN VĂN CAO HỌC ĐỀ TÀI: PHƯƠNG PHÁP ĐIỀU KHIỂN TRỰC TIẾP TẢI ĐIỀU HỊA GVHD: PGS-TS Phan Thị Thanh Bình HVTH: Nguyễn Hồng Khoa LỜI CAM ĐOAN Học viên xin cam đoan sản phẩm riêng học viên Các số liệu sử dụng luận văn trung thực, có nguồn gốc rõ ràng Những kết trình bày luận văn lấy từ chương trình mơ học viên chưa xuất công trình nghiên cứu khác Contents  DANH MỤC CÁC CHỮ VIẾT TẮT TRONG LUẬN VĂN 3  DANH SÁCH CÁC HÌNH VẼ 4  DANH SÁCH CÁC BẢNG BIỂU 6  LỜI MỞ ĐẦU 7  * Tính cấp thiết đề tài: 7  * Ý nghĩa thực tiễn: 8  * Mục tiêu luận văn: 9  CHƯƠNG 1: TỔNG QUAN VỀ DLC VÀ ĐIỀU HỊA KHƠNG KHÍ 10  I/ Khái niệm DLC: 10  II/ Nhóm điều khiển tải trước (sớm): 11  III/ Nhóm điều khiển tải thời gian thực: 12  Page 1    Học viên thực hiện: Nguyễn Hoàng Khoa IV/ Hiện trạng sử dụng lượng nước ta: 13  V/ Các chương trình quản lý phụ tải Việt Nam: 13  VI/ Quản lý tải điều hòa: 16  CHƯƠNG 2: MỘT CÁCH TIẾP CẬN TRONG ĐIỀU KHIỂN TẢI ĐIỀU HÒA 19  I/ Giới thiệu: 19  II/ Điều khiển độ thoải mái nhiệt độ hệ thống điều hịa khơng khí: 20  III/Ảnh hưởng nhiệt độ trời đến tải điều hòa: 22  IV/ Nhiệt độ nước ướp lạnh thay đổi thời gian nhàn rỗi (off shift): 23  V/ Hiện tượng trả lại tải trình DLC: 24  VI/ Những hạn chế tải sa thải: 25  VII/ Thuật toán DLC với TCL cho hệ thống điều hịa khơng khí: 27  VIII/ Thời gian ON bắt buộc: 31  IX/ Phương án nhóm điều khiển trực tiếp tải (g-DLC): 33  CHƯƠNG 3: ĐIỀU KHIỂN TRỰC TIẾP TẢI THEO NHÓM 35  I/Mô tả ý tưởng, công thức cần thiết thuật toán: 35  1/ Mô tả ý tưởng: 35  2/ Các công thức cần thiết: 35  3/ Thuật toán: 43  II/Chạy chương trình cho số liệu khác: 53  IV/ Nhận xét: 62  CHƯƠNG 4: KẾT LUẬN VÀ ĐỀ XUẤT 63  I/ Những kết thu được: 63  II/ Kết luận hướng phát triển: 63  *Các tài liệu tham khảo: 65  PHỤ LỤC 66  I/ Chương trình chính: 66  II/ Các hàm phụ: 91    Page 2    Học viên thực hiện: Nguyễn Hoàng Khoa DANH MỤC CÁC CHỮ VIẾT TẮT TRONG LUẬN VĂN Tên viết tắt Tên đầy đủ Nghĩa DLC Direct Load Control Điều khiển tải trực tiếp DSM Demand Side Management Quản lý nhu cầu phụ tải LEE Least Enthalpy Estimation Ước lượng enthalpy tối thiểu TCL Thermal Comfort Level Mức độ thoải mái nhiệt LMP Load Managemnet Chương trình quản lý tải Program BTC Based Thermal Comfort Dựa vào thoải mái nhiệt FC Fuzzy Control Điều khiển mờ SCADA Supervisory Control and Giám sát kiểm soát thu thập AC Data Acquisition liệu Air Conditioning Điều hịa khơng khí – máy lạnh HVAC ĐHKK Heating, Ventilating and Hệ thống sưởi, thơng gió Air Conditioning điều hịa khơng khí Điều hịa khơng khí Điều hịa khơng khí Page 3    Học viên thực hiện: Nguyễn Hồng Khoa DANH SÁCH CÁC HÌNH VẼ Thứ tự Hình Tên hình vẽ Trang Cơ cấu tiêu thụ điện tịa nhà 17 Hình 2.1 Điều kiện nhiệt độ mơi trường hộ tiêu thụ 21 Hình 2.2 Hai phòng điều hòa dùng chung hệ thống AC 22 Hình 2.3 Biểu đồ so sánh lượng tiêu tốn sử dụng kĩ 23 thuật FC CV Hình 2.4 So sánh nhiệt độ nước làm lạnh phương pháp 24 FC CV Hình 2.5 Hiện tượng trả lại tải điều hịa khơng khí 25 Hình 2.6 Thuật tốn cho phương pháp điều khiển trực tiếp LEE 27 Hình 2.7 Biểu đồ kết thực nghiệm (thời gian off 15 29 phút/giờ) Hình 2.8 Biểu đồ kết thực nghiệm (thời gian off 25 30 phút/giờ) Hình 2.9 Biểu đồ group-DLC 32 Hình 3.1 Chuỗi thời gian thể đặc tính thời gian 35 chu kỳ DLC Hình 3.2 Ma trận chưa xác định Ma trận xác định 38 Hình 3.3 Biểu đồ cắt tải theo thời gian sau thực 39 phương pháp g-DLC Hình 3.4 Giao diện điều khiển mơ lập trình GUI 41 Hình 3.5 Ma trận xác định cho liệu nhập từ Hình 3.4 46 Hình 3.6 Kết thu từ mơ máy lạnh chu kỳ 49 [12 stages – 5min per stage] Page 4    Học viên thực hiện: Nguyễn Hồng Khoa Hình 3.7 Kết thu từ mô máy lạnh chu kỳ 50 [12 stages – 5min per stage] Hình 3.8 Kết thu từ mô máy lạnh chu kỳ 51 [14 stages – per stage] Hình 3.9 Kết thu từ mô máy lạnh chu kỳ 52 [10 stages – per stage] Hình 3.10 Kết thu từ mơ máy lạnh chu kỳ 53 [10 stages – per stage] Hình 3.11 Kết thu từ mô máy lạnh chu kỳ 54 [10 stages – per stage] Hình 3.12 Kết thu từ mô máy lạnh chu kỳ 56 [10 stages – per stage] Hình 3.13 Kết thu từ mơ máy lạnh chu kỳ 57 [12 stages – per stage] – lỗi Hình 3.14 Kết thu từ mô máy lạnh chu kỳ [12 stages – per stage] – hợp lý Page 5    59 Học viên thực hiện: Nguyễn Hoàng Khoa DANH SÁCH CÁC BẢNG BIỂU Thứ tự Tên bảng biểu Trang Table I Bảng thời gian ON bắt buộc sau OFF tương ứng với 31 thời lượng OFF Table II Bảng thời gian ON bắt buộc sau OFF tương ứng với thời lượng OFF giá trị công suất cắt từ 10 – 20 kW Page 6    39 Học viên thực hiện: Nguyễn Hoàng Khoa LỜI MỞ ĐẦU * Tính cấp thiết đề tài: Hiện nước ta, tình trạng thiếu điện chưa khắc phục, khơng muốn nói cịn phụ thuộc nhiều vào Thủy điện Nhiệt điện Hiện tượng khan điện vùng nông thôn phổ biến, đồng thời thành phố phải trì lịch cắt điện định kỳ, đặc biệt vào mùa khô Chưa kể cố phát sinh khiến việc cắt điện trở nên thường xuyên số tỉnh, thành nước Song song với việc cắt điện - thiếu điện, ngành điện phải đối mặt với vấn đề đáp ứng nhu cầu phụ tải cao điểm Chính ngành điện chi khoản lớn để xây dựng nguồn điện nhằm đáp ứng nhu cầu trên, vấn đề khác lại nảy sinh: thấp điểm, nguồn điện khơng phát huy hiệu quả, gây lãng phí đầu tư Tỷ lệ Pmin / Pmax nước ta khoảng 0.4 – 0.7 [7], tỷ lệ thấp so với nước giới Do đó, Chính phủ đề chương trình tiết kiệm điện “Đảm bảo sử dụng điện tiết kiệm hiệu quả” giai đoạn 2007 – 2015 [8] Mục tiêu giảm 5% phụ tải đỉnh Một giải pháp hồn thiện chương trình quản lý nhu cầu phụ tải (Demand Side Management - DSM) điều khiển trực tiếp tải (Direct Load Control - DLC) Tại Việt Nam, DLC thử nghiệm tải tiêu thụ lớn mang tính chất báo hiệu hệ thống thiếu hụt công suất chưa tùy Page 7    Học viên thực hiện: Nguyễn Hoàng Khoa end %Step2: Tinh Vkm val_vkm = step2_7u(val_dkm(1,i),val_stage,val_minstage,val_p ow01,val_pow02,val_pow03,val_pow04 ,val_pow05,val_pow06,val_pow07); val_vkmplus(1,i) = val_vkm; %Lay gia tri cua Wkm val_wkm01=step0(val_pow01,(round(val_vkm/5)*5)); val_wkm02=step0(val_pow02,(round(val_vkm/5)*5)); val_wkm03=step0(val_pow03,(round(val_vkm/5)*5)); val_wkm04=step0(val_pow04,(round(val_vkm/5)*5)); val_wkm05=step0(val_pow05,(round(val_vkm/5)*5)); val_wkm06=step0(val_pow06,(round(val_vkm/5)*5)); val_wkm07=step0(val_pow07,(round(val_vkm/5)*5)); %Step3: Ma tran Ckm %buoc du, ko can lam cung dc vi da tich hop o step4 %val_ckm = step3(val_n,val_stage,val_lkm01,val_lkm02,val_lkm0 3,val_lkm04,val_lkm05,val_lkm06); %Step4: Ma tran Ckm moi va gia tri least square error val_ckmao = step4_7u(val_stage,val_minstage,val_vkm,val_lkm01, val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_pow01, val_pow02,val_pow03,val_pow04,val_pow05,val_pow06, val_pow07,val_dkm(1,i)); val_kqm = step4_7u_kq(val_stage,val_minstage,val_vkm,val_lkm 01,val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_pow01, val_pow02,val_pow03,val_pow04,val_pow05,val_pow06, Page 84    Học viên thực hiện: Nguyễn Hoàng Khoa val_pow07,val_dkm(1,i)); val_kqmplus(1,i) = val_kqm; val_dif = step4_7u_d(val_stage,val_minstage,val_vkm,val_lkm0 1,val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_pow01, val_pow02,val_pow03,val_pow04,val_pow05,val_pow06, val_pow07,val_dkm(1,i)); val_difplus(1,i) = val_dif; val_per = step4_7u_p(val_stage,val_minstage,val_vkm,val_lkm0 1,val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_pow01, val_pow02,val_pow03,val_pow04,val_pow05,val_pow06, val_pow07,val_dkm(1,i)); val_perplus(1,i) = val_per; %Step5: Tinh Ukm cho tung tai val_ukm01 = step5(val_ckmao,val_stage,val_minstage,1); val_ukm02 = step5(val_ckmao,val_stage,val_minstage,2); val_ukm03 = step5(val_ckmao,val_stage,val_minstage,3); val_ukm04 = step5(val_ckmao,val_stage,val_minstage,4); val_ukm05 = step5(val_ckmao,val_stage,val_minstage,5); val_ukm06 = step5(val_ckmao,val_stage,val_minstage,6); val_ukm07 = step5(val_ckmao,val_stage,val_minstage,7); %Step6: Tinh Lkm cho buoc tiep theo vong lap val_lkmnext01 = step6(val_stage,val_minstage,val_ukm01,val_vkm,val _wkm01); Page 85    Học viên thực hiện: Nguyễn Hoàng Khoa val_lkmnext02 = step6(val_stage,val_minstage,val_ukm02,val_vkm,val _wkm02); val_lkmnext03 = step6(val_stage,val_minstage,val_ukm03,val_vkm,val _wkm03); val_lkmnext04 = step6(val_stage,val_minstage,val_ukm04,val_vkm,val _wkm04); val_lkmnext05 = step6(val_stage,val_minstage,val_ukm05,val_vkm,val _wkm05); val_lkmnext06 = step6(val_stage,val_minstage,val_ukm06,val_vkm,val _wkm06); val_lkmnext07 = step6(val_stage,val_minstage,val_ukm07,val_vkm,val _wkm07); %end %StepPlus: Ve thi val_powtb_1 = stepplus_7u(val_stage,val_ckmao,val_pow01,val_pow0 2,val_pow03,val_pow04, val_pow05,val_pow06,val_pow07); val_dkmplus_1 = val_dkm(1,i) * ones(1,val_stage); val_powtb=[val_powtb;val_powtb_1]; val_dkmplus=[val_dkmplus;val_dkmplus_1]; end %set(handles.sta01,'string',val_powtb(3,:)); val_powtb_ve_tam=0; val_dkmplus_ve_tam=0; for i=1:val_numcycle val_powtb_ve_tam=[val_powtb_ve_tam val_powtb(i+1,:)]; val_dkmplus_ve_tam=[val_dkmplus_ve_tam val_dkmplus(i+1,:)]; end val_powtb_ve = val_powtb_ve_tam(2:length(val_powtb_ve_tam)); Page 86    Học viên thực hiện: Nguyễn Hoàng Khoa val_dkmplus_ve = val_dkmplus_ve_tam(2:length(val_dkmplus_ve_tam)); timeset = (val_stage * val_minstage) * val_numcycle - val_minstage; time = 0:(val_minstage):(timeset); set(handles.sta01,'string','Done'); set(handles.stavkm, 'String', val_vkmplus); set(handles.stakqm, 'String', val_kqmplus); axes(handles.axes1); stairs(handles.axes1,time,val_powtb_ve); grid on hold on; stairs(handles.axes1,time,val_dkmplus_ve,'g '); hold off; set(gca,'XTick',[0:(val_stage * val_minstage):(val_stage * val_minstage * 10)]); end %////////////// %Khi numload=8 if val_n == for i=1:val_numcycle if i>1 val_lkm01=val_lkmnext01; val_lkm02=val_lkmnext02; val_lkm03=val_lkmnext03; val_lkm04=val_lkmnext04; val_lkm05=val_lkmnext05; val_lkm06=val_lkmnext06; val_lkm07=val_lkmnext07; val_lkm08=val_lkmnext08; end %Step2: Tinh Vkm val_vkm = step2_8u(val_dkm(1,i),val_stage,val_minstage,val_p ow01,val_pow02, Page 87    Học viên thực hiện: Nguyễn Hoàng Khoa val_pow03,val_pow04,val_pow05,val_pow06,val_pow07, val_pow08); val_vkmplus(1,i) = val_vkm; %Lay gia tri cua Wkm val_wkm01=step0(val_pow01,(round(val_vkm/5)*5)); val_wkm02=step0(val_pow02,(round(val_vkm/5)*5)); val_wkm03=step0(val_pow03,(round(val_vkm/5)*5)); val_wkm04=step0(val_pow04,(round(val_vkm/5)*5)); val_wkm05=step0(val_pow05,(round(val_vkm/5)*5)); val_wkm06=step0(val_pow06,(round(val_vkm/5)*5)); val_wkm07=step0(val_pow07,(round(val_vkm/5)*5)); val_wkm08=step0(val_pow08,(round(val_vkm/5)*5)); %Step3: Ma tran Ckm %buoc du, ko can lam cung dc vi da tich hop o step4 %val_ckm = step3(val_n,val_stage,val_lkm01,val_lkm02,val_lkm0 3,val_lkm04,val_lkm05,val_lkm06); %Step4: Ma tran Ckm moi va gia tri least square error val_ckmao = step4_8u(val_stage,val_minstage,val_vkm,val_lkm01, val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_lkm08, val_pow01,val_pow02,val_pow03,val_pow04,val_pow05, val_pow06, val_pow07,val_pow08,val_dkm(1,i)); val_kqm = step4_8u_kq(val_stage,val_minstage,val_vkm,val_lkm 01,val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_lkm08, val_pow01,val_pow02,val_pow03,val_pow04,val_pow05, val_pow06, val_pow07,val_pow08,val_dkm(1,i)); val_kqmplus(1,i) = val_kqm; Page 88    Học viên thực hiện: Nguyễn Hoàng Khoa val_dif = step4_8u_d(val_stage,val_minstage,val_vkm,val_lkm0 1,val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_lkm08, val_pow01,val_pow02,val_pow03,val_pow04,val_pow05, val_pow06, val_pow07,val_pow08,val_dkm(1,i)); val_difplus(1,i) = val_dif; val_per = step4_8u_p(val_stage,val_minstage,val_vkm,val_lkm0 1,val_lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06,val_lkm07,val_lkm08, val_pow01,val_pow02,val_pow03,val_pow04,val_pow05, val_pow06, val_pow07,val_pow08,val_dkm(1,i)); val_perplus(1,i) = val_per; %Step5: Tinh Ukm cho tung tai val_ukm01 = step5(val_ckmao,val_stage,val_minstage,1); val_ukm02 = step5(val_ckmao,val_stage,val_minstage,2); val_ukm03 = step5(val_ckmao,val_stage,val_minstage,3); val_ukm04 = step5(val_ckmao,val_stage,val_minstage,4); val_ukm05 = step5(val_ckmao,val_stage,val_minstage,5); val_ukm06 = step5(val_ckmao,val_stage,val_minstage,6); val_ukm07 = step5(val_ckmao,val_stage,val_minstage,7); val_ukm08 = step5(val_ckmao,val_stage,val_minstage,8); %Step6: Tinh Lkm cho buoc tiep theo vong lap val_lkmnext01 = step6(val_stage,val_minstage,val_ukm01,val_vkm,val _wkm01); Page 89    Học viên thực hiện: Nguyễn Hoàng Khoa val_lkmnext02 = step6(val_stage,val_minstage,val_ukm02,val_vkm,val _wkm02); val_lkmnext03 = step6(val_stage,val_minstage,val_ukm03,val_vkm,val _wkm03); val_lkmnext04 = step6(val_stage,val_minstage,val_ukm04,val_vkm,val _wkm04); val_lkmnext05 = step6(val_stage,val_minstage,val_ukm05,val_vkm,val _wkm05); val_lkmnext06 = step6(val_stage,val_minstage,val_ukm06,val_vkm,val _wkm06); val_lkmnext07 = step6(val_stage,val_minstage,val_ukm07,val_vkm,val _wkm07); val_lkmnext08 = step6(val_stage,val_minstage,val_ukm08,val_vkm,val _wkm08); %end %StepPlus: Ve thi val_powtb_1 = stepplus_8u(val_stage,val_ckmao,val_pow01,val_pow0 2,val_pow03,val_pow04,val_pow05,val_pow06,val_pow0 7,val_pow08); val_dkmplus_1 = val_dkm(1,i) * ones(1,val_stage); val_powtb=[val_powtb;val_powtb_1]; val_dkmplus=[val_dkmplus;val_dkmplus_1]; end %set(handles.sta01,'string',val_powtb(3,:)); val_powtb_ve_tam=0; val_dkmplus_ve_tam=0; for i=1:val_numcycle val_powtb_ve_tam=[val_powtb_ve_tam val_powtb(i+1,:)]; val_dkmplus_ve_tam=[val_dkmplus_ve_tam val_dkmplus(i+1,:)]; end Page 90    Học viên thực hiện: Nguyễn Hoàng Khoa val_powtb_ve = val_powtb_ve_tam(2:length(val_powtb_ve_tam)); val_dkmplus_ve = val_dkmplus_ve_tam(2:length(val_dkmplus_ve_tam)); timeset = (val_stage * val_minstage) * val_numcycle - val_minstage; time = 0:(val_minstage):(timeset); set(handles.sta01,'string','Done'); set(handles.stavkm, 'String', val_vkmplus); set(handles.stakqm, 'String', val_kqmplus); axes(handles.axes1); stairs(handles.axes1,time,val_powtb_ve); grid on hold on; stairs(handles.axes1,time,val_dkmplus_ve,'g '); hold off; set(gca,'XTick',[0:(val_stage * val_minstage):(val_stage * val_minstage * 10)]); end II/ Các hàm phụ: Step 0: Xác định thời gian ON bắt buộc sau OFF máy lạnh dựa vào bảng Table II: function wkm=step0(tai,vkm) mat=[15 20 25 30;20 20 25 30; 20 20 25 30;20 25 30 35;20 25 30 35;25 25 30 35; 25 25 30 35;25 30 30 35;25 30 35 40;30 30 35 40;30 30 35 40]; if vkm = 25 for i=0:10 if tai == (10+i) wkm = mat(i+1,4); end end end end Step 2: Xác định thời gian OFF chu kỳ xét dựa vào công thức (1): function val_vkm= step2(val_dkm,val_stage,val_minstage,val_tai01 ,val_tai02, val_tai03,val_tai04,val_tai05,val_tai06) val_vkm= val_dkm * val_stage * val_minstage / (val_tai01 + val_tai02 + val_tai03+ val_tai04+ val_tai05+ val_tai06); end Page 92    Học viên thực hiện: Nguyễn Hoàng Khoa Step 3: Xây dựng ma trận chưa xác định dựa vào giá trị biết: [lấy code phần lập trình máy lạnh group, số lượng khác tương tự] function val_ckm=step3_6u(val_n,val_stage,val_lkm01,val _lkm02,val_lkm03, val_lkm04,val_lkm05,val_lkm06) %val_ckm = * ones(val_n,val_stage); val_ckm = zeros(val_n,val_stage); for i = : (val_lkm01) val_ckm(1,i) = 1; end for i = : (val_lkm02) val_ckm(2,i) = 1; end for i = : (val_lkm03) val_ckm(3,i) = 1; end for i = : (val_lkm04) val_ckm(4,i) = 1; end for i = : (val_lkm05) val_ckm(5,i) = 1; end for i = : (val_lkm06) val_ckm(6,i) = 1; end end Step 4: Tìm ma trận xác định tính tốn độ chênh lệch cơng suất cắt thực tế yêu cầu: [lấy code phần lập trình máy lạnh group, số lượng khác tương tự] function [mat_ckm kqua kqua_min]=step4_6u_full_tn(val_stage,val_minst age, Page 93    Học viên thực hiện: Nguyễn Hoàng Khoa val_v,val_lkm1,val_lkm2,val_lkm3,val_lkm4,val_ lkm5,val_lkm6, tai1,tai2,tai3,tai4,tai5,tai6,dkm) dk1=round(val_v/val_minstage); ctam= ones(6,val_stage); %Ma tran Ckm ctam_nd= zeros(6,val_stage); %Ma tran phu de tinh cong suat off-shift dem=1; tai=zeros(1,6); %Chuoi chua gia tri cong suat cat cua cac tai val_lkm1 val_lkm2 val_lkm3 val_lkm4 val_lkm5 val_lkm6 val_lkm1a hien vong val_lkm2a val_lkm3a val_lkm4a val_lkm5a val_lkm6a = = = = = = round(val_lkm1); round(val_lkm2); round(val_lkm3); round(val_lkm4); round(val_lkm5); round(val_lkm6); = round(val_lkm1); %Bien tam de thuc lap = round(val_lkm2); = round(val_lkm3); = round(val_lkm4); = round(val_lkm5); = round(val_lkm6); tai(1,1)=tai1; tai(1,2)=tai2; tai(1,3)=tai3; tai(1,4)=tai4; tai(1,5)=tai5; tai(1,6)=tai6; skm_b=10^9; %skm before de so sanh Page 94    Học viên thực hiện: Nguyễn Hoàng Khoa d_f = zeros(6,val_stage); %ma tran d final de xuat ctama=ctam; ctamb=ctam_nd; for for for for for for for n1=1:(13-dk1-val_lkm6) m1=1:(13-dk1-val_lkm5) l1=1:(13-dk1-val_lkm4) k1=1:(13-dk1-val_lkm3) j1=1:(13-dk1-val_lkm2) i1=1:(13-dk1-val_lkm1) i=1:dk1 ctama(1,val_lkm1a+i)=0; ctamb(1,val_lkm1a+i)=1; ctama(2,val_lkm2a+i)=0; ctamb(2,val_lkm2a+i)=1; ctama(3,val_lkm3a+i)=0; ctamb(3,val_lkm3a+i)=1; ctama(4,val_lkm4a+i)=0; ctamb(4,val_lkm4a+i)=1; ctama(5,val_lkm5a+i)=0; ctamb(5,val_lkm5a+i)=1; ctama(6,val_lkm6a+i)=0; ctamb(6,val_lkm6a+i)=1; end %mat_ckm=ctam; d=zeros(1,val_stage); for j=1:val_stage for i=1:6 d(j)=d(j) + (tai(1,i) * ctamb(i,j)); end end skm=0; for j=1:val_stage skm=skm + (d(j) - dkm)^2; end kqua(dem) = skm/val_stage; dem = dem +1; if skm < skm_b %so sanh Skm voi Skm o lan tinh gan nhat, cai nao nho hon thi giu lai skm_b=skm; Page 95    Học viên thực hiện: Nguyễn Hoàng Khoa d_f = ctama; dif = 0; for j=1:val_stage dif = dif + d(j) - dkm; end end val_lkm1a=val_lkm1a+1; ctama=ctam; ctamb=ctam_nd; end val_lkm2a=val_lkm2a+1; val_lkm1a=val_lkm1; ctama=ctam; ctamb=ctam_nd; end val_lkm3a=val_lkm3a+1; val_lkm2a=val_lkm2; val_lkm1a=val_lkm1; ctama=ctam; ctamb=ctam_nd; end val_lkm4a=val_lkm4a+1; val_lkm3a=val_lkm3; val_lkm2a=val_lkm2; val_lkm1a=val_lkm1; ctama=ctam; ctamb=ctam_nd; end val_lkm5a=val_lkm5a+1; val_lkm4a=val_lkm4; val_lkm3a=val_lkm3; val_lkm2a=val_lkm2; val_lkm1a=val_lkm1; ctama=ctam; ctamb=ctam_nd; end val_lkm6a=val_lkm6a+1; val_lkm5a=val_lkm5; val_lkm4a=val_lkm4; val_lkm3a=val_lkm3; val_lkm2a=val_lkm2; val_lkm1a=val_lkm1; Page 96    Học viên thực hiện: Nguyễn Hoàng Khoa ctama=ctam; ctamb=ctam_nd; end val_lkm_max=max([val_lkm1 val_lkm2 val_lkm3 val_lkm4 val_lkm5 val_lkm6]); if val_stage+1-dk1-val_lkm_max