电能仿真基于matlab电能质量谐波波形仿真【含Matlab源码 2263期】.zip

% Power Quality Harmonic Interactive Demo % This interactive demo generates power quality harmonic waveform up to % 15th orders based on numerical modeling approach. % This demo also visualize the harmonic in frequency domain and computes % the its Total Harmonic Distortion. % It allows user to set the magnitude and phase for each harmonic order. % The harmonic waveform can then save as .mat file format. function PQHarmonicDemoMain figure('NumberTitle','off','Toolbar','none',... 'Name','Power Quality Harmonic Interactive Demo'); hTD = zeros(1,320); h = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; hlast = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; ph = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; phlast = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; subplot(3,1,1); subplot(3,1,2); uicontrol('Style', 'text', 'String', '幅值 (pu)',... 'Position', [80 90 80 20],'horizontalAlignment', 'left'); uicontrol('Style', 'text', 'String', '相位 (deg)',... 'Position', [80 60 80 20],'horizontalAlignment', 'left'); uicontrol('Style', 'text', 'String', '谐波阶数',... 'Position', [80 30 80 20],'horizontalAlignment', 'left'); hMagTxt = uicontrol('Style','edit','string','1','Enable','inactive',... 'Position', [170 90 40 20],'BackgroundColor','w',... 'horizontalAlignment','left'); hMagSldr = uicontrol('Style','slider','Min',0,'Max',1 ,'Value',1,... 'SliderStep', [0.01, 0.1],... 'Position',[205 90 15 20],'Callback',@harmonic); hPhTxt = uicontrol('Style','edit','string','0','Enable','inactive',... 'Position', [170 60 40 20],'BackgroundColor','w',... 'horizontalAlignment','left'); hPhSldr = uicontrol('Style','slider','Min',0,'Max',180 ,'Value',0,... 'SliderStep', [1/18, 0.1],... 'Position',[205 60 15 20],'Callback',@harmonic); hOrder=uicontrol('Style', 'popup', 'String',1:16,... 'Position', [170 30 50 20],'Callback',@OrderSelect,... 'horizontalAlignment', 'left'); uicontrol('Style', 'text', 'String', '系统频率 (Hz)',... 'Position', [250 90 100 20],'horizontalAlignment', 'left'); uicontrol('Style', 'text', 'String', '采样频率 (Hz)',... 'Position', [250 60 100 20],'horizontalAlignment', 'left'); hFreqSelect=uicontrol('Style', 'popup', 'String',{'50','60'},'Value',1,... 'Position', [350 90 50 20],'Callback',@harmonic,... 'horizontalAlignment', 'left'); hFsTxt = uicontrol('Style','edit','string','1600','Enable','inactive',... 'Position', [350 60 50 20],'BackgroundColor','w',... 'horizontalAlignment','left'); hPreset=uicontrol('Style', 'popup', 'String',... {'纯正弦波','1ph整流器谐波'},'Value',1,... 'Position', [250 20 150 30],'Callback',@Preset,... 'horizontalAlignment', 'left'); uicontrol('Style', 'pushbutton', 'String', '导出数据',... 'Position', [430 30 80 80],'Callback',@Export); harmonic; function Export(~,~) uisave('hTD','harmonic'); end function Preset(~,~) switch get(hPreset,'Value') case 1 h = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; hlast = [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; ph = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; phlast = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; case 2 h = [1,0,0.8,0,0.6,0,0.4,0,0.2,0,0.1,0,0,0,0]; hlast = [1,0,0.8,0,0.6,0,0.4,0,0.2,0,0.1,0,0,0,0]; ph = [0,0,180,0,0,0,180,0,0,0,180,0,0,0,0]; phlast = [0,0,180,0,0,0,180,0,0,0,180,0,0,0,0]; end OrderSelect; harmonic; end function OrderSelect(~,~) switch get(hOrder,'Value') case 1 set(hMagSldr,'Value',hlast(1)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(1)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 2 set(hMagSldr,'Value',hlast(2)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(2)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 3 set(hMagSldr,'Value',hlast(3)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(3)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 4 set(hMagSldr,'Value',hlast(4)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(4)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 5 set(hMagSldr,'Value',hlast(5)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(5)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 6 set(hMagSldr,'Value',hlast(6)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(6)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 7 set(hMagSldr,'Value',hlast(7)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(7)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 8 set(hMagSldr,'Value',hlast(8)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(8)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 9 set(hMagSldr,'Value',hlast(9)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(9)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 10 set(hMagSldr,'Value',hlast(10)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(10)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 11 set(hMagSldr,'Value',hlast(11)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(11)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 12 set(hMagSldr,'Value',hlast(12)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(12)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 13 set(hMagSldr,'Value',hlast(13)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(13)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 14 set(hMagSldr,'Value',hlast(14)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(14)); set(hPhTxt,'String',get(hPhSldr,'Value')); case 15 set(hMagSldr,'Value',hlast(15)); set(hMagTxt,'String',get(hMagSldr,'Value')); set(hPhSldr,'Value',phlast(15)); set(hPhTxt,'String',get(hPhSldr,'Value')); end end function harmonic(~,~) if get(hFreqSelect,'Value')==1 f = 50; fs = 1600; t = 1/fs:1/fs:0.2; set(hFsTxt,'String','1600'); else f =