现在的位置: 首页 > 综合 > 正文

Matlab:任意矩阵计算分布密度(海明距离的分布密度)

2013年07月31日 ⁄ 综合 ⁄ 共 518字 ⁄ 字号 评论关闭

如何画分布密度?如下图所示:


x轴是海明距离,y轴是分布密度。我写一个函数,输入为海明距离hd。

function out  = YanMiDu(hd)

xMin = min(hd);

xMax = max(hd);

x = linspace(xMin, xMax, 20); %将最大值 最小值区间分成19个小段

y = hist(hd,x);  %计算各个区间,hd出现的次数

y = y/length(hd); %hd是个行向量,length计算个数.换算成百分比

如果想画bar壮概率密度,则用下面语句:

bar(x,y)

如果画分布函数:

s = 0;

for i=2:length(x)
 s=[s,trapz(x([1:i]),y([1:i]))];
end
figure;
plot(x,s,x,s,'*')

 

如果曲线需要拟合:

a= polyfit(x,y,2);
yy= polyval(a, x);

plot(x, yy); %拟合后的曲线

 

注意:有时候曲线不闭合,就是在x轴最大值的地方,没有闭合。在x = linspace(xMin, xMax, 20);
可以将其改为:x = linspace(xMin, xMax+0.01, 20);

xMax加的那个0.01就是为了让曲线闭合,更好看!当然参数根据自己的需要设。

 

抱歉!评论已关闭.