rlc_analyzer.m

function varargout = rlc_analyzer(varargin)
%instantiate gui singleton
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @rlc_analyzer_OpeningFcn, ...
                   'gui_OutputFcn',  @rlc_analyzer_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(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

function rlc_analyzer_OpeningFcn(hObject, eventdata, handles, varargin)
%profile on;

%%initialize application variables
handles.output = hObject;

%current_circuit assumes a state of series or parallel
handles.series = 0;
handles.parallel = 1;
handles.current_circuit = handles.series;

%metrics used
handles.neper_frequency = 0;
handles.resonant_radian_frequency = 0;
handles.damped_frequency = 0;
handles.s1 = 0;
handles.s2 = 0;

%damp_type assumes a state of underdamped, overdamped or critdamped
handles.damp_type = -1;
handles.underdamped = 0;
handles.overdamped = 1;
handles.critdamped = 2;

%initial conditions
handles.v0 = 0;
handles.il = 0;

handles.i0 = 0;
handles.vc = 0;

handles.xf = 0;

%RLC parameters
handles.resistance =str2double(get(handles.R_text_edit, 'String'));
handles.inductance = str2double(get(handles.L_text_edit, 'String'))/10^3;
handles.capacitance = str2double(get(handles.C_text_edit, 'String'))/10^6;

%Build GUI that GUIDE doesn't support
handles.tabs = uitabgroup(handles.tab_container);
handles.results_tab = uitab(handles.tabs,'Title','Results');
handles.graph_tab = uitab(handles.tabs,'Title','Graph', 'ButtonDownFcn', @(hObject,eventdata)update_view(hObject,eventdata));

handles.results_table =  uitable(handles.results_tab,'units', 'normalized', 'Position',[0 0 1 1]);

handles.results_table.ColumnName = {'Metrics','Value','Units'};
handles.graph = axes('Parent',handles.graph_tab);

%Evaluate variables to models and draw views
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);

draw_circuit(hObject, eventdata, handles);


function update_model(hObject, eventdata, handles)
%calculate frequencies and roots
handles.resonant_radian_frequency = 1/(handles.inductance*handles.capacitance)^.5;
if handles.current_circuit == handles.series
    handles.neper_frequency = handles.resistance /(2*handles.inductance);
else
    handles.neper_frequency = 1/(2*handles.resistance*handles.capacitance); 
end

handles.damped_frequency = (handles.resonant_radian_frequency^2 - handles.neper_frequency^2)^.5;
handles.s1 = -handles.neper_frequency + (handles.neper_frequency^2 - handles.resonant_radian_frequency^2)^.5;
handles.s2 = -handles.neper_frequency - (handles.neper_frequency^2 - handles.resonant_radian_frequency^2)^.5;

%determine damp type
if isreal(handles.s1) && isreal(handles.s2)
    if handles.s1 == handles.s2
        handles.damp_type = handles.critdamped;
    else
        handles.damp_type = handles.overdamped;
    end
else
    handles.damp_type = handles.underdamped;
end

%determine damp type
switch handles.current_circuit;
    case handles.series
        handles.i0 = str2double(get(handles.ic_edit_text_1, 'String'));
        handles.vc = str2double(get(handles.ic_edit_text_2, 'String'));
    case handles.parallel
        handles.v0 = str2double(get(handles.ic_edit_text_1, 'String'));
        handles.il = str2double(get(handles.ic_edit_text_2, 'String'));
        
end

handles.xf = str2double(get(handles.ic_edit_text_3, 'String'));

guidata(hObject, handles);


function update_view(hObject, eventdata, handles)
handles = guidata(hObject);

%render IC text
switch handles.current_circuit;
    case handles.series
        set(handles.ic_name_1, 'String', 'I0:')
        set(handles.ic_units_1, 'String', 'A')
        set(handles.ic_name_2, 'String', 'VC:')
        set(handles.ic_units_2, 'String', 'V')
        set(handles.ic_name_3, 'String', 'IF:')
        set(handles.ic_units_3, 'String', 'A')
    case handles.parallel
        set(handles.ic_name_1, 'String', 'V0:')
        set(handles.ic_units_1, 'String', 'V')
        set(handles.ic_name_2, 'String', 'IL:')
        set(handles.ic_units_2, 'String', 'A')
        set(handles.ic_name_3, 'String', 'VF:')
        set(handles.ic_units_3, 'String', 'V')
end

%update states
draw_results(hObject, eventdata, handles);
draw_graph(hObject, eventdata, handles);


function draw_results(hObject, eventdata, handles)
%draws model to table view
type = '';
switch handles.current_circuit;
    case handles.series
        type = 'Series';
    case handles.parallel
        type = 'Parallel';
end

%Determine text for ciruit type cell
damping = '';
switch handles.damp_type;
case handles.overdamped 
   damping = 'Overdamped';
case handles.underdamped 
   damping = 'Underdamped';
case handles.critdamped 
   damping = 'Critically Damped';
end
handles.results_data = {
                    'R', handles.resistance,'Ohms';...
                    'L', handles.inductance, 'Henries';...
                    'C', handles.capacitance, 'Farads'
                    'RLC Type', type, '';...
                    'Neper Frequency', handles.neper_frequency, 'rad/s';...
                    'Resonant Radian Frequency', handles.resonant_radian_frequency, 'rad/s';...
                    'Damped Frequency', handles.damped_frequency, 'rad/s';...
                    's1', handles.s2, 'rad/s';...
                    's2',handles.s1, 'rad/s';...
                    'Damping', damping, ''};
                
switch handles.current_circuit;
    case handles.series
        handles.results_data = cat(1,handles.results_data, {'Initial Current ', handles.i0, 'A';'Capacitor Voltage',handles.vc,'V'; 'Final Current', handles.xf, 'A'});
    case handles.parallel
        handles.results_data = cat(1,handles.results_data, {'Initial Voltage ', handles.v0, 'V';'Inductor Current',handles.il,'A';'Final Voltage', handles.xf, 'V'});
end       
handles.results_table.Data =  handles.results_data;
guidata(hObject, handles);


function save_state_csv(filename, pathname, table, col)
myTable =  cell2table(table, 'VariableNames', col);
writetable(myTable,[strcat(pathname,filename)],'WriteRowNames',true);

function draw_axes(hObject, eventdata, handles)
cla(handles.graph);
syms v(t)
C=handles.capacitance;
R=handles.resistance;
L=handles.inductance;
switch handles.current_circuit;
    case handles.series
        x0 = handles.i0;
        dx0 = (-handles.vc-handles.i0*R)/L;
        eqn = -L*diff(v,t,2)+ handles.xf/(C)== R*diff(v)+1./(C).*v;
        Dv = diff(v,t);
        ySol(t) = dsolve(eqn, [v(0)==x0 Dv(0)==dx0]);
        t = 'i(t)=';
        u = 'Current (A)';
    case handles.parallel
        x0 = handles.v0;
        dx0 = (-handles.il-handles.v0/R)/C;
        eqn = -C.*diff(v,t,2) + + handles.xf/(L)==1./R.*diff(v)+1./L.*v;
        Dv = diff(v,t);
        ySol(t) = dsolve(eqn, [v(0)==x0 Dv(0)==dx0]);
        t = 'v(t)=';
        u = 'Voltage (V)';
end
tmax = 20/handles.resonant_radian_frequency;
x = 0:tmax/100:tmax; 


ySol = vpa(ySol);
A= double(coeffs(simplify(ySol)));

y = double(ySol(x));

indexmin = find(min(y) == y);
xmin = x(indexmin);
ymin = y(indexmin);

ymin = min(y);
ymax = max(y);
[pks_max,locs_max] = findpeaks(y) ;
[pks_min, locs_min] =findpeaks(ymax-y); 

axes(handles.graph);
hold(handles.graph,'on');
switch handles.damp_type;
case handles.overdamped % User selects peaks.
   plot(handles.graph, x, y, 'linewidth', 2);
case handles.underdamped % User selects peaks.
   %plot(handles.graph, x, A*exp(-handles.neper_frequency*x), '--r');
   %plot(handles.graph, x, -A*exp(-handles.neper_frequency*x), '--r');
   plot(x, y, 'linewidth', 2);
case handles.critdamped % User selects peaks.
   plot(handles.graph, x, y, 'linewidth', 2);
end

xlim(handles.graph, [0 tmax]);
incs = (ymax-ymin)/15;
ylim(handles.graph, [ymin-incs ymax+incs]);

for idx = 1:numel(locs_max)
    loc = locs_max(idx);
    strmax = ['',num2str(y(loc))];
    text(x(loc),incs/2+y(loc),strmax);
    plot(x(loc), y(loc), 'o', 'LineWidth',1,'MarkerEdgeColor',...
    [ 0    0.4470    0.7410],'MarkerFaceColor',[ 0    0.4470    0.7410],'MarkerSize',6);
    
end

for idx = 1:numel(locs_min)
    loc = locs_min(idx);
    strmax = ['',num2str(y(loc))];
    text(x(loc),-incs/2+y(loc),strmax);
    plot(x(loc), y(loc), 'o', 'LineWidth',1,'MarkerEdgeColor',...
    [ 0    0.4470    0.7410],'MarkerFaceColor',[ 0    0.4470    0.7410],'MarkerSize',6);
    
end

n = vpa((ySol), 3);
d1 = digits(4);
title(handles.graph, { ['$' t latex(simplify(n)) '$']}, 'Interpreter', 'latex');
digits(d1);
xlabel(handles.graph, 'Time (s)', 'Interpreter', 'latex');
ylabel(handles.graph, u, 'Interpreter', 'latex');
hold(handles.graph,'off');


function draw_graph(hObject, eventdata, handles)
if(handles.tabs.SelectedTab == handles.graph_tab)
    draw_axes(hObject, eventdata, handles)
end 


function varargout = rlc_analyzer_OutputFcn(hObject, eventdata, handles) 
varargout{1} = handles.output;


function menu_file_Callback(hObject, eventdata, handles)


function menu_save_Callback(hObject, eventdata, handles)
[FileName,PathName] = uiputfile('*.csv','Save RLC As');
if(~(FileName==0))
save_state_csv(FileName, PathName, handles.results_data, handles.results_table.ColumnName);
end


function menu_load_Callback(hObject, eventdata, handles)
[FileName,PathName] = uigetfile('*.csv','Select an RCL csv');
if(~(FileName==0))
temp = readtable(strcat(PathName, FileName));
handles.resistance = str2double(temp.Value{1});
handles.capacitance = str2double(temp.Value{3});
handles.inductance = str2double(temp.Value{2});
switch temp.Value{4};
    case 'Series'
        handles.current_circuit = handles.series;
    case 'Parallel'
        handles.current_circuit = handles.parallel;
end
set(handles.select_menu, 'Value', handles.current_circuit+1);

set(handles.R_text_edit,'string',num2str(handles.resistance));
set(handles.L_text_edit,'string',num2str(handles.inductance));
set(handles.C_text_edit,'string',num2str(handles.capacitance));

if(height(temp)>10)
    switch handles.current_circuit;
    case handles.series
        handles.i0 = str2double(temp.Value{11});
        handles.vc = str2double(temp.Value{12});
        
        set(handles.ic_edit_text_1,'string',num2str(handles.i0));
        set(handles.ic_edit_text_2,'string',num2str(handles.vc));
        
    case handles.parallel
        handles.v0 = str2double(temp.Value{11});
        handles.il = str2double(temp.Value{12});
        
        set(handles.ic_edit_text_1,'string',num2str(handles.v0));
        set(handles.ic_edit_text_2,'string',num2str(handles.il));
    end 
if(height(temp)>12)
handles.xf = str2double(temp.Value{13});
set(handles.ic_edit_text_3,'string',num2str(handles.xf));
end
end

guidata(hObject, handles);
handles = guidata(hObject);
draw_circuit(hObject, eventdata, handles);
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);
end


function select_menu_Callback(hObject, eventdata, handles)
str = get(hObject, 'String');
val = get(hObject,'Value');
switch str{val};
case 'Series' % User selects peaks.
   handles.current_circuit = handles.series;
case 'Parallel' % User selects membrane.
   handles.current_circuit = handles.parallel;
end
% Save the handles structure.
draw_circuit(hObject, eventdata, handles);
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function draw_circuit(hObject, eventdata, handles)
axes(handles.type_axes)
if handles.current_circuit == handles.series
    matlabImage = imread('./img/series.png');
else
    matlabImage = imread('./img/parallel.png');
end
image(matlabImage)
axis off
axis image


function select_menu_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


function calculate_button_Callback(hObject, eventdata, handles)
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function R_text_edit_Callback(hObject, eventdata, handles)
resistance = str2double(get(hObject, 'String'));
if isnan(resistance)
    set(hObject, 'String', 0);
    errordlg('Input must be a number','Error');
end

handles.resistance = resistance;
guidata(hObject,handles);
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function R_text_edit_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


function L_text_edit_Callback(hObject, eventdata, handles)
inductance = str2double(get(hObject, 'String'));
if isnan(inductance)
    set(hObject, 'String', 0);
    errordlg('Input must be a number','Error');
end

handles.inductance = inductance/10^3;
guidata(hObject,handles)
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function L_text_edit_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


function C_text_edit_Callback(hObject, eventdata, handles)
capacitance = str2double(get(hObject, 'String'));
if isnan(capacitance)
    set(hObject, 'String', 0);
    errordlg('Input must be a number','Error');
end

handles.capacitance = capacitance/10^6;
guidata(hObject,handles)
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function C_text_edit_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


function menu_about_Callback(hObject, eventdata, handles)
uiwait(msgbox({'Made by Team StEthan for ENGR 222 Winter 2017';'';'Testing and Logistics:                       Stephen Poanessa';'Programming and Documentation:   Ethan Lew';'';'My software never has bugs. It just develops random features.';'';'""Legacy code" often differs from its suggested alternative by actually working and scaling."-Bjarne Stroustrup, C++ Creator' } ,'About','modal'));


function menu_help_Callback(hObject, eventdata, handles)


function menu_help_bug_Callback(hObject, eventdata, handles)
web('https://github.com/EthanJamesLew/RLC_Analyzer/issues', '-browser');


function menu_help_git_Callback(hObject, eventdata, handles)
web('https://github.com/EthanJamesLew/RLC_Analyzer', '-browser');


function menu_file_exit_Callback(hObject, eventdata, handles)
delete(handles.figure1);
%profile viewer;


function ic_edit_text_1_Callback(hObject, eventdata, handles)
val = str2double(get(hObject, 'String'));
if isnan(val)
    set(hObject, 'String', 0);
    errordlg('Input must be a number','Error');
end
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function ic_edit_text_1_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


function ic_edit_text_2_Callback(hObject, eventdata, handles)
val = str2double(get(hObject, 'String'));
if isnan(val)
    set(hObject, 'String', 0);
    errordlg('Input must be a number','Error');
end
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);


function ic_edit_text_2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


function menu_file_print_Callback(hObject, eventdata, handles)
set(0,'DefaultFigureVisible','off');
draw_axes(hObject, eventdata, handles);
h = figure;
copyobj(handles.graph,h);
set(0,'DefaultFigureVisible','on');
printpreview(h);

function left_panel_SizeChangedFcn(hObject, eventdata, handles)

function figure1_SizeChangedFcn(hObject, eventdata, handles)
handles = guidata(hObject);
set(handles.results_tab, 'units', 'pixels');
t = get(handles.results_tab, 'Position');
t_width = t(3)/3;
set(handles.results_table, 'ColumnWidth', {t_width,t_width,t_width});


set(handles.figure1, 'units', 'pixels');
set(handles.left_panel, 'units', 'pixels');
set(handles.tab_container, 'units', 'pixels');
t = get(handles.figure1, 'Position');
if(t(4) > 585 & t(3) > 201)
set(handles.left_panel, 'Position', [0 t(4)-585 201 585]);
set(handles.tab_container, 'Position', [201 0 t(3)-201 t(4)]);
end
guidata(hObject, handles);



function ic_edit_text_3_Callback(hObject, eventdata, handles)
val = str2double(get(hObject, 'String'));
if isnan(val)
    set(hObject, 'String', 0);
    errordlg('Input must be a number','Error');
end
update_model(hObject, eventdata, handles);
update_view(hObject, eventdata, handles);

function ic_edit_text_3_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end