AENC/switchchain Changeset - ea18a7bd1a0d · Centrum Wiskunde & Informatica (CWI)

Changeset - ea18a7bd1a0d

Parent rev.

Child rev.

[Not reviewed]

0 2 0

Tom Bannink - 8 years ago 2017-07-03 17:22:59
tombannink@gmail.com

Plot log version of exponential fit for mixing time

2 files changed with 15 insertions and 3 deletions:

triangle_analysis.m

triangle_exponent_plots.m

0 comments (0 inline, 0 general)

triangle_analysis.m

➞

Show inline comments

 (* Use fit["BestFitParameters"] to get parameters *)
 (* Use fit[t] to get fit value *)
 fitFuncsT=Map[#[t]&,fitList];
 tmixList=Map[tmix/.#["BestFitParameters"]&,fitList];
 timeplot1=ListPlot[coarseData,Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime},PlotStyle->Opacity[0.5]];
 Show[timeplot1,Plot[fitFuncsT,{t,1,maxTime},PlotRange->All]]
 (* Log version of exponential fits *)
 fitAverages=Map[c/.#["BestFitParameters"]&,fitList];
 shiftedFitFuncsT=MapIndexed[#1[t]-fitAverages[[#2[[1]]]]&,fitList];
 shiftedCoarseData=MapIndexed[MovingAverage[#1[[2]],1000][[1;;-1;;skipPts]]-fitAverages[[#2[[1]]]]&,selectedData];
 (* Plot log version *)
 timeplot2=ListLogPlot[shiftedCoarseData[[1;;5]],Joined->True,PlotRange->{0*minCount+0.1,maxCount},DataRange->{0,measureSkip*maxTime},PlotStyle->Opacity[0.5]];
 Show[timeplot2,LogPlot[Evaluate[shiftedFitFuncsT[[1;;5]]],{t,1,maxTime},PlotRange->All,PlotStyle->Dotted]]
 (* ::Subsection:: *)
 (*Plot success rate over "time"*)

triangle_exponent_plots.m

➞

Show inline comments

 (* For visual, shift the tau values slightly left or right to distinguish the two datasets *)
 tauValues=averagesGrouped[[All,1,1,1,2]];
 averagesExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
-Transpose[{tauValues-0.005,averagesFitsExtra}]];
+Transpose[{tauValues-0.001,averagesFitsExtra}]];
 mediansExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues+0.005,mediansFitsExtra}]];
 plot4=Show[ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars},Joined->True,PlotMarkers->Automatic,Frame->True,FrameLabel->{"tau","triangle exponent"},PlotRange->{{2,3},{0,1.6}},ImageSize->300],Plot[3/2(3-tau),{tau,2,3},PlotStyle->{Dashed}]]
 Transpose[{tauValues+0.001,mediansFitsExtra}]];
 plot2=Show[ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars},Joined->True,PlotMarkers->Automatic,Frame->True,FrameLabel->{"tau","triangle exponent"},PlotRange->{{2,3},{0,1.6}},ImageSize->300],Plot[3/2(3-tau),{tau,2,3},PlotStyle->{Dashed}]]
 Export[NotebookDirectory[]<>"plots/triangle_exponent.pdf",plot2]

0 comments (0 inline, 0 general)