performprony.m

matlab 非平稳信号分析

legend('Prony Sub-Approximate','Measured');
ylabel('Signal'),xlabel('Time');


% --------------------------------------------------------------------
function plotnew_frequency_modesfit_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_frequency_modesfit (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
IIGuiData = prepare_guidata(handles);
figure('Color','White');
plot(IIGuiData.f2,IIGuiData.fft_iapp2,'r',IIGuiData.f2,IIGuiData.fft_data2,'b');
title(sprintf('Prony Approximation of Model Order %d, Picking %d Residues in Frequency Domain',IIGuiData.NewVal,IIGuiData.SUB_N));
legend('Prony Sub-Approximate','Measured');% place it at lower right corner of the axis
xlabel('Frequency');
ylabel('Magnitude (dB)');

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


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


% --------------------------------------------------------------------
function plotnew_serror_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_serror (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
    IIGuiData = prepare_guidata(handles);
    figure('Color','White');
    plot(IIGuiData.test_time,IIGuiData.seerror,'b');xlabel('Time'), ylabel('Squared Error')
    title(sprintf('Squared Error for Model Order %d and Picking %d Residues',IIGuiData.NewVal,size(IIGuiData.SUB_IND,2)))
     
% --------------------------------------------------------------------
function plotnew_polezero_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_polezero (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
    IIGuiData = prepare_guidata(handles);
    figure('Color','White');
    zplane(IIGuiData.inummz,IIGuiData.idennz,'y');axis([-1.1 1.1 -1.1 1.1]),xlabel('Real Part'), ylabel('Imaginary Part')
    title('Poles Zero Plot')
    

% --------------------------------------------------------------------
function plotnew_sortedres_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_sortedres (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
    IIGuiData = prepare_guidata(handles);
    figure('Color','White');
    max_aa=max(IIGuiData.aa_list)+1;
    Index=1:size(IIGuiData.SUB_IND,2);
    plot(Index,IIGuiData.aa_list,'b'), axis([0.9 max(Index)+0.1 0 max_aa]),xlabel('Index'), ylabel('Residue Magnitude')
    title(sprintf('Prony Residue Magnitude "Decay" for %d Modes',max(Index)))



% --------------------------------------------------------------------
function plotnew_poles_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_poles (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
    IIGuiData = prepare_guidata(handles);
    figure('Color','White');
    zplane(0,IIGuiData.idennz,'y');axis([-1.1 1.1 -1.1 1.1]);xlabel('Real Part')
    ylabel('Imaginary Part'),title('Poles Plot');

% --------------------------------------------------------------------
function plotnew_allres_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_allres (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
    IIGuiData = prepare_guidata(handles);
    figure('Color','White');
    max_aa=max(IIGuiData.all_res)+1;
    Index=1:size(IIGuiData.res_SUBIND,2);
    plot(Index,IIGuiData.all_res,'b'), axis([0.9 max(Index)+0.1 0 max_aa]),xlabel('Index'), ylabel('Residue Magnitude')
    title(sprintf('Prony Residue Magnitude "Decay" for %d Modes',max(Index)))

% --------------------------------------------------------------------
function plotnew_allenergy_Callback(hObject, eventdata, handles)
% hObject    handle to plotnew_allenergy (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
    IIGuiData = prepare_guidata(handles);
    figure('Color','White');
    Index=1:IIGuiData.NewVal;
    semilogy(Index,IIGuiData.energy(IIGuiData.SUB_IND'),'b');
    axis([0 max(Index)+1 0 1.1*max(IIGuiData.energy)]),xlabel('Index'),ylabel('Energy')
    title('Energy of Prony Modes');


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


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


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


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



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

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

% Hint: place code in OpeningFcn to populate axes_pronyfit


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

% Hint: place code in OpeningFcn to populate axes_modesfit


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

% Hint: place code in OpeningFcn to populate axes_seerror


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

% Hint: edit 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



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

% Hints: get(hObject,'String') returns contents of lstbx_numeric as text
%        str2double(get(hObject,'String')) returns contents of lstbx_numeric as a double

%mcwHndl=handles.lstbx_numeric;
%set(gcf,'UserData',mcwHndl);




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

% Hint: edit 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




% --------------------------------------------------------------------
function menu_gridon_Callback(hObject, eventdata, handles)
% hObject    handle to menu_gridon (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
axes(handles.axes_pronyfit);
grid on;
axes(handles.axes_seerror);
grid on;
axes(handles.axes_modesfit);
grid on;
% --------------------------------------------------------------------
function menu_gridoff_Callback(hObject, eventdata, handles)
% hObject    handle to menu_gridoff (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
axes(handles.axes_pronyfit);
grid off;
axes(handles.axes_seerror);
grid off;
axes(handles.axes_modesfit);
grid off;

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



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


% --------------------------------------------------------------------
function menu_compare_Callback(hObject, eventdata, handles)
% hObject    handle to menu_compare (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% Call CompareSessions GUI if any session is saved else give a message to
% the user to save the session

if (handles.menu_saveflag==1)
    comparesessions;
    IIGuiData = prepare_guidata(handles);
    assignin('base','IIGuiData',IIGuiData);

else
    % Message to the user that data is saved
    warndlg('Please Save the Session to Use the Compare Sessions Tool !!',' Warning Save Session');
    whitebg('white');    
end



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

% Hint: edit 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



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

% Hints: get(hObject,'String') returns contents of edit_morder as text
%        str2double(get(hObject,'String')) returns contents of edit_morder as a double


% --- Executes during object creation, after setting all properties.
function ppmenu_pronyplots_CreateFcn(hObject, eventdata, handles)
% hObject    handle to ppmenu_pronyplots (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_pronyplots.
function ppmenu_pronyplots_Callback(hObject, eventdata, handles)
% hObject    handle to ppmenu_pronyplots (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_pronyplots contents as cell array
%        contents{get(hObject,'Value')} returns selected item from ppmenu_pronyplots
    


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

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



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

% Hints: get(hObject,'String') returns contents of edit6 as text
%        str2double(get(hObject,'String')) returns contents of edit6 as a double