我已实现了以下功能来估算Parzen Density of a matrix, parzen.m function [retval] = parzen (matrix, dataPoint, variance)
[r c] = size(matrix);
A = ones(r, c)*dataPoint;
sub = matrix - A;
up = sub
我试图从具有给定密度的分布中抽取1000个数字,这个分布的密度是f(x)=(1/3)* X^2 -1 < x < 2.到目前为止,直方图,这显然是错误的,但我认为我正朝着正确的方向前进。任何帮助,我如何解决这个问题将不胜感激。 X=rand(0,1000); Y=rand(1,10000);
a=-1;
b=2;
c=2;
f = (1/3).*X.^2;
for i = 1:leng