Thiết kế, lắp đặt mô hình điều hòa không khí và ứng dụng phần mềm matlab vào mô phỏng

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH ĐỒ ÁN TỐT NGHIỆP NGÀNH CNKT Ô TÔ THIẾT KẾ, LẮP ĐẶT MƠ HÌNH ĐIỀU HỊA KHƠNG KHÍ VÀ ỨNG DỤNG PHẦN MỀM MATLAB VÀO MÔ PHỎNG GVHD: THS NGUYỄN THÀNH TUYÊN SVTH: NGUYỄN KHẮC HÙNG KHUẤT ĐỨC DUY SKL008262 Tp.Hồ Chí Minh, tháng 8/2021 TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA CƠ KHÍ ĐỘNG LỰC ĐỒ ÁN TỐT NGHIỆP THIẾT KẾ, LẮP ĐẶT MƠ HÌNH ĐIỀU HỊA KHƠNG KHÍ VÀ ỨNG DỤNG PHẦN MỀM MATLAB VÀO MÔ PHỎNG SVTH: NGUYỄN KHẮC HÙNG MSSV: 17145304 SVTH: KHUẤT ĐỨC DUY MSSV: 17145268 GVHD: Th.S NGUYỄN THÀNH TUYÊN Tp Hồ Chí Minh, ngày 25 tháng 08 năm 2021 TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA CƠ KHÍ ĐỘNG LỰC ĐỒ ÁN TỐT NGHIỆP Chuyên ngành: Công nghệ Kỹ thuật ô tô Tên đề tài THIẾT KẾ, LẮP ĐẶT MƠ HÌNH ĐIỀU HỊA KHƠNG KHÍ VÀ ỨNG DỤNG PHẦN MỀM MATLAB VÀO MƠ PHỎNG SVTH: NGUYỄN KHẮC HÙNG MSSV: 17145304 SVTH: KHUẤT ĐỨC DUY MSSV: 17145268 GVHD: Th.S NGUYỄN THÀNH TUYÊN Tp Hồ Chí Minh, ngày 25 tháng năm 2021 TRƯỜNG ĐH SƯ PHẠM KỸTHUẬT CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM TP HỒ CHÍ MINH Độc Lập – Tự Do – Hạnh Phúc KHOA CƠ KHÍ ĐỘNG LỰC ──────── ──────── TP Hồ Chí Minh, ngày 25 tháng năm 2021 NHIỆM VỤ ĐỒ ÁN TỐT NGHIỆP Họ tên sinh viên: Nguyễn Khắc Hùng MSSV: 17145304 Email: 17145304@student.hcmute.edu.vn Số điện thoại: 0393716337 Khuất Đức Duy MSSV: 17145268 Email: 17145268@student.hcmute.edu.vn Số điện thoại: 0976215311 Ngành: Công nghệ Kỹ Thuật ô tơ Khóa: 2017-2021 Lớp: 171452A Tên đề tài: “Thiết kế, lắp đặt mơ hình điều hịa khơng khí ứng dụng phần mềm Matlab vào mô phỏng.” Nhiệm vụ đề tài • Ơn lại kiến thức tảng hệ thống điều hịa khơng khí tơ • Tra cứu mở rộng tài liệu có liên quan đề tài • Nghiên cứu sâu rộng giải thuật, lập trình thiết kế điện • Thi cơng cải tiến hồn thiện mơ hình điều hịa thực tế dựa vật tư xưởng • Mơ hệ thống điều hồ Sản phẩm đề tài • Mơ hình điều hịa thi cơng hồn thiện đặt xưởng điện khoa sẵn sàng cho việc học tập giảng dạy • Biên soạn hồn chỉnh sách thiết minh kèm mơ hình • Mơ hệ thống thành công MATLAB/SIMULINK Ngày giao nhiệm vụ đề tài: 10/03/2021 Ngày hoàn thành nhiệm vụ: Theo kế hoạch khoa CKĐ TRƯỞNG BỘ MÔN GIẢNG VIÊN HƯỚNG DẪN TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP.HCM CỘNG HỊA XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Độc Lập – Tự Do – Hạnh Phúc ──────── ──────── Bộ môn : Điện – Điện tử ô tô PHIẾU NHẬN XÉT ĐỒ ÁN TỐT NGHIỆP (Dành cho giảng viên hướng dẫn) Họ tên Sinh viên: Nguyễn Khắc Hùng MSSV: 17145304 Khuất Đức Duy MSSV: 17145268 Tên đề tài: “Thiết kế, lắp đặt mơ hình điều hịa khơng khí ứng dụng phần mềm Matlab vào mơ phỏng.” Ngành: Công nghệ Kỹ thuật ô tô Họ tên GV hướng dẫn: Th.S Nguyễn Thành Tuyên Ý KIẾN NHẬN XÉT Nhận xét tinh thần, thái độ làm việc sinh viên Nhận xét kết thực ĐATN 2.1 Kết cấu, cách thức trình bày ĐATN: 2.2 Nội dung đồ án: 2.3 Kết đạt được: 2.4 Những tồn (Nếu có): Đánh giá: TT Mục đánh giá Điểm tối đa Hình thức kết cấu ĐATN 30 Đúng format với đầy đủ hình thức nội dung 10 Mục tiêu, nhiệm vụ, tổng quan đề tài 10 Tính cấp thiết đề tài 10 Nội dung ĐATN 50 Khả ứng dụng kiến thức toán học, khoa học kỹ thuật, khoa học xã hội,… Khả thực hiện/phân tích/tổng hợp/đánh giá 10 Khả thiết kế chế tạo hệ thống, thành phần, quy trình đáp ứng yêu cầu đưa với ràng buộc thực tế 15 Khả cải tiến phát triển 15 Khả sử dụng công cụ kỹ thuật, phần mềm chuyên ngành… Đánh giá khả ứng dụng đề tài 10 Sản phẩm cụ thể ĐATN 10 Tổng điểm Điểm đạt 100 Kết luận  Được phép bảo vệ  Khơng phép bảo vệ Tp Hồ Chí Minh, ngày 25 tháng năm 2021 Giảng viên hướng dẫn TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP.HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA CƠ KHÍ ĐỘNG LỰC Độc Lập – Tự Do – Hạnh Phúc ──────── ──────── Bộ môn : Điện – Điện tử ô tô PHIẾU NHẬN XÉT ĐỒ ÁN TỐT NGHIỆP (Dành cho giảng viên phản biện) Họ tên Sinh viên: Nguyễn Khắc Hùng MSSV: 17145304 Khuất Đức Duy MSSV: 17145268 Tên đề tài: “Thiết kế, lắp đặt mơ hình điều hịa khơng khí ứng dụng phần mềm Matlab vào mơ phỏng.” Ngành: Công nghệ Kỹ thuật ô tô Họ tên GV phản biện: (Mã GV): Ý KIẾN NHẬN XÉT Kết cấu, cách thức trình bày ĐATN: Nội dung đồ án: Kết đạt được: Những kiến thức cịn thiếu sót Câu hỏi Đánh giá TT Mục đánh giá Điểm tối đa Hình thức kết cấu ĐATN 30 Đúng format với đầy đủ hình thức nội dung 10 Mục tiêu, nhiệm vụ, tổng quan đề tài 10 Tính cấp thiết đề tài 10 Nội dung ĐATN 50 Khả ứng dụng kiến thức toán học, khoa học kỹ thuật, khoa học xã hội,… Khả thực hiện/phân tích/tổng hợp/đánh giá 10 Khả thiết kế chế tạo hệ thống, thành phần, quy 15 Điểm đạt trình đáp ứng yêu cầu đưa với ràng buộc thực tế Khả cải tiến phát triển 15 Khả sử dụng công cụ kỹ thuật, phần mềm chuyên ngành… Đánh giá khả ứng dụng đề tài 10 Sản phẩm cụ thể ĐATN 10 Tổng điểm 100 Kết luận  Được phép bảo vệ  Không phép bảo vệ Tp Hồ Chí Minh, ngày 25 tháng 08 năm 2021 Giảng viên phản biện TRƯỜNG ĐH SƯ PHẠM KỸTHUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA CƠ KHÍ ĐỘNG LỰC ──────────── XÁC NHẬN HOÀN THÀNH ĐỒ ÁN Tên đề tài: Thiết kế, lắp đặt mơ hình điều hịa khơng khí ứng dụng phần mềm Matlab vào mô Họ tên Sinh viên: Nguyễn Khắc Hùng MSSV: 17145304 Khuất Đức Duy MSSV: 17145268 Ngành: Công nghệ Kỹ thuật ô tô Sau tiếp thu điều chỉnh theo góp ý Giảng viên hướng dẫn, Giảng viên phản biện thành viên Hội đồng bảo Đồ án tốt nghiệp hoàn chỉnh theo yêu cầu nội dung hình thức Chủ tịch Hội đồng: Giảng viên hướng dẫn: Giảng viên phản biện: Tp Hồ Chí Minh, ngày 25 tháng năm 2021 LỜI CẢM ƠN Lời đầu tiên, chúng em xin chân thành cảm ơn tồn thể q thầy trường Đại Học Sư Phạm Kỹ Thuật Thành phố Hồ Chí Minh nhiệt tình giảng dạy truyền đạt kiến thức quý báu cho chúng em suốt thời gian học tập trường Đặc biệt với giúp đỡ q thầy Khoa Cơ Khí Động Lực tạo điều kiện cho chúng em hoàn thành đồ án tốt nghiệp thời gian quy định Nhóm chúng em xin gửi lời cảm ơn chân thành đến: • Giảng viên Th.S Nguyễn Thành Tuyên, thầy hướng dẫn tận tình theo sát đơn đốc chúng em suốt q trình thực đồ án • Chúng em xin gửi lời cảm ơn đến thầy khoa Cơ khí động lực – trường Đại học Sư Phạm Kỹ Thuật hỗ trợ kiến thức, tư vấn kỹ thuật cho chúng em mượn dùng thiết bị xưởng thực hành • Cảm ơn gia đình bạn bè hỗ trợ, động viên, khuyến khích chúng em tự tin sống cố gắng vươn lên học tập Một lần nữa, chúng em xin chân thành cảm ơn kính chúc quý thầy cô Trường Đại Học Sư Phạm Kỹ Thuật Thành phố Hồ Chí Minh, đặc biệt quý thầy khoa Cơ Khí Động Lực dồi sức khỏe, niềm vui nhiệt huyết với nghề giáo để góp phần vào nghiệp trăm năm trồng người Chúng em xin chân thành cảm ơn! Nhóm thực đề tài Nguyễn Khắc Hùng Khuất Đức Duy i Hình 5.27: Đồ thị so sánh nhiệt độ môi chất lạnh van giãn nở với cơng suất khác 91 Hình 5.28: Đồ thị so sánh nhiệt độ môi chất lạnh dàn lạnh với công suất khác 5.9 Kết luận Sau hồn thành mơ hệ thống điều hoà MATLAB/SIMULINK , chúng em đạt mục tiêu đề trì nhiệt độ định cabin xe Biết nhiệt độ môi chất làm lạnh qua thành phần hệ thống điều hồ Việc xác thực cho tính thực nghiệm mơ hình xác minh việc thay đổi thông số đầu vào “công suất máy nén”,”khối lượng khơng khí bên cabin”,”số lượng hành khách xe” Từ kết có từ đồ thị cho thấy tính hợp lí mơ hình 92 CHƯƠNG 6: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN CỦA ĐỀ TÀI 6.1 Kết luận Sau thời gian nghiên cứu tài liệu với hướng dẫn theo sát thầy hướng dẫn nhóm chúng em thi cơng hồn chỉnh đề tài “Thiết kế, lắp đặt mơ hình điều hịa khơng khí ứng dụng phần mềm Matlab vào mô phỏng” đồng thời cải tiến thành công số tính giúp mơ hình dễ dàng tiếp cận đưa vào thực tiễn giảng dạy Với việc mơ hệ thống điều hồ tơ ứng dụng MATLAB/SIMULINK giúp chúng em hiểu rõ hệ thống Thấy yếu tố tác động đến hệ thống Xây dựng nên mơ hình mơ hệ thống điều hoà hỗ trợ cho bạn sinh viên dễ dàng tiếp cận hiểu nguyên lí hoạt động hệ thống Giúp bạn có nhìn trực quan Về phần thuyết minh, thuyết minh nhóm chúng em biên soạn rõ ràng dựa tài liệu chuyên ngành kiến thức chúng em tự rút được, đảm bảo cung cấp đầy đủ kiến thức hệ thống điều hịa cách vận hành mơ hình 6.2 Hướng phát triển Như biết hệ thống điều hịa tơ hệ thống tiện nghi thiết yếu hàng đầu ô tô đại Vì vậy, việc trang bị cho sinh viên ngành tô kiến thức lý thuyết hết mơ hình thực tiễn hệ thống điều hịa tô vấn đề cần thiết Song, để trang bị tơ nhằm mục đích giảng dạy học tập hệ thống điều hịa chi phí q lớn Do đó, mơ hình hệ thống điều hịa tự động tơ chúng em câu trả lời cho vấn đề Qua việc mơ hệ thống điều hồ MATLAB/SIMUINK với việc sử dụng định luật cân lượng để xây dựng mơ hình mơ Từ biết thơng số ảnh hưởng đến mơ hình để phát triển hệ thống điều hồ thực tế giúp hệ thống làm việc với hiệu suất cao tiết kiệm lượng 93 TÀI LIỆU THAM KHẢO Tiếng Việt [1] Nguyễn Oanh (2008), “ Ô tô hệ Điện lạnh ô tô”, nhà xuất Giao Thông Vận Tải, 202 trang [2] “Hệ thống sưởi ấm điều hịa tơ”, TS Lê Thanh Phúc [3] Diễn đàn Axeoto.com : https://axeoto.com/ [4] Diễn đàn Otohui.com : https://www.oto-hui.com/ [5] Tài liệu đào tạo kỹ thuật viên chuẩn đoán Toyota, “ Điều hịa khơng khí” [6] “Giáo trình phương pháp tính”, Th.S Lê Thị Thanh Hải Tiếng anh [1] 2005 CAMRY ELECTRICAL WIRING DIAGRAM, cardiagn.com, 251 Pages [2] Steven Daly (2006), “Automotive Air Conditioning and Climate Control Systems”, 382 pages [3] A State Space Approach for the Dynamic Analysis of Automotive Air Conditioning System [4] Boufadene, M (2018) Modeling and Control of AC Machine using MATLAB®/SIMULINK CRC Press 94 PHỤ LỤC clear all close all clc %% cabinspace mr = 9000.50; % mass of air inside cabin [g] cpr = 1.008; % specificheat at constant pressure at room in kJ/kg K Tin = 27; % Intial Temperature of cabin space Cpe = 1.008; % Constant pressure at evaprator in kJ/kg K P = 1.2 ; % density of air in kg/m³ Uo = 4; % heat transfer co-efficient of vechile wall Ao = 30; % surface area of cabin Ta = 21; % Temperature of ambient air Mf = 0.001 ; % mass flow rate to cabin g/s Cpa = 1.008; % Sp.heat at constant pressure at ambient point v = 8; % volume of cabin Qs = 50 ; % heat load due to solar radiation Qps = 70 ; % sensible heat load per passengers N=1 ; Kf = 1.1 ; % no of passanger % Heat gain co-efficient of fan V = 0.115 ; % volumetricair flow % % Wm = 0.058; % moisture content of the mixture air 95 Ws = 0.009; % moisture content of the supply air Wa = 0.0158; % Moisture content of the ambient Air Wr = 0.0113; % Moisture content of the air inside the vehicle %% compressor parameters n = 6; % number of cyclinders in compressor D = 0.14; % diameter of the cylinder sp = 0.49; % relative piston stroke sp1 = 10; % piston stroke Nc =11.50; % speed of compressor (rps) nv = (0.3596+(1.1072*sp)-(0.4025*sp*sp)+(0.0001175*Nc)-(2.449*1.0000e08*Nc*Nc)); % volumetric efficiency of compressor vs = 0.087; % Refrigerant specific volume f1 = 0.281; % flow rate1 f2 = 0.115; % flow rate2 cpr1 = 1.436; % Specific heat of Refrigirant cpr2 = 0.955; % Specific heat of Refrigirant %% Evaporator parameters Cp_E = 1.366; % Specific heat of Refrigirant in Evaprator hfg = 2450; % Latent heat of vaporization of water in kJ/kg Wm = 0.0113; % moisture content of the mixture air Wg = 0.009; % moisture content of the mixture air a1 = 0.125; % Heat transfer coefficient between air and evaporator wall a2 = 0.027; % Heat transfer coefficient between air and evaporator wall A1 = 0.53; % Heat transfer area of dry – cooling region in m A2 = 3.265; % Heat transfer area of wet – cooling region in m hr2 = 3938; % Enthalpy of refrigerant at the evaporator core inlet hr1 = 2335; % Enthalpy of refrigerant at the evaporator core outlet 96 Vh1 = 0.004; % Air side volume of the evaporator core Vh2 = 0.016; % Air side volume of the evaporator core %% Thermal Expansion valve Mfv = 0.01; % mass flow of vapour Mfl = 0.001; % mass flow of liquid cpv = 0.906; % Specific heat of Refrigirant in vapour phase cpl = 1.166; % Specific heat of Refrigirant in liquid phase %% condensor parameters Cp_C = 1.147; % Specific heat of Refrigirant in condensor hfg = 2450; % Latent heat of vaporization of water in kJ/kg Wm = 0.0113; % moiture content Wg = 0.009; % moiture content a1 = 0.125; % Heat transfer coefficient between air and evaporator wall a2 = 0.027; % Heat transfer coefficient between air and evaporator wall A1 = 2.83; % Heat transfer area of dry – cooling region in m A2 = 0.55; % Heat transfer area of wet – cooling region in m hr2 = 4288; % Enthalpy of refrigerant at the evaporator core inlet hr1 = 2335; % Enthalpy of refrigerant at the evaporator core outlet Vh1 = 0.004; % Air side volume of the evaporator core Vh2 = 0.016; % Air side volume of the evaporator core %% Connecting model sim_time = 1000; sim('AC_Vfinal') %% plotting graph results %% Ploting results Cabin figure('Name','Đầu hệ thống') % Naming The figure 97 set(gcf,'color','cyan') % -plot for Cabin Temperature % plot(T(:,1), T(:,2),'-r','linewidth',1) % DE of of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') ylim([21 28]) % setting axis limit grid on title('Đồ thị nhiệt độ cabin xe') legend('Nhiệt độ') %% Ploting results Evaporator figure('Name','Đầu dàn lạnh') % Naming The figure set(gcf,'color','cyan') % plot for Evaporator temperature % plot(Evp(:,1), Evp(:,2),'-r','linewidth',1) % DE of of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') grid on title('Đồ thị nhiệt độ môi chất lạnh dàn lạnh') legend('Nhiệt độ') %% Ploting results Compressor figure('Name','Đầu máy nén') % Naming The figure set(gcf,'color','cyan') % plot for Compressor Temperature % plot(Comp(:,1), Comp(:,2),'-r','linewidth',1) % DE of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') 98 ylim([0 70]) grid on title('Đồ thị nhiệt độ môi chất lạnh máy nén') legend('Nhiệt độ') %% Ploting results Condenser figure('Name','Đầu dàn nóng') % Naming The figure set(gcf,'color','cyan') % -plot for Condensor Temperature % plot(Con(:,1), Con(:,2),'-r','linewidth',1) % DE of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') grid on title('Đồ thị nhiệt độ mơi chất lạnh dàn nóng') legend('Nhiệt độ') %% Ploting results TEV figure('Name','Đầu van giãn nở') % Naming The figure set(gcf,'color','cyan') % -plot for Tev Temperature % plot(Tev(:,1), Tev(:,2),'-r','linewidth',1) % DE of Temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') grid on title('Đồ thị nhiệt độ môi chất lạnh van giãn nở') legend('Nhiệt độ') 99 %% Sensitivity analysis % -system response-analysis % % -For differrent cabin space i.e with different quantity of air % for mr = [6000 9000 12000] sim_time = 1000; sim('AC_Vfinal') figure(6) % Naming The figure set(gcf,'color','cyan') % plot for Cabin Temperature -% % nexttile plot(T(:,1), T(:,2)) % DE of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') ylim([21 28]) grid on title('Đồ thị thể ảnh hưởng không gian đến nhiệt độ cabin') end figure(6) Legend=cell(3,1); Legend{1}=' Hatchback'; Legend{2}=' Sedan'; Legend{3}=' SUV'; legend(Legend); % resetting the value to intial condition mr = 9000; % -% %% % Response for different no of passenger inside the car -% for N = [1 3] sim_time = 1000; 100 sim('AC_Vfinal') figure(7) % Naming The figure set(gcf,'color','cyan') % plot for Cabin Temperature -% plot(T(:,1), T(:,2)) % DE of of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') grid on ylim([21 28]) title('Đồ thị thể ảnh hưởng số lượng người xe đến nhiệt độ cabin') end figure(7) Legend=cell(3,1); Legend{1}=' người'; Legend{2}=' người'; Legend{3}=' người'; legend(Legend); % resetting the value N = 1; % -Components response-analysis % % -For different power input -% %% for Nc = [8.6 11.75 16.00] sim_time = 1000; sim('AC_Vfinal') figure(8) % Naming The figure % subplot(2,2,1) plot(Evp(:,1), Evp(:,2)) % DE of temp hold on 101 xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') title('Đồ thị so sánh nhiệt độ môi chất lạnh dàn lạnh với công suất khác nhau') figure(9) plot(Con(:,1), Con(:,2)) % DE of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') grid on title('Đồ thị so sánh nhiệt độ môi chất lạnh dàn nóng với cơng suất khác nhau') figure(10) plot(Comp(:,1), Comp(:,2)) % DE of temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') title('Đồ thị so sánh nhiệt độ môi chất lạnh máy nén với công suất khác nhau') figure(11) plot(Tev(:,1), Tev(:,2)) % DE of Temp hold on xlabel('Thời gian [s]') % Labeling axis ylabel('Nhiệt độ [°C]') grid on title('Đồ thị so sánh nhiệt độ môi chất lạnh van giãn nở với công suất khác nhau') end figure(8) Legend=cell(3,1); Legend{1}=' Công suất đầu vào nhỏ hơn'; Legend{2}=' Công suất đầu vào giữ nguyên'; Legend{3}=' Công suất đầu vào lớn hơn'; legend(Legend); figure(9) 102 Legend=cell(3,1); Legend{1}=' Công suất đầu vào nhỏ hơn'; Legend{2}=' Công suất đầu vào giữ nguyên'; Legend{3}=' Công suất đầu vào lớn hơn'; legend(Legend); figure(10) Legend=cell(3,1); Legend{1}=' Công suất đầu vào nhỏ hơn'; Legend{2}=' Công suất đầu vào giữ nguyên'; Legend{3}=' Công suất đầu vào lớn hơn'; legend(Legend); figure(11) Legend=cell(3,1); Legend{1}=' Công suất đầu vào nhỏ hơn'; Legend{2}=' Công suất đầu vào giữ nguyên'; Legend{3}=' Công suất đầu vào lớn hơn'; legend(Legend); Nc = 11.75; % presenting graphs % figure(1) movegui('northwest') figure(2) movegui('west') figure(3) movegui('southwest') figure(4) movegui('north') figure(5) movegui('center') 103 figure(6) movegui('south') figure(7) movegui('northeast') figure(8) movegui('east') figure(9) movegui('southeast') figure(10) movegui('southeast') figure(11) movegui('southeast') 104 S K L 0