- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2022-03-17 12:12:51+13:00
- Desc:
- Move the folders to BG; Add open channel testing
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/64f/rawfile/f60338fb964a811a3c12944d752e3b0f5165ce39/BG/scripts_BG/Na_gating_m.m
function Na_gating_m(T,phi,V)
rng default % For reproducibility
%% Steady-state gating parameters and time constants of HH
[alpha,beta]=HH_gate_m(-(V+65),phi);
tau = calc_tau(alpha,beta);
g_ss = calc_gss(alpha,beta);
%% Fit new kinetic parameters for the bond graph of state models
% params: [kf, zf, kr, zr];
error_func_alpha = @(params) square_error(alpha - calc_alpha(params,V/1000,T));
error_func_beta = @(params) square_error(beta - calc_beta(params,V/1000,T));
error_func_gss = @(params) square_error(g_ss - p2gss(params,V,T));
error_func_tau = @(params) square_error(tau- p2tau(params,V,T));
%error_func = @(params) error_func_alpha(params) + error_func_beta(params) + error_func_gss(params) + error_func_tau(params);
error_func = @(params) 10.*error_func_alpha(params) + error_func_beta(params) + 1000000.* error_func_gss(params) + 100000.* error_func_tau(params);
lb = [0; -10; 0; -10];
ub = [Inf; 10; Inf; 10];
options_ps = optimoptions('particleswarm','UseParallel',true,'HybridFcn',@fmincon,'SwarmSize',200,'FunctionTolerance',1e-8);
[params_vec,fval,exitflag,output] = particleswarm(error_func,4,lb,ub,options_ps);
save(['Na_m_parameters.mat'],'params_vec');
MSE = (square_error(alpha - calc_alpha(params_vec,V/1000,T)) + square_error(beta - calc_beta(params_vec,V/1000,T))...
+ square_error(g_ss - p2gss(params_vec,V,T)) +square_error(tau- p2tau(params_vec,V,T)))/length(V);
%% Plot for comparison
h = figure;
x0=100;
y0=50;
height=350*2;
width=400*2;
set(gcf,'position',[x0,y0,width,height])
gate_plot(alpha,beta,g_ss,tau,params_vec,V,T,'m',MSE);
filename='m_gate';
savefig(h,sprintf('%s.fig',filename))
exportgraphics(h,sprintf('%s.eps',filename));
