How to create a histogram without using the matlab function
6 views (last 30 days)
Show older comments
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 Comments
Guillaume
on 1 Mar 2020
@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.
Accepted Answer
Guillaume
on 1 Mar 2020
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');
More Answers (1)
Steven Lord
on 28 Feb 2020
Consider using discretize to bin the data then passing that grouping information into groupsummary or splitapply.
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!