function varargout = GuiTimeResponse(varargin)
% GUITIMERESPONSE M-file for GuiTimeResponse.fig
% GUITIMERESPONSE, by itself, creates a new GUITIMERESPONSE or raises the existing
% singleton*.
%
% H = GUITIMERESPONSE returns the handle to a new GUITIMERESPONSE or the handle to
% the existing singleton*.
%
% GUITIMERESPONSE('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in GUITIMERESPONSE.M with the given input arguments.
%
% GUITIMERESPONSE('Property','Value',...) creates a new GUITIMERESPONSE or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before GuiTimeResponse_OpeningFunction gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to GuiTimeResponse_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help GuiTimeResponse
% Last Modified by GUIDE v2.5 27-Dec-2003 17:24:24
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @GuiTimeResponse_OpeningFcn, ...
'gui_OutputFcn', @GuiTimeResponse_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 GuiTimeResponse is made visible.
function GuiTimeResponse_OpeningFcn(hObject, eventdata, handles, varargin)
% Choose default command line output for GuiTimeResponse
handles.output = hObject;
% Load data
global strFilterObject
if length(varargin) == 1
DoWhat = varargin{1}; % Impulse, Step, Ramp, or Parabola
strFilterObject=Utility_zpk(strFilterObject); % find poles, zeros
else
temp=load('matlab');
disp([mfilename ' called in debug mode using matlab.mat datafile'])
strFilterObject = temp.strFilterObject;
DoWhat = 'Impulse';
end
switch DoWhat
case 'Impulse', set(handles.uitxTitle,'String','Impulse Response')
case 'Step', set(handles.uitxTitle,'String','Step Response')
case 'Ramp', set(handles.uitxTitle,'String','Ramp Response')
case 'Parabola', set(handles.uitxTitle,'String','Parabola Response')
end
handles.DoWhat = DoWhat;
set(handles.uiTimeResponse,'Name',strFilterObject.sTitle)
% save the basic variables
handles.TMin = 0;
handles.TMax = 3 * max(-1./real(strFilterObject.vPoles)); % 3 time constants
handles.DrawPeak = 0;
handles.DrawSettle = 0;
handles.fDrawSettle = 10; % 10%
handles.DrawFinal = 0;
if isempty(strFilterObject.vPoles1)
handles.DrawStandard = 0;
set(handles.uitxStandard,'Enable','off')
set(handles.uicbStandard,'Enable','off')
set(handles.uifrBox,'BackgroundColor',[.9 .9 .9]);
else
handles.DrawStandard = 1;
set(handles.uitxStandard,'Enable','on')
set(handles.uicbStandard,'Enable','on')
set(handles.uicbStandard,'Value',1)
end
% Calculate y and plot it
guidata(hObject, handles);
Recalculate(handles);
function varargout = GuiTimeResponse_OutputFcn(hObject, eventdata, handles)
varargout{1} = handles.output; % default output
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Limit Callbacks
% uitxTMin, uitxTMax
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function uitxTMin_Callback(hObject, eventdata, handles)
sTMin = get(hObject,'String');
if findstr('default',lower(sTMin))
sTMin = '0';
set(hObject,'String','default ')
end
TMin = str2num(sTMin);
if isempty(TMin)
errordlg('Enter the simulation start time, or "default"','Error')
set(hObject,'String',num2str(handles.TMin))
elseif TMin >= handles.TMax
errordlg('Simulation start time must be at least zero but less than the maximum simulation time','Error')
set(hObject,'String',num2str(handles.TMin));
else
handles.TMin = TMin;
guidata(hObject, handles);
Recalculate(handles)
end
function uitxTMax_Callback(hObject, eventdata, handles)
global strFilterObject
sTMax = get(hObject,'String');
if findstr('default',lower(sTMax)) % 3 time constants
sTMax = num2str(3 * max(-1./real(strFilterObject.vPoles)));
set(hObject,'String','default ')
end
TMax = str2num(sTMax);
if isempty(TMax)
errordlg('Enter the simulation end time, or "default"','Error')
set(hObject,'String',num2str(handles.TMax))
elseif TMax <= handles.TMin
errordlg('Simulation end time must be greater than the simulation start time','Error')
set(hObject,'String',num2str(handles.TMax));
else
handles.TMax = TMax;
guidata(hObject, handles);
Recalculate(handles)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Analysis Callbacks
% uicbPeak, uitxPeakY, uitxPeakT
% uicbSettle, uipmSettleY, uipmSettleT
% uicbFinal, uitxFinal, uitxFinal
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function uicbPeak_Callback(hObject, eventdata, handles)
handles.DrawPeak = get(hObject,'Value');
guidata(hObject, handles);
Recalculate(handles)
function uicbSettle_Callback(hObject, eventdata, handles)
handles.DrawSettle = get(hObject,'Value');
guidata(hObject, handles);
Recalculate(handles)
function uipmSettleY_Callback(hObject, eventdata, handles)
vcell = get(hObject,'String');
handles.fDrawSettle = str2num(vcell{get(hObject,'Value')});
guidata(hObject, handles);
Recalculate(handles)
function uicbFinal_Callback(hObject, eventdata, handles)
handles.DrawFinal = get(hObject,'Value');
guidata(hObject, handles);
Recalculate(handles)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Component Value Callback
% uicbStandard
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function uicbStandard_Callback(hObject, eventdata, handles)
handles.DrawStandard = get(hObject,'Value');
guidata(hObject, handles);
Recalculate(handles)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Helper functions
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Recalculate(handles)
global strFilterObject
%---------------------------Setup---------------------------
DoWhat=handles.DoWhat;
TMin = handles.TMin;
TMax = handles.TMax;
z = strFilterObject.vZeros;
p = strFilterObject.vPoles;
k = strFilterObject.fK;
tiny = 1e-40;
tol = 1e-10;
switch DoWhat
case 'Impulse'
yFinal = real(k*prod([tiny;tiny-z])/prod(tiny-p));
case 'Step'
yFinal=real(k*prod([tiny-z])/prod(tiny-p));
case 'Ramp'
yFinal=real(k*prod([tiny-z])/prod([tiny;tiny-p]));
case 'Parabola'
yFinal=real(k*prod([tiny-z])/prod([tiny;tiny;tiny-p]));
end
if abs(yFinal) < tol
yFinal = 0;
elseif abs(yFinal) > 1/tol
yFinal = inf;
end
%-----------------------DrawResponse-------------------------
[y,t] = Utility_CalculateTime(z,p,k,TMin,TMax,DoWhat);
if TMin>0 % eliminate result before TMin
keep=find(t>=TMin);
ydraw=y(keep);
tdraw=t(keep);
else
ydraw = y; tdraw=t;
end
plot(tdraw,ydraw,'k-')
axlimits=axis;
hold on
xlabel('time (s)')
%----------------------DrawImpulse-------------------------
if isequal(DoWhat,'Impulse') & isequal(strFilterObject.sPurpose,'Highpass')
xtext=axlimits(1)+(axlimits(2)-axlimits(1))*0.25;
ytext=mean(axlimits(3:4));
stext=sprintf('plus an impulse at t=0 of area %g',k);
text(xtext,ytext,stext)
end
%-----------------------DrawPeak---------------------------
if handles.DrawPeak
[yPeak,index]=max(y);
tPeak=t(index);
vx = [TMin tPeak tPeak]; vy = [yPeak yPeak axlimits(3)];
if yFinal < inf
plot(vx,vy)
set(handles.uitxPeakT,'String',num2str(tPeak))
set(handles.uitxPeakY,'String',num2str(yPeak))
else
set(handles.uitxPeakT,'String','Infinity')
set(handles.uitxPeakY,'String','Infinity')
end
else
set(handles.uitxPeakT,'String','')
set(handles.uitxPeakY,'String','')
end
%-----------------------DrawSettle---------------------------
if handles.DrawSettle
% settles when y is within max-min% of the final value
yRange = max(y) - min(y); % find in case DrawSettle is called with Tmin > o
yLow = yFinal - handles.fDrawSettle*yRange/100;
yHigh = yFinal + handles.fDrawSettle*yRange/100;
if y(end) > yHigh | y(end) < yLow
if yFinal == inf
sTSettle = 'Infinity';
else
sTSettle = '?'; % increase the time plotted to find it
end
else
index1=find(y<yHigh & y>yLow);
if length(index1)==1
tSettle = t(end);
else
index2=find(diff(index1)>1);
if isempty(index2) % as soon as inbounds, it stays there
tSettle = t(index1(1));
else
index2 = index2(end)+1;
if index2>length(index1), index2=index1(end); end
tSettle = t(index1(index2));
end
end
sTSettle=num2str(tSettle);
vx1 = [TMin TMax NaN TMin TMax NaN tSettle tSettle];
vy1 = [yHigh yHigh NaN yLow yLow NaN axlimits(3) yHigh];
plot(vx1,vy1,'g-')
end
set(handles.uitxSettleT,'String',sTSettle)
set(handles.uipmSettleY,'Enable','on')
else
set(handles.uitxSettleT,'String','')
set(handles.uipmSettleY,'Enable','off')
end
%-----------------------DrawFinal---------------------------
if handles.DrawFinal
vx = [TMin TMax]; vy=[yFinal yFinal];
plot(vx,vy,'r-')
set(handles.uitxFinalT,'String','Infinity')
if yFinal == inf
set(handles.uitxFinalY,'String','Infinity')
else
set(handles.uitxFinalY,'String',num2str(yFinal))
end
else
set(handles.uitxFinalT,'String','')
set(handles.uitxFinalY,'String','')
end
%-----------------------DrawStandard---------------------------
if handles.DrawStandard % blue
z1 = strFilterObject.vZeros1; p1 = strFilterObject.vPoles1; k1 = strFilterObject.fK1;
[yStandard,tStandard] = Utility_CalculateTime(z1,p1,k1,TMin,TMax,DoWhat);
if TMin>0 % eliminate result before TMin
keep=find(t>=TMin);
yStandard=yStandard(keep);
tStandard=tStandard(keep);
end
plot(tStandard,yStandard,'b-')
end
%-----------------------Cleanup---------------------------
axlimits=axis;
axlimits(1) = TMin; axlimits(2) = TMax;
axis(axlimits)
hold off
