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

Changeset - 68c5df4ebb15

Parent rev.

Child rev.

[Not reviewed]

0 1 2

Tom Bannink - 8 years ago 2017-05-29 12:44:42
tombannink@gmail.com

Add plots of triangle exponent

3 files changed with 10 insertions and 4 deletions:

plots/avgtris_n.pdf

bin+new file

plots/triangle_exponent.pdf

bin+new file

triangle_exponent_plots.m

0 comments (0 inline, 0 general)

plots/avgtris_n.pdf

➞

Show inline comments

		new file 100644
		binary diff not shown

plots/triangle_exponent.pdf

➞

Show inline comments

		new file 100644
		binary diff not shown

triangle_exponent_plots.m

➞

Show inline comments

@@ @@ -16,7 +16,7 @@ Needs["ErrorBarPlots`"] @@
 (* When importing from exponent-only-data file *)
-gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent.m"];
+gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_mix32.m"];
 gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *)
 averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
@@ @@ -51,7 +51,10 @@ fitsExtra[[1]]["ParameterErrors"] @@
 fitsExtra[[1]]["ParameterConfidenceIntervals"]
 Show[ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->True,PlotMarkers->Automatic,AxesLabel->{"n","\[LeftAngleBracket]triangles\[RightAngleBracket]"},PlotLegends->taulabels],Plot[fits,{logn,1,2000}]]
 plot1=Show[ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->False,PlotMarkers->Automatic,AxesLabel->{"n","\[LeftAngleBracket]triangles\[RightAngleBracket]"},PlotLegends->taulabels],Plot[fits,{logn,1,2000}]]
 Export[NotebookDirectory[]<>"plots/avgtris_n.pdf",plot1]
 (* ::Subsection:: *)
 tauValues=averagesGrouped[[All,1,1,1,2]];
 exponents=Transpose[{tauValues,fits[[All,2,1]]}];
-Show[ListPlot[exponents,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","T(\[Tau])"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
+Show[ListPlot[exponents,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
 (* ::Subsection:: *)
 tauValues=averagesGrouped[[All,1,1,1,2]];
 exponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues,fitsExtra}]];
 Show[ErrorListPlot[exponentsErrorBars,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","T(\[Tau])"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
 plot2=Show[ErrorListPlot[exponentsErrorBars,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)