Pearson Johnson , pearsrnd johnsrnd ( , )
:
%
data = load('data.dat');
%
moments = {mean(data),std(data),skewness(data),kurtosis(data)};
newData = pearsrnd(moments{:}, [4096-499 1]);
%
extendedData = [data; newData];
%
subplot(121), hist(data), xlabel('x'), ylabel('Frequency')
subplot(122), hist(extendedData), xlabel('x'), ylabel('Frequency')
johnsrnd:
%
quantiles = quantile(data, normcdf([-1.5 -0.5 0.5 1.5]));
newData = johnsrnd(quantiles, [4096-499 1]);
, , ecdf ksdensity.
, ( !).