FX4_ANALYZE_Lp_Lc

분석된 힘과 길이를 이용하여 persistence length 및 contour length 결정

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%  Top level section 4: FX4_ANALYZE_Lp_Lc %%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 

if 1

    

    %%% This is the forth section in a series of routines to analyze

    %%% force-extension data taken with the new multi bead code for MT

    %%% (J. P. Cnossen, D. Dulin and N. H. Dekker Rev Sci Instrum 85, (2014).)

    %%%

    %%% This section assumes that you have carried out fits of the WLC model for

    %%% all beads and all magnet heights/forces

    %%%

    %%% Author: Jan Lipfert

    %%% Date:  2014-07-17

    

    

    

    

    %%% Collect the results of the WLC fits, for the various force

    %%% determination methods

    PSD_Lps = [];

    PSD_Lcs = [];

    AV_Lps = [];

    AV_Lcs = [];

    SA_Lps = [];

    SA_Lcs = [];

    

    for i=1:Nbeads

        

        PSD_Lps = [PSD_Lps bead(i).PSD_Lp];

        PSD_Lcs = [PSD_Lcs bead(i).PSD_Lc];

        AV_Lps = [AV_Lps bead(i).AV_Lp];

        AV_Lcs = [AV_Lcs bead(i).AV_Lc];

        SA_Lps = [SA_Lps bead(i).SA_Lp];

        SA_Lcs = [SA_Lcs bead(i).SA_Lc];

        

    end

    

    

    %%% Histogram the results

    if 1

        centersVals = [0:5:200];

        

        figure(8); clf; box on; hold on;

        subplot(1,2,1); box on; hold on;

        [counts centers] = hist(PSD_Lps*1000, centersVals);

        plot(centers, counts, 'bo-','linewidth', 2)

        

        [counts centers] = hist(AV_Lps*1000, centersVals);

        plot(centers, counts, 'ro-','linewidth', 2)

        

        [counts centers] = hist(SA_Lps*1000, centersVals);

        plot(centers, counts, 'go-','linewidth', 2)

        

        set(gca,'LineWidth', 1,'FontSize', 15, 'FontWeight', 'bold', 'TickLength',[0.02 0.02])

        xlabel('Lp (nm)'); ylabel('Counts')

        title(['Lp(PSD) = ' num2str(median(PSD_Lps*1000),3) ' nm'  ])

        

        

        centersVals = [0:0.2:8];

        

        subplot(1,2,2); box on; hold on;

        [counts centers] = hist(PSD_Lcs, centersVals);

        plot(centers, counts, 'bo-','linewidth', 2)

        

        [counts centers] = hist(AV_Lcs, centersVals);

        plot(centers, counts, 'ro-','linewidth', 2)

        

        [counts centers] = hist(SA_Lcs, centersVals);

        plot(centers, counts, 'go-','linewidth', 2)

        

        set(gca,'LineWidth', 1,'FontSize', 15, 'FontWeight', 'bold', 'TickLength',[0.02 0.02])

        xlabel('Lc (um)'); ylabel('Counts')

        title(['Lp(PSD) = ' num2str(median(PSD_Lcs),3) ' um' ])

        legend('PSD', 'Allan', 'Aartjan', 'location', 'northwest')

        

    end

    

    

end