comparesessions.m

Prony分析是获取系统振荡模式特征的一种非常有效的方法

function varargout = comparesessions(varargin)
% Implements Compare Sessions GUI
% See also: GUIDE, GUIDATA, GUIHANDLES

gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @comparesessions_OpeningFcn, ...
                   'gui_OutputFcn',  @comparesessions_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin & isstr(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before comparesessions is made visible.
function comparesessions_OpeningFcn(hObject, eventdata, handles, varargin)

% Choose default command line output for comparesessions
handles.output = hObject;
% Set the background of the figure to White
set(hObject,'Color','white');
% Set the toolbar for the figure
pronytoolbar(hObject,'on');
% Flag to check if the GUI is picking up data from file or from workspace
handles.flag_menuopen=0;
handles.idennz_mat=0;
handles.Index_Selected=0;
% Update handles structure
guidata(hObject, handles);

% UIWAIT makes comparesessions wait for user response (see UIRESUME)
% uiwait(handles.figure1);



% --- Outputs from this function are returned to the command line.
function varargout = comparesessions_OutputFcn(hObject, eventdata, handles)
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;


% --- Executes during object creation, after setting all properties.
function lstbx_session_CreateFcn(hObject, eventdata, handles)
% hObject    handle to lstbx_session (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: listbox controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc
    set(hObject,'BackgroundColor','white');
else
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end

% --- Executes on selection change in lstbx_session.
function lstbx_session_Callback(hObject, eventdata, handles)
% hObject    handle to lstbx_session (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: contents = get(hObject,'String') returns lstbx_session contents as cell array
%        contents{get(hObject,'Value')} returns selected item from lstbx_session


% --- Executes on button press in push_close.
function push_close_Callback(hObject, eventdata, handles)
% hObject    handle to push_close (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
close(gcf);

% --------------------------------------------------------------------
function menu_compareset_Callback(hObject, eventdata, handles)
% hObject    handle to menu_compareset (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)


% --------------------------------------------------------------------
function menu_open_Callback(hObject, eventdata, handles)
% hObject    handle to menu_open (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% Flag to check if the GUI is picking up data from file or from workspace

% Update handles structure

[filename, pathname] = uigetfile( ...
    {'*.cmp', 'All CMP-Files (*.cmp)'},...
    'Open Prony Comparator File');
% If "Cancel" is selected then return
if isequal([filename,pathname],[0,0])
    handles.flag_menuopen=0;
    guidata(hObject, handles);
    return
    % Otherwise load the file.
else
    cd (pathname);
    CmpData=load(eval('filename'),'-mat');
    assignin('base','CmpData',CmpData);
    handles.flag_menuopen=1;
    guidata(hObject, handles);
end
setstatus(gcbf,'Select Sessions from Listbox and Choose an Option from "Plots"');
set(handles.ppmenu_plots,'enable','on');
maxNum=evalin('base','CmpData.maxNum');
SaveData=evalin('base','CmpData.SaveData');
[promptstr]=  prepare_comparelist(maxNum,SaveData);
set(handles.lstbx_session,'String',promptstr)

% --------------------------------------------------------------------
function menu_close_Callback(hObject, eventdata, handles)
% hObject    handle to menu_close (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
close(gcf);

% --------------------------------------------------------------------
function menu_view_Callback(hObject, eventdata, handles)
% hObject    handle to menu_view (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)


% --------------------------------------------------------------------
function menu_load_Callback(hObject, eventdata, handles)
% hObject    handle to menu_load (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

setstatus(gcbf,'Select Sessions from Listbox and Choose an Option from "Plots"');
set(handles.ppmenu_plots,'enable','on');
maxNum=evalin('base','maxNum');
SaveData=evalin('base','SaveData');
[promptstr]=  prepare_comparelist(maxNum,SaveData);
set(handles.lstbx_session,'ForegroundColor','black');
set(handles.lstbx_session,'String',promptstr);



% --- Executes during object creation, after setting all properties.
function ppmenu_plots_CreateFcn(hObject, eventdata, handles)
% hObject    handle to ppmenu_plots (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: popupmenu controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc
    set(hObject,'BackgroundColor','white');
else
    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end


% --- Executes on selection change in ppmenu_plots.
function ppmenu_plots_Callback(hObject, eventdata, handles)
% hObject    handle to ppmenu_plots (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: contents = get(hObject,'String') returns ppmenu_plots contents as cell array
%        contents{get(hObject,'Value')} returns selected item from ppmenu_plots

% --- Executes on button press in push_compareplots.
function push_compareplots_Callback(hObject, eventdata, handles)
% hObject    handle to push_compareplots (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

    %evalin('base','clc');
    if(handles.flag_menuopen==1)
        SaveData=evalin('base','CmpData.SaveData');
    else
        SaveData=evalin('base','SaveData');
    end
    Index_Selected=get(handles.lstbx_session,'Value');
    handles.Index_Selected=Index_Selected;
    guidata(hObject, handles);
    % Get the column size of Index_Selected
    Run_limit=size(Index_Selected,2);
    % if nothing is selected
    if (isempty(Index_Selected))
        errordlg('Please Select the Sessions to be Compared','No Session Selected !!');
        whitebg('white'); 
    end
    
    for index=1:Run_limit
        Run=Index_Selected(index);
        mode_sizes(index)=SaveData(Run).modes;
    end
        max_modes=max(mode_sizes);
        
    
val=get(handles.ppmenu_plots,'Value');
switch val
    case 1 % poles plot
        
        for index=1:Run_limit
            Run=Index_Selected(index);
            [dp_x,dp_y,ai,a_list,SUB_IND,p,energy]= preparecomparedata(SaveData,Run);
            [inummz,idennz]=residuez(a_list(SUB_IND),p(SUB_IND),1);
            [z,p,k] = tf2zp(inummz,idennz);
            p_size=size(p,1); % Get the column size of p
            diff_size=max_modes-p_size;
            if (p_size<max_modes)
                dummy(1:diff_size,1)=-9;   % Create dummy poles of -9
                dummy_rowsize=size(dummy,1);
                if(dummy_rowsize>diff_size)
                    dummy_new=dummy(1:diff_size);
                else
                    dummy_new=dummy(1:diff_size);
                end
                p1=vertcat(p,dummy_new);
            else
                p1=p;
            end
            idennz_mat(:,index)=p1;
        end
        zplane(-9,idennz_mat,{'w^','x'});  % dummy zero at -9
        axis([-1.1 1.1 -1.1 1.1]);xlabel('Real Part'), ylabel('Imaginary Part');
        title('Poles Plot');
        legend_poles(handles,Index_Selected,Run_limit);
        handles.idennz_mat=idennz_mat;
        % Update handles structure
        guidata(hObject, handles);
        
        
    case 2 % Squared Error
        S=['m:';'b-';'y:';'r-';'g:'];
        row_S=size(S,1);
        
        for index=1:Run_limit
            colorindex=mod(index,row_S);% limit colorindex to the row size of S
            % avoid zero value of colorindex, it will happen when row_
            if (colorindex==0)
                colorindex=row_S;
            end
            Run=Index_Selected(index);
            [dp_x,dp_y,ai,a_list,SUB_IND,p,energy]= preparecomparedata(SaveData,Run);
            seerror=(dp_y(:)-ai(:)).^2;
            axes(handles.axes_pole);
            plot(dp_x,seerror,S(colorindex));
            hold on;
        end
         xlabel('Time'), ylabel('Squared Error')
         title('Squared Error Plot');
         legend_seerror(handles,Index_Selected,Run_limit);
         hold off;
            
     case 3 % Residue Plot
         S=['m:';'b-';'y:';'r-';'g:'];
         row_S=size(S,1);
         for index=1:Run_limit
            colorindex=mod(index,row_S);% limit colorindex to the row size of S
            % avoid zero value of colorindex, it will happen when row_
            if (colorindex==0)
                colorindex=row_S;
            end
            Run=Index_Selected(index);
            [dp_x,dp_y,ai,a_list,SUB_IND,p,energy]= preparecomparedata(SaveData,Run);
            SUB_N=size(SUB_IND,2);
            aa_list=abs(2*a_list(SUB_IND'));
            axes(handles.axes_pole);
            Index=1:SUB_N;
            semilogy(Index,aa_list,S(colorindex))
            if(Run_limit==1)
                axis([0 max(SUB_N)+1 0 1.1*max(aa_list)]);
            end
            hold on;
        end
        xlabel('Index'), ylabel('Residue Magnitude')
        title('Prony Residue Magnitude "Decay"')
        legend_seerror(handles,Index_Selected,Run_limit);
        hold off;
        
   case 4 % Energy Plot
         S=['m:';'b-';'y:';'r-';'g:'];
         row_S=size(S,1);