% DoOnceThrough executes one curvf fit iteration of the DoMuMatchFitLSQRd % routine and generates figures, %Initialization tic delete fig*.fig close all; clear all; %Set the biconcave shape template, enter material constants and cell dimensions: INPUTS = make_INPUTS_FitMud; % INPUTS{1} is the AREA _TEMPLATE % INPUTS{2} is PROPS % INPUTS{3} is dims % INPUTS{4} is INDEX (# of points in different cell regions) % INPUTS{5} is rstart % INPUTS{6} is lam0 dims = INPUTS{3}; Lp = dims{1}; Rp = dims{3}; nvals = size(Lp); PROPS = INPUTS{2}; % The cell data to be fit are loaded in Make inputs via a call to load cell % data. The following could be used for plotting, but for now are commented % out to avoid confusion. There are three different experimental days. % These are data points from 6/18, 6/19 and 6/20/90, corrected data (See MUTEST.xls) % For 620, Dp = 1.7 um; For 6/19 and 6/18, Dp = 1.1 um % PRpdata620 = [10.3618 12.4342 14.5065 16.5789 20.7236 24.8684 20.7236 16.5789 14.5065 12.4342 10.3618 8.2895 6.2171]; % LoRpdata620 = [1.2742 1.3762 1.4911 1.5989 1.7500 1.9875 1.8435 1.6672 1.6036 1.4543 1.2899 1.1992 0.9612]; % LoRpErr620 = [0.2405 0.2087 0.2001 0.1843 0.1693 0.1699 0.1479 0.1577 0.0586 0.1501 0.1197 0.1343 0.1211]; % PRpdata620 = [10.3618 12.4342 14.5065 16.5789 20.7236 24.8684 ]; % LoRpdata620 = [1.2742 1.3762 1.4911 1.5989 1.7500 1.9875 ]; % LoRpErr620 = [0.2405 0.2087 0.2001 0.1843 0.1693 0.1699 ]; % For 6/19, Dp = 1.13 um % clo % PRpdata619 = [13.734 16.4808 19.2276 21.9744 27.468 21.9744 16.4808 10.9872 8.2404]; % LoRpdata619 = [1.22679 1.42679 1.55536 1.68929 1.90357 1.72143 1.52321 1.19107 0.97024]; % LoRpErr619 = [0.06469 0.08685 0.10762 0.113 0.11269 0.135 0.14037 0.1243 0.1171]; % PRpdata619 = [13.734 16.4808 19.2276 21.9744 27.468 ]; % LoRpdata619 = [1.22679 1.42679 1.55536 1.68929 1.90357 ]; % LoRpErr619 = [0.06469 0.08685 0.10762 0.113 0.11269 ]; % % % For 6/18, Dp = 1.12 um % % PRpdata618 = [8.31397 11.0853 13.85662 16.6280 19.3993 22.1706 27.7133 22.1706 16.6280 11.0853]; % LoRpdata618 = [1.16056 1.32491 1.49178 1.6182 1.75727 1.9115 2.03894 1.86549 1.66119 1.33755]; % LoRpErr618 = [0.11911 0.1365 0.11593 0.13279 0.10204 0.07782 0.12415 0.08072 0.09389 0.16409]; PRpdata618 = [8.31397 11.0853 13.85662 16.6280 19.3993 22.1706 27.7133 ]; LoRpdata618 = [1.16056 1.32491 1.49178 1.6182 1.75727 1.9115 2.03894 ]; LoRpErr618 = [0.11911 0.1365 0.11593 0.13279 0.10204 0.07782 0.12415 ]; % %varparam = [bovers0, sigap]; LoRp = Lp./Rp; %Set values for bovers0, sigap: varparams1 = [25.0, 2, 160.0]; [sumsqerr, PxRp_th1] = SumSqrErr_plots(LoRp,varparams1,INPUTS); varparams2 = [25.0, 2, 180.0]; [sumsqerr, PxRp_th2] = SumSqrErr_plots(LoRp,varparams2,INPUTS); varparams3 = [25.0, 2, 200.0]; [sumsqerr, PxRp_th3] = SumSqrErr_plots(LoRp,varparams3,INPUTS); varparams4 = [25.0, 2, 140.0]; [sumsqerr, PxRp_th4] = SumSqrErr_plots(LoRp,varparams4,INPUTS); %% figure(52); % errorbar(PRpdata618,LoRpdata618,LoRpErr618,'bs','MarkerSize',14,'linewidth',4); % hold on; % Lmax = PROPS{2}; x = PxRp_th4; y = LoRp; p=varparams4(1); sigap = varparams4(2); Lmax =varparams4(3); plot(x,y,'LineWidth',4.0, 'DisplayName',['Model: s_0=[20 ... 97] nm,',' s_{max}=',num2str(Lmax,'%1.0f'),' nm, p=',num2str(p,'%1.0f'),' nm']); hold on x = PxRp_th1; y = LoRp; p=varparams1(1); sigap = varparams1(2); Lmax =varparams1(3); plot(x,y,'LineWidth',4.0, 'DisplayName',['Model: s_0=[20 ... 97] nm,',' s_{max}=',num2str(Lmax,'%1.0f'),' nm, p=',num2str(p,'%1.0f'),' nm']); hold on x = PxRp_th2; y = LoRp; p=varparams2(1); sigap = varparams2(2); Lmax =varparams2(3); plot(x,y,'LineWidth',4.0, 'DisplayName',['Model: s_0=[20 ... 97] nm,',' s_{max}=',num2str(Lmax,'%1.0f'),' nm, p=',num2str(p,'%1.0f'),' nm']); hold on x = PxRp_th3; y = LoRp; p=varparams3(1); sigap = varparams3(2); Lmax =varparams3(3); plot(x,y,'LineWidth',4.0, 'DisplayName',['Model: s_0=[20 ... 97] nm,',' s_{max}=',num2str(Lmax,'%1.0f'),' nm, p=',num2str(p,'%1.0f'),' nm']); hold on y2=(x-x(1)).*(y(2)-y(1))/(x(2)-x(1))+y(1); % plot(x,y2,'k--','LineWidth',4.0); % errorbar(PRpdata620,LoRpdata620,LoRpErr620,'rs', ... % 'MarkerEdgeColor',[.5,1,.5], ... % 'MarkerFaceColor',[.5,1,.5]); % errorbar(PRpdata619,LoRpdata619,LoRpErr619,'bs') % plot set set(gca,'FontSize',36) %set labels %ylabel('L_p/R_p','FontSize',36) % x-axis label %xlabel('\DeltaP R_p (pN/\mum)','FontSize',36) % y-axis label xlim([0 45]); ylabel('{\boldmath $L_p/R_p$}','FontSize',40,'Interpreter','latex') % x-axis label xlabel('{\boldmath $\Delta P R_p (pN/\mu m)$}','FontSize',40,'Interpreter','latex') % y-axis label % modelx=sprintf('Model mu_0=%0.2f (pN/um), l_0=%0.2f (nm), L_{max}=%0.2f (nm), p=%0.2f (nm)',mu0,x(1),x(2),x(3)); % modelx2=sprintf('Model mu_0=%0.2f (pN/um), alpha =%0.3f',mu0,alph); % lgnd=legend({'Experiment Jun 19',modelx}); % lgnd=legend({'Mohandas, N., & Evans, E','Model \mu_0=%0.2f (pN/um)',mu0,'m','Model \mu_0=%0.2f (pN/um)',mu0,'m'}) % set(lgnd,'Interpreter','latex'); %legends %lgnd=legend({['Svetina, et al. 2016'], ... %['Model \mu_0=', num2str(mu1, '%1.0f'),'pN/\mum, l_0=',num2str(x(1),'%1.0f'),'nm',', L_{c}=',num2str(x(2),'%1.0f'),'nm, p=',num2str(x(3),'%1.0f'),'nm, \alpha_{0}=0'], ... %['Model \mu_0=',num2str(mu2,'%1.0f'),'pN/\mum, l_0=',num2str(xa(1),'%1.0f'),'nm',', L_{c}=',num2str(xa(2),'%1.0f'),'nm, p=',num2str(xa(3),'%1.0f'),'nm, \alpha_{0}=',num2str(alph,'%1.4f')]}); % lgnd=legend({['June 18, 1990, R_{p}=0.56\mum'], ... % ['Model \mu_0=', num2str(mu3, '%1.0f'),'pN/\mum, L_0_{eff}=',num2str(sd(1),'%1.0f'),'nm',', L_{c}=',num2str(xq(1),'%1.0f'),'nm, p=',num2str(xq(2),'%1.0f'),'nm, \alpha_{0}=0'], ... % ['Model \mu_0=',num2str(mu2,'%1.0f'),'pN/\mum, l_0=',num2str(xa(1),'%1.0f'),'nm',', L_{c}=',num2str(xa(2),'%1.0f'),'nm, p=',num2str(xa(3),'%1.0f'),'nm, \alpha_{0}=',num2str(alph,'%1.4f')]}); lgnd=legend; %({['Experiment, R_{p}=0.56 \mum'], ... %['Model, \mu_0=', num2str(mu3, '%1.0f'),' pN/\mum, s_0=[20 ... 97] nm,',char(10),'s_{max}=',num2str(Lmax,'%1.0f'),' nm, p=',num2str(varparams(1),'%1.0f'),' nm, \sigma_{ap}=',num2str(sigap,'%1.0f'), ' pN/\mum'],['Linear model without hardening']}); %lgnd=legend({['June 18, 1990, R_{p}=0.56\mum'], ... % lgnd=legend({['June 18, 1990, R_{p}=0.56\mum, \sigma_{ap}= ',num2str(sigma_ap, '%1.0f')], ... % ['Model \mu_0=', num2str(mu1, '%1.0f'),'pN/\mum, l_0=',num2str(xa3(1),'%1.0f'),'nm',', L_{c}=',num2str(xa3(2),'%1.0f'),'nm, p=',num2str(xa3(3),'%1.0f'),'nm, \alpha=0'], ... % ['Model \mu_0=',num2str(mu2,'%1.0f'),'pN/\mum, l_0=',num2str(xa4(1),'%1.0f'),'nm',', L_{c}=',num2str(xa4(2),'%1.0f'),'nm, p=',num2str(xa4(3),'%1.0f'),'nm']}); set(lgnd,'FontSize',28); set(lgnd,'FontName','Arial','Location','southeast'); %% % poolobj = gcp('nocreate'); % delete(poolobj); spreadfigures; toc