Tính quan sát được trong bài toán đánh giá trạng thái lưới phân phối

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA - NGUYỄN HUY THÔNG TÍNH QUAN SÁT ĐƯC TRONG BÀI TOÁN ĐÁNH GIÁ TRẠNG THÁI LƯỚI PHÂN PHỐI CHUYÊN NGÀNH: MẠNG VÀ HỆ THỐNG ĐIỆN MÃ SỐ NGÀNH: 2.06.07 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, THÁNG 06 NĂM 2005 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 TP HỒ CHÍ MINH Cán hướng dẫn khoa học: TS PHAN THỊ THANH BÌNH 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ĩ 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 2005 Có thể tham khảo luận văn tại: THƯ VIỆN TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SĐH ĐỘC LẬP – TỰ DO – HẠNH PHÚC - Tp.HCM, ngaøy thaùng năm 2005 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN HUY THÔNG Phái: Nam Ngày sinh: 14 – 05 – 1965 Nơi sinh : TIỀN GIANG Chuyên ngành: KỸ THUẬT ĐIỆN MSHV: KTĐI13.012 I TÊN ĐỀ TÀI: TÍNH QUAN SÁT ĐƯC TRONG BÀI TOÁN ĐÁNH GIÁ TRẠNG THÁI LƯỚI PHÂN PHỐI II NHIỆM VỤ VÀ NỘI DUNG: - Giới thiệu tổng quan toán đánh giá trạng thái lưới phân phối - Phân tích phương pháp giải toán đánh giá trạng thái lưới phân phối, phân tích tính quan sát xử lý số liệu đo lường xấu - Giới thiệu thuật toán phù hợp - Thực viết chương trình áp dụng vào ví dụ III NGÀY GIAO NHIỆM VUÏ: IV NGÀY HOÀN THÀNH NHIỆM VỤ: V HỌ VÀ TÊN CB HƯỚNG DẪN: TS PHAN THỊ THANH BÌNH CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH Nội dung đề cương luận văn thạc só Hội đồng Chuyên Ngành thông qua Ngày PHÒNG ĐÀO TẠO SĐH tháng năm 2005 KHOA QUẢN LÝ NGÀNH Lời cảm ơn 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, đặc biệt CôPhan Thị Thanh Bình, người hướng dẫn thực luận án cao học Cám ơn tất bạn bè thân thuộc, bạn đồng nghiệp bên cạnh động viên giúp đỡ nhiều Đặc biệt cảm ơn vợ tôi, Phương Viên, người gánh vác việc gia đình, chăm sóc Bé Tí để hoàn thành Luận văn TP Hồ Chí Minh, ngày 29/06/2005 Nguyễn Huy Thông TĨM TẮT Mục tiêu hệ thống tự động hoá lưới phân phối hướng đến việc quản lý điều khiển lưới phân phối hiệu Điều kiện tiên để thực thành công hệ thống tự động lưới phân phối giải hiệu tốn phân tích tính quan sát được, phát liệu xấu đánh giá trạng thái lưới phân phối Trong Luận án giới thiệu giải thuật hiệu để đánh giá trạng thái áp dụng cho lưới phân phối ba pha hình tia Phương pháp mở rộng cho tính quan sát được, q trình xử lý liệu xấu nghiên cứu chi tiết Luận án Kỹ thuật “forward, backward” áp dụng để đánh giá, ước tính cơng suất nhánh, điện áp nút, phụ tải nút dựa số liệu đo lường Nếu khơng có đầy đủ số liệu đo lường phù hợp toán đánh giá trạng thái thực Phương pháp đề nghị Luận án áp dụng để phân tích ví dụ lưới phân phối Chương trình viết phần mềm Matlab 6.0 ABSTRACT A distribution automation system (DAS) aims for better management and control of the distribution networks A efficient network observability, bad data detection and state estimation solution technique is a prerequisite for the success of DAS This thesis presents an efficient and robust three phase state estimation (SE) algorithm for application to radial distribution networks The extension of the method to the network observability analysis and bad data processing is discussed in details in thesis Forward propagation and backward propagation technique will be applied to estimate the line flows, node voltages ang loads at each node, based on the measured quantities The SE cannot be excuted without adequate number of measurements The proposed method in this theis has been tested to analyze several examples of the distribution networks The program is written using Matlab 6.0 TÓM TẮT LÝ LỊCH TRÍCH NGANG Họ tên: Ngày, tháng, năm sinh: Nơi sinh: Địa liên lạc: Điện thoại quan: Nhà: Nguyễn Huy Thông 14/5/1965 Tiền Giang 242/47/47 Nguyễn Thiện Thuật P.3, Q.3, TP.HCM 307 945 399 831 – 0908 227 565 QUAÙ TRÌNH ĐÀO TẠO: Tháng 7/1996: Tốt nghiệp Đại Học Kỹ Thuật – ĐHQG TPHCM Năm 2002 – 2004: Học viên cao học khóa 13 Kỹ Thuật Điện – ĐH Bách Khoa TP.HCM QUÁ TRÌNH CÔNG TÁC: Năm 1996 – 1999: công tác Xí Nghiệp Thiết Kế Điện – Công ty Điện Lực TP.HCM Năm 1999 – 2005: công tác Công ty Tư Vấn Xây Dựng Điện –Tổng Công ty Điện Lực Việt Nam Trang Bảng danh mục từ viết tắt Luận văn A – ampere (đơn vị đo dòng điện) DMS – Distribution Management System DS – Disconnected Switch EMS – Energy Management System kV – kilovolt (đơn vị đo điện áp- 1kV = 103 V ) KVL – Kirchoff Voltage Law LAV – Least Absolute Value LBS – Load Breaker Switch LMS – Least Median of Squares LTD – Line Tension Disconnector MatLab – Matrix Laboratory: tên phần mềm pu – per unit (hệ đơn vị tương đối) RTU – Remote Telemetry Unit SCADA – System Control and Data Acquisition VAr – đơn vị đo công suất phản kháng W – watt (đơn vị đo công suất tác dụng) WLAV – Weighted Least Absolute Value WLS – Weighted Least Squares Tính quan sát toán đánh giá trạng thái lưới phân phối Trang Mục lục Mục lục Bảng danh mục từ viết tắt Luận văn Chương - Đặt vấn đề 10 1.1 Tổng quan phát triển ngành điện Việt Nam 10 1.2 Giới thiệu EMS/DMS 11 1.2.1 Hoàn cảnh đời SCADA/EMS 11 1.2.2 EMS (Energy Management System) 13 1.2.3 Các chức EMS 14 1.2.4 Các yêu cầu cần có EMS 14 1.2.5 Hệ thống EMS Trung tâm Điều Độ Hệ Thống Điện Miền Việt Nam 15 1.2.6 DMS (Distribution Management System) khác biệt với EMS 18 1.3 Đánh giá trạng thái hệ thống tính quan sát toán đánh giá trạng thái hệ thống 19 1.4 Phạm vi luận văn 20 Chương - Bài toán đánh giá trạng thái hệ thống 22 2.1 Tổng quan toán đánh giá trạng thái hệ thống 22 2.1.1 Lịch sử đời – phát triển toán đánh giá trạng thái 22 2.1.2 Đánh giá trạng thái hệ thống xây dựng mơ hình hệ thống 22 2.2 Những đặc điểm toán đánh giá trạng thái 24 2.2.1 Nhiệm vụ toán đánh giá trạng thái 24 2.2.2 Mơ tả tốn đánh giá trạng thái 24 2.2.3 Tiêu chuẩn xác suất cực đại 26 2.2.4 Cấu trúc ma trận Jacobi H 29 2.2.5 Nguồn liệu cho toán đánh giá trạng thái 31 2.2.6 Dữ liệu đầu vào toán đánh giá trạng thái 32 2.2.7 Dữ liệu đầu toán đánh giá trạng thái 33 2.2.8 Những khó khăn tốn đánh giá trạng thái 34 2.3 Tổng quan giải pháp giải toán đánh giá trạng thái 35 2.3.1 Bài toán đánh giá trạng thái động 36 2.3.2 Bài toán đánh giá trạng thái tĩnh 36 2.4 Bài toán đánh giá trạng thái lưới phân phối 42 2.4.1 Đặc điểm cấu trúc hệ thống phân phối 42 2.4.2 Đặc điểm giám sát tự động điều khiển hệ thống phân phối 43 2.4.3 Những khó khăn, thách thức tốn đánh giá trạng thái lưới phân phối 44 Tính quan sát toán đánh giá trạng thái lưới phân phối Trang 2.5 Kết luận 46 Chương - Tính quan sát toán đánh giá trạng thái lưới phân phối 48 3.1 Khái quát tính quan sát (observability) 48 3.1.1 Tổng quan 48 3.1.2 Mô hình tốn để nghiên cứu tính quan sát 48 3.1.3 Định nghĩa tính quan sát theo quan điểm toán học 49 3.1.4 Ví dụ tính quan sát 49 3.2 Tính quan sát tốn đánh giá trạng thái lưới phân phối 50 3.2.1 Giới thiệu 50 3.2.2 Khái niệm 50 3.3 Những phương pháp dùng để phân tích tính quan sát tốn đánh giá trạng thái lưới phân phối 51 3.3.1 Tổng quan 51 3.3.2 Phương pháp phân tích tính quan sát dựa vào việc phân tích số 51 3.3.3 Phương pháp phân tích tính quan sát dựa vào kỹ thuật “forward, backward” 53 3.4 Ảnh hưởng vấn đề đo lường đến tính quan sát toán đánh giá trạng thái lưới phân phối 54 3.5 Ảnh hưởng đo lường Ampere đến tính quan sát toán đánh giá trạng thái lưới phân phối 57 3.5.1 Giới thiệu 57 3.5.2 Đặt vấn đề 57 3.5.3 Đo lường dòng điện tới hạn phần dư đo lường mở rộng 58 3.5.4 Xử lý đo lường xấu có xuất đo lường Ampere 61 3.5.5 Áp dụng thuật tốn khảo sát tính quan sát vào lưới phân phối 62 3.6 Vai trị tính quan sát tốn đánh giá trạng thái 65 3.7 Kết luận 67 Chương - Thuật toán phân tích tính quan sát được, xử lý liệu đo lường xấu giải toán đánh giá trạng thái lưới phân phối 68 4.1 Mơ hình lưới phân phối trung 68 4.1.1 Mơ hình nút 68 4.1.2 Mơ hình nhánh 69 4.1.3 Nhận xét mơ hình phần tử lưới phân phối 70 4.2 Thuật tốn phân tích tính quan sát 70 4.3 Thuật toán xử lý số liệu xấu 74 4.4 Thuật toán giải toán đánh giá trạng thái lưới phân phối 74 4.4.1 Tổng quan giải thuật 74 Tính quan sát tốn đánh giá trạng thái lưới phân phối Trang 4.4.2 Lưu đồ giải thuật 76 4.4.3 Khởi gán ban đầu cho giá trị giải thuật 82 4.4.4 Điều kiện hội tụ 82 4.5 Nhận xét giải thuật đề nghị 82 Chương - Cấu trúc chương trình ví dụ áp dụng giải thuật 83 5.1 Tổng quan 83 5.2 Dữ liệu ban đầu 84 5.2.1 File “branch.xls” 84 5.2.2 File “bus.xls” 85 5.2.3 File “TableParentChild.xls” 85 5.3 Cấu trúc chương trình 86 5.3.1 Chương trình giải tốn đánh giá trạng thái lưới phân phối 86 5.3.2 Chương trình xử lý số liệu đo lường xấu 86 5.3.3 Chương trình phân tích tính quan sát 86 5.4 Dữ liệu kết 87 5.5 Ví dụ áp dụng 87 5.5.1 Các thông số mạng 87 5.5.2 Kết sau chạy chương trình phân tích tính quan sát lưới phân phối nút 88 5.5.3 Kết sau chạy chương trình đánh giá trạng thái lưới phân phối nút 89 5.6 Ví dụ áp dụng 92 5.6.1 Các thông số mạng 93 5.6.2 Kết sau chạy chương trình phân tích tính quan sát lưới phân phối 19 nút 95 5.6.3 Kết sau chạy chương trình đánh giá trạng thái lưới phân phối 19 nút 96 Chương - Kết luận 100 TÀI LIỆU KHAM KHẢO 101 Phụ lục – Kết thực chương trình 102 Phụ lục - Mã nguồn chương trình MatLab 110 Tính quan sát toán đánh giá trạng thái lưới phân phối Trang 136 % Lay tong cong suat cac nhanh NodeChild-TableNodeChild [Pa,Qa,Pb,Qb,Pc,Qc]=GetPQBranchs(NodeChild,TableNodeChild,branch); % Lay cong suat tai nut NodeChild [A1,A2,B1,B2,C1,C2]=GetPQBus(NodeChild,bus); % Gan tong cong suat cho nhanh NodeParent-NodeChild branch=SetPQBranch(NodeParent,NodeChild,branch,Pa+A1,Qa+A2,Pb+B1,Qb+B2,Pc+C 1,Qc+C2); % Dat cac chi so tai nhanh tuong ung Jpqf la Jpqf=SetJpqf(NodeParent,NodeChild,2,Jpqf); case % Nut NodeChild khong co tri so luong otherwise % Mot cac nhanh NodeChild-TableNodeChild khong co tri so luong % Lay tong cong suat cac nhanh NodeChild-TableNodeChild [Pa,Qa,Pb,Qb,Pc,Qc]=GetPQBranchs(NodeParent,NodeChild,branch); % Lay cong suat tai nut NodeChild [A1,A2,B1,B2,C1,C2]=GetPQBus(NodeChild,bus); Pa=Pa-A1; Qa=Qa-A2; Pb=Pb-B1; Qb=Qb-B2; Pc=Pc-C1; Qc=Qc-C2; % Lay tong cong suat cac nhanh NodeChild-TableNodeChild [A1,A2,B1,B2,C1,C2]=GetPQBranchs(NodeChild,TableNodeChild,branch); branch=SetPQBranch(NodeChild,TableNodeChild(location-2),branch,Pa-A1,QaA2,Pb-B1,Qb-B2,Pc-C1,Qc-C2); % Dat cac chi so tai nhanh tuong ung Jpqf la Jpqf=SetJpqf(NodeChild,TableNodeChild(location-2),2,Jpqf); end end i=i+1; end return; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % STEP SEVEN: TINH TOAN CONG SUAT TAI CAC NHANH CHUA CO TRI SO DO LUONG (FORWARD) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [Jpqd,bus]=Step07(TableParentChild,StatusDS,branch,bus,Jpqd,Jpqf) % Tao duong dan di nguoc tu cac nut den nut goc Path=Backward(TableParentChild,StatusDS); i=1; while (i Nut %d\n',TableParentChild); % In bang mo ta trang thai dao cach ly if isempty(StatusDS) fprintf(fid,'\n%s\n','Luoi phan phoi khong co dao cach ly phan doan'); else StatusDS=StatusDS'; fprintf(fid,'\n%s\n','Trang thai cac dao cach ly (1-dong; 0-mo )'); fprintf(fid,'Nut %d -> nut %d: %d\n',StatusDS); end Tính quan sát toán đánh giá trạng thái lưới phân phối Trang 140 % In ma tran tong tro nhanh fprintf(fid,'\n%s\n','Ma tran tong tro nhanh'); fprintf(fid,'%s\n','Zaa Zab Zac'); fprintf(fid,'%s\n','Zba Zbb Zbc'); fprintf(fid,'%s\n','Zca Zcb Zcc'); [row,col]=size(Zbranch); for r=1:row for c=1:col switch c case fprintf(fid,'Tu nut %d ',Zbranch(r,c)); case fprintf(fid,'den nut %d:\n',Zbranch(r,c)); otherwise if (c==5)|(c==8)|(c==11) if (imag(Zbranch(r,c))>=0) fprintf(fid,' %.5f+%.5fi\n',real(Zbranch(r,c)),imag(Zbranch(r,c))); else fprintf(fid,' %.5f-%.5fi\n',real(Zbranch(r,c)),abs(imag(Zbranch(r,c)))); end else if (imag(Zbranch(r,c))>=0) fprintf(fid,' %.5f+%.5fi ',real(Zbranch(r,c)),imag(Zbranch(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(Zbranch(r,c)),abs(imag(Zbranch(r,c)))); end end end end end % In ma tran cong suat luong tren nhanh fprintf(fid,'\n%s\n','Cong suat luong tren nhanh'); fprintf(fid,'%s\n',' Sa Sb [row,col]=size(PQMbranch); for r=1:row for c=1:col switch c case fprintf(fid,'Tu nut %d ',PQMbranch(r,c)); case fprintf(fid,'den nut %d: ',PQMbranch(r,c)); otherwise if (c==5) if imag(PQMbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQMbranch(r,c)),imag(PQMbranch(r,c))); Tính quan sát tốn đánh giá trạng thái lưới phân phối Sc'); Trang 141 else fprintf(fid,' %.5f%.5fi\n',real(PQMbranch(r,c)),abs(imag(PQMbranch(r,c)))); end else if imag(PQMbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQMbranch(r,c)),imag(PQMbranch(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQMbranch(r,c)),abs(imag(PQMbranch(r,c)))); end end end end end % In ma tran cong suat tinh toan tren nhanh fprintf(fid,'\n%s\n','Cong suat tinh toan tren nhanh'); fprintf(fid,'%s\n',' Sa Sb Sc'); [row,col]=size(PQCbranch); for r=1:row for c=1:col switch c case fprintf(fid,'Tu nut %d ',PQCbranch(r,c)); case fprintf(fid,'den nut %d: ',PQCbranch(r,c)); otherwise if (c==5) if imag(PQCbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQCbranch(r,c)),imag(PQCbranch(r,c))); else fprintf(fid,' %.5f%.5fi\n',real(PQCbranch(r,c)),abs(imag(PQCbranch(r,c)))); end else if imag(PQCbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQCbranch(r,c)),imag(PQCbranch(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQCbranch(r,c)),abs(imag(PQCbranch(r,c)))); end end end end end Tính quan sát toán đánh giá trạng thái lưới phân phối Trang 142 % In ma tran so tren nhanh fprintf(fid,'\n%s\n','Trong so tren nhanh'); fprintf(fid,'%s\n','Tu nut den nut Jpf Jqf W'); [row,col]=size(Jpqf); for r=1:row for c=1:col switch c case fprintf(fid,' %d ',Jpqf(r,c)); fprintf(fid,'%.3f\n',WBranch(Jpqf(r,1),Jpqf(r,2),Jpqf)); otherwise fprintf(fid,' %d ',Jpqf(r,c)); end end end % In ma tran dien ap nut fprintf(fid,'\n%s\n','Dien ap tai cac nut '); fprintf(fid,'%s\n',' Va Vb Vc Fia [row,col]=size(Vabc); for r=1:row for c=1:col switch c case fprintf(fid,'Nut %d: ',Vabc(r,c)); case % Khong lam gi ca otherwise if (c==8) fprintf(fid,'%.5f\n',Vabc(r,c)); else fprintf(fid,'%.5f ',Vabc(r,c)); end end end end Fib % In ma tran cong suat luong cua phu tai tai nut fprintf(fid,'\n%s\n','Phu tai luong tai cac nut'); fprintf(fid,'%s\n','Nut Sa Sb Sc'); [row,col]=size(PQML); for r=1:row for c=1:col switch c case fprintf(fid,' %d ',PQML(r,c)); case Tính quan sát toán đánh giá trạng thái lưới phân phối Fic'); Trang 143 % Khong lam gi ca otherwise if (c==5) if imag(PQML(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQML(r,c)),imag(PQML(r,c))); else fprintf(fid,' %.5f+%.5fi\n',real(PQML(r,c)),abs(imag(PQML(r,c)))); end else if imag(PQML(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQML(r,c)),imag(PQML(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQML(r,c)),abs(imag(PQML(r,c)))); end end end end end % In ma tran cong suat luong cua phu tai tai nut fprintf(fid,'\n%s\n','Phu tai tinh toan tai cac nut'); fprintf(fid,'%s\n','Nut Sa Sb Sc'); [row,col]=size(PQCL); for r=1:row for c=1:col switch c case fprintf(fid,' %d ',PQCL(r,c)); case % Khong lam gi ca otherwise if (c==5) if imag(PQCL(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQCL(r,c)),imag(PQCL(r,c))); else fprintf(fid,' %.5f-%.5fi\n',real(PQCL(r,c)),abs(imag(PQCL(r,c)))); end else if imag(PQCL(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQCL(r,c)),imag(PQCL(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQCL(r,c)),abs(imag(PQCL(r,c)))); end end end end end Tính quan sát tốn đánh giá trạng thái lưới phân phối Trang 144 % In ma tran so tai nut fprintf(fid,'\n%s\n','Trong so tai nut'); fprintf(fid,'%s\n','Nut Jpd Jqd W'); [row,col]=size(Jpqd); for r=1:row for c=1:col switch c case fprintf(fid,' %d ',Jpqd(r,c)); fprintf(fid,'%.3f\n',WNode(Jpqd(r,1),Jpqd)); otherwise fprintf(fid,' %d ',Jpqd(r,c)); end end end else % Ghi cac ket qua tinh toan permission='a'; fid = fopen(filename,permission); if (Iteration==1) fprintf(fid,'\n%s','********************************************'); fprintf(fid,'\n%s\n','KET QUA TINH TOAN'); fprintf(fid,'%s\n','********************************************'); end fprintf(fid,'\n************* LAN LAP THU %d *************',Iteration); % In ma tran dien ap nut fprintf(fid,'\n%s\n','Dien ap tai cac nut '); fprintf(fid,'%s\n',' Va Vb Vc Fia [row,col]=size(Vabc); for r=1:row for c=1:col switch c case fprintf(fid,'Nut %d: ',Vabc(r,c)); case % Khong lam gi ca otherwise if (c==8) fprintf(fid,'%.5f\n',Vabc(r,c)); else fprintf(fid,'%.5f ',Vabc(r,c)); end end end end Fib % In ma tran cong suat luong tren nhanh Tính quan sát tốn đánh giá trạng thái lưới phân phối Fic'); Trang 145 fprintf(fid,'\n%s\n','Cong suat luong tren nhanh'); fprintf(fid,'%s\n',' Sa Sb [row,col]=size(PQMbranch); for r=1:row for c=1:col switch c case fprintf(fid,'Tu nut %d ',PQMbranch(r,c)); case fprintf(fid,'den nut %d: ',PQMbranch(r,c)); otherwise if (c==5) if imag(PQMbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQMbranch(r,c)),imag(PQMbranch(r,c))); else fprintf(fid,' %.5f%.5fi\n',real(PQMbranch(r,c)),abs(imag(PQMbranch(r,c)))); end else if imag(PQMbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQMbranch(r,c)),imag(PQMbranch(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQMbranch(r,c)),abs(imag(PQMbranch(r,c)))); end end end end end % In ma tran cong suat tinh toan tren nhanh fprintf(fid,'\n%s\n','Cong suat tinh toan tren nhanh'); fprintf(fid,'%s\n',' Sa Sb [row,col]=size(PQCbranch); for r=1:row for c=1:col switch c case fprintf(fid,'Tu nut %d ',PQCbranch(r,c)); case fprintf(fid,'den nut %d: ',PQCbranch(r,c)); otherwise if (c==5) if imag(PQCbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQCbranch(r,c)),imag(PQCbranch(r,c))); Tính quan sát toán đánh giá trạng thái lưới phân phối Sc'); Sc'); Trang 146 else fprintf(fid,' %.5f%.5fi\n',real(PQCbranch(r,c)),abs(imag(PQCbranch(r,c)))); end else if imag(PQCbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQCbranch(r,c)),imag(PQCbranch(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQCbranch(r,c)),abs(imag(PQCbranch(r,c)))); end end end end end % In ma tran cong suat trung binh so tren nhanh fprintf(fid,'\n%s\n','Cong suat trung binh co so tren nhanh'); fprintf(fid,'%s\n',' Sa Sb Sc'); [row,col]=size(PQWbranch); for r=1:row for c=1:col switch c case fprintf(fid,'Tu nut %d ',PQWbranch(r,c)); case fprintf(fid,'den nut %d: ',PQWbranch(r,c)); otherwise if (c==5) if imag(PQWbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi\n',real(PQWbranch(r,c)),imag(PQWbranch(r,c))); else fprintf(fid,' %.5f%.5fi\n',real(PQWbranch(r,c)),abs(imag(PQWbranch(r,c)))); end else if imag(PQWbranch(r,c))>=0 fprintf(fid,' %.5f+%.5fi ',real(PQWbranch(r,c)),imag(PQWbranch(r,c))); else fprintf(fid,' %.5f-%.5fi ',real(PQWbranch(r,c)),abs(imag(PQWbranch(r,c)))); end end end end end Tính quan sát toán đánh giá trạng thái lưới phân phối Trang 147 end fclose(fid); return; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function W=WNode(Node,Jpqd) % So hang cua ma tran dong dien nhanh tinh toan [a,b]=find(Jpqd(:,1)==Node); switch Jpqd(a,3) case W=0.5; case W=0.45; case W=0.4; end return; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % WDATOb2FILE: GHI CAC KET QUA CHUONG TRINH OBSERPRO.M VAO FILE %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function WDatOb2File(filename,permission,bus,branch,Jpqd,Jpqf,ObIsGood) switch permission case 'w' fid = fopen(filename,permission); fprintf(fid,'%s\n','********************************************'); fprintf(fid,'%s\n',' KET QUA CUA CHUONG TRINH PHAN TICH TINH QUAN SAT DUOC'); fprintf(fid,'%s\n','********************************************'); case 'a' fid = fopen(filename,permission); switch ObIsGood case fprintf(fid,'\n%s\n','***LUOI PHAN PHOI DANG XET KHONG CO TINH QUAN SAT DUOC***'); case fprintf(fid,'\n%s\n','***SAU KHI BO SUNG CAC SO LIEU DO LUONG GIA***'); case fprintf(fid,'\n%s\n','***LUOI PHAN PHOI DANG XET CO TINH QUAN SAT DUOC***'); end % In cong suat luong cua phu tai tai nut fprintf(fid,'\n%s\n','Phu tai luong tai cac nut'); fprintf(fid,'%s\n','Nut Pa Qa Pb Qb Pc Qc'); Tính quan sát toán đánh giá trạng thái lưới phân phối Trang 148 [row,col]=size(bus); for r=1:row for c=1:col switch c case fprintf(fid,' %d ',bus(r,c)); case % Khong lam gi ca otherwise if (c