Hcmute xác định vị trí chống sét van trên cơ sở giảm rủi ro hư hỏ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 BÁO CÁO TỔNG KẾT ÐỀ TÀI KH&CN CẤP TRUỜNG TRỌNG ÐIỂM XÁC ÐỊNH VỊ TRÍ CHỐNG SÉT VAN TRÊN CƠ SỞ GIẢM RỦI RO HƯ HỎNG MÃ SỐ: T2013-01TÐ S K C0 Tp Hồ Chí Minh, tháng 12 - 2013 Luan van 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 BÁO CÁO TỔNG KẾT ĐỀ TÀI KH&CN CẤP TRƯỜNG TRỌNG ĐIỂM XÁC ĐỊNH VỊ TRÍ CHỐNG SÉT VAN TRÊN CƠ SỞ GIẢM RỦI RO HƯ HỎNG MÃ SỐ: T2013-01TĐ Chủ nhiệm đề tài: PGS.TS Quyền Huy Ánh Tp Hồ Chí Minh, tháng 12 năm 2013 Luan van 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 -o0o - BÁO CÁO TỔNG KẾT ĐỀ TÀI KH&CN CẤP TRƯỜNG TRỌNG ĐIỂM XÁC ĐỊNH VỊ TRÍ CHỐNG SÉT VAN TRÊN CƠ SỞ GIẢM RỦI RO HƯ HỎNG Mã số: T2013-01TĐ Chủ nhiệm đề tài: PGS-TS Quyền Huy Ánh TP HCM, 12/2013 Luan van DANH SÁCH THÀNH VIÊN VÀ ĐƠN VỊ PHỐI HỢP CHÍNH PGS.TS Quyền Huy Ánh, Đại học Sư phạm Kỹ thuật Tp HCM, CN đề tài ThS Lê Quang Trung, Cao đẳng nghề Lilama, đơn vị phối hợp i Luan van MỤC LỤC Trang Trang bìa phụ i Danh sách thành viên tham gia NCĐT đơn vị phối hợp ii Mục lục iii Danh mục bảng biểu Danh mục chữ viết tắt Danh mục hình vẽ .3 Thông tin kết nghiên cứu tiếng Việt tiếng Anh MỞ ĐẦU 10 0.1 Tổng quan tình hình nghiên cứu ngồi nước 10 0.2 Tính cấp thiết đề tài 10 0.3 Mục tiêu đề tài 11 0.4 Cách tiếp cận 11 0.5 Phương pháp nghiên cứu 11 0.6 Đối tượng phạm vi nghiên cứu 11 0.7 Nội dung nghiên cứu 11 CHƯƠNG 1: TỔNG QUAN VỀ CÁC PHƯƠNG PHÁP LẮP ĐẶT CHỐNG SÉT VAN 1.1 Các phương pháp lựa chọn vị trí lắp đặt chống sét van 13 1.2 Nhận xét 23 CHƯƠNG 2: CƠ SỞ TÍNH TỐN 24 2.1 Rủi ro hư hỏng 24 2.2 Logic mờ 29 2.3 Mô SimPowerSystem 31 ii Luan van CHƯƠNG 3: CẤU HÌNH MẠNG ĐIỆN 33 3.1 Các cấu hình trạm biến áp .33 3.2 Mơ hình phần tử mạng 34 3.3 Mơ hình mạng điện cần mơ 39 CHƯƠNG 4: PHƯƠNG PHÁP & CHƯƠNG TRÌNH TÍNH TỐN 42 4.1 Đặc tính ngẫu nhiên xung sét 42 4.2 Tiêu chí bảo vệ 43 4.3 Sơ đồ khối chương trình 44 4.4 Chương trình tính tốn FUPOSAR 54 CHƯƠNG 5: KẾT LUẬN 61 TÀI LIỆU THAM KHẢO 62 PHỤ LỤC ……………………………………………………………………… 63 Phụ lục A Phần code Matlab giao diện chương trình FUPOSA dạng file.m Phụ lục B Phần code Malab chương trình tính rủi ro cho cấu hình máy biến áp dạng file.m Phụ lục C Phần code Malab chương trình tính rủi ro cho cấu hình hai máy biến áp dạng file.m Phụ lục D Hợp đồng triển khai nhiệm vụ khoa học công nghệ cấp trường trọng điểm năm 2013, số T2013-01TĐ/KHCN-GV Phụ lục E Bài báo đăng tạp chí Phát triển Khoa học Công nghệ, Đại học Quốc gia Tp Hồ Chí Minh, tập 15, số K2/2012 iii Luan van Đề tài nghiên cứu khoa học T2013-01TĐ DANH MỤC BẢNG BIỂU Trang Bảng 1.1 Điện áp cực đại máy biến áp bảo vệ chống sét van 16 Bảng 1.2 Mức cách điện (BIL) mức điện áp bảo vệ chống sét van (Up) đại với Up = 4pu 20 Bảng 1.3 Độ dốc giá trị điện áp phóng điện cách điện đường dây không 20 Bảng 1.4 Chiều dài cho phép lớn Dk đoạn cáp với bên bảo vệ thiết bị chống sét 21 Bảng 1.5 Khoảng cách phân cách tối đa cho phép chống sét cáp máy biến áp Hình 1.8 với b = Hai chống sét van nối hai đầu cáp, MBA khơng có chống sét van 22 Bảng 3.1: Các dạng sóng tiêu chuẩn 38 Bảng 4.1 : Luật mờ định mờ xác định vị trí lắp đặt chống sét van 53 PGS-TS Quyền Huy Ánh Luan van Đề tài nghiên cứu khoa học T2013-01TĐ DANH MỤC CÁC CHỮ VIẾT TẮT BIL: Basic Insulation Level MOV: Metal Oxide Varistor AI: Artificial Intelligence FDCL: Fuzzy Dependency and Command Language MIQ: Machine Intelligence quotient MBA : Máy biến áp CSV : Chống sét van PGS-TS Quyền Huy Ánh Luan van Đề tài nghiên cứu khoa học T2013-01TĐ DANH MỤC HÌNH VẼ HÌNH TRANG Hình 1.1 Mơ hình Petersen bảo vệ trạm 13 Hình 1.2 Sơ đồ bảo vệ trạm có xét đến khoảng cách phân cách 14 Hình 1.3 Mơ hình Benoit de Metz-Noblat đường dây trạm biến áp 15 Hình 1.4 Sóng tới sóng phản xạ trạm biến áp có chống sét van 16 Hình 1.5 Điện áp trên đường dây bảo vệ chống sét van theo thời gian (trường hợp Bảng 1.1) 17 Hình 1.6 Quá áp cực đại đầu cực máy biến áp (B) theo thời gian  (trường hợp Bảng 1.1) 18 Hình 1.7 Quá điện áp cuối đường dây 19 Hình 1.8 Trạm biến áp kết nối với cáp ngầm 21 Hình 2.1: Dạng sóng hàm mật độ xác suất xuất áp, hàm xác suất phóng điện đánh thủng rủi ro hư hỏng 24 Hình 2.2: Lưu đồ tính tốn rủi ro hư hỏng 26 Hình 2.3 Đặc tuyến hệ số phối hợp thống kê rủi ro hư hỏng 27 Hình 2.4: Đặc tuyến hàm xác suất phóng điện đánh thủng P(V) theo số lượng lớp cách điện đồng tâm 28 Hình 2.5: Lưu đồ tính toán rủi ro hư hỏng phương pháp đồ thị 28 Hình 3.1: Cấu hình trạm máy biến áp 33 Hình 3.2: Cấu hình trạm hai máy biến áp 34 Hình 3.3: Mơ hình hình PI cho đường dây 35 Hình 3.4 : Các tham số lựa chọn mơ hình hình PI 35 PGS-TS Quyền Huy Ánh Luan van Đề tài nghiên cứu khoa học T2013-01TĐ Hình 3.5 : Mơ hình tổng trở đặc tính 36 Hình 3.6 : Mơ hình tổng trở đặc tính đơn giản 36 Hình 3.7 : Đặc tuyến V-I chống sét van 37 Hình 3.8 : Các thơng số lựa chọn chống sét van 37 Hình 3.9: Dạng xung sét tiêu chuẩn 38 Hình 3.10: Thơng số mơ hình nguồn phát xung sét 39 Hình 3.11 Mơ hình cho trạm máy biến áp 39 Hình 3.12 Mơ hình trạm hai máy biến áp chống sét van nằm đoạn cáp ngầm khác 40 Hình 3.13 Mơ hình trạm hai máy biến áp sử dụng hai chống sét van 41 Hình 4.1: Xác suất xuất xung sét có đại lượng ngẫu nhiên lớn giá trị xác định 43 Hình 4.2: Phân phối giá trị biến ngẫu nhiên 43 Hình 4.3: Lưu đồ tính tốn bảo vệ theo rủi ro trung bình cho cấu hinh máy biến áp 46 Hình 4.4: Lưu đồ sử dụng chống sét van để bảo vệ cho trạm hai máy biến áp 48 Hình 4.5: Lưu đồ sử dụng hai chống sét van để bảo vệ cho trạm hai máy biến áp 49 Hình 4.6 Sơ đồ khối chương trình máy tính xác định vị trí lắp đặt chống sét van sử dụng logic mờ theo tiêu chuẩn rủi ro cho phép 50 Hình 4.7 Lưu đồ chương trình tính tốn 51 Hình 4.8:Cấu tạo định mờ xác định vị trí chống sét van 52 Hình 4.9: Các biến ngôn ngữ ngõ vào sai số rủi ro hư hỏng lần lặp thứ (k) 53 PGS-TS Quyền Huy Ánh Luan van % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - Executes on button press in checkbox_eco function checkbox_eco_Callback(hObject, eventdata, handles) global cri cri=1; maxVal = get(handles.checkbox_eco,'Min') set(handles.checkbox_per,'Value',maxVal); set(handles.edit_adrisk,'Enable','off'); set(handles.edit_sl,'Enable','on'); % hObject handle to checkbox_eco (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of checkbox_eco % - Executes on button press in checkbox_per function checkbox_per_Callback(hObject, eventdata, handles) global cri cri=2; maxVal = get(handles.checkbox_per,'Min'); set(handles.checkbox_eco,'Value',maxVal); set(handles.edit_sl,'Enable','off'); set(handles.edit_adrisk,'Enable','on'); % hObject handle to checkbox_per (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of checkbox_per function edit_sl_Callback(hObject, eventdata, handles) % hObject handle to edit_sl (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_sl as text % str2double(get(hObject,'String')) returns contents of edit_sl as a double % - Executes during object creation, after setting all properties function edit_sl_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_sl (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_adrisk_Callback(hObject, eventdata, handles) Luan van % hObject handle to edit_adrisk (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_adrisk as text % str2double(get(hObject,'String')) returns contents of edit_adrisk as a double % - Executes during object creation, after setting all properties function edit_adrisk_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_adrisk (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - Executes on button press in pushbutton3 function pushbutton3_Callback(hObject, eventdata, handles) closereq open('front_open.fig') % hObject handle to pushbutton3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % - Executes on button press in checkbox_Idata function checkbox_Idata_Callback(hObject, eventdata, handles) global d; d=2; v=get(handles.checkbox_Idata,'Min') set(handles.checkbox_Cdata,'Value',v); set(handles.edit_PV,'Enable','off'); set(handles.edit_SV,'Enable','off'); set(handles.edit_NV,'Enable','off'); set(handles.edit_im,'Enable','on'); % hObject handle to checkbox_Idata (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of checkbox_Idata % - Executes on button press in checkbox_Cdata function checkbox_Cdata_Callback(hObject, eventdata, handles) global d; d=1; v=get(handles.checkbox_Cdata,'Min') set(handles.checkbox_Idata,'Value',v); set(handles.edit_PV,'Enable','on'); set(handles.edit_SV,'Enable','on'); set(handles.edit_NV,'Enable','on'); set(handles.edit_im,'Enable','off'); % hObject handle to checkbox_Cdata (see GCBO) Luan van % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of checkbox_Cdata function edit_PV_Callback(hObject, eventdata, handles) % hObject handle to edit_PV (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_PV as text % str2double(get(hObject,'String')) returns contents of edit_PV as a double % - Executes during object creation, after setting all properties function edit_PV_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_PV (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_SV_Callback(hObject, eventdata, handles) % hObject handle to edit_SV (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_SV as text % str2double(get(hObject,'String')) returns contents of edit_SV as a double % - Executes during object creation, after setting all properties function edit_SV_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_SV (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function edit_im_Callback(hObject, eventdata, handles) % hObject handle to edit_nv (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_nv as text % str2double(get(hObject,'String')) returns contents of edit_nv as a double Luan van % - Executes during object creation, after setting all properties function edit_im_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_nv (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - Executes on button press in pushbutton1 function pushbutton1_Callback(hObject, eventdata, handles) global filename pathname [filename,pathname]=uigetfile('*.xls','nhan file du lieu duoc chon'); set(handles.edit_im,'string', [pathname,filename]); % hObject handle to pushbutton1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) function edit_NV_Callback(hObject, eventdata, handles) % hObject handle to edit_NV (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of edit_NV as text % str2double(get(hObject,'String')) returns contents of edit_NV as a double % - Executes during object creation, after setting all properties function edit_NV_CreateFcn(hObject, eventdata, handles) % hObject handle to edit_NV (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - Executes during object creation, after setting all properties function ax2_CreateFcn(hObject, eventdata, handles) I=imread('ch2.JPG'); imshow(I); % hObject handle to ax2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: place code in OpeningFcn to populate ax2 % - Executes on button press in pushbutton_sol function pushbutton_sol_Callback(hObject, eventdata, handles) % lay gia tri thong so mo hinh % l=get(handles.edit_L,'string'); % lay gia tri chieu dai duong day Luan van % l= str2num(l); % chuyen doi LOA=get(handles.edit_LOA,'string'); % Chieu dai doan OA la 50m LOA=str2num(LOA); LABt=get(handles.edit_AB,'string'); LABt=str2num(LABt); LACt=get(handles.edit_AC,'string'); LACt=str2num(LACt); LAB=min(LABt,LACt); %Chieu dai doan AB=min LAC=max(LABt,LACt); %Chieu dai doan AC=max sl=get(handles.edit_sl,'string'); % sl la so luong arester sl=str2num(sl); dl=LOA/15; l1=[0:dl:LOA]; l3=[0:dl:LAC]; m=length(l1); C=get(handles.edit_C,'string'); C=str2num(C); V=get(handles.edit_V,'string'); V=str2num(V); Re=get(handles.edit_R,'string'); Re=str2num(Re); % - Lay gia tri thong so mo hinh -% - tao du lieu xung set global d cri pathname filename if d == aa=get(handles.edit_PV,'string'); % lay gia tri chieu dai duong day [peak1 peak2]= strtok(aa); peakmu=str2num(peak1); % chuyen doi peaksima=str2num(peak2); bb=get(handles.edit_SV,'string'); % lay gia tri chieu dai duong day [slope1 slope2]= strtok(bb); slopemu=str2num(slope1); % chuyen doi slopesima=str2num(slope2); sl1=get(handles.edit_NV,'string'); sl1=str2num(sl1); A=normrnd(peakmu,peaksima,sl1,1); S=normrnd(slopemu,slopesima,sl1,1);A=A';S=S'; A=abs(A);S=abs(S); % A=exprnd(pu1,sl1,1); S=exprnd(pu2,sl1,1);A=A';S=S'; end if d==2 path=[pathname filename]; data=xlsread(path); A=data(:,2); S=data(:,3); A=A'; S=S'; end n=length(A); Tds=[]; Ts=[]; for i=1:n tds=A(i)/S(i); Tds=[Tds tds]; ts= tds*2.5; Ts=[Ts ts]; end % -tao du lieu xung set Luan van % - tinh toan voi bao ve bang arester tieu chuan trung binh be %nhat -if (cri==1)&(sl==1) LAB=min(LABt,LACt); %Chieu dai doan AB=min LAC=max(LABt,LACt); %Chieu dai doan AC=max dl=LOA/15; l1=[0:dl:LOA]; l2=[0:dl:LAB]; m=length(l1); R1=[]; R2=[]; R3=[]; R4=[];R=[]; MU1=[]; MU2=[]; MU3=[]; MU4=[]; SIGMA1=[]; SIGMA2=[]; SIGMA3=[]; SIGMA4=[]; L1=LOA/1000; L2=LAB/1000; L3=LAC/1000; for d=1:m if (d==1) la=(l1(d)+0.01)/1000; %chuyen sang he don vi km elseif (d==m) la=(l1(d)-0.01)/1000; else la=l1(d)/1000; end lb=LOA/1000-la; NODE1=[]; NODE2=[]; NODE3=[]; NODE4=[]; for i=1:n tds= Tds(i); ts= Ts(i); t=tds*1.5*1e-6; Imax= A(i); options = simset('SrcWorkspace','current'); sim('cauhinh2a1',[],options); max1=max(node1); max2=max(node2); max3=max(node3); max4=max(node4); NODE1=[NODE1 max1]; NODE2=[NODE2 max2]; NODE3=[NODE3 max3]; NODE4=[NODE4 max4]; end [mu1,sigma1] = normfit(NODE1); MU1=[MU1, mu1] ; SIGMA1=[SIGMA1, sigma1]; [mu2,sigma2] = normfit(NODE2); MU2=[MU2, mu2]; SIGMA2=[SIGMA2, sigma2]; [mu3,sigma3] = normfit(NODE3); MU3=[MU3, mu3] ; SIGMA3=[SIGMA3, sigma3]; [mu4,sigma4] = normfit(NODE4); MU4=[MU4, mu4] ; SIGMA4=[SIGMA4, sigma4]; r1=risk(mu1,sigma1,1); r2=risk(mu2,sigma2,1); r3=risk(mu3,sigma3,1); Luan van r4=risk(mu4,sigma4,1); R1=[R1,r1]; R2=[R2,r2]; R3=[R3,r3]; R4=[R4,r4]; r=(r3+r4+r1*0.1+r2*0.1)/(1+1+0.1+0.1); R=[R,r]; end l2=[0:dl:LAB]; assignin('base', 'LAB', LAB) assignin('base', 'l2', l2) assignin('base', 'l1', l1) m=length(l2); for d=1:m if (d==1) la1=(l2(d)+0.01)/1000; %chuyen sang he don vi km elseif (d==m) la1=(l2(d)-0.01)/1000; else la1=l2(d)/1000; end lb1=LAB/1000-la1; NODE1=[]; NODE2=[]; NODE3=[]; NODE4=[]; for i=1:n tds= Tds(i); ts= Ts(i); t=tds*1.5*1e-6; Imax= A(i); options = simset('SrcWorkspace','current'); sim('cauhinh2b1',[],options); max1=max(node1); max2=max(node2); max3=max(node3); max4=max(node4); NODE1=[NODE1 max1]; NODE2=[NODE2 max2]; NODE3=[NODE3 max3]; NODE4=[NODE4 max4]; end [mu1,sigma1] = normfit(NODE1); MU1=[MU1, mu1] ; SIGMA1=[SIGMA1, sigma1]; [mu2,sigma2] = normfit(NODE2); MU2=[MU2, mu2]; SIGMA2=[SIGMA2, sigma2]; [mu3,sigma3] = normfit(NODE3); MU3=[MU3, mu3] ; SIGMA3=[SIGMA3, sigma3]; [mu4,sigma4] = normfit(NODE4); MU4=[MU4, mu4] ; SIGMA4=[SIGMA4, sigma4]; r1=risk(mu1,sigma1,1); r2=risk(mu2,sigma2,1); r3=risk(mu3,sigma3,1); r4=risk(mu4,sigma4,1); R1=[R1,r1]; R2=[R2,r2]; R3=[R3,r3]; R4=[R4,r4]; r=(r3+r4+r1*0.1+r2*0.1)/(1+1+0.1+0.1); R=[R,r]; Luan van end RI=[R1;R2;R3;R4]; % xanh da troi, xanh la, do, xanh nhat l2=l2+LOA; L=[l1,l2]; %plot(L,RI); k=length(L) [mini,index]=min(R); if index==0 posi=0; elseif index==k posi=LOA+LAB; else posi= L(index); end assignin('base', 'L', L) assignin('base', 'RI', RI) assignin('base', 'index', index) assignin('base', 'posi', posi) set(handles.edit_result1,'string',posi); end % bao ve bang chong set van if (cri==1)&(sl==2) LAB=min(LABt,LACt); %Chieu dai doan AB=min LAC=max(LABt,LACt); %Chieu dai doan AC=max dl=LOA/15; l1=[0:dl:LOA]; l3=[0:dl:LAC]; m=length(l1); R1=[]; R2=[]; R3=[]; R4=[];R=[]; MU1=[]; MU2=[]; MU3=[]; MU4=[]; SIGMA1=[]; SIGMA2=[]; SIGMA3=[]; SIGMA4=[]; L1=LOA/1000; L2=LAB/1000; L3=LAC/1000; for d=1:m if (d==1) la=(l1(d)+0.01)/1000; %chuyen sang he don vi km elseif (d==m) la=(l1(d)-0.01)/1000; else la=l1(d)/1000; end lb=LOA/1000-la; NODE1=[]; NODE2=[]; NODE3=[]; NODE4=[]; for i=1:n tds= Tds(i); ts= Ts(i); t=tds*1.5*1e-6; Imax= A(i); options = simset('SrcWorkspace','current'); sim('cauhinh2a1',[],options); max1=max(node1); max2=max(node2); max3=max(node3); max4=max(node4); NODE1=[NODE1 max1]; Luan van NODE2=[NODE2 max2]; NODE3=[NODE3 max3]; NODE4=[NODE4 max4]; end [mu1,sigma1] = normfit(NODE1); MU1=[MU1, mu1] ; SIGMA1=[SIGMA1, sigma1]; [mu2,sigma2] = normfit(NODE2); MU2=[MU2, mu2]; SIGMA2=[SIGMA2, sigma2]; [mu3,sigma3] = normfit(NODE3); MU3=[MU3, mu3] ; SIGMA3=[SIGMA3, sigma3]; [mu4,sigma4] = normfit(NODE4); MU4=[MU4, mu4] ; SIGMA4=[SIGMA4, sigma4]; r1=risk(mu1,sigma1,1); r2=risk(mu2,sigma2,1); r3=risk(mu3,sigma3,1); r4=risk(mu4,sigma4,1); R1=[R1,r1]; R2=[R2,r2]; R3=[R3,r3]; R4=[R4,r4]; r=(r3+r4+r1*0.1+r2*0.1)/(1+1+0.1+0.1); R=[R,r]; end [mini,index]=min(R); if index==0 posi=0; elseif index==m posi=LOA; else posi= l1(index); end assignin('base', 'index', index) assignin('base', 'posi', posi) set(handles.edit_result1,'string',posi); % -R=[]; l3=[0:dl:LAC]; m=length(l3); for d=1:m if (d==1) la1=(l3(d)+0.01)/1000; %chuyen sang he don vi km elseif (d==m) la1=(l3(d)-0.01)/1000; else la1=l3(d)/1000; end lb1=LAC/1000-la1; NODE1=[]; NODE2=[]; NODE3=[]; NODE4=[]; for i=1:n tds= Tds(i); ts= Ts(i); t=tds*1.5*1e-6; Imax= A(i); Luan van options = simset('SrcWorkspace','current'); sim('cauhinh2d1',[],options); max1=max(node1); max2=max(node2); max3=max(node3); max4=max(node4); NODE1=[NODE1 max1]; NODE2=[NODE2 max2]; NODE3=[NODE3 max3]; NODE4=[NODE4 max4]; end [mu1,sigma1] = normfit(NODE1); MU1=[MU1, mu1] ; SIGMA1=[SIGMA1, sigma1]; [mu2,sigma2] = normfit(NODE2); MU2=[MU2, mu2]; SIGMA2=[SIGMA2, sigma2]; [mu3,sigma3] = normfit(NODE3); MU3=[MU3, mu3] ; SIGMA3=[SIGMA3, sigma3]; [mu4,sigma4] = normfit(NODE4); MU4=[MU4, mu4] ; SIGMA4=[SIGMA4, sigma4]; r1=risk(mu1,sigma1,1); r2=risk(mu2,sigma2,1); r3=risk(mu3,sigma3,1); r4=risk(mu4,sigma4,1); R1=[R1,r1]; R2=[R2,r2]; R3=[R3,r3]; R4=[R4,r4]; r=(r3+r4+r1*0.1+r2*0.1)/(1+1+0.1+0.1); R=[R,r]; end % RI=[R1;R2;R3;R4]; % xanh da troi, xanh la, do, xanh nhat l3=l3+LOA; L=[l1,l3]; %plot(L,RI); % k=length(L) [mini,index]=min(R); if index==0 posi=0; elseif index==m posi=LAC; else posi= l3(index); end assignin('base', 'index', index) assignin('base', 'posi', posi) set(handles.edit_result2,'string',posi); end if cri==2 Re=get(handles.edit_R,'string'); Re=str2num(Re); % hObject handle to pushbutton_sol (see GCBO) b=readfis('F.fis'); acce= get(handles.edit_adrisk,'string'); % lay gia tri acce acce= str2num(acce); % chuyen doi acce Luan van % kiem tra thu lieu co the dung arrester hay khong %kiem tra voi truong hop chay tren doan OA % dieu kien can la tai o, A, B : (r1o,r2A,r4A,r3B)