Phân tích ảnh hưởng của hiện tượng trồi sụt điện áp (flicker) tới chất lượng điện năng

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA -o0o - NGÔ THÀNH MẠNH PHÂN TÍCH ẢNH HƯỞNG CỦA HIỆN TƯNG TRỒI SỤT ĐIỆN ÁP (FLICKER) TỚI CHẤT LƯNG ĐIỆN NĂNG CHUYÊN NGÀNH : THIẾT BỊ, MẠNG & NHÀ MÁY ĐIỆN MÃ SỐ NGÀNH : 2.06.07 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH,Tháng 07 năm 2005 ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA o0o Cộng Hòa Xã Hội Chủ Nghĩa Việt Nam Độc Lập – Tự Do – Hạnh Phúc o0o NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên : Ngô Thành Mạnh Ngày tháng năm sinh : 25 - 02 - 1980 Phái : Nam Nơi Sinh : Bắc Ninh Chuyên Ngành : Thiết Bị, Mạng Và Nhà Máy Điện Mã Ngành : 2.06.07 I Tên đề tài: PHÂN TÍCH ẢNH HƯỞNG CỦA HIỆN TƯƠNG TRỒI SỤT ĐIỆN ÁP (FLICKER) ĐẾN CHẤT LƯỢNG ĐIỆN NĂNG II Nhiệm vụ nội dung: Tìm hiểu ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện Tìm hiểu lý thuyết Wavelet ứng dụng Wavelet hệ thống điện Ứng dụng phép biến đổi Wavelet vào tốn nhận dạng tín hiệu điện áp Sử dụng phần mềm ATP-EMTP để mơ chuyển dạng sóng qua Matlab Chương trình nhận dạng tượng trồi sụt điện áp III Ngày giao nhiệm vụ : 01 – 02 – 2005 IV Ngày hoàn thành : 01 – 07 – 2005 V Họ tên CB hướng dẫn : Tiến Sĩ NGUYỄN HỮU PHÚC Cán Bộ Hướng Dẫn Chủ Nhiệm Ngành Bộ Môn Quản Lý Ngành - - - Nội dung đề cương luận văn thạc sĩ hội đồng chun ngành thơng qua Phịng Đào Tạo SĐH Ngày …… tháng …… Năm 2005 Khoa Quản Lý Ngành LỜI CẢM ƠN Trước hết, xin trân trọng bày tỏ lòng biết ơn sâu sắc đến Thầy TS Nguyễn Hữu Phúc, người tận tình hướng dẫn thực hoàn thành luận văn tốt nghiệp Với kiến thức có ngày hôm nay, kết trình học tập rèn luyện lâu dài; tất công ơn tất q Thầy, Cô Trường Đại Học Bách Khoa đem đến hành trang kiến thức cho vào đời Đặc biệt xin chân thành biết ơn tới Thầy Cô Bộ môn Hệ Thống Điện tạo điều kiện thuận lợi hỗ trợ cho nhiều trình học tập, thời gian làm luận văn tốt nghiệp Cuối cùng, xin gởi lời cảm ơn chân thành đến bạn bè, đồng nghiệp gia đình tạo cho niềm tin nỗ lực cố gắng để hoàn thành luận văn TP.Hồ Chí Minh, tháng 07 năm 2005 Người thực Ngô Thành Mạnh BK TPHCM ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA -o0o - PHÂN TÍCH ẢNH HƯỞNG CỦA HIỆN TƯNG TRỒI SỤT ĐIỆN ÁP (FLICKER) TỚI CHẤT LƯNG ĐIỆN NĂNG GVHD:TS NGUYỄN HỮU PHÚC HVTH : NGƠ THÀNH MẠNH MSHV : 01803470 I ĐẶT VẤN ĐỀ • Việc nghiên cứu, phân tích nguyên nhân gây dao động điện áp hệ thống điện đóng vai trị quan trọng • Các dao động điện áp, gây ảnh hưởng đến tình trạng làm việc bình thường thiết bị hệ thống điện • Tìm hiểu đặc điểm chúng để từ nhận biết đưa phương án khắc phục LÝ DO CHỌN ĐỀ TÀI • Hiện tượng trồi sụt điện áp làm giảm chất lượng điện năng, làm cho điện áp khơng ổn định • Hiện EVN có quan tâm đến tượng trồi sụt điện áp, cụ thể nhà máy thép có cơng suất lớn đưa vào hoạt động EVN yêu cầu độ trồi sụt nhỏ 10% A1 KẾT QUẢ & MỤC ĐÍCH CỦA ĐỀ TÀI • Mục đích Luận Văn nhằm giới thiệu cơng cụ có khả mơ loại sóng điện áp tồn hệ thống điện • Nhờ việc chuyển liệu mô từ ATP-EMTP qua Matlab để thực phân tích • Wavelet vốn tính tốn mẫu liệu – cho kết Dạng sóng điện áp A3 A4 Dạng sóng điện áp • Giá trị hiệu dụng điện áp thời điểm t1: t1+TF / U (t)dt • UF(t1) = TF ∫ t1−TF / • Với TF = 20 ms ω F = 2π /TF giá trị hiệu dụng UF là: T • Ufm = Lim [ ∫ U F ( t ) dt ] T →∞ T = U* Các nguyên nhân gây tượng(Flicker) • Đóng cắt tải lớn • Khởi động động có cơng suất lớn (đặc biệt động hoạt động theo chu kỳ tuần hồn) • Dao động bất thường tải • Các thiết bị hàn • Lị hồ quang A5 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện % n_miscella = 'DWAn1d_Miscella'; ind_graph_area = 1; ind_view_mode = 2; ind_savepath = 3; nb4_stored = 3; % tag_m_savesyn = 'Save_Syn'; tag_m_savecfs = 'Save_Cfs'; tag_m_savedec = 'Save_Dec'; tag_pus_anal = 'Pus_Anal'; tag_pus_comp = 'Pus_Comp'; tag_pus_stat = 'Pus_Stat'; tag_pop_viewm = 'View_Mode'; tag_pus_dispopt = 'Pus_Options'; tag_valapp_scr = 'ValApp_Scr'; tag_txtapp_scr = 'TxtApp_Scr'; tag_valdet_scr = 'ValDet_Scr'; tag_txtdet_scr = 'TxtDet_Scr'; tag_declev = 'Pop_DecLev'; tag_txtdeclev = 'Txt_DecLev'; tag_Def_DispM = 'Default_DispM'; switch option case 'create' % Khai Bao bien toan cuc % -[Def_Btn_Height,Def_Btn_Width,Pop_Min_Width, Def_Txt_Height,X_Spacing,Y_Spacing,Def_FraBkColor ] = mextglob('get', 'Def_Btn_Height','Def_Btn_Width','Pop_Min_Width', 'Def_Txt_Height','X_Spacing','Y_Spacing','Def_FraBkColor' ); % Man hinh cua so % win_title = 'Chuong trinh nhan dang song dien ap'; [win_dw1dtool,pos_win,win_units,str_numwin, frame0,pos_frame0,Pos_Graphic_Area,pus_close] = wfigmngr('create',win_title,winAttrb,'ExtFig_Tool',mfilename,1,1,0); varargout{1} = win_dw1dtool; % Tao menu % -[m_files,m_load,m_save,m_view] = wfigmngr('getmenus',win_dw1dtool,'file','load','save','view'); m_loadsig = uimenu(m_load, 'Label','&Signal', 'Position',1, 'Callback', ['dw1dmngr(''load_sig'',' str_numwin ');'] ); m_savesyn = uimenu(m_save, 'Label','Signal', 'Position',1, Phần Phụ Lục Trang PL.2 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện 'Enable','Off', 'Tag',tag_m_savesyn, 'Callback', ['dw1dmngr(''save_synt'',' str_numwin ');'] ); % Them cac lua chon hien thi % viewmode_Names = strvcat( 'Show and Scroll', 'Full Decomposition', 'Separate Mode', 'Superimpose Mode', 'Tree Mode', 'Show and Scroll (Stem Cfs)' ); viewmode_Names = viewmode_Names(1:end,:); DW1D_Display_Mode = 2; m_disp = uimenu(m_view, 'Label','Default Display &Mode', 'Separator','On' ); nbm = size(viewmode_Names,1); dm = zeros(1,nbm); for k=1:nbm dm(k) = uimenu(m_disp, 'Label',viewmode_Names(k,:), 'Checked','Off' ); end dmstr = num2mstr(dm); for k = 1:nbm set(dm(k), 'Userdata',k, 'CallBack', ['dw1dvdrv(''default'',' str_numwin ',' dmstr ',' sprintf('%.0f',k) ');']); end dw1dvdrv('default',win_dw1dtool,dm,DW1D_Display_Mode); % Sau nhap dang song, cho xu ly % wwaiting('msg',win_dw1dtool,'Xin cho may tinh xu ly dang song vua nhap '); % Tao vi tri % -dx = X_Spacing; dx2 = 2*dx; dy = Y_Spacing; dy2 = 2*dy; d_txt = (Def_Btn_Height-Def_Txt_Height); % Vi tri cua cac khoi tren cua so % -x_frame0 = pos_frame0(1); cmd_width = pos_frame0(3); xlocINI = pos_frame0([1 3]); ybottomINI = pos_win(4)-3.5*Def_Btn_Height-dy2; y_low w_uic h_uic Phần Phụ Lục = ybottomINI-4.5*Def_Btn_Height-dy2; = 2.5*Def_Btn_Width; = 2*Def_Btn_Height; Trang PL.3 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện bdx x_left pos_anal = (cmd_width-w_uic)/2; = x_frame0+bdx; = [x_left, y_low, w_uic, h_uic]; x_left y_low push_width pos_stat pos_comp pos_comp(1) y_low = = = = = = = x_frame0+dx2; pos_anal(2)-h_uic-2*dy2; (cmd_width-3*dx2)/2; [x_left, y_low, push_width, h_uic]; pos_stat; pos_stat(1)+pos_stat(3)+dx2; 10*y_low-h_uic-dy; pos_viewm_txt Def_Txt_Height]; pop_width x_left y_low pos_pop_viewm = [x_left, y_low, 2*Def_Btn_Width, = = = = cmd_width-4*dx2; x_frame0+(cmd_width-pop_width)/2; pos_viewm_txt(2)-1.5*Def_Btn_Height-dy; [x_left, y_low, pop_width, Def_Btn_Height]; pos_pus_dispopt = [x_left, y_low, pop_width, Def_Btn_Height]; y_low x_left w_txt pos_txtapp_scr xl pos_valapp_scr xl pos_txtdet_scr xl pos_valdet_scr y_low pos_txtlev_dec Def_Txt_Height]; xl pos_declev = = = = = = = = = = pos_pop_viewm(2)-Def_Btn_Height-dy2; pos_frame0(1)+dx2; Pop_Min_Width/2+dx/4; [x_left,y_low+d_txt/2, w_txt, Def_Txt_Height]; pos_txtapp_scr(1)+pos_txtapp_scr(3)+dx; [xl, y_low, Pop_Min_Width, Def_Btn_Height]; pos_valapp_scr(1)+pos_valapp_scr(3)+dx2; [xl, y_low+d_txt/2, w_txt, Def_Txt_Height]; pos_txtdet_scr(1)+pos_txtdet_scr(3)+dx; [xl, y_low, Pop_Min_Width, Def_Btn_Height]; = pos_pop_viewm(2)-Def_Btn_Height-dy2; = [x_left, y_low+d_txt/2, Def_Btn_Width, = pos_txtlev_dec(1)+Def_Btn_Width+dx2; = [xl, y_low, Pop_Min_Width, Def_Btn_Height]; % -str_levels = int2str([1:max_lev_anal]'); str_anal = 'Analyze'; str_stat = 'Statistics'; str_comp = 'Compress'; str_viewm_txt = 'Display mode : '; str_pop_viewm = viewmode_Names; str_pus_dispopt = 'More Display Options'; str_txtapp_scr = 'App.'; str_valapp_scr = str_levels; str_txtdet_scr = 'Det.'; str_valdet_scr = str_levels; str_txtlev_dec = 'at level'; str_declev = str_levels; % -cba_anal = ['dw1dmngr(''anal'',' str_numwin ');']; cba_stat = ['dw1dmngr(''stat'',' str_numwin ');']; cba_comp = ['dw1dmngr(''comp'',' str_numwin ');']; cba_viewm = ['dw1dvmod(''ch_vm'',' str_numwin ');']; Phần Phụ Lục Trang PL.4 Luận văn tốt nghiệp cba_mdisp cba_app_scr cba_det_scr cba_declev Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện = = = = ['dw1ddisp(''create'',' str_numwin ');']; ['dw1dtool(''app'',' str_numwin ');']; ['dw1dtool(''det'',' str_numwin ');']; ['dw1dtool(''dec'',' str_numwin ');']; % Du lieu, bien doi Wavelet % -utanapar('create',win_dw1dtool, 'xloc',xlocINI,'bottom',ybottomINI, 'enable','off', 'wtype','dwt' ); commonProp = { 'Parent',win_dw1dtool, 'Style','Pushbutton', 'Unit',win_units, 'Enable','off' }; pus_anal = uicontrol(commonProp{:}, 'Position',pos_anal, 'String',xlate(str_anal), 'Tag',tag_pus_anal, 'Callback',cba_anal, 'Interruptible','On' ); pus_stat = uicontrol(commonProp{:}, 'Position',pos_stat, 'String',xlate(str_stat), 'Tag',tag_pus_stat, 'Callback',cba_stat ); pus_comp = uicontrol(commonProp{:}, 'Position',pos_comp, 'String',xlate(str_comp), 'Tag',tag_pus_comp, 'Callback',cba_comp ); txt_viewm = uicontrol( 'Parent',win_dw1dtool, 'Style','text', 'Unit',win_units, 'HorizontalAlignment','left', 'Position',pos_viewm_txt, 'String',str_viewm_txt, 'Backgroundcolor',Def_FraBkColor ); view_mode = 1; pop_viewm = uicontrol( 'Parent',win_dw1dtool, 'Style','Popup', 'Unit',win_units, 'Position',pos_pop_viewm, 'String',str_pop_viewm, 'Value',view_mode, 'UserData',view_mode, 'Enable','off', 'Callback',cba_viewm, Phần Phụ Lục Trang PL.5 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện 'Tag',tag_pop_viewm ); pus_dispopt = uicontrol( 'Parent',win_dw1dtool, 'Style','Pushbutton', 'Unit',win_units, 'Position',pos_pus_dispopt, 'String',xlate(str_pus_dispopt), 'Visible','off', 'Enable','off', 'Tag',tag_pus_dispopt, 'Interruptible','On', 'Callback',cba_mdisp ); txt_txtapp_scr = uicontrol( 'Parent',win_dw1dtool, 'Style','Text', 'Unit',win_units, 'Position',pos_txtapp_scr, 'Visible','Off', 'String',str_txtapp_scr, 'Backgroundcolor',Def_FraBkColor, 'Tag',tag_txtapp_scr ); pop_valapp_scr = uicontrol( 'Parent',win_dw1dtool, 'Style','Popup', 'Unit',win_units, 'Position',pos_valapp_scr, 'Visible','Off', 'String',str_valapp_scr, 'Enable','off', 'Tag',tag_valapp_scr, 'Callback',cba_app_scr ); txt_txtdet_scr = uicontrol( 'Parent',win_dw1dtool, 'Style','Text', Phần Phụ Lục Trang PL.6 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện 'Unit',win_units, 'Position',pos_txtdet_scr, 'Visible','Off', 'String',str_txtdet_scr, 'Backgroundcolor',Def_FraBkColor, 'Tag',tag_txtdet_scr ); pop_valdet_scr = uicontrol( 'Parent',win_dw1dtool, 'Style','Popup', 'Unit',win_units, 'Position',pos_valdet_scr, 'Visible','Off', 'String',str_valdet_scr, 'Enable','off', 'Tag',tag_valdet_scr, 'Callback',cba_det_scr ); txt_txtlev_dec = uicontrol( 'Parent',win_dw1dtool, 'Style','Text', 'Unit',win_units, 'Position',pos_txtlev_dec, 'Visible','Off', 'String',str_txtlev_dec, 'Backgroundcolor',Def_FraBkColor, 'Tag',tag_txtdeclev ); pop_declev = uicontrol( 'Parent',win_dw1dtool, 'Style','Popup', 'Unit',win_units, 'Position',pos_declev, 'Visible','Off', 'String',str_declev, 'Enable','off', 'Tag',tag_declev, Phần Phụ Lục Trang PL.7 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện 'Callback',cba_declev ); % utcolmap('create',win_dw1dtool, 'xloc',xlocINI, 'enable','off', 'briflag',0, 'bkcolor',Def_FraBkColor); % -pos_txt1 = get(txt_txtapp_scr,'Position'); pos_opt = get(pus_dispopt,'Position'); pos_opt(2) = pos_txt1(2)-Def_Btn_Height-2*dy2; set(pus_dispopt,'Position',pos_opt); set([txt_txtapp_scr, txt_txtdet_scr, pop_valapp_scr, pop_valdet_scr, pus_dispopt],'Visible','on'); drawnow % Giup (CSHelp) % hdl_DW1D_DISPMODE = [txt_viewm,pop_viewm]; hdl_DW1D_MOREDISP = [pus_dispopt]; wfighelp('add_ContextMenu',win_dw1dtool,hdl_DW1D_DISPMODE,'DW1D_DISPMODE' ); wfighelp('add_ContextMenu',win_dw1dtool,hdl_DW1D_MOREDISP,'DW1D_MOREDISP' ); % - % -Pos_Graphic_Area = wfigmngr('normalize',win_dw1dtool,Pos_Graphic_Area); % Cac du lieu duoc cap nhat % -utanapar('set_cba_num',win_dw1dtool,[m_files;pus_anal]); % Cac gia tri luu giu % -wmemtool('ini',win_dw1dtool,n_InfoInit,nb0_stored); wmemtool('ini',win_dw1dtool,n_param_anal,nb1_stored); wmemtool('ini',win_dw1dtool,n_coefs_longs,nb2_stored); wmemtool('ini',win_dw1dtool,n_synt_sig,nb3_stored); wmemtool('ini',win_dw1dtool,n_miscella,nb4_stored); wmemtool('wmb',win_dw1dtool,n_miscella, ind_graph_area,Pos_Graphic_Area); % Tao cac dinh dang mau % -cbcolmap('set',win_dw1dtool,'pal',{'pink',default_nbcolors}); % dw1dvmod('ini_vm',win_dw1dtool); Phần Phụ Lục Trang PL.8 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện % Ket thuc viec cho doi % wwaiting('off',win_dw1dtool); case {'app','det'} fig = varargin{1}; pop_viewm = findobj(get(fig,'Children'),'flat','tag',tag_pop_viewm); view_mode = get(pop_viewm,'Value'); switch view_mode case , dw1dscrm(option,fig); case , dw1dcfsm(option,fig); end case 'dec' fig = varargin{1}; dw1ddecm('dec',fig); case 'close' % -fig = varargin{1}; called_win = wfindobj('figure','Userdata',fig); delete(called_win); set(fig,'Visible','Off') drawnow feval('dw1dfile','del',fig); otherwise errargt(mfilename,'Unknown Option','msg'); error('*'); end Giao diện chương trình function WindowOption(Extra,State) if nargin < Extra = 1; end if nargin < State = 'mucnangluong'; end scrsz = get(0,'ScreenSize' ); Phase_Test = figure( 'menubar', 'none','NumberTitle','off','name', 'Main Program', 'color' , [0.7 0.7 0.7], 'Position' ,[scrsz(1) scrsz(2) scrsz(3) scrsz(4)-30],'Resize','on' ); File = uimenu('label', 'File'); New = uimenu(File, 'label' , 'New','Callback' ,'mucnangluong(''userdata'',''initialise'')'); Quit = uimenu(File, 'label', 'Quit','callback', 'close all'); if Extra == ExtraMenu = uimenu('label','CHUONG TRINH'); ExtraFull = uimenu(ExtraMenu, 'label', 'NHAP DANG SONG','callback' , 'mucnangluong'); end Phần Phụ Lục Trang PL.9 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện Help = uimenu('label', 'GIUP DO'); P_Help = uimenu(Help, 'label', 'Phase Test Help','callback', 'helpwin Voltage'); clear all CHƯƠNG TRÌNH NHẬN DẠNG VÀ PHÂN BIỆT HIỆN TƯỢNG TRỒI SỤT ĐIỆN ÁP VỚI CÁC HIỆN TƯỢNG DAO ĐỘNG ĐIỆN ÁP KHÁC [filename, pathname]=uigetfile('D:\LVTN\Nhan dang\'); if filename~=0 t=filename; load(t) s=vU; u=1 x=abs(max(s)); y=abs(min(s)); end [C,L] = wavedec(s,13,'db4'); for i=1:13 P{i} = wrcoef('d',C,L,'db4',i); E(i)=((sum(P{i}.^2))/length(P{i}))^(1/2); end figure(1); subplot(4,1,1);plot(s);axis([1 5001 min(s) max(s)]); title('CHƯƠNG TRÌNH NHẬN DẠNG HIỆN TƯNG FLICKER','FontName','VNITimes'); subplot(4,1,2);plot(abs(P{1}));axis([1 5001 max(abs(P{1}))]); pause subplot(4,1,3);plot(abs(P{2}));axis([1 5001 max(abs(P{2}))]); pause subplot(4,1,4);plot(abs(P{3}));axis([1 5001 max(abs(P{3}))]); pause figure(2);bar(E);axis([1 13 ceil(max(E*10))/10]); pause title('CÁC MỨC NĂNG LƯNG','FontName','VNI-Times'); plot (E, 'DisplayName', 'song dua vao', 'YDataSource', 'dang song dua vao'); figure(gcf) pause m=(E(9)+E(10)+E(11)); n=(E(4)+E(5)+E(6)); if m>n & x1.04*u; String='DANG SONG BAN DUA VAO CO HIEN TUONG FLICKER !'; else String='DANG SONG BAN DUA VAO KHONG PHAI HIEN TUONG FLICKER ! !'; end Title='Identify Message'; uiwait(msgbox(String,Title,'warn')); break close Phần Phụ Lục Trang PL.10 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện B MÃ NGUỒN MATLAB CÁC VÍ DỤ MINH HỌA TRONG CHƯƠNG Chương trình tính FT cho tín hiệu dừng hình 2.1 t = 0:0.000127:0.127; x = sin(2*pi*200*t) + sin(2*pi*300*t)+ sin(2*pi*400*t)+ sin(2*pi*600*t); y = fft(x,2048); m = abs(y); subplot(2,1,1); plot((0:1000),x);title('Orgin Signal x(t)') ylabel('Amplitude'); grid on xlabel('Time [Samples]') f=(0:255)*(0.000127)^(-1)/2048; subplot(2,1,2); plot(f,m(1:256));title('Fourier Transform of x(t)') ylabel('Amplitude'); grid on xlabel('Frequency [Hertz]') pause close clear all Chương trình tính FT cho tín hiệu không dừng hình 2.2 deltat=0.0001; t1=0:deltat:0.04-deltat; t2=0.04:deltat:0.08-deltat; t3=0.08:deltat:0.12-deltat; t4=0.12:deltat:0.16-deltat; x1 = sin(2*pi*200*t1); x2 = sin(2*pi*300*t2); x3 = sin(2*pi*400*t3); x4 = sin(2*pi*600*t4); x = [x1 x2 x3 x4]; length(x) y = fft(x,2048); m = abs(y); subplot(2,1,1); plot((0:1599),x);title('Orgin Signal x(t)'); ylabel('Amplitude'); grid on xlabel('Time [Samples]') f=(0:184)*(deltat)^(-1)/2048; subplot(2,1,2); plot(f,m(1:185));title('Fourier Transform of x(t)'); ylabel('Amplitude'); grid on xlabel('Frequency [Hertz]') pause close clear all Chương trình tính FT cho tín hiệu xung Dirac deltat=0.001; k=1; for t=-1:deltat:1 if t==0 x(k) = 1; else x(k)=0; end k=k+1; end Phần Phụ Lục Trang PL.11 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện y = fft(x,2048); m = abs(y); subplot(2,1,1); plot((-1000:1000),x,'r');title('Signal \delta(t)'); axis([-1000 1000 1.3]) ylabel('Amplitude'); grid on xlabel('Time [Samples]') f=(-200:200)*(deltat)^(-1)/2048; subplot(2,1,2); plot(f*2*pi,m(1:401),'r');title('Fourier Transform of \delta(t)');axis([-600 600 1.3]) ylabel('Amplitude'); grid on xlabel('Frequency [Radians]') pause close clear all Chương trình tính STFT cho tín hiệu không dừng hình 2.7 deltat=0.0006; t1=0:deltat:0.05-deltat; t2=0.05:deltat:0.1-deltat; t3=0.1:deltat:0.15-deltat; t4=0.15:deltat:0.2-deltat; x1 = sin(2*pi*200*t1); x2 = sin(2*pi*300*t2); x3 = sin(2*pi*400*t3); x4 = sin(2*pi*600*t4); t = [t1 t2 t3 t4]; x = [x1 x2 x3 x4]; subplot(2,2,1); plot(t,x);hold on; length_window=20; window=kaiser(length_window,5); plot(t(1:length_window),window,'r'); title('Orgin Signal & Narrow Window'); axis([0 0.2 -1 1]) ylabel('Amplitude'); xlabel('Time [seconds]') subplot(2,2,2);specgram(x,512,(deltat)^(1),kaiser(length_window,5),16);colorbar title('Short-Time Fourier Transform') subplot(2,2,3); plot(t,x);hold on; length_window=100; window=kaiser(length_window,5); plot(t(1:length_window),window,'r'); title('Orgin Signal & Wide Window'); axis([0 0.2 -1 1]) ylabel('Amplitude'); xlabel('Time [seconds]') subplot(2,2,4);specgram(x,512,(deltat)^(1),kaiser(length_window,5),90);colorbar title('Short-Time Fourier Transform') pause close clear all Chương trình tính CWT cho tín hiệu không dừng hình 2.10 % Set wavelet name,interval and number of samples wname = 'morl'; A = 0; B = 64; P = 500000; clc Phần Phụ Lục Trang PL.12 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện % Compute the sampling period and the sampled function, and the true frequencies tmin=0; tmax=0.2;P_new=round((tmax-tmin)*P/(B-A)); delta = (B-A)/(P-1); t=linspace(tmin,tmax,P_new-1); tt=[1 round(P_new/4) round(P_new/2) round(3/4*P_new) P_new]; tab_OMEGA = 2*pi*[200,300,400,600]; tab_FREQ = tab_OMEGA/(2*pi); tab_COEFS = [1,1,1,1]; x=[]; for k = 1:length(tab_OMEGA); x = [x tab_COEFS(k)*sin(tab_OMEGA(k)*t(tt(k):(tt(k+1)-1)))]; end subplot(221); plot(t,x);title('Orgin Signal'); axis([0 0.2 -1 1]) ylabel('Amplitude'); xlabel('Time [seconds]') % Set scales and use scal2frq to compute the array of pseudo-frequencies scales = [60:-0.5:1]; tab_PF = scal2frq(scales,wname,delta); % Compute the nearest pseudo-frequencies and the corresponding scales for k=1:length(tab_OMEGA) [dummy,ind] = min(abs(tab_PF-tab_FREQ(k))); PF_app(k) = tab_PF(ind); SC_app(k) = scales(ind); end subplot(222);cwt(x,scales,wname,'plot');colorbar title('Abs Values of Ca,b Coefficients'); xlabel('Space [samples]') ylabel('Inverse Scales ') pause close % Continuous analysis and plot str1 = strvcat(' Signal x(t) = sin(2\pi.200t)+sin(2\pi.200t)+sin(2\pi.300t)+sin(2\pi.400t)+sin(2\pi.600t )', [' True frequencies (in Hz): = [' num2str(tab_FREQ,4) ']' ]); str2 = ['Array of pseudo-frequencies and scales: ']; str3 = [num2str([tab_PF',scales'],4)]; str4 = ['Pseudo-frequencies = ' num2str(PF_app,4)]; str5 = ['Corresponding scales = ' num2str(SC_app,4)]; figure; cwt(x,scales,wname,'plot'); ax = gca; colorbar axTITL = get(ax,'title'); axXLAB = get(ax,'xlabel'); set(axTITL,'String',str1) set(axXLAB,'String',strvcat(str4, str5)) pause clc; disp(strvcat(' ',str1,' ',str2,str3,' ',str4,str5)) close clear all Chương trình tính FT, STFT CWT cho tín hiệu Dirac hình 2.12 % Set wavelet name,interval and number of samples wname = 'morl';A = 0; B = 64; P = 500000; % Compute the sampling period and the sampled function, and the true frequencies Phần Phụ Lục Trang PL.13 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện tmin=0; tmax=0.2;P_new=round((tmax-tmin)*P/(B-A));delta = (B-A)/(P-1); t=linspace(tmin,tmax,P_new); k=1; for tt=1:P_new if tt==round(P_new/2) x(k) = 1; else x(k)=0; end k=k+1; end y = fft(x,2048); m = abs(y); subplot(221); plot(t,x); colorbar title('Signal \delta(t) with 1563 samples'); axis([tmin tmax 1.3]) ylabel('Amplitude'); xlabel('Time [Seconds]') f=(-200:200)*(delta)^(-1)/2048; subplot(222); plot(f*2*pi,m(1:401));colorbar title('Fourier Transform of \delta(t)'); axis([-600 600 1.3]) ylabel('Amplitude'); xlabel('Frequency [Radians]') length_window=20; window=kaiser(length_window,5); subplot(223);specgram(x,512,(delta)^(1),kaiser(length_window,5),16);colorbar title('Short-Time Fourier Transform') scales = [30:-5:1]; clc;%tab_PF = scal2frq(scales,wname,delta); subplot(224);cwt(x,scales,wname,'plot');colorbar title('Abs Values of Ca,b Coefficients'); xlabel('Space [samples]') ylabel('Inverse Scales ') pause close clear all 10 Chương trình tính FT, STFT CWT cho tín hiệu không dừng hình 13 % Set wavelet name,interval and number of samples wname = 'morl';A = 0; B = 64; P = 500000; % Compute the sampling period and the sampled function, and the true frequencies tmin=0; tmax=0.2;P_new=round((tmax-tmin)*P/(B-A));delta = (B-A)/(P-1); t=linspace(tmin,tmax,P_new-1); tt=[1 round(P_new/4) round(P_new/2) round(3/4*P_new) P_new]; tab_OMEGA = 2*pi*[200,300,400,600]; tab_FREQ = tab_OMEGA/(2*pi); tab_COEFS = [1,1,1,1]; x=[]; for k = 1:length(tab_OMEGA); x = [x tab_COEFS(k)*sin(tab_OMEGA(k)*t(tt(k):(tt(k+1)-1)))]; end Phần Phụ Lục Trang PL.14 Luận văn tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện subplot(221); plot(t,x); colorbar; axis([tmin tmax -1 1]) title('Sine Signal includes freq - 1563 samples'); ylabel('Amplitude'); xlabel('Time [Seconds]') y = fft(x,2048); m = abs(y); f=(0:180)*(delta)^(-1)/2048; subplot(222); plot(f,m(1:181));colorbar title('Fourier Transform of \delta(t)'); axis([0 800 max(m)]) ylabel('Amplitude'); xlabel('Frequency [Radians]') length_window=50; window=kaiser(length_window,5); subplot(223);specgram(x,512,(delta)^(1),kaiser(length_window,5),40);colorbar axis([0 tmax 800]) title('Short-Time Fourier Transform') scales = [35:-0.5:10]; tab_PF = scal2frq(scales,wname,delta); for k=1:length(tab_OMEGA) [dummy,ind] = min(abs(tab_PF-tab_FREQ(k))); PF_app(k) = tab_PF(ind); SC_app(k) = scales(ind); end subplot(224);cwt(x,scales,wname,'plot');colorbar title('Abs Values of Ca,b Coefficients'); xlabel('Space [samples]') ylabel('Inverse Scales ') pause close str1 = strvcat(' Signal x(t) = sin(2\pi.200t)+sin(2\pi.200t)+sin(2\pi.300t)+sin(2\pi.400t)+sin(2\pi.600t )', [' True frequencies (in Hz): = [' num2str(tab_FREQ,4) ']' ]); str2 = ['Array of pseudo-frequencies and scales: ']; str3 = [num2str([tab_PF',scales'],4)]; str4 = ['Pseudo-frequencies = ' num2str(PF_app,4)]; str5 = ['Corresponding scales = ' num2str(SC_app,4)]; figure; cwt(x,scales,wname,'plot'); ax = gca; colorbar axTITL = get(ax,'title'); axXLAB = get(ax,'xlabel'); set(axTITL,'String',str1) set(axXLAB,'String',strvcat(str4, str5)) clc; disp(strvcat(' ',str1,' ',str2,str3,' ',str4,str5)) pause close clear all Phần Phụ Lục Trang PL.15 TĨM TẮT LÝ LỊCH TRÍCH NGANG Họ tên : Ngô Thành Mạnh Ngày tháng năm sinh : 25 - 02– 1980 Địa liên lạc Vấp Phái : Nam Nơi Sinh : Bắc Ninh : C34, Căn 26, Đường 26-3, Phường 17, Quận Gị Q trình đào tạo : - Từ 1997 – 2002, Học trường ĐH Bách Khoa Tp.HCM - Từ 2003 đến nay, theo học cao học ngành Thiết Bị, Mạng Nhà Máy Điện Khóa 14, Trường ĐH Bách Khoa Tp.HCM Quá trình cơng tác : - Từ 2002, cơng tác Ban Quản Lý Dự Án Khí Điện Đạm Cà Mau - Từ 2004, công tác Trung Tâm Điều Độ Hệ Thống Điện Miền Nam ... ẢNH HƯỞNG CỦA HIỆN TƯỢNG TRỒI SỤT ĐIỆN ÁP TỚI CHẤT LƯỢNG ĐIỆN NĂNG I TÌM HIỂU ẢNH HƯỞNG CỦA HIỆN TƯỢNG TRỒI SỤT ĐIỆN ÁP TỚI CHẤT LƯỢNG ĐIỆN NĂNG 1.1 1.1 Khái niệm tượng trồi sụt điện áp. .. tốt nghiệp Phân tích ảnh hưởng tượng trồi sụt điện áp tới chất lượng điện Chương ẢNH HƯỞNG CỦA HIỆN TƯỢNG TRỒI SỤT ĐIỆN ÁP TỚI CHẤT LƯỢNG ĐIỆN NĂNG Có thể nói, việc nghiên cứu, phân tích nguyên... Máy Điện Mã Ngành : 2.06.07 I Tên đề tài: PHÂN TÍCH ẢNH HƯỞNG CỦA HIỆN TƯƠNG TRỒI SỤT ĐIỆN ÁP (FLICKER) ĐẾN CHẤT LƯỢNG ĐIỆN NĂNG II Nhiệm vụ nội dung: Tìm hiểu ảnh hưởng tượng trồi sụt điện áp tới