How to create a histogram without using the matlab function

Question

Ole 2020-2-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/508164-how-to-create-a-histogram-without-using-the-matlab-function

评论： Ole 2020-3-1

采纳的回答： Guillaume

How to create a histogram without using the matlab hist function.

Given scattered data x = randn(1,100); y = randn(1,100);

with coresponding phase p = randn(1,100), having phase phase = exp(1i*p);

would like to create uniform grid, add the phase for the data points that are inside each bin.

This is to create a na intensity plot or coherent sum of the scattered data.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Guillaume 2020-3-1

@Ole, can you give a formal mathematical definition of what it is you want to calculate for a bin, because as you can see we're a bit confused.

Ole 2020-3-1

There is scatted data (points) in space given as vectors x, y. Let say we discretize the data and in each bin (pixel) fall a set of points {x(k), y(k)}. Each point x(k), y(k) is associated with phase exp(1i*p(k)). In each bin (n,m) I would like to sum the phases for the points that fall inside the bin {x(k), y(k)} -> sum( exp(1i*p(k)) ). So the out put is a matrix M(n,m) that is with reduced size because of the binning and the values of the matrix are M(n,m) = sum( exp(1i*p(k)) ). And would like also to have the number of points x(k), y(k) that are in n,m bin (what histogram does).

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Guillaume 2020-3-1

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/508164-how-to-create-a-histogram-without-using-the-matlab-function#answer_417967

在 MATLAB Online 中打开

As Steven said, use discretize to find the bin indices and then one of the many aggregation functions in matlab. With 3 vector inputs, I'd use the older accumarray:

%demo data and bin definitions
x = randn(1, 100);
y = randn(1, 100);
p = randn(1, 100);
phase = exp(1i*p);
xedges = [-Inf, -3:3, +Inf]; %there will be one less bin that there are edges
yedges = [-Inf, -3:3, +Inf]; %see documentation of discretize
%histogram building
destrow = discretize(y, ybins);
destcol = discretize(x, xbins);
phasehistogram = accummaray([destrow, destcol], phase, [numel(ybins), numel(xbins)] - 1);

Or you could put the vectors in a table, and call groupsummary which would do the binning and summing for you:

%demo data and bin definitions
x = randn(1, 100);
y = randn(1, 100);
p = randn(1, 100);
phase = exp(1i*p);
xedges = [-Inf, -3:3, +Inf]; %there will be one less bin that there are edges
yedges = [-Inf, -3:3, +Inf]; %see documentation of discretize
%table construction and histogram:
phasetable = table(x, y, phase);
phasehistogram = groupsummary(phasetable, {'x', 'y'}, {xedges, yedges}, 'sum', 'phase');