这是迭代的二维最小熵反褶积实现根据迭代方法在原始的纸:
R.A.韦根,最小熵反褶积,地质勘探,16卷,爱思唯尔科学出版,阿姆斯特丹,1978年。21-35页。
更新到2015年包括卷积调整修复。提出了
在第二个引用论文。这是重要的修复地中海
deconvolving的平凡解卷积不连续面:
overlapMode = '有效':使用卷积调整修复;
overlapMode =“满”:使用原始R.A.•威金斯提出的卷积定义。我不推荐这个。
读者注意,你可能想要把我的一些其他MED-based提交:
此外:
这是一个类似的反褶积问题造成卡洛斯Cabrelli基于地中海
有价值的解决方案直接为最优反褶积过滤器。我建议
检查出来如果你对医学感兴趣,然而往往不是最佳的解决方案
预期达成的地中海是一个更好的结果,因为此外能够更好的实现
一个冲动的最优解。
MOMEDA:
这是最佳的解决方案为旋转周期性的冲动和建议
机器故障而不是地中海或此外。因为它是价值,它可以
快速生成光谱诊断机器的健康。

算法参考:
R.A.韦根,最小熵反褶积,地质勘探,卷。
16日,爱思唯尔科学出版,阿姆斯特丹,1978年。21-35页。

卷积的调整:
基准线麦当劳,<人>,多点最优最小熵反褶积和卷积
解决办法:应用振动故障检测,未出版

输入:
x:
信号进行最小熵反褶积。如果一个单一的
列/行指定的数据,1 d滤波器设计
最小熵产生的信号。如果一个2 d数据
指定矩阵,一个一维滤波器设计
每一列的平均熵最小化的过滤
数据。

filterSize:
这是有限的长度inpulse过滤器过滤
设计。使用价值大约30是适当的根据
数据。调查使用性能差异
不同的值。

白蚁:(可选)
这是终止的迭代次数。如果
迭代的数量超过这个数字,地中海的过程
将完成。指定要使用[]30的默认值。

termDelta:(可选)
这是终止条件。如果峰态的变化
低于这个阈值,迭代之间的迭代
过程将终止。使用默认值指定[]
0.01。您可以指定一个值为0的只有终止
白蚁的条件,即。执行的确切数字
迭代。

overlapMode:(可选)
你应该总是使用这个参数包括“有效”
卷积修复纠正错误deconvolving地中海
虚假的冲动。看到算法参考部分
卷积的细节调整。您可以使用“全部”如果你想要的
复制原始的医疗效果,但不建议
上面的原因。

plotMode:
如果这个值是> 0,情节将生成的迭代
性能和产生的信号。

输出:
y_final:
输入信号(s) x,由地中海造成过滤器过滤。
这是获得简单:y_final =过滤器(f_final 1 x);

f_final:
最后1 d地中海有限脉冲响应滤波器的格式。

kurtIter:
地中海峰度根据迭代。kurtIter(结束)
最终峰度,ie。每个y_final总结的峰度
y_final列。总和(峰度(每一列y_final))

例子:
% - - - - - - - - - - - - - - - - 1 d反褶积例子
n = 0:999;
x =[罪(n / 30) + 0.2 * (mod (n, 21) = = 0));
[y_final f_final kurt] = med2d (x ', 30100年,[],'有效',1);

% - - - - - - - - - - - - - - - - 2 d反褶积例子
%这将主要提取impulse-like
%扰动引起的0.2 * (mod (n, 21) = = 0)
%和阴谋的结果。
n = 0:999;
x =[罪(n / 30) + 0.2 * (mod (n, 21) = = 0);
罪(n / 13) + 0.2 * (mod (n, 21) = = 0)];
[y_final f_final kurt] = med2d (x ', 30100年,[],'有效',1);

注意:
解决方案是不能保证最优解决方案
熵minimizataion问题,解决方案只是一个地方
最小的熵,因此一个好的选择。