% 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 =
