How to generate vector with non-linear spacing?

17 次查看(过去 30 天)
Hello,
I am attempting to generate an array with non-linear spacing between the elements but struggling to figure out a way to impliment this. I need an array x of points spaced from a real, negative value a to a real, positive value b. Ideally the spacing should follow a Gaussian distribution such that many points are clustered about x=0 in the vicinity of some finite width d. Conversely the spacing should be large at the boundaries. More generally it doesn't have to be Gaussian but some guidance on a scheme where I can specify a length d<<|b-a| across which many points should be distributed would be really appreciated.
I tried to work out an algorithm of sorts from pen and paper but it didn't really work and resulted in something that did not appear to be immediately solveable.
Thanks

采纳的回答

John D'Errico
John D'Errico 2024-8-14
编辑:John D'Errico 2024-8-14
It sounds like you want a deterministic, non-random set of points. You can use the inverse gaiussian CDF to do what you want.
a = -2.5;
b = 3;
N = 40;
x = norminv(linspace(normcdf(a),normcdf(b),N));
xline(x,'r')
So most dense at zero, and a density that drops off away from there. I chose limits that are not centered around zero, but had they been centered, the set would be perfectly symmetric.
Of course, you can adjust the speed of fall off, by use of a different variance on the implicit Gaussian.

更多回答(2 个)

Matt J
Matt J 2024-8-14
编辑:Matt J 2024-8-14
a=-1;
b=+1;
d=1/(b-a)^2;
y=d*randn(1,500);
x=sort(y+(a+b)/2);
stem(x,x.^0); xlim([a,b])
  2 个评论
Ryan
Ryan 2024-8-14
Thanks for the suggestion, this is a possible option. Any way to get the spacing more regular?
Ryan
Ryan 2024-8-14
Ah I see what I meant by Gaussian is not in the sense of a random distribution but of a spacing that followed a pattern akin to a contiunous Gaussian function. I had envisaged an algorithm that produces regular values
Where Δis some constant. I don't think this algorithm "works" as such but I think it illustrates what I'm after. For what I am working on a random distribution may also work.

请先登录,再进行评论。


Aquatris
Aquatris 2024-8-14
The answer from another question here:
A = 1;
B = 5;
std = .9; % need to play with this to get desired shape within desired range,
% if too high, it does not become a normal distribution
% if too low, it does not cover the whole A-B range
mean = 3;
step1Result = mean + std * randn(1e3,1); % create random data with mean and std
step2Result = min(max(step1Result,A),B); % make sure they end up in A-B range
hist(step2Result)
ylabel('# of occurence')
xlabel('Data Value')

类别

Help CenterFile Exchange 中查找有关 Descriptive Statistics 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by