Dự báo phụ tải bằng phần mềm matlab đồ án tốt nghiệp khoa công nghệ điện

TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP TP.HCM KHOA CÔNG NGHỆ ĐIỆN KHÓA LUẬN TỐT NGHIỆP DỰ BÁO PHỤ TẢI BẰNG PHẦN MỀM MATLAB Sinh viên : Phạm Tiến Đạt MSSV : 1406 7251 – DHDI10D Sinh viên : Nguyễn Hoàng Huy MSSV : 1404 8821 – DHDI10C GVHD : TS Trần Thanh Ngọc THÀNH PHỐ HỒ CHÍ MINH, THÁNG 05 NĂM 2018 PHIẾU GIAO ĐỀ TÀI KHÓA LUẬN TỐT NGHIỆP Họ tên sinh viên Phạm Tiến Đạt MSSV:14067251 Nguyễn Hoàng Huy MSSV:14048821 Đề tài DỰ BÁO PHỤ TẢI BẰNG PHẦN MỀM MATLAB Nội dung - Nghiên cứu phương pháp dự báo - Tìm hiểu phần mềm MATLAB - Áp dụng dự báo phụ tải Kết Lập trình chương trình phần mềm MATLAB, từ chương trình phục vụ cho việc dự báo điện cho tương lai với liệu qua khứ Giảng viên hướng dẫn Tp HCM, ngày tháng năm Sinh viên i NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN ii MỤC LỤC PHIẾU GIAO ĐỀ TÀI KHÓA LUẬN TỐT NGHIỆP i NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN ii MỤC LỤC iii DANH SÁCH CÁC HÌNH VẼ vi DANH SÁCH CÁC BẢNG viii LỜI MỞ ĐẦU CHƯƠNG 1: HỒI QUY ĐƠN BIẾN 1.1 Khái niệm: 1.2 Cơ sở lý thuyết hồi quy đơn biến: .5 1.2.1 Các hệ số: .5 1.2.1.1 Hệ số 𝛽1, 𝛽2: 1.2.1.2 Độ lệch chuẩn (SE): .8 1.2.1.3 Kiểm định giả thuyết (T- statistic): 1.2.1.4 P-value: 10 1.2.1.5 MSE – RMSE: .10 1.2.1.6 Hệ số xác định (R-squared) 11 1.2.1.7 Hệ số xác định hiệu chỉnh( Adj R-squared) 13 1.2.2 Hàm Fitlm .13 1.2.2.1 Cú pháp 13 1.2.2.2 Ví dụ: .14 1.2.3 Hồi quy mở rộng: 15 1.2.3.1 Hồi quy qua góc toạ độ: 16 1.2.3.2 Mơ hình tuyến tính logarit: 17 CHƯƠNG 2: HỒI QUY ĐA BIẾN .19 2.1 Phân tích hồi quy đa biến: .19 2.2 Cơ sở lý thuyết: .19 2.3 Áp dụng 21 2.3.1 Dữ liệu 21 2.3.3 Kết 24 iii 2.3.4 So sánh nhận xét kết 27 CHƯƠNG 3: DỰ BÁO NHU CẦU ĐIỆN NĂNG GIAI ĐOẠN 2006-2020 30 3.1 Dân số (POP) : 41 3.2 Tốc độ tăng dân số (POPGR) (%): 42 3.3 Số người/hộ (HH): 43 3.4 Tỷ lệ điện khí hóa (ELEC) (%): 44 3.5 Lượng khách hàng khu vực dân dụng văn phòng (CUST): 45 3.6 Tỉ trọng GDP ngành: 46 3.6.1 Dự báo tỷ trọng GDP nông nghiệp (SHAG) (%) 46 3.6.2 Dự báo tỷ trọng GDP công nghiệp (SHIN) (%) 48 3.6.3 Dự báo tỷ trọng GDP dịch vụ thương mại (SHCMDB) (%) 49 3.7 Tăng trưởng kinh tế: 50 3.7.1 Dự báo GDP tổng (GDPGRDB) (Tỷ VND): .50 3.7.2 Dự báo GDP nông nghiệp (GDPAGDB) (Tỷ VND): 51 3.7.3 Dự báo GDP công nghiệp (GDPINDB) (Tỷ VND): 52 3.7.4 Dự báo GDP dịch vụ thương mại (GDPCMDB) (Tỷ VND): 53 3.8 Tốc độ tăng trưởng kinh tế: .54 3.8.1 Dự báo tốc độ tăng GDP nông nghiệp (GRAGDB) (%): .54 3.8.2 Dự báo tốc độ tăng GDP công nghiệp (GRINDB) (%): .55 3.8.3 Dự báo tốc độ tăng GDP dịch vụ thương mại (GRCMDB) (%): 56 3.9 Dự báo nhu cầu điện năng: 57 3.9.1 Nhu cầu điện công nghiệp (INELDB)(GWh): 57 3.9.2 Nhu cầu điện dân dụng (REELDB) (GWh): .58 3.9.3 Nhu cầu điện thương mại (CMELDB) (GWh): 60 3.9.4 Nhu cầu điện nông nghiệp (AGELDB) (GWh): 61 3.9.5 Nhu cầu điện khác (OTELDB) (GWh): 62 3.9.5 Tổng nhu cầu điện (FNELDB) (Gwh): 63 3.10 Dự báo tỉ lệ tổn thất (LOSSDB) (%) .64 3.11 Dự báo điện sản xuất (ELELDB) (GWh): 66 3.12 So sánh kết quả: .67 3.13 Đánh giá mơ hình dự báo 70 iv CHƯƠNG 4: KẾT LUẬN .71 TÀI LIỆU THAM KHẢO .72 PHỤ LỤC 73 LỜI CẢM ƠN 79 v DANH SÁCH CÁC HÌNH VẼ Số hình Tên hình Trang 3.1 Biểu đồ dự báo dân số Matlab Simple E giai đoạn 2006 - 2020 29 3.2 Biểu đồ dự báo tốc độ tăng dân số Matlab Simple E giai đoạn 2006 - 2020 30 3.3 Biểu đồ dự báo người/hộ Matlab Simple E giai đoạn 2006 - 2020 31 3.4 Biểu đồ dự báo điện khí hóa Matlab Simple E giai đoạn 2006 – 2020 32 3.5 Biểu đồ dự báo khách hàng Matlab Simple E giai đoạn 2006 – 2020 33 3.6 Biểu đồ dự báo tỉ trọng GDP nông nghiệp Matlab Simple E giai đoạn 2006 – 2020 35 3.7 Biểu đồ dự báo tỉ trọng GDP công nghiệp Matlab Simple E giai đoạn 2006 – 2020 36 3.8 Biểu đồ dự báo tỉ trọng GDP TM DV Matlab Simple E giai đoạn 2006 – 2020 37 3.9 Biểu đồ dự báo GDP tổng Matlab Simple E giai đoạn 2006 – 2020 38 3.10 Biểu đồ dự báo GDP nông nghiệp Matlab Simple E giai đoạn 2006 – 2020 39 3.11 Biểu đồ dự báo GDP công nghiệp Matlab Simple E giai đoạn 2006 – 2020 40 3.12 Biểu đồ dự báo GDP DVTM Matlab Simple E giai đoạn 2006 – 2020 41 3.13 Biểu đồ dự báo tốc độ tăng GDP nông nghiệp Matlab Simple E giai đoạn 2006 – 2020 42 vi 3.14 Biểu đồ dự báo tốc độ tăng GDP công nghiệp Matlab Simple E giai đoạn 2006 – 2020 43 3.15 Biểu đồ dự báo tốc độ tăng TMDV Matlab Simple E giai đoạn 2006 – 2020 44 3.16 Biểu đồ dự báo nhu cầu điện công nghiệp Matlab Simple E giai đoạn 2006 – 2020 46 3.17 Biểu đồ dự báo nhu cầu điện dân dụng Matlab Simple E giai đoạn 2006 – 2020 47 3.18 Biểu đồ dự báo nhu cầu điện thương mại Matlab Simple E giai đoạn 2006 – 2020 48 3.19 Biểu đồ dự báo nhu cầu điện nông nghiệp Matlab Simple E giai đoạn 2006 – 2020 49 3.20 Biểu đồ dự báo nhu cầu điện khác Matlab Simple E giai đoạn 2006 – 2020 50 3.21 Biểu đồ dự báo tổng nhu cầu điệnbằng Matlab Simple E giai đoạn 2006 – 2020 51 3.22 Biểu đồ dự báo tỷ lệ tổn thất Matlab Simple E giai đoạn 2006 – 2020 53 3.23 Biểu đồ dự báo điện sản xuất Matlab Simple E giai đoạn 2006 – 2020 54 vii DANH SÁCH CÁC BẢNG Số bảng Tên bảng Trang 1.1 Mẫu liệu thu nhập (X) – chi tiêu (Y) (USD) tháng 10 hộ gia đình 2.1 Dữ liệu đầu vào giai đoạn 1982-1995 18 2.2 Dữ liệu đầu vào giai đoạn 1982-1995 (tiếp theo) 18 2.3 So sánh kết Regression 23 2.4 So sánh kết Regression 24 2.5 So sánh kết Regression 24 2.6 So sánh kết Regression 25 2.7 So sánh kết Regression 25 viii Phạm Tiến Đạt – Nguyễn Hồng Huy Khóa luận tốt nghiệp LỜI MỞ ĐẦU Sự cần thiết đề tài Ngành điện ngành công nghiệp mũi nhọn, ảnh hưởng trực tiếp đến tất ngành kinh tế kinh tế Với phát triển mạnh mẽ nay, ngành điện vừa có vai trị cung cấp lượng thúc đẩy phát triển kinh tế ngành vừa trực tiếp tham gia phục vụ đời sống xã hội sinh hoạt người Nhất nước ta thúc đẩy nhanh trình cơng nghiệp hố đại hố để hội nhập với kinh tế khu vực giới nhu cầu điện ngày tăng nhanh đòi hỏi dự báo xác cần thiết Nếu cung khơng đáp ứng đủ cầu gây thiệt hại kinh tế, cầu thấp cung gây lãng phí thiệt hại kinh tế khơng ngành điện mà cịn ảnh hưởng đến ngành khác Nhằm đáp ứng nhu cầu phát triển kinh tế - xã hội, việc tiến hành dự báo nhu cầu điện để phục vụ công tác phát triển nguồn lưới điện, kịp thời phục vụ cung cấp ổn định hệ thống điện đảm bảo chất lượng điện Trong thông tư Quyết định số 07/QĐ-ĐTĐL ngày 14 tháng 03 năm 2013 Cục trưởng Cục Điều tiết điện lực ban hành quy định nguyên tắc, trình tự, thủ tục phương pháp dự báo nhu cầu hệ thống điện quốc gia: - Đối tượng áp dụng: o Đơn vị vận hành hệ thống điện thị trường điện (Trung tâm Điều độ hệ thống điện quốc gia) o Đơn vị bán buôn điện (Công ty Mua bán điện) o Đơn vị phân phối điện (các Tổng công ty điện lực) o Đơn vị truyền tải điện (Tổng công ty Truyền tải điện Quốc gia) o Khách hàng sử dụng điện nhận điện trực tiếp từ hệ thống điện truyền tải - Nguyên tắc chung: o Dự báo nhu cầu phụ tải điện hệ thống điện quốc gia dự báo nhu cầu cho toàn phụ tải điện cung cấp điện từ hệ thống điện truyền tải, Phạm Tiến Đạt – Nguyễn Hoàng Huy Khóa luận tốt nghiệp c Kết tính tốn đồ thị: Số liệu tỷ lệ tổn thất giai đoạn 1990 – 2005: Năm LOSS Năm LOSS 1990 28.73 1998 18.27 1991 28.51 1999 17.04 1992 28.58 2000 15.68 1993 26.47 2001 15.57 1994 24.45 2002 15.55 1995 23.55 2003 15.67 1996 21.07 2004 14.56 1997 20.01 2005 14.30 Dự báo tỷ lệ tổn thất giai đoạn 2006 – 2020: Năm 2006 2007 2008 2009 2010 MATLAB 14.97 14.713 14.47 14.24 14.03 Simple E 14.04 2011 Năm MATLAB 13.83 13.78 2012 13.64 13.52 2013 13.47 13.26 2014 13.30 13 2015 13.14 Simple E 12.9 12.79 12.72 12.55 12.48 Năm 2016 2017 2018 2019 2020 MATLAB 12.99 12.84 12.71 12.57 12.45 Simple E 12.35 12.31 12.27 12.22 12.4 Đồ thị dự báo: Hình 3.22: Biểu đồ dự báo tỷ lệ tổn thất Matlab Simple E giai đoạn 2006 – 2020 65 Phạm Tiến Đạt – Nguyễn Hoàng Huy Khóa luận tốt nghiệp 3.11 Dự báo điện sản xuất (ELELDB) (GWh): a Hàm dự báo: ELELDB= FNELDB / (1- LOSSDB/100) b Kết tính tốn đồ thị: Số liệu tỷ lệ tổn thất giai đoạn 1990 – 2005: Năm 1990 1991 1992 1993 1994 1995 1996 1997 19131 ELEL 8679.3 9209 9704.3 10661.4 12288.3 14647.71 16944.66 58 Năm 1998 1999 2000 2001 2002 2003 2004 2005 ELEL 21688.6 23558.5 26561.27 30607.7 35795.75 40924.89 46755.62 53567 Dự báo điện sản xuất giai đoạn 2006 – 2020: Năm 2006 2007 2008 2009 2010 MATLAB 61332.83 69781.54 79359.82 90226.64 102562.49 Simple E 60667.7 69023.3 78471.8 89163.4 101267 Năm 2011 2012 2013 2014 2015 MATLAB 116572.64 132490.82 150583.40 171154.06 194549.25 Simple E 115187 130970.7 148955.8 169192.3 192379.2 Năm 2016 2017 2018 2019 2020 MATLAB 221164.36 251450.83 285924.27 325173.85 369873.05 Simple E 218705.4 248712.3 282868.6 312723.8 365891.8 Đồ thị dự báo: Hình 3.23: Biểu đồ dự báo điện sản xuất Matlab Simple E giai đoạn 2006 – 2020 66 Phạm Tiến Đạt – Nguyễn Hồng Huy Khóa luận tốt nghiệp 3.12 So sánh kết quả: Chọn số liệu dự báo cho tổng nhu cầu điện số liệu dự báo điện sản xuất, để so sánh Để kiểm chứng kết dự báo nhu cầu điện Việt Nam giai đoạn 2006 – 2020 phần mềm MATLAB, thực phép so sanh với kết dự báo nhu cầu điện Việt Nam từ “Luận văn thạc sỹ kỹ thuật Trần Thị Ngoạt” Số liệu để so sánh với kết dự báo MATLAB trích từ đề tài “DỰ BÁO NHU CẦU ĐIỆN NĂNG CHO VIỆT NAM GIAI ĐOẠN 2006 – 2020 TRÊN PHẦN MỀM SIMPLE E” thuộc “Luận văn thạc sỹ kỹ thuật Trần Thị Ngoạt”, số liệu dự báo từ luận văn so sánh kiểm chứng từ số liệu Viện lượng Dưới bảng số liệu đồ thị dự báo: Dự báo nhu cầu điện tổng giai đoạn 2006 – 2020: Năm 2006 2007 2008 2009 2010 MATLAB 52149.14 59514.07 67875.07 77373.33 88169.32 Simple E 52149.969 59511.871 67861.44 77340.3 88102.63 Năm 2011 2012 2013 2014 2015 MATLAB 100445.71 114410.71 130303.12 148389.87 168985.03 Simple E 100327.99 114219.5 130008.6 147958.7 168370.3 Năm 2016 2017 2018 2019 2020 MATLAB 192440.18 219155.86 249597.42 284288.49 323831.94 Simple E 191585.906 217996.351 248047.506 282248.308 321179.83 67 Phạm Tiến Đạt – Nguyễn Hồng Huy Khóa luận tốt nghiệp Đồ thị dự báo: Biểu đồ dự báo tổng nhu cầu điện Matlab Simple E giai đoạn 2006 – 2020 Dự báo điện sản xuất giai đoạn 2006 – 2020: Năm 2006 2007 2008 2009 2010 MATLAB 61332.83 69781.54 79359.82 90226.64 102562.49 Simple E 60667.7 69023.3 78471.8 89163.4 101267 Năm 2011 2012 2013 2014 2015 MATLAB 116572.64 132490.82 150583.40 171154.06 194549.25 Simple E 115187 130970.7 148955.8 169192.3 192379.2 Năm 2016 2017 2018 2019 2020 MATLAB 221164.36 251450.83 285924.27 325173.85 369873.05 Simple E 218705.4 248712.3 282868.6 312723.8 365891.8 68 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy Đồ thị dự báo: Biểu đồ dự báo điện sản xuất Matlab Simple E giai đoạn 2006 – 2020 Nhận xét so sánh kết quả: Lấy hai số liệu quan trọng kết dự báo cho tổng nhu cầu điện kết dự báo cho điện sản xuất giai đoạn 2006 – 2020, ta thấy hai có sai số chênh lệch không lớn 10 năm đầu, giai đoạn 2016 – 2020 có chênh lệch nhiều Do việc dự báo cho giai đoạn dài nên không khỏi tránh khỏi sai số kết mang ý nghĩa hoạch định trước cho tương lai xa Nhìn chung nhu cầu điện điện sản xuất có xu hướng tăng, trung bình hàng năm tăng khoảng 13.8 % giai đoạn 2005 – 2020 Số liệu dự báo công cụ MATLAB lớn số liệu dự báo công cụ Simple E, chênh lệch không đáng kể việc áp dụng thuật toán hai phần mềm khác nhau, tác động từ khứ ảnh hưởng đến kết dự báo 69 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hồng Huy 3.13 Đánh giá mơ hình dự báo Qua đồ thị, dể dàng nhận thấy có sai số dự báo MATLAB Simple E, sai số khơng đáng kể Sai số đến từ khác biệt thuật toán phần mềm đến từ biến động lớn số liệu khứ Ta sử dụng MATLAB để báo điện năng, việc dự báo xác tạo chủ động cho dự án, chiến lược thực tương lai Qua đó, giúp lựa chọn phương án tối ưu công tác xây dựng, quản lý, vận hành, điều tiết hệ thống lưới điện 70 Phạm Tiến Đạt – Nguyễn Hồng Huy Khóa luận tốt nghiệp CHƯƠNG 4: KẾT LUẬN Trong công đại hóa đất nước bước hội nhập tồn cầu Các ngành kinh tế trọng điểm nói chung ngành điện nói riêng ngành cơng nghiệp mũi nhọn, ảnh hưởng trực tiếp đến tất ngành kinh tế Với phát triển vũ bảo, ngành điện vừa có vai trị cung cấp lượng phục vụ đời sống xã hội sinh hoạt người góp phần thúc đẩy phát triển kinh tế Giai đoạn nước ta thúc đẩy nhanh trình cơng nghiệp hố đại hố để hội nhập với kinh tế khu vực giới nhu cầu điện ngày tăng nhanh đòi hỏi cần thiết dự báo đặt biệt dự báo xác Điều quan trọng việc lựa chọn phương pháp dự báo công cụ dự báo có phù hợp với thực tế hay khơng Qua phân tích tình hình nhu cầu điện Việt Nam giai đoạn 1990 – 2006 tất ngành, nhận thấy rõ tất ngành có xu tiêu thụ điện tăng nhanh theo năm, đặt biệt nhu cầu điện công nghiệp nhu cầu điện dân dụng, kèm với gia tăng tốc độ tăng GDP nơng nghiệp công nghiệp Tiến hành dự báo nhu cầu điện cho ngành tiến đến dự báo điện sản xuất, sau so sánh kết dự báo từ MATLAB dự báo từ Simple E Qua bước so sánh, mức độ chênh lệch không lớn 71 Phạm Tiến Đạt – Nguyễn Hồng Huy Khóa luận tốt nghiệp TÀI LIỆU THAM KHẢO [1] Trần Thị Ngoạt, Dự báo nhu cầu điện cho Việt Nam giai đoạn 2006 – 2020 phần mềm Simple E, Luận văn thạc sĩ kỹ thuật, Trường đại học Nông nghiệp I, Năm 2007 [2] Dr Kathleen Ann Cullen, Forecasting Electricity Demand using Rgresssion and Monte Carlo Simulation Under Conditions of Insuffcient Data, Luận văn thạc sỹ, Đại học West Virginia, Năm 1999 [3] Quy trình dự báo nhu cầu phụ tải điện hệ thống điện quốc gia – Ban hành kèm theo Quyết định số 07/QĐ-ĐTĐL ngày 14 tháng 03 năm 2013 Cục trưởng Cục Điều tiết điện lực [4] T.S Nguyễn Quang Dong, Bài giảng kinh tế lượng, NXB Thống Kê, Năm 2003 72 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy PHỤ LỤC Code MATLAB clear all; clc; format short; % Du lieu dau vao: X=[1990;1991;1992;1993;1994;1995;1996;1997;1998;1999;2000; 2001;2002;2003;2004;2005]; % Dan so POP=[66.02;67.24;68.45;69.64;70.82;72;73.16;74.31;75.46;76.6; 77.64;78.69;79.73;80.9;82.07;83.21]; % Nguoi/Ho HH=[4.97;4.94;4.91;4.87;4.83;4.376;4.74;4.72;4.69;4.67;4.64; 4.62;4.59;4.56;4.53;4.4955]; % Luong khach hang CUST=[7890000;8975000;9123000;9256000;9546000;9763000;9850000; 10486000;10490000;11220000;13236090;13250000;13450000; 13590000;15645275;16367040]; % GDP tong GDP=[131968;139634;151782;164043;178534;195567;213833;231264; 244596;256272;273666;292535;313247;335989;361877;389018]; % Ti NN SHAG=[31.83;30.74;30.22;28.88;27.43;26.24;25.06;24.17;23.66; 23.76;23.28;22.43;21.82;21;20.16;19.43]; % Ti CN SHIN=[25.17;25.63;26.59;27.71;28.87;29.94;31.34;32.64;33.43; 34.36;35.41;36.57;37.39;38.47;39.39;40.56]; % Ti TMDV SHCM=[43;43.64;43.19;43.41;43.7;43.82;43.6;43.2;42.91;41.88; 41.3;41;40.79;40.53;40.45;40.01]; % GDP NN GDPAG=[42003;42917;45869;47373;48968;51319;53577;55895;57866; 60896;63717;65618;68352;70574;72939;75592]; % GDP CN GDPIN=[33221;35783;40359;45454;51540;58550;67016;75474;81764; 88047;96913;106986;117126;129247;142559;157778]; % GDP TMDV GDPCM=[56744;60934;65554;71216;78026;85698;93240;99895;104966; 107330;113036;119931;127770;136168;146379;155648]; % Dien dan dung REEL=[2778;2957.5;3155.4;3708.7;4576.5;5369.4;6426.7;7516.5; 9103;9918.4;10985.6;12651.1;14333.2;16339.7;18914.5; 21839.7]; % Nhu cau Dien cong nghiep INEL=[2875.6;3080.2;3196.9;3477.4;3944.1;4619.4;5503.2;6162.7; 6781.1;7568.4;9088.4;10503.2;12681.2;14613.8;16916.6; 19177.4]; % Nhu cau dien thuong mai 73 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy CMEL=[128.9;132.9;141.8;164.4;197;547.5;660.5;751.1;861.5;933; 1083.7;1251.3;1373;1618.7;1873.8;2096.2]; % Nhu cau dien nong nghiep AGEL=[189.8;193.4;202;225.4;263;241.3;283.4;315.4;338.1;403.8; 428.3;465.2;505.6;586.9;679.4;766.7]; % Nhu cau dien khac OTEL=[213.1;219.8;234.7;263.4;302.9;420.5;500.9;556.9;641.7; 720.5;811.6;972.3;1335.2;1350.9;1563.8;1774]; % Ti le dien hoa ELEC=[0.56;0.56;0.56;0.57;0.57;0.58;0.62;0.67;0.71;0.75;0.8; 0.81;0.82;0.84;0.86;0.8843]; % Ti le ton that LOSS= [28.73;28.51;28.58;26.47;24.45;23.55;21.07;20.01;18.27; 17.04;15.68;15.57;15.55;15.67;14.56;14.3]; n=length(X); x=[2006;2007;2008;2009;2010;2011;2012;2013;2014;2015;2016; 2017;2018;2019;2020]; n1=length(x); GDPGR=7.5; e=2.7183; % Du bao dan so mdl0=fitlm(X,POP); % a0=-2196.4 b0=1.137; % constant adjusted a0= POP(n)-b0*X(n); for i=1:n1 POPDB(i)=a0+b0*x(i); end POPDB; % Du bao gia tang dan so POPGRDB(1)=((POPDB(1)/POP(n))-1)*100; for i=2:n1 POPGRDB(i)=((POPDB(i)/POPDB(i-1))-1)*100; end POPGRDB; % Du bao GDP tong nGDP= length(GDP); GDPDB=zeros(1,n1); GDPDB(1)=GDP(n)*(1+(GDPGR/100)); for i=2:n1 GDPDB(i)=GDPDB(i-1)*(1+(GDPGR/100)); end GDPDB; % Du bao ti ngành nong nghiep, cong nghiep thuong mai LSHAG=log(SHAG); mdl1=fitlm(X,LSHAG); B1= -0.031769; LSHIN=log(SHIN); mdl2=fitlm(X,LSHIN); 74 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy B2= 0.032789; LSHCM=log(SHCM); mdl3=fitlm(X,LSHCM); B3=-0.0061727; % constant adjusted b1=e^B1; b2=e^B2; b3=e^B3; a1=SHAG(n)/(b1^2005); a2=SHIN(n)/(b2^2005); a3=SHCM(n)/(b3^2005); for i=1:1:n1 SHAGDB(i)=a1*(b1^(x(i))); SHINDB(i)=a2*(b2^(x(i))); SHCMDB(i)=a3*(b3^(x(i))); end SHAGDB; SHINDB; SHCMDB; SHTL=SHAGDB+SHINDB+SHCMDB; % Du bao tang truong kinh te nganh for i=1:n1 GDPAGDB(i)= GDPDB(i)*(SHAGDB(i)/SHTL(i)); GDPINDB(i)= GDPDB(i)*(SHINDB(i)/SHTL(i)); GDPCMDB(i)= GDPDB(i)*(SHCMDB(i)/SHTL(i)); end GDPAGDB; GDPINDB; GDPCMDB; % Du bao toc tang truong nganh GRAGDB(1)=((GDPAGDB(1)/GDPAG(n))-1)*100; GRINDB(1)=((GDPINDB(1)/GDPIN(n))-1)*100; GRCMDB(1)=((GDPCMDB(1)/GDPCM(n))-1)*100; for i=2:n1 GRAGDB(i)=((GDPAGDB(i)/GDPAGDB(i-1))-1)*100; GRINDB(i)=((GDPINDB(i)/GDPINDB(i-1))-1)*100; GRCMDB(i)=((GDPCMDB(i)/GDPCMDB(i-1))-1)*100; end GRAGDB; GRINDB; GRCMDB; % Du bao so nguoi/ho LHH=log(HH); mdl4=fitlm(X,LHH); b4=-0.0059517; a4=e^13.435; for i=1:n1 HHDB(i)=a4*(e^(b4*x(i))); end HHDB; % Ti le dien hoa 75 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy for i=1:n abc(i)=GDP(i)/POP(i); LELEC(i)=log(ELEC(i)/(1-ELEC(i))); end LELEC; mdl5=fitlm(abc,LELEC); % a5=-1.4654 b5=0.00074845; % constant adjusted a5=log(ELEC(n)/(1-ELEC(n)))-b5*(GDP(n)/POP(n)); for i=1:n1 ELECDB(i)=e^(a5+b5*(GDPDB(i)/POPDB(i)))/(1+e^(a5+b5*(GDPDB(i)/ POPDB(i)))); end ELECDB; % Ti le ton that t=1:1:n; T=17:1:31; LLOSS=log(LOSS); LX=log(t); mdl6= fitlm(LX,LLOSS); a6=e^3.5778; b6=-0.3076; for i=1:n1 LOSSDB(i)=a6*((T(i)^b6)); end LOSSDB; % Du bao luong khach hang for i=1:n1 DS(i)=((POPDB(i)*10^6)/HHDB(i))*ELECDB(i); end CUSTDB=DS'; % Du bao nhu cau Dien cong nghiep LINEL=log(INEL); LGDPIN=log(GDPIN); mdl7=fitlm(LGDPIN,LINEL); % a7=-5.2882 b7=1.2576; % constant adjust a7=INEL(n)/(GDPIN(n))^b7; for i=1:n1 INELDB(i)=a7*(GDPINDB(i)^b7); end INELDB; % Nhu cau Dien dan dung LREEL=log(REEL); GC=zeros(n,1); for i=1:n GC(i)=GDP(i)/CUST(i); end 76 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hồng Huy LGC=log(GC); X2=zeros(n1,1); for i=1:n1 X2(i)=GDPDB(i)/CUSTDB(i); end LX2=log(X2); LCUST=log(CUST); X1=[LGC,LCUST]; mdl8=fitlm(X1,LREEL); % a8= -15.372 b8=2.0118; c8=1.9799; % constant adjust a8=(REEL(16)/(((GDP(16)/CUST(16))^b8)*(CUST(16)^c8))); for i=1:n1 REELDB(i)=a8*(X2(i)^b8)*(CUSTDB(i)^c8); end REELDB; % Nhu cau dien thuong mai LAG1GDPCM=[56744;60934;65554;71216;78026;85698;93240;99895; 104966;107330;113036;119931;127770;136168;146379]; GDPCMY=[60934;65554;71216;78026;85698;93240;99895;104966; 107330;113036;119931;127770;136168;146379;155648]; LAG1CMEL=[128.9;132.9;141.8;164.4;197;547.5;660.5;751.1;861.5; 933;1083.7;1251.3;1373;1618.7;1873.8]; CMELY=[132.9;141.8;164.4;197;547.5;660.5;751.1;861.5;933; 1083.7;1251.3;1373;1618.7;1873.8;2096.2]; X3=[GDPCMY,LAG1GDPCM,LAG1CMEL]; mdl9=fitlm(X3,CMELY); % a9=-708.71; b9=0.025259; c9=-0.01466; d9=0.53511; % constant adjust a9=CMELY(15)-b9*GDPCMY(15)-c9*LAG1GDPCM(15)-d9*LAG1CMEL(15); CMELDB(1)=a9+b9*(GDPCMDB(1))+c9*GDPCM(16)+d9*CMEL(16); for i=2:n1 CMELDB(i)=a9+b9*(GDPCMDB(i))+c9*GDPCMDB(i-1)+d9*CMELDB(i-1); end CMELDB; % Nhu cau dien nong nghiep LAG1GDPAG=[42003;42917;45869;47373;48968;51319;53577;55895; 57866;60896;63717;65618;68352;70574;72939]; GDPAGY=[42917;45869;47373;48968;51319;53577;55895;57866;60896; 63717;65618;68352;70574;72939;75592]; LAG1AGEL=[189.8;193.4;202;225.4;263;241.3;283.4;315.4;338.1; 403.8;428.3;465.2;505.6;586.9;679.4]; AGELY=[193.4;202;225.4;263;241.3;283.4;315.4;338.1;403.8; 428.3;465.2;505.6;586.9;679.4;766.7]; X4=[GDPAGY,LAG1GDPAG,LAG1AGEL]; 77 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy mdl10=fitlm(X4,AGELY); % a10=-80.736; b10=-0.0036233; c10=0.0056315; d10=0.906; % constant adjust a10=AGELY(15)-b10*GDPAGY(15)-c10*LAG1GDPAG(15)d10*LAG1AGEL(15); AGELDB(1)=a10+b10*(GDPAGDB(1))+c10*GDPAG(16)+d10*AGEL(16); for i=2:n1 AGELDB(i)=a10+b10*GDPAGDB(i)+c10*GDPAGDB(i1)+d10*AGELDB(i-1); end AGELDB; % Nhu cau dien khac LAGOTEL=[213.1;219.8;234.7;263.4;302.9;420.5;500.9;556.9; 641.7;720.5;811.6;972.3;1335.2;1350.9;1563.8]; OTELY=[219.8;234.7;263.4;302.9;420.5;500.9;556.9;641.7;720.5; 811.6;972.3;1335.2;1350.9;1563.8;1774]; GDPY=[139634;151782;164043;178534;195567;213833;231264;244596; 256272;273666;292535;313247;335989;361877;389018]; LAGPOP=[66.02;67.24;68.45;69.64;70.82;72;73.16;74.31;75.46; 76.6;77.64;78.69;79.73;80.9;82.07]; X5=zeros(n1,1); for i=1:n1 X5(i)=GDPY(i)/LAGPOP(i); end X6=[X5,LAGOTEL]; mdl11=fitlm(X6,OTELY); % a11=-417.89 b11=0.20651; c11=0.75997; % constant adjust a11=OTELY(15)-b11*(GDPY(15)/LAGPOP(15))-c11*(LAGOTEL(15)); OTELDB(1)=a11+b11*(GDPDB(1)/POP(16))+c11*OTELY(15); for i=2:n1 OTELDB(i)=a11+b11*(GDPDB(i)/POPDB(i-1))+c11*OTELDB(i-1); end OTELDB; % Nhu cau dien nang tieu thu toan quoc for i=1:n1 FNELDB(i)= INELDB(i)+REELDB(i)+CMELDB(i)+AGELDB(i)+OTELDB(i); end FNELDB; %luong dien nang san xuat for i=1:n1 ELELDB(i)=FNELDB(i)/(1-(LOSSDB(i)/100)); end ELELDB; 78 Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hoàng Huy LỜI CẢM ƠN Sau năm tháng học tập, để hồn thành khóa luận tốt nghiệp trước hết em xin gửi đến quý thầy, cô giáo khoa Công Nghệ Điện trường Đại Học Công Nghiệp Thành Phố Hồ Chí Minh lời cảm ơn chân thành Đặc biệt, em xin gửi đến thầy TS Trần Thanh Ngọc, người tận tình hướng dẫn, giúp đỡ em hồn thành khóa luận tốt nghiệp lời cảm ơn sâu sắc Đồng thời em xin cảm ơn khoa Công Nghệ Điện nhà trường tạo điều kiện cho em làm khóa luận, bên cạnh cho chúng em hội nghiên cứu tìm hiểu kỹ phần mềm Matlab vận dụng lý thuyết trao dồi năm học tập trường Vì kiến thức thân cịn hạn chế q trình thực hiện, hịa thiện báo cáo em không tránh khỏi sai sót, kính mong nhận ý kiến đóng góp từ q thầy Và cuối em xin kính chúc quý thầy, cô dồi sức khỏe thành công nghiệp cao quý 79 ... o Dự báo nhu cầu phụ tải điện hệ thống điện quốc gia dự báo nhu cầu cho toàn phụ tải điện cung cấp điện từ hệ thống điện truyền tải, Khóa luận tốt nghiệp Phạm Tiến Đạt – Nguyễn Hồng Huy trừ phụ. .. trừ phụ tải điện có nguồn cung cấp điện riêng không nhận điện từ hệ thống điện quốc gia o Dự báo nhu cầu phụ tải điện hệ thống điện quốc gia bao gồm dự báo nhu cầu phụ tải điện năm, tháng, tuần,... 3.10 Biểu đồ dự báo GDP nông nghiệp Matlab Simple E giai đoạn 2006 – 2020 39 3.11 Biểu đồ dự báo GDP công nghiệp Matlab Simple E giai đoạn 2006 – 2020 40 3.12 Biểu đồ dự báo GDP DVTM Matlab Simple