|
用MATLAB求分布
Matlab可以求分布的数字解,假设数据为一个向量A,假设样本量为n。
首先确定数据的取值范围,按照一定的规则扩展原数据的上下界,例如:X_min=min(A), X_max=max(A)。
然后,确定分点向量;例如:Number_Division=100; Vector_Division = rand(Number_Division,1); Vector_D = X_min + Vector_Division*(X_max -X_min );
最后,确定原数据小于各个分点值的频率。这里有个加速运算的技巧。
Matrix_X = repmat(X,1,n);
%% 将数据向量复制,每个列数据相同。
Matrix_D = repmat(transpose(Vector_D),n,1);
%% 将分点向量复制,每个行数据相同。
Matrix_Delta = Matrix_X - Matrix_D;
%% 每个列的Delta都是原始数据与某个分点的差值。
Matrix_Signs = [Matrix_Delta <0];
Vector_Signs_Negative = sum(Matrix_Signs);
%% 每个列的Signs值都是原始数据小于某个分点的个数,刚好能求出频率。
Vector_Distribution = rdivide(Vector_Signs_Negative,n);
%%每个列表示小于某个分点的频率,结果为行向量。 |
|