这是一个使用Nic Price实施Hartigan’s Dip测试来识别单峰分布的脚本.棘手的一点是计算xpdf,它不是概率密度函数,而是一个有序的样本.

p_value是获得测试统计量的概率,至少与实际观察到的一样极端,假设零假设为真.在这种情况下,零假设是分布是单峰的.

close all; clear all;

function [x2, n, b] = compute_xpdf(x)

x2 = reshape(x, 1, prod(size(x)));

[n, b] = hist(x2, 40);

% This is definitely not probability density function

x2 = sort(x2);

% downsampling to speed up computations

x2 = interp1 (1:length(x2), x2, 1:1000:length(x2));

end

nboot = 500;

sample_size = [256 256];

% Unimodal

sample2d = normrnd(0.0, 10.0, sample_size);

[xpdf, n, b] = compute_xpdf(sample2d);

[dip, p_value, xlow, xup] = HartigansDipSignifTest(xpdf, nboot);

figure;

subplot(1,2,1);

bar(n, b)

title(sprintf('Probability of unimodal %.2f', p_value))

% Bimodal

sample2d = sign(sample2d) .* (abs(sample2d) .^ 0.5);

[xpdf, n, b] = compute_xpdf(sample2d);

[dip, p_value, xlow, xup] = HartigansDipSignifTest(xpdf, nboot);

subplot(1,2,2);

bar(n, b)

title(sprintf('Probability of unimodal %.2f', p_value))

print -dpng modality.png

matlab双峰滤波,MATLAB中的单峰或双峰分布相关推荐

  1. matlab实现中值滤波程序,中值滤波流程(matlab平滑滤波和中值滤波程序)

    matlab平滑滤波和中值滤波程序 glRasterPos2i(100,100); //定位当前光标 glutBitmapCharacter(GLUT_BITMAP_9_BY_15,'H'); //写 ...

  2. matlab soble滤波,MATLAB 图像滤波去噪分析及其应用

    MATLAB 图像滤波去噪分析及其应用 第1章 图像颜色空间相互转换与MATLAB实现 1.1 图像颜色空间原理 1.1.1 RGB颜色空间 1.1.2 YCbCr颜色空间 1.1.3 YUV颜色空间 ...

  3. kalman滤波的matlab,kalman滤波matlab实现

    <kalman滤波matlab实现>由会员分享,可在线阅读,更多相关<kalman滤波matlab实现(8页珍藏版)>请在人人文库网上搜索. 1.function S = do ...

  4. matlab高斯滤波跟中值滤波区别,matlab图像滤波处理代码讲解 均值滤波 中值滤波 高斯滤波...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 (1)均值滤波 function imres = averageFilter(im,N) % imres为均值滤波后的矩阵,参数im为图像的路径,N为滤波 ...

  5. matlab均值滤波,中值滤波

    滤波核3×3的均值滤波. F=imread('D:\Desktop\66.png'); F=rgb2gray(F); f=double(F); [row,col]=size(f); f1=zeros( ...

  6. matlab空域滤波,MATLAB数字图像处理基本操作及空域滤波

    MATLAB数字图像处理基本操作及空域滤波 实验一 MATLAB 数字图像处理基本操作及空域滤波实验目的1.了解有关数字图像处理的基本概念,熟悉Matlab软件中关于数字图像处理的基本命令,掌握利用M ...

  7. matlab加权滤波,matlab实现七种滤波方法

    对于信号二,小波的去噪效果非常不错,虽然得到波形不是很平滑,但是上升沿和下降沿保持的非常高,基本可以看到棱角. %******************************************* ...

  8. 数字图像处理 中值滤波 MATLAB实验

    一.原理_中值滤波 中值滤波的基本思想是将图像中每个像素的灰度值用其邻域内像素灰度的中值代替,它是一种非线性平滑滤波算法. 设加噪图像为 f(x,y) ,经中值滤波处理后的图像为g(x,y) ,则: ...

  9. 均值滤波、中值滤波【MATLAB】【图像处理】

    1. 均值滤波:即线性滤波,有:  a.邻域平均, b. 加权平均两种. a. 邻域平均:用一个像素邻域平均值作为滤波结果,此时滤波器模板的所系数都取1,为保证输出图仍在原来的灰度值范围,在算得R后要 ...

  10. 中值滤波matlab

    摘  要 中值滤波技术是一种在去除噪声的同时能较好保护图像边缘细节的非线性技术,在图像增强和恢复等领域中得到了广泛的应用.文章阐述了中值滤波的原理和特点,并使用软件工具MATLAB实现了图像的中值滤波 ...

最新文章

  1. 端口保护:如果你不把我当回事,我就会让你好看
  2. 新款iPhone,比绿更绿
  3. 设计模式6---代理模式(Proxy Pattern)
  4. springboot报错---@RunWith(SpringRunner.class)
  5. 【SqlServer】Sqlserver中的DOS命令操作
  6. 三星手机 java_如何在三星手机上安装Java ME应用程序?
  7. 起止时间控制,显示格式控制
  8. spark sql建表的异常
  9. html5之Canvas坐标变换应用-时钟实例
  10. docker下nginx反向代理和负载均衡配置
  11. 对《致加西亚的信》的异议
  12. 配置Atari Gym环境
  13. ABP框架—后台:后台程序发布至IIS(6)
  14. 以开发之名|线上家装新美学——梦想之家,由你来定
  15. 微信小程序实战(仿小米商城)
  16. java工具包Lombok
  17. Python爬虫|Python爬虫入门:请求
  18. python绘制语谱图(详细注释)
  19. 用户管理和文件权限部分练习题
  20. RIoTBoard开发板系列笔记(三)—— 移植Gstreamer

热门文章

  1. 900行c语言贪吃蛇,挑战CMD版贪吃蛇,能超过我的分算你赢
  2. r语言boxcox异方差_R教程-15:线性回归中的异方差
  3. 数字单位转换(例:元->万元)
  4. 大学计算机在线作业答案,大学计算机第一次在线作业及答案
  5. 云服务器功能检验方法
  6. 淘宝APP用户行为数据分析案例——Python
  7. 彻底解决乱码问题(一):为何会出现乱码
  8. 复联4里的钢铁侠原型是谁?(剧透警告慎入)
  9. ETL数据处理平台,快速实现数据集成
  10. 财务记账微信小程序开发