小编Mar*_*lig的帖子

我经常不得不同时可视化多个数据集,一般在ListPlot或它的日志 - 同伴.由于数据集的数量通常大于容易区分的线条样式的数量,并且创建大型情节传说仍然有点不明显,我仍然在寻找一种在我的图中注释不同线条/集合的好方法.提示在屏幕上工作时不错,但如果我需要pritn情节他们不帮助.

最近,我使用Mesh选项来枚举我的数据集,并发现了一些奇怪的东西

GraphicsGrid[Partition[Table[ListPlot[
Transpose@
 Table[{Sin[x], Cos[x], Tan[x], Cot[x]}, {x, 0.01, 10, 0.1}], 
PlotMarkers -> {"1", "2", "3", "4"}, Mesh -> i, Joined -> True, 
PlotLabel -> "Mesh\[Rule]" <> ToString[i], ImageSize -> 180], {i, 
1, 30}], 4]]

结果在我的机器上看起来像这样(Windows 7 x64,Mathematica 8.0.1):

有趣的是,对于Mesh-> 2,8和10,结果看起来像我预期的,其余的没有.要么我不理解Mesh选项,要么它不理解我.

这是我的问题:

1. ListPLot中的Mesh是错误的还是我错误地使用它？
2. 我怎样才能对连续集的网格点进行x移位以避免叠印？
3. 你对如何在一个情节中注释/枚举多个数据集有任何其他建议吗？

有没有人有将数据导出为FITS具有自定义元数据(FITS标题)信息的文件的经验？到目前为止,我只能FITS使用标准的Mathematica FITS标头模板生成文件.该文档没有提供是否支持自定义元数据导出以及如何完成它的提示.

comp.soft-sys.math.mathematica的以下建议不起作用:

header=Import[<some FITS file>, "Metadata"]; 
Export<"test.fits",data ,"Metadata"->header] 

要么

Export["test.fits",{"Data"->data,"Metadata"->header}]

将我自己的元数据导出到FITS文件的正确方法是什么？

干杯,

马库斯

更新:来自Wolfram支持的回复:"Mathematica尚不支持导出 FITS文件的元数据.该示例指的是导入此数据.我们计划在未来支持此...""还有计划包括二进制表到FITS导入功能."

我会尝试提出一些解决方法.

我试着整合以下功能:

test[a_] := 
 Module[{dnda, cg, atg, acg, alphag, betag, sig, b1, b2, dndavsg, a01,
 a02, bc5},
 alphag = -1.96;
 betag = -0.813;
 atg = 6.93*^-6;
 acg = 0.000348;
 cg = 2.95*^-13;
 dnda = (cg/a) (a/atg)^alphag;
 If[betag >= 0, 
       dnda = dnda (1 + betag a/atg), 
       dnda = dnda/(1 - betag a/atg)
 ];
 If[a > atg, dnda = dnda Exp[((atg - a)/acg)^3]];
 a01 = 3.5 10^-8;
 a02 = 3 10^-7;
 sig = 0.4;
 b1 = 2.0496 10^-7;
 b2 = 9.6005 10^-11; …