Matlab - 浓度与距离图（2D）

Question

我有几组数据。每个集合都是一个数字列表，它是一个粒子从 0 开始行进的距离。每个集合都与有限时间相关联，因此集合 1 是 T=0 时的距离；第 2 组是 T=1 时的距离，依此类推。每组大小为粒子总数，每组大小相同。

我想绘制一条浓度与距离的线。

例如，如果有 1000 个粒子（集合的大小）；在时间 T=0 时，绘图将只是一条直线 x=0，因为所有粒子都为 0（该集合包含 1000 个零）。所以 x=0 处的浓度 =100%，其他距离处的浓度为 0%

在 T=1 和 T=2 以此类推，距离会增加（通常），所以我可能有这样的集合：（只是一个例子）

T1 = (1.1,2.2,3.0,1.2,3.2,2.3,1.4...) 等 T2 = (2.9,3.2,2.6,4.5,4.3,1.4,5.8...) 等

很可能每个集合中的每个数字在该集合中都是唯一的

我们的目标是绘制几张图（我最终可以将它们绘制在一张图上），在 y 轴上显示浓度，在 x 轴上显示距离。我想随着 T 增加 T0、T1、T2，那么图会变平，直到各处的浓度大致相同。

x 轴（距离）有一个固定的最大值，每个绘图都相同。因此，例如，某些集合将有一条曲线在 y 轴（浓度）上以较低的 x（距离）值达到零，但随着时间的增加，我设想一条几乎平坦的线，该线不与x轴（浓度处处非零）

我已经用直方图进行了尝试，但它并没有真正给出我想要的结果。我想要一个线图，但必须尝试将距离放入常识大小的箱子中。

谢谢你

一些粗略的数据

Y1 = 1.0e-09 * [0.3358, 0.3316, 0.3312, 0.3223, 0.2888, 0.2789, 0.2702,...
    0.2114, 0.1919, 0.1743, 0.1738, 0.1702, 0.0599, 0.0003, 0, 0, 0, 0, 0, 0];

Y2 = 1.0e-08 * [0.4566, 0.4130, 0.3439, 0.3160, 0.3138, 0.2507, 0.2483,...
    0.1714, 0.1371, 0.1039, 0.0918, 0.0636, 0.0502, 0.0399, 0.0350, 0.0182,...
    0.0010, 0, 0, 0];

Y3 = 1.0e-07 * [0.2698, 0.2671, 0.2358, 0.2250, 0.2232, 0.1836, 0.1784,...
    0.1690, 0.1616, 0.1567, 0.1104, 0.0949, 0.0834, 0.0798, 0.0479, 0.0296,...
    0.0197, 0.0188, 0.0173, 0.0029];

这些数据集仅包含 20 个粒子的距离。 Y0 集为零。我要处理数千个，所以数据集太大了。

谢谢

Answer 1

嗯，基本上，你只是错过了hold 命令。但首先，将所有数据放在一个矩阵中，如下所示：

Y = [1.0e-09 * [0.3358, 0.3316, 0.3312, 0.3223, 0.2888, 0.2789, 0.2702,...
    0.2114, 0.1919, 0.1743, 0.1738, 0.1702, 0.0599, 0.0003, 0, 0, 0, 0, 0, 0];
    1.0e-08 * [0.4566, 0.4130, 0.3439, 0.3160, 0.3138, 0.2507, 0.2483,...
    0.1714, 0.1371, 0.1039, 0.0918, 0.0636, 0.0502, 0.0399, 0.0350, 0.0182,...
    0.0010, 0, 0, 0];
    1.0e-07 * [0.2698, 0.2671, 0.2358, 0.2250, 0.2232, 0.1836, 0.1784,...
    0.1690, 0.1616, 0.1567, 0.1104, 0.0949, 0.0834, 0.0798, 0.0479, 0.0296,...
    0.0197, 0.0188, 0.0173, 0.0029]];

然后你需要分别绘制每个时间步，并使用hold on将它们粘贴到相同的轴上：

hold on
for r = size(Y,1):-1:1
    histogram(Y(r,:));
end
hold off
T_names = [repmat('T',size(Y,1),1) num2str((size(Y,1):-1:1).')];
legend(T_names)

这将为您提供（使用示例数据）：

请注意，在循环中我向后迭代行 - 这只是为了让较窄的直方图绘制在较宽的位置上，这样您就可以清楚地看到所有这些。

编辑

如果您想要连续线而不是 bin，您必须首先通过 histcounts 获取直方图值，然后将它们绘制成一条线：

hold on
for r = 1:size(Y,1)
    [H,E] = histcounts(Y(r,:));
    plot(E,[H(1) H])
end
hold off
T_names = [repmat('T',size(Y,1),1) num2str((1:size(Y,1)).')];
legend(T_names)

虽然您的小示例数据看起来并不那么令人印象深刻：