转载来源
特别感谢
对于该参数的含义询问了ZAH_zah,十分感谢!
函数
- [modes its]=ceemdan(x,Nstd,NR,MaxIter)
输入参数及含义
-
Nstd = 0.2;
- 信噪比,作用为了抑制混叠效应。
- 噪声权重(范围:0-inf,一般选择0-1)
-
NR = 500;
- 噪声添加次数,进行500次的噪声添加
- 一般选择50-100
-
MaxIter = 5000;
- 对emd 内部最大包络次数设定。
- 分解出分量的最大个数(假设ceemdan本来分解出10个分量,你将该值设置为5:那么最终就会获得5个分量。如果你将该值设置为20:最终还是会获得10个分量)
- 该值一般设置为inf或一个很大的数。
- 如果你要求CEEMDAN每次只获得一个分量,那么就可以将该参数设置为1。不过我们一般不会认为去阻止CEEMDAN的分解过程,所以将该参数设置为inf就行。
(该部分实验后有点小问题,可能解读不太准确,大家可以自行设置值进行实验)
输出参数
[modes its]=ceemdan(ecg,0.2,500,5000);
- 其中modes 为IMF分量,最后一列为残余剩量
- its为内部包络迭代的次数,因此最大不会超过MaxIter(这里为5000)
盒图介绍
-
盒图是在1977年由美国的统计学家约翰·图基(John Tukey)发明的。
-
它由五个数值点组成:最小值(min),下四分位数(Q1),中位数(median),上四分位数(Q3),最大值(max)。也可以往盒图里面加入平均值(mean)。
-
如上图。下四分位数、中位数、上四分位数组成一个“带有隔间的盒子”。
-
上四分位数到最大值之间建立一条延伸线,这个延伸线成为**“胡须(whisker)”**。
-
由于现实数据中总是存在各式各样地“脏数据”,也成为“离群点”,于是为了不因这些少数的离群数据导致整体特征的偏移,将这些离群点单独汇出,而盒图中的胡须的两级修改成最小观测值与最大观测值。
-
这里有个经验,就是最大(最小)观测值设置为与四分位数值间距离为1.5个IQR(中间四分位数极差)。即:
- IQR = Q3-Q1,即上四分位数与下四分位数之间的差,也就是盒子的长度。
- 最小观测值为min = Q1 -1.5*IQR,如果存在离群点小于最小观测值,则胡须下限为最小观测值,离群点单独以点汇出。如果没有比最小观测值小的数,则胡须下限为最小值。
- 最大观测值为max = Q3 +1.5*IQR,如果存在离群点大于最大观测值,则胡须上限为最大观测值,离群点单独以点汇出。如果没有比最大观测值大的数,则胡须上限为最大值。
盒图的优势:
- 盒图能够有效地帮助我们识别数据的特征:
- 直观地识别数据集中的异常值(查看离群点)。
- 判断数据集的数据离散程度和偏向(观察盒子的长度,上下隔间的形状,以及胡须的长度)。
matlab盒图函数
-
boxplot(x) :创建 x 中数据的箱线图
- 如果 x 是向量,boxplot 绘制一个箱子。如果 x 是矩阵,boxplot 为 x的每列绘制一个箱子。在每个箱子上,中心标记表示中位数,箱子的底边和顶边分别表示第 25 个和 75 个百分位数。须线会延伸到不是离群值的最远端数据点,离群值会以 ‘+’ 符号单独绘制。
-
boxplot(x,g):
- 使用 g 中包含的一个或多个分组变量创建箱线图。boxplot 为具有相同的一个或多个 g 值的各组 x 值创建一个单独的箱子。
-
boxplot(ax,___) :
- 使用坐标区图形对象 ax 指定的坐标区和任何上述语法创建箱线图。
-
boxplot(___,Name,Value) :
- 使用由一个或多个 Name,Value对组参数指定的附加选项创建箱线图。例如,您可以指定箱子样式或顺序。
参考文献
MATLAB中基于CEEMDAN与FastICA算法的单通道信号规律噪声的滤除方法
盒图(boxplot)
matlab boxplot
ceemd的参数怎么设置,有什么要求吗