ldynabif.m

一种求解目标函数最小化的MATLAB仿真程序.

% Dynamic Bifurcation Analysis Program (dynabif.m)
% LDYNABIF  up GUI for dynamic bifurcation analysis.
% This M-file is designed to obtain the lower part of the nose curve
% It allows users to specify search direction, tolerance, number of steps etc..

% ****************************************************************************
if exist('dlgfig')
    set(dlgfig,'Visible','off');
end

% Check whether data is loaded or not
%************************************

if ~exist('DataFlag'),DataFlag=0;end
if DataFlag==0,getdata;end;
if DataFlag==0,return;end;
if (strcmp(computer,'SOL2' ))
  mexwild = '*.mexsol';
  lensub=7;
elseif (strcmp(computer, 'PCWIN'))
   mexwild = '*.dll';
    lensub=4;
end

[LFfilename,LFpathname]=uigetfile(mexwild,'choose system for analysis');


% Check the correct .dll file
%********************************

if LFfilename~=0,
    CurrentSystem=LFfilename(1:length(LFfilename)-lensub);
    if isempty(findstr(CurrentFileName,CurrentSystem))
        titleStr=[' Inconsistent Names '];
        textStr=['                                              '
                 ' Names of the data file and executable must   '
                 ' be consistent. Do you want to choose another '
                 ' executable now ?                             '
                ];
        dlgfig=mydlg(titleStr,textStr,'dynabif | rtnnow')
        return;
    end
   
    % Define a figure window and assign a handle to it.
    % **************************************************
    
    position=get(0,'DefaultFigurePosition');
    position=position-[70 70 0 0];
    dynamic_fig=figure(...
        'NumberTitle','off',...
        'Name','Voltage Stability Toolbox - Dynamic Bifurcation Analysis',...
        'Resize','off',...
        'Position',position,...
        'Color',[0.7 0.8 0.9]);

%********************************************************
% Help menu
% *******************************************************

helpmenu=uimenu(dynamic_fig,...
    'Label','Help');

    hlp=uimenu(helpmenu,...
        'Label','Help',...
        'CallBack','helpfun(ldb_hlpTitle,hlpStr5)');

ldb_hlpTitle=['Voltage Stability Toolbox - Dynamic Bifurcation Analysis Help Window'];
            
    % Add the LFpathname to the current path
    % *****************************************
    
    path(path,LFpathname);
    
    % Load the CurrentSystem (IEEE#.mat) mat file
    %*********************************************
     
    eval(['load ',CurrentSystem,'.mat']);
    
    % Add the title to the dynamic bifurcation analysis window (i.e.,Current System:IEEE#)
    %************************************************************************************
    
    titlepanel;
         
    % Present editable data: states & params
    % ***************************************
    
    k_states=length(x);		% x's are the states(i.e.,angle and voltage magnitude) 
    k_params=length(param);	% param is the parameters(i.e., net injections at each bus)
    
    %Label the states (i.e., delta1, V2, etc...)
    %********************************************
    
    lfmklabl; %Put the initial data and states values from the input


    % Default Search Direction
    %***************************
    
    p=zeros(size(param));
    
    % Add the lower panel to the current simulation window
    % *****************************************************
    
    lowerpanel;

    ParamSlider=uicontrol(gcf,...
        'Style','popupmenu',...
        'Position',[display_length*45-30,106,150,16],...
        'BackgroundColor','Y',...
        'String','Parameter Values|Search Direction',...
        'Callback','lfpardis');

    %*********************************
    % Load Flow Controls
    %*********************************
    % Tolerence
    
    LFAbsTol=.000001;
    LFRelTol=.0001;

    % *************************
    % NR and NRS Step Control	
    % *************************
    % NR and NRS Steps
    
    NR_steps=100;
    NRS_Steps=100;
    % *************************
    % pop-up algortithm selector
    % *************************
    % Choose either Newton-Raphson-Seydel or Newton-Raphson algorithm
    % ****************************************************************
    %lhpfsrf implements continuation mehod and lfcomp runs load flow once for the given parameter
    LF_Alg_List=['lhpfsrf';'lfcomp '];
    LF_Alg=LF_Alg_List(1,1:7);
    
    % set up the upper panel in the current sumilation window
    % *********************************************************
    upperpanel1;

    % Set up the algorithm display
    % ********************************
    
    alg_popup=uicontrol(gcf,...
        'Style','popupmenu',...
        'Position',[30,310,223,18],...
        'HorizontalAlignment','Center',...
        'String','Newton-Raphson-Seydel|Newton-Raphson',...
        'Callback','LF_Alg=LF_Alg_List(get(alg_popup,''Value''),1:7);');
end;