[已答复] 直方图拟合

[复制链接]
tkl1050 发表于 2010-5-3 17:03:48
x_chuli.txt (98.11 KB, 下载次数: 53239)

17 条回复


ctgu000 发表于 2010-5-3 17:32:38

为什么要用拟合呢?

plot就可以画直方图的曲线图
clear all
clc
I=imread('11112.jpg');
I=rgb2gray(I);
I=double(I)
Value=unique(I(:))
Count=[hist(I(:),Value)]'
plot(Value,Count) %曲线图
figure,imhist(uint8(I))%直方

[ 本帖最后由 ctgu000 于 2010-5-3 17:34 编辑 ]

tkl1050 发表于 2010-5-3 18:50:49
谢谢你的回复,我拟合的目的是为了找出其概率密度函数!

谢中华 发表于 2010-5-3 19:39:17
你的数据服从尖峰厚尾的分布,常用分布中没有合适的分布能够拟合你的数据,建议进行核密度估计,请查阅相关资料。绘制频率直方图、进行正态分布拟合、进行核密度估计的matlab代码如下:
  1. % 导入数据
  2. load x_chuli.txt
  3. % 调用ecdf函数计算xc处的经验分布函数值f_ecdf
  4. [f_ecdf, xc] = ecdf(x_chuli);
  5. % 新建图形窗口,然后绘制频率直方图,直方图对应7个小区间
  6. figure;
  7. ecdfhist(f_ecdf, xc, 97);
  8. hold on;
  9. xlabel('x');  % 为X轴加标签
  10. ylabel('f(x)');  % 为Y轴加标签

  11. % 调用ksdensity函数进行核密度估计
  12. [f_ks1,xi1,u1] = ksdensity(x_chuli);
  13. % 绘制核密度估计图,并设置线条为黑色实线,线宽为3
  14. plot(xi1,f_ks1,'r','linewidth',2)

  15. ms = mean(x_chuli);  % 计算x的平均值
  16. ss = std(x_chuli);  % 计算x的标准差
  17. % 计算xi1处的正态分布密度函数值,正态分布的均值为ms,标准差为ss
  18. f_norm = normpdf(xi1,ms,ss);
  19. % 绘制正态分布密度函数图,并设置线条为红色点划线,线宽为3
  20. plot(xi1,f_norm,'k-.','linewidth',2)

  21. % 为图形加标注框,标注框的位置在坐标系的左上角
  22. legend('频率直方图','核密度估计图', '正态分布密度图', 'Location','NorthWest')
复制代码

绘制的图形如下:
untitled.jpg
可见正态分布无法拟合你的数据,核密度函数还是不错的。

tkl1050 发表于 2010-5-4 14:19:03
谢谢,谢老师的回复,还想请教一个问题,核密度估计能得到具体的函数表达式吗?

谢中华 发表于 2010-5-4 19:46:26
核密度估计能得到函数表达式,是一个求和的形式,自己查阅资料吧。

hit218 发表于 2012-5-17 20:46:17
很好的回复,QQQ

venmy 发表于 2012-6-6 16:20:59
这样的帖子值得一看

venmy 发表于 2012-6-6 16:21:36
这样的帖子值得一看

lixiangtnt 发表于 2012-8-17 10:11:48
谢谢 老师函数指教 很受用

小晴子182 发表于 2013-11-22 17:12:59
看不懂啊!

iprize 发表于 2015-12-29 17:26:44
频率直方图还能大于1?

非法智慧 发表于 2016-3-15 21:09:05
谢中华 发表于 2010-5-3 19:39
你的数据服从尖峰厚尾的分布,常用分布中没有合适的分布能够拟合你的数据,建议进行核密度估计,请查阅相关 ...

大神,能不能帮我看看EM算法求混合高斯模型的算法?
我有两个1*N(N>80000)的矩阵,需要拟合在一起~求高斯混合模型

孚平之水 发表于 2016-10-12 11:23:59
hit218 发表于 2012-5-17 20:46
很好的回复,QQQ

matlab可以得出概率密度函数吗?

XCc0900720130 发表于 2016-11-26 16:22:41
表示感谢

章鱼小丸子 发表于 2016-11-28 17:14:04
请问老师,用ecdfhist将多组数据的频率直方图画在同一axes时,如何让不同数据组有不同颜色?

XCc0900720130 发表于 2016-12-5 13:40:34
谢老师的回复很受用。

ershuai 发表于 2017-8-16 19:45:30
谢谢楼主,很受用,学到了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表