我想使用网格线在2d图上创建
millimeter graphing paper的效果,以显示多变量函数如何依赖于1变量.不同变量的尺度差异很大,所以我的天真方法(我之前使用过)似乎不起作用.
我现在所拥有的例子：
<< ErrorBarPlots`
Cmb[x_,y_,ex_,ey_] := {{N[x],N[y]},ErrorBar[ex,ey]};
SetAttributes[Cmb,Listable];
ELP[x_,ey_,name_] :=
ErrorListPlot[
Cmb[x,y,ex,ey],PlotRange -> FromTo[x,y],PlotLabel -> name,Joined -> True,Frame -> True,GridLines -> GetGrid,ImageSize -> {600}
]
FromTo(我想在框架中留下5％的边距)和GetGrid不能完全按照我的要求工作.
在某些轴上,变量会有很多10的阶数.而且我不希望,一个轴有10个网格线的订单多于其他轴.最重要的是,我希望网格线与刻度线对齐.
样本数据：
ELP[
{4124961/25000000,27573001/100000000,9162729/25000000,44635761/
100000000,15737089/25000000,829921/1562500,4405801/4000000,23068809/25000000,329386201/100000000,58079641/100000000},{1/10,1/5,3/10,2/5,3/5,1/2,1/2},{2031/(250000 Sqrt[10]),5251/(500000 Sqrt[10]),3027/(
250000 Sqrt[10]),6681/(500000 Sqrt[10]),3967/(250000 Sqrt[10]),911/(62500 Sqrt[10]),2099/(100000 Sqrt[10]),4803/(
250000 Sqrt[10]),18149/(500000 Sqrt[10]),7621/(500000 Sqrt[10])},{1/2000,1/1000,3/2000,1/500,3/1000,1/400,1/400},"T2,m"
]
会导致：
而我天真的GetGrid,在某种程度上有效：
FromTo[x_,y_] := Module[{dx,dy},dx = (Max[x] - Min[x])*0.1;
dy = (Max[y] - Min[y])*0.1;
{{Min[x] - dx,Max[x] + dx},{Min[y] - dy,Max[y] + dy}}];
GetGrid[min_,max_] := Module[{step,i},step = (max - min)/100;
Table[
{min + i*step,If[Equal[Mod[i,10],0],Directive[Gray,Thick,Opacity[0.5]],5],Directive[LightGray,Opacity[0.5]]
]]},{i,1,100}]
]
题
如何使GridLines与ticks对齐？
编辑：有
GetTicks[x_,dx = (Max[x] - Min[x])*0.1;
dy = (Max[y] - Min[y])*0.1;
{
Min[x] - dx + Table[i*dx*1.2,9}],Min[y] - dy + Table[i*dy*1.2,9}]
}];
ELP[x_,FrameTicks -> GetTicks[x,ImageSize -> {600},AspectRatio -> 1
]
我可以得到：
而且这要好得多.但我想转移网格而不是刻度线.
编辑：@Sjoerd C. de Vries
您的解决方案可以实现我想要归档和运行的功能.我还注意到,如果我采用样本数据的前5个元素,那么绘图将是(元素被排序并且回归线被添加).
请注意,最左边的元素就像离网格一样.
解决方法
不要使用FrameTicks,而是正确地移动网格.这是第一种方法.晚餐等待.
getGrid[min_,max_] :=
Module[{step,Print[{min,max}];
step = 1/100;
Table[
{
Floor[min,0.1] + i*step,Opacity[0.5]]
]
]
},(Ceiling[max,0.1] - Floor[min,0.1])/step // Round}
]
]
使用适合网格的AspectRatio(可能是x和y范围的比率)
晚餐后更新
为了使它对于不同的值范围(根据您的注释)更加健壮,我生成将由ListPlot选择的刻度并基于我的步骤：
getGrid[min_,i,j},i = Cases[(Ticks /.
AbsoluteOptions[ListPlot[{{min,min},{max,max}}],Ticks])[[1]],{a_,___,{_,AbsoluteThickness[0.25`]}} :> a];
step = i[[2]] - i[[1]];
Table[
{
i[[1]] + j*step/10,If[Equal[Mod[j,{j,10 Length[i]}
]
]
并获得产生方形栅格的宽高比
getAspect[{{minX_,maxX_},{minY_,maxY_}}] :=
Module[{stepx,stepy,rx,ry},i = (Ticks /.AbsoluteOptions[ListPlot[{{minX,minY},{maxX,maxY}}],Ticks]);
rx = Cases[i[[1]],AbsoluteThickness[0.25`]}} :> a];
stepx = rx[[2]] - rx[[1]];
ry = Cases[i[[2]],AbsoluteThickness[0.25`]}} :> a];
stepy = ry[[2]] - ry[[1]];
((maxY - minY)/stepy)/((maxX - minX)/stepx)
]
测试
ELP[x_,name_] :=
ErrorListPlot[Cmb[x,GridLines -> getGrid,PlotRangePadding -> 0,AspectRatio -> getAspect[FromTo[x,y]],y]]
ELP[{4124961/25000000,44635761/100000000,3027/(250000 Sqrt[10]),1/100000 6681/(500000 Sqrt[10]),4803/(250000 Sqrt[10]),m"]
在这里,我将y值除以20,并将x值乘以10000,以显示网格仍然良好：
最后更新(我希望)
这使用FindDivisions作为suggested by belisarius.但是,我按照Margus的要求使用了毫米纸的三级线结构标准：
getGrid[x_,y_] :=
FindDivisions[{x,y},{10,2,5}] /. {r_,s_,t_} :>
Join[
{#,Opacity[0.5]]} & /@ r,{#,Opacity[0.5]]} & /@ Union[Flatten[s]],Opacity[0.5]]} & /@ Union[Flatten[t]]
]
和
getAspect[{{minX_,stepy},stepx = (#[[2]] - #[[1]]) &@FindDivisions[{minX,maxX},10];
stepy = (#[[2]] - #[[1]]) &@FindDivisions[{minY,maxY},10];
((maxY - minY)/stepy)/((maxX - minX)/stepx)
]
警告！！！
我刚注意到,如果你在MMA中有这个：
然后你把它复制到SO(只是ctrl-c ctrl-v),你得到这个：
(maxY - minY)/stepy/(maxX - minX)/stepx
这在数学上是等价的.它应该是这样的：
((maxY - minY)*stepx)/((maxX - minX)*stepy)
