Ứng dụng phần mềm matlabsimulink mô phỏng hệ phát điện gió sử dụng máy phát điện đồng bộ nam châm vĩnh cửu

HỌC VIỆN NÔNG NGHIỆP VIỆT NAM KHOA CƠ ĐIỆN  KHÓA LUẬN TỐT NGHIỆP ỨNG DỤNG PHẦN MỀM MATLAB/SIMULINK MÔ PHỎNG HỆ PHÁT ĐIỆN GIÓ SỬ DỤNG MÁY PHÁT ĐIỆN ĐỒNG BỘ NAM CHÂM VĨNH CỬU Hà Nội - 2021 HỌC VIỆN NÔNG NGHIỆP VIỆT NAM KHOA CƠ ĐIỆN  KHÓA LUẬN TỐT NGHIỆP ỨNG DỤNG PHẦN MỀM MATLAB/SIMULINK MƠ PHỎNG HỆ PHÁT ĐIỆN GIĨ SỬ DỤNG MÁY PHÁT ĐIỆN ĐỒNG BỘ NAM CHÂM VĨNH CỬU Sinh viên thực Mã sinh viên Lớp Bộ môn Chuyên ngành Người hướng dẫn : PHAN THÀNH LỘC : 597873 : K59 – HTĐB : HỆ THỐNG ĐIỆN : HỆ THỐNG ĐIỆN : ThS ĐÀO XUÂN TIẾN Hà Nội - 2021 LỜI CÁM ƠN Năng lượng tái tạo, lượng xanh nghiên cứu ứng dụng rộng rãi nước Những tác động người vào thiên nhiên làm thiên nhiên biến đổi Thực mục tiêu ngành lượng tìm cho bước tiến phát triển Do em khảo nghiệm nghiên cứu đề tài sử dụng lượng gió Trong q trình nghiên cứu đề tài, em thu nhận học tập nhiều kinh nghiệm quý báu Em xin trân trọng gửi lời cảm ơn chân thành tới thầy Ths Đào Xuân Tiến, người tận tình hướng dẫn, giúp đỡ cho em ý kiến, đóng góp quý báu để em hoàn thành đồ án Em xin gửi lời cảm ơn tới quý Thầy, Cô khoa Cơ Điện trang bị cho em khối lượng kiến thức bổ ích q báu q trình học tập nghiên cứu Những kiến thức tạo tảng vững giúp em hoàn thành tốt đồ án Tuy em cố gắng, chắn đồ án cịn nhiều thiếu sót cần điều chỉnh bổ sung Rất mong nhận ý kiến đóng góp từ q Thầy Cơ để tri thức em hoàn thiện trường Hà Nội, ngày tháng năm 202 Sinh viên thực i MỤC LỤC LỜI CÁM ƠN MỤC LỤC ii DANH MỤC BẢNG iv DANH MỤC HÌNH v CÁC TỪ VIẾT TẮT TRONG ĐỒ ÁN viii MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN VỀ NĂNG LƯỢNG GIÓ 1.1 Sự hình thành gió 1.2 Sử dụng lượng gió để tạo điện 1.3 Lợi ích lượng gió mang lại 1.4 Tình hình sử dụng lượng điện gió số nước giới 11 1.5 Tình hình sử dụng điện gió Việt Nam 13 CHƯƠNG 2: TỔNG QUAN VỀ HỆ THỐNG PHÁT ĐIỆN GIÓ 15 2.1 Cấu trúc hệ thống phát điện gió 15 2.2 Khí động học gió 16 2.3 Các loại cấu trúc turbine gió 17 2.4 Các phận turbine gió 18 2.4.1 Cánh quạt 18 2.4.2 Bộ truyền động 23 2.4.3 Máy phát 26 CHƯƠNG 3: MÔ PHỎNG HỆ PHÁT ĐIỆN GIÓ SỬ DỤNG MÁY PHÁT ĐIỆN ĐỒNG BỘ NAM CHÂM VĨNH CỬU BẰNG PHẦN MỀM MATLAB/SIMULINK 32 3.1 Mơ hình hóa phần tử hệ phát điện gió 32 3.2 Giới thiệu Matlab/Simulink 35 3.2.1 Khái quát Matlab 35 ii CHƯƠNG 4: THIẾT LẬP MƠ PHỎNG HỆ PHÁT ĐIỆN GIĨ VỚI MÁY PHÁT ĐIỆN ĐỒNG BỘ NAM CHÂM VĨNH CỬU 39 4.1 Mô hệ thống Matlab/Simulink 39 4.2 Thiết lập mô hệ thống phát điện gió 42 4.3 Mơ hình chi tiết 44 4.3.1 Khối Tuabin gió 44 4.3.2 Khối máy phát điện đồng nam châm vĩnh cửu 47 4.3.3 Khối chỉnh lưu 48 4.3.4 Khối P&O MPPT 48 4.3.5 Mơ hình mơ hệ phát điện gió sử dụng MPPT 50 4.4 Kết mô 50 4.4.1 Hệ mô không sử dụng MPPT 51 4.4.2 Hệ mô sử dụng MPPT 56 4.4.3 Hệ mô thay đổi tải sử dụng MPPT 61 CHƯƠNG 5: THỰC NGHIỆM XÂY DỰNG ĐẶC TÍNH LÀM VIỆC CỦA TUABIN PHÁT ĐIỆN GIÓ ĐỒNG BỘ NAM CHÂM VĨNH CỬU 64 5.1 THIẾT LẬP MƠ HÌNH THỰC NGHIỆM 64 5.1.1 Mơ hình thực nghiệm 64 5.1.2 Phương pháp thu thập số liệu 66 5.1.3 Kết nhận 66 KẾT LUẬN VÀ KIẾN NGHỊ 68 TÀI LIỆU THAM KHẢO 70 iii DANH MỤC BẢNG Bảng 1.1- Bảng phân bố lượng điện gió số nước giới năm 2005 [10] 12 Bảng 4.1 Các Block thường dùng để mô đồ án 39 Bảng 5.1 Các thiết bị sử dụng mơ hình 64 Bảng 5.2 Bảng số liệu thu thập mơi trường mơ hình hoạt động 67 iv DANH MỤC HÌNH Hình 1.1 Sự hình thành gió khí Hình 1.2 Turbine gió 12 KW Hình 1.3 Turbine ngồi nước Mỹ Hình 1.4 Cánh đồng phát điện gió Bình Thuận Hình 2.1 Cấu trúc hệ thống phát điện gió 15 Hình 2.2 turbine trục ngang (trái), turbine trục dọc (phải) 17 Hình 2.3 Turbine ngược chiều gió (trái) turbine thuận chiều gió (phải) 18 Hình 2.4 Lực nâng lực kéo tác động vào cánh quạt 19 Hình 2.5 Mơ hình BEM- blade element method 19 Hình 2.6 Actuator model 20 Hình 2.7 Sự khác biệt số lượng cánh quạt 21 Hình 2.8 Cp loại cánh quạt khác turbine gió 23 Hình 2.9 Cấu tạo hộp số (Gearbox) 25 Hình 2.10 Sơ đồ khối minh họa cho hệ thống truyền động trực tiếp; cánh quạt nối trực tiếp vào rotor máy phát; có biến đối dùng linh kiện bán dẫn để tương thích tần số phát tải 25 Hình 2.11 – Hệ thống điện gió sử dụng động lồng sóc tốc độ khơng đổi kết hợp với khởi động mềm 27 Hình 2.12 - Máy phát lồng sóc điều chỉnh tốc độ, trang bị biến đổi điện áp AC/DC/AC 27 Hình 2.13 - Máy phát dây quấn đồng ứng dụng turbine gió 28 Hình 2.14 - Turbine gió sử dụng máy phát dây quấn khơng đồng 29 Hình 2.15 - Sơ đồ nguyên lý máy phát không đồng nguồn kép 29 Hình 2.16 - Turbine gió sử dụng máy phát kích từ vĩnh cửu 30 Hình 2.17 - Turbine gió dùng máy phát DC 31 Hình 3.1 Mơ hình mạch điện tương hệ tọa độ dq PMSG 34 v Hình 3.2 Cơng ty MatWorks 36 Hình 3.3 Phần mềm Matlab-Simulink 36 Hình 4.1 Khối Tuabin gió mơ Matlab 44 Hình 4.2 Mơ hình Tuabin thích họa 44 Hình 4.3 Tốc độ Tuabin(p.u) cung cấp với tốc độ gió tối đa 12m/s 46 Hình 4.4 Các thơng số khối Turbine 46 Hình 4.5 Cấu hình máy phát điện đồng nam châm vĩnh cửu 47 Hình 4.6 Tham số máy phát điện 47 Hình 4.7 Khối chỉnh lưu, thông số khối chỉnh lưu 48 Hình 4.8 Khối mơ MPPT 48 Hình 4.9 Khối IO đặc trưng cho khối IO1 IO2 hình 4.8 49 Hình 4.11 Bộ chuyển đổi Boost Converter 49 Hình 4.12 Mơ hình mơ hệ thống Tuabin phát điện gió 50 Hình 4.13 Hệ mô không sử dụng MPPT với tốc độ gió thay đổi 6m/s, 8m/s, 10m/s, 12m/s Hệ mơ với R=2000(Ω) 51 Hình 4.14 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 51 Hình 4.15 Đồ thị đặc tính điện áp(V), dòng điện(A) đầu máy phát 51 Hình 4.16 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 52 Hình 4.17 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 52 Hình 4.18 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 53 Hình 4.19 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 53 Hình 4.20 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 54 Hình 4.21 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 54 Hình 4.22 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 54 Hình 4.23 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 55 Hình 4.24 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 55 Hình 4.25 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 56 vi Hình 4.26 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 56 Hình 4.27 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 57 Hình 4.28 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 57 Hình 4.29 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 58 Hình 4.30 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 58 Hình 4.31 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 58 Hình 4.32 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 59 Hình 4.33 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 59 Hình 4.34 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 60 Hình 4.35 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 60 Hình 4.36 Đồ thị đặc tính điện áp(V), dòng điện(A) đầu máy phát 61 Hình 4.37 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 61 Hình 4.38 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 62 Hình 4.39 Đồ thị đặc tính điện áp(V), dòng điện(A) đầu máy phát 62 Hình 4.40 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 62 Hình 5.1 Sơ đồ thực nghiệm tuabin phát điện gió 64 Hình 5.2 Đồng hồ vạn ZT – S1(bên trái), Đồng hồ vạn XL830(ở giữa), Điện trở thực nghiệm(bên phải) 65 Hình 5.3 Mơ hình thực nghiệm 66 Hình 5.4 Tuabin lắp thực nghiệm 67 vii CÁC TỪ VIẾT TẮT TRONG ĐÒ ÁN P&O Perturb & observe MPP Maximum Power Point MPPT Maximum Power Point Tracker PMG Permanent Magnet Generator PMSG Permanet Magnet Synchronous Generator WRSG Wound Rotor Synchronous Generator WRIG Wound Rotor Induction Generator SCIG Squirrel Cage Induction Generator DFIG Doubly Fed Induction Generator DCG DC Generator VRG Variable Reluctance Generator INCond Incremental Conductance FLC Fuzzy logic controller IEA International Energy Agency HAWT Horizontal Axis Wind Turbine VAWT Vertical Axis Wind Turbine BEM blade element method viii Hình 4.27 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát Hình 4.28 Đồ thị đặc tính điện áp(V) tải qua khối IGBT Nhận xét: - Trong trường hợp nhiều khác biệt so với khơng có MPPT tốc độ gió nhỏ, tốc độ rotor mức 128.5 rad/s (1227 vòng/phút) sau 0.115(s) momen động đạt giá trị 0.1N.m - Trong 0.144(s) để điện áp đầu myas phát xác lập giá trị 59.5(V) 0.09(A) - Đối với điện áp phía tải khơng thay đổi đáng kể mức 177(V) tải qua MBA 8.7(V) tải không qua MBA 57 4.4.2.2 Mô hệ tốc độ gió 8m/s Hình 4.29 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát Hình 4.30 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát Hình 4.31 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 58 Nhận xét: - Với tốc độ gió 8m/s ta thu giá trị q trình mơ khơng khác nhiều so với khơng có MPPT - Tốc độ rotor thu 174.5 rad/s (1666 vòng/phút) sau 0.105(s) để đạt tới chế độ xác lập momen động 0.15N.m - Điện áp đầu máy phát 124(V) dòng điện thu giá trị 0.23(A) - Về phần nối tải qua MBA tuyến tính giá trị điện áp vào khoảng 177(V) tải bình thường 8.7(V) 4.4.2.3 Mơ hệ tốc độ gió 10m/s Hình 4.32 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát Hình 4.33 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát 59 Hình 4.34 Đồ thị đặc tính điện áp(V) tải qua khối IGBT Nhận xét: - Trong lần mô sử dụng MPPT ta thu tốc độ rotor mức 222 rad/s (2120 vòng/phút) sau 0.09(s) momen 0.2N.m - Đối với điện áp đâu máy phát 0.1(s) để tiến tới giá trị xác lập 96(V), 0.2(A) - Với phía tải ta thu giá trị điện áp qua MBA 180(V) tải không qua MBA 8.8(V) 4.4.2.4 Mô hệ tốc độ gió 12m/s Hình 4.35 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát 60 Hình 4.36 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát Hình 4.37 Đồ thị đặc tính điện áp(V) tải qua khối IGBT Nhận xét: - Trong lần mô rotor 0.08(s) để đạt 277 rad/s (2645 vòng/phút) Ở đầu máy phát 0.1(s) để máy phát đạt giá trị xác lập 123(V) 0.25(A) - Về phía tải MBA ta thu giá trị 184(V) tải không qua MBA 9.2(V) 4.4.3 Hệ mô thay đổi tải sử dụng MPPT Sau mô thay đổi giá trị tốc độ gió khác nhau, giá trị điện áp, dòng điện tốc độ rotor thay đổi có thay đổi số tác động đầu vào Tuabin Trong trường hợp em mô hệ thống thay đổi 61 thông số tải để so sánh với kết thu phía Lần mơ em thay đổi giá trị tải mức gió 10m/s giá trị R=1000(Ω) Khi thay đổi tải ta thu giá trị đặc tính điện áp, dịng điện, tốc độ rotor điện áp phía tải sau: Hình 4.38 Đồ thị đặc tính tốc độ rotor(rad/s) momen(N.m) máy phát Hình 4.39 Đồ thị đặc tính điện áp(V), dịng điện(A) đầu máy phát Hình 4.40 Đồ thị đặc tính điện áp(V) tải qua khối IGBT 62 Nhận xét: - Với trường hợp thay đổi tải lần rotor 0.092(s) để đạt giá trị xác lập 227.5 rad/s (2645 vịng/phút) 0.15N.m - Phía đầu điện áp máy phát ngưỡng 100(V) sau 0.115(s) Phần điện áp phía tải qua MBA tuyến tính thay đổi tải giá trị điện áp cao tăng lên tới 412(V) tải bình thường mức 20.6(V) Nhận xét chung: Hệ thống mô vận hành tốt, đưa chế độ mô để so sánh thông số qua thời gian Cũng từ kết mô cho thấy việc sử dụng MPPT giúp cho hệ thống có giá trị điện áp, dịng điện tốc độ rotor đạt giá trị xác lập nhanh so với khơng có MPPT thay đổi tốc độ gió Tuy nhiên thời gian độ hệ thống xuất thay đổi tốc độ gió hay phụ tải ngắn Các kết đạt có giao động khác giao động nhỏ, không đáng kể Từ mơ cho thấy hệ thống có sử dụng MPPT hoạt động tối ưu khơng có MPPT Qua giúp em hiểu thêm phần vận hành hệ thống sử dụng trợ giúp 63 CHƯƠNG THỰC NGHIỆM XÂY DỰNG ĐẶC TÍNH LÀM VIỆC CỦA TUABIN PHÁT ĐIỆN GIÓ ĐỒNG BỘ NAM CHÂM VĨNH CỬU 5.1 THIẾT LẬP MƠ HÌNH THỰC NGHIỆM 5.1.1 Mơ hình thực nghiệm Hình 5.1 Sơ đồ thực nghiệm tuabin phát điện gió Ở sơ đồ thực nghiệm sinh viên thực thay giá trị điện trở khác để đo thơng số dịng điện điện áp tương ứng với giá trị điện trở khác Bảng 5.1 Các thiết bị sử dụng mô hình Tên thiết bị Thơng số Số lượng Tuabin phát điện gió 50W – AC 12 pha Đồng hồ vạn – 20V – 10 A Điện trở 10W – 10R Điện trở 5W – 30R Điện trở 5W – 1kR 64 Hình 5.2 Đồng hồ vạn ZT – S1(bên trái), Đồng hồ vạn XL830(ở giữa), Điện trở thực nghiệm(bên phải) - Đồng hồ vạn XL830:  Kiểu hiển thị: Kỹ thuật số  Kích thước: 143mm x 69mm  Sử dụng pin: 9V  Khoảng đo điện áp chiều: 200mV -> 1000V  Khoảng đo điện áp xoay chiều: 200mV -> 750V  Khoảng đo cường độ dòng điện: 20mA ->20A  Khoảng đo điện trở: 200Ω - 2MΩ  Khoảng đo tụ điện: 20nF ->200uF - Đồng hồ vạn ZT – S1:  Kiểu hiển thị: Kỹ thuật số  Kích thước: 118mm x 59mm  Khoảng đo điện áp chiều: 4V -> 600V  Khoảng đo điện áp xoay chiều: 4V -> 600V  Khoảng đo cường độ dòng điện: 1000mA -> 10A 65  Khoảng đo điện trở: 4k -> 40MΩ 5.1.2 Phương pháp thu thập số liệu Xây dựng mơ hình lắp đặt Đo thơng số tác động đến mơ hình thực nghiệm tốc độ gió ngồi mơi trường, thay đổi giá trị điện trở Do việc đo phụ thuộc vào thay đổi tốc độ gió mơ trường xung quanh nên kết đo liên tục theo thời gian ghi kết lại 5.1.3 Kết nhận Hình 5.3 Mơ hình thực nghiệm 66 Hình 5.4 Tuabin lắp thực nghiệm Bảng 5.2 Bảng số liệu thu thập mơi trường mơ hình hoạt động R(Ω) 5W – 1kR 5W – 30R 10W – 10R Điện áp 6.3 7.7 8.24 7.1 7.97 8.8 Dòng điện 1.26 1.09 1.37 1.24 1.15 1.08 8.37 8.42 7.92 8.53 8.22 1.31 1.34 1.28 1.37 1.35 9.35 9.49 9.68 9.13 8.92 1.43 1.47 1.54 1.52 1.49 67 KẾT LUẬN VÀ KIẾN NGHỊ Kết đạt Thông qua đề tài: “ỨNG DỤNG PHẦN MỀM MATLAB/SIMULINK MÔ PHỎNG HỆ PHÁT ĐIỆN GIÓ” Đồ án em thu thập nội dung sau: - Xây dựng mô thành cơng hệ phát điện gió phần mềm Matlab/Simulink - Sử dụng thuật tốn để áp dụng vào mơ hệ phát điện gió giúp so sánh giá trị điện áp, dịng điện,tốc độ rotor mơ - Thông qua phần mềm Matlab/Simulink giúp em biết thêm thơng tin máy phát điện, tuabin gió, hệ thống phát điện gió Hiểu biết thêm thơng tin thuật toán MPPT - Hệ thống phát điện sử dụng MPPT có thời gian xác lập nhanh không sử dụng MPPT Cải thiện hệ thống thêm khối chức vào hệ thống Những mặt hạn chế Mặc dù có cố gắng điều kiện vừa học vừa làm nên đồ án cịn nhiều mặt hạn chế thiếu xót Đa số thông số hay tham số mô dựa phương pháp thử sai nên việc tối ưu hệ thống chưa đảm bảo tính xác thực hệ thống chưa cao Nhiều mơ cịn sử dụng toolbox có sẵn Matlab/Simulink khối Tuabin gió, Máy phát điện, khối chỉnh lưu Các số liệu thu thập hệ mô phần mềm hệ thực nghiệm có phần chủ quản, em chưa đem so sánh tiếp cận tốt với mơ hình nghiên cứu hay hoạt động ngồi mơi trường Tất tham số thay đổi đặt hệ mô em hoàn toàn dựa phương pháp thử để đạt giá trị điện áp đầu tải để giá trị hình Sin Nên giá trị khơng xác cao Mong thầy thơng 68 cảm cho em góp ý sai xót mà hệ thống em gặp phải Em xin chân thành cảm ơn tới giảng viên ThS Đào Xuân Tiến tạo điều kiện giúp đỡ em trình thực đồ án Em xin chân thành cảm ơn thầy cô môn tham gia đanh giá cho em! 69 TÀI LIỆU THAM KHẢO MATLAB Simulink dành cho kỹ sư điều khiển tự động, Nguyễn Phùng Quang (2006), NXB Khoa Học Và Kỹ Thuật Hà Nội Bài giảng Điện tử công suất, Nguyễn Văn Đạt, NXB Đại Học Nông Nghiệp Hà Nội Giáo trình lý thuyết thiết kế ứng dụng lượng gió, McGowan, J G, Manwell, J F NXB Trường Đại Học Thủy Lợi Cơ sở lượng tái tạo, Đặng Đình Thống – Lê Danh Liên, NXB Khoa Học Và Kỹ Thuật (2006) Tìm điểm cơng suất cực đại máy phát điện gió PMSG 200W, Nguyễn Thái Sơn, Tạp Chí Khoa Học Trường Đại Học Cần Thơ (2019) Tìm điểm cơng suất cực đại cho máy phát điện gió – Máy đồng nam châmvĩnh cửu phương pháp P&O, Bùi Hữu Hân, Luận Văn Thạc Sĩ (2017) Wikimedia – Năng Lượng Gió: https://vi.wikipedia.org/wiki/N%C4%83ng_l%C6%B0%E1%BB%A3ng_gi %C3%B3 Trần Thanh Lâm & Nguyễn Linh Chi (2015) Phân tích chi phí - lợi ích khai thác điện gió khu vực Bạc Liêu Nguyễn Trung Dũng Nguyễn Tuấn Anh Ứng dụng lượng gió –xu chung đánh giá góc độ kinh tế môi trường Trường Đại Học Thủy Lợi 10 Lê Hồi Chí (2008) Nghiên cứu tổng quan lượng gió nhà máy điện gió Phương Mai - Việt Nam Trường Đại học Kỹ Thuật Công Nghệ TP Hồ Chí Minh, 2008 11 Hồ Phạm Huy Ánh, Nguyễn Hữu Phúc, Nguyễn Văn Tài, Phạm Đình Trực, Nguyễn Quang Nam, Trần Công Binh Phan Quang Ấn (2013) Kỹ thuật hệ thống lượng tái tạo, Đại Học Bách Khoa TP.HCM 70 Một số tài liệu mạng internet: https://www.technologymag.net/ http://nangluongvietnam.vn/news/vn/trang-chu http://www.dientuvietnam.net/forums/ http://doc.edu.vn/default.aspx Tham khảo số mơ hình Matlab/Simulink số Youtube: Ấn Độ, Isarel, https://www.youtube.com/watch?v=5EVT_GaGZ40&t=971s https://www.youtube.com/watch?v=Tjm-Oa4GKW0&t=338s https://www.youtube.com/watch?v=beXGDBiCJTk&t=671s https://www.youtube.com/watch?v=Y3kXvRZAmDI&t=256s 71