正态云公式matlab代码,云模型代码 - MATLAB等数学软件专版 - 经管之家(原人大经济论坛)...

本人用两种方法生成正态分布随机数，并产生云滴绘图。但一种离散程度很高，一种基本不离散，希望找出原因！！！源代码如下：

load G1.txt;    %%% 矩阵G1为一个评语的取值范围(一列)。

[m,n]=size(G1);

E=ones(1,1);   %%% 矩阵E用于存放熵En

S=ones(1,1);   %%% 矩阵S用于存放样本方差(二阶平方距)

H=ones(1,1);   %%% 矩阵H用于存放超熵He

F=mean(G1,1);   %%% 矩阵F用于存放各项指标评价结果的期望值Ex

E=(pi/2)^(1/2)*(1/m);

S=(1/(m-1));                %%%%%  二阶平方距

H=(abs(S-E^2))^(1/2);

F                  %%%%%%% 某个评语的期望值Ex

E                  %%%%%%% 某个评语的熵En

H                  %%%%%%% 某个评语的超熵He

n=1000;        %%% 设定随机数数目

x=zeros(1,n);

y=zeros(1,n);

f=F*ones(1,n);

for i=1:n

x(i)=normrnd(E,H);

Enn(i)=x(i);

end

hold on        %%% 绘制此评语的云图

for i=1:n

x(i)=normrnd(F,Enn(i));

y(i)=exp(-(x(i)-f).^2/(2*Enn.^2));

plot(x(i),y(i),'*')

end

load G1.txt;    %%% 矩阵G1为一个评语的取值范围(一列)。

[m,n]=size(G1);

E=ones(1,1);   %%% 矩阵E用于存放熵En

S=ones(1,1);   %%% 矩阵S用于存放样本方差(二阶平方距)

H=ones(1,1);   %%% 矩阵H用于存放超熵He

F=mean(G1,1);   %%% 矩阵F用于存放各项指标评价结果的期望值Ex

E=(pi/2)^(1/2)*(1/m);

S=(1/(m-1));                %%%%%  二阶平方距

H=(abs(S-E^2))^(1/2);

F                  %%%%%%% 某个评语的期望值Ex

E                  %%%%%%% 某个评语的熵En

H                  %%%%%%% 某个评语的超熵He

hold on

for i=1:1000

Enn=randn(1)*He+En;

x(i)=randn(1)*Enn+Ex;

y(i)=exp(-(x(i)-Ex)^2/(2*Enn^2));

plot(x(i),y(i),'*')

end