vst_main.m

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

% VST_MAIN	Main Program of Voltage Stability Toolbox (VST).
% 		VST_MAIN displays welcome message and sets up menus. 
%

% There are four menus 1) Model, 2) Analysis, 3) Edit, and 4) Help

% =================================================================
% Clear all windows, variables, etc.
clc
close all


% =================================================================
% Set data-loaded flag: 0 - no data loaded
DataFlag=0;


% =================================================================
% Create main window, menus, and menuitems

% -------------------------------------------------
% Initialize default figure variables
position=get(0,'DefaultFigurePosition');
position=position-[20 20 0 0];
Main_Fig=figure(...
    'NumberTitle','off',...
    'Name','Voltage Stability Toolbox',...
    'Color',[0.7 0.8 0.9],...
    'Pointer','watch',...
    'Position',position,...
    'Resize','off');

% -------------------------------------------------
% MODEL MENU
Model_Menu=uimenu(Main_Fig,...
    'Label','Model');

    % Import menuitem
    import=uimenu(Model_Menu,...
        'Label','Import');

        statieee=uimenu(import,...
            'Label','IEEE Static',...
            'CallBack','imptieee');

        dynaieee=uimenu(import,...
            'Label','IEEE Dynamic',...
            'Enable','off');

    
    % Load menuitem
    open_db=uimenu(Model_Menu,...
        'Label','Load',...
        'Separator','on',...
        'Callback','getdata');


    % Build menuitem
    build_model=uimenu(Model_Menu,...
        'Label','Build',...
        'Separator','on');
% this was off before, I made it on
        nmc=uimenu(build_model,...
            'Label','Network Model C',...
            'CallBack','C_Flag=1;eqmodel(bus_type,bus_condc,bus_suscp,NumBus,tap_bus,z_bus,brch_r,brch_x,NumBranch,C_Flag);',...
            'Enable','on');
% this was off before, I made it on
        nmm=uimenu(build_model,...
            'Label','Network Model Maple',...
            'CallBack','C_Flag=0;eqmodel(bus_type,bus_condc,bus_suscp,NumBus,tap_bus,z_bus,brch_r,brch_x,NumBranch,C_Flag);',...
            'Enable','on');

        pocm=uimenu(build_model,...
           'Label','Classical PoC',...
           'CallBack','mex_poc');
    

            


    % Compile menuitem
    compile_model=uimenu(Model_Menu,...
        'Label','Compile',...
        'CallBack','mex_cmpl');

% -------------------------------------------------
% ANALYSIS MENU
Analysis_Menu=uimenu(Main_Fig,...
    'Label','Analysis');

    % Load Flow menuitem
    load_flow=uimenu(Analysis_Menu,...
        'Label','Load Flow');

        snrm=uimenu(load_flow,...
            'Label','Standard NR',...
            'CallBack','loadflow');

        %convm=uimenu(load_flow,...
           %'Label','Convergent (NRS)',...
        %'Callback','conv_loadflow');

	% Simulation menuitem
	simulate=uimenu(Analysis_Menu,...
      'Label','Simulation',...
      'CallBack','sim_0');
     
     % Static Bifurcation menuitem
    static_bif=uimenu(Analysis_Menu,...
        'Label','Static Bifurcation',...
        'Separator','on',...
        'Callback','statbif');

    % Dynamic Bifurcation menuitem
    	dynamic_bif=uimenu(Analysis_Menu,...
        'Label','Dynamic Bifurcation Analysis');
        
        z_bif=uimenu(dynamic_bif,...
        'Label','Zoom around nose point ',...
        'Callback','dynabif');
        
        low_bif=uimenu(dynamic_bif,...
        'Label','Zoom Lower Part',...
        'Callback','ldynabif');
                
        % Eigenvalues of the system matrix, Asys (which is a reduced system)
       eig_loc=uimenu(Analysis_Menu,...
          'Label','Eigenvalue of system matrix',...);
          'Callback','chlpe');
       
    % Sensitivity information obtain from left and right eigencevtor at the point of collapse point
    sensitivity=uimenu(Analysis_Menu,...
       'Label','Sensitivity around saddle node bifurcation',...);
          'Callback','eigen_vec');
    % Singular point calculations menuitems
    	sing_point=uimenu(Analysis_Menu,...
        'Label','Singular Point Analysis');
        
       % sing_NRS=uimenu(sing_point,...
          % 'Label','NRS only (upper)',...
          % 'CallBack','singpo_NRS');

        
         sing_NR_NRS=uimenu(sing_point,...
        'Label','NR and NRS ',...
        'Callback','singpo_NR_NRS');
        
        

    % -------------------------------------------------
% EDIT MENU
% Edit menuitem
    Edit_Menu=uimenu(Main_Fig,...
        'Label','Edit');
	
	% Edit VST Data Menuitem
    edit_db=uimenu(Edit_Menu,...
    	'Label','Edit VST Data',...
    	'Callback','vst_edit');
    	
%----------------------------------------------------
    	
% Help menu
helpmenu=uimenu(Main_Fig,...
    'Label','Help');

    hlp=uimenu(helpmenu,...
        'Label','Help',...
        'Separator','on',...
        'CallBack','helpfun(hlpTitle,hlpStr1,hlpStr2,hlpStr3,hlpStr4,hlpStr5)');

	abTitle=['Voltage Stability Toolbox - About VST'];
	aboutStr=[
        '                                               '
        ' Voltage Stability Toolbox  Version 2.0        '
        '                                               '
        ' Copyright (c) 1999, CEPE, Drexel University   '
        '                                               '
        ' Center for Electric Power Engineering         '
        ' Drexel University                             '
        ' Philadelphia, PA 19104                        '
        ' Phone:[phone]                          '
        ' Fax  :[phone]                          '
        '                                               '
        ' Contact Information:                          '
        ' Saffet Ayasun                                 '
        ' Email: [e-mail]             '
        '                                               '
        ' Chika Nwankpa                                 '
        ' Email: [e-mail]           '
       ];