你好,请问斑竹有关于EMD边缘效应问题的MATLAB源程序么,特别是镜象法.我的邮箱是tjf741@sina.com,非常感谢!

rilling的EMD程序里面就有啊 就是用镜像法处理的边界问题
http://zhlong.ys168.com 这是zhlong版主收集的一些资料，2005版的程序里面的emd程序里已经有了镜像法的内容

[ 本帖最后由 zhangnan3509 于 2007-12-13 20:16 编辑 ]

谁有能处理图象这种二维数据的EMD程序啊?
mihaohao@126.com

Image analysis by bidimensional empirical mode decomposition


Bidimensional Empirical Mode Decomposition Modified for Texture Analysis



哪位大侠有的话麻烦传我邮箱 mihaohao@126.com 或给个连接.谢谢!

在文献中有看到对imf分量进行阀值处理，请有做过次方面的朋友，发个程序参考一下，我的邮箱zhuokangzhi@126.com，十分感谢！

我用小波包对一个非平稳信号进行分解，找到了最优树结构，我画出了在小波基下，和小波包最优基下的时频面，得到的图，横坐标都不是时域的时间？纵坐标又说明了什么呢？请高手指点！！！！！！

程序和文件如下：

数据文件：

   因为文件很大，所以分两部分。下载后把1和2 的数据合起来就好了

麻烦哪位前辈发给我一份齿轮故障数据（格式是txt的）。最好这个数据已能成功在MATLAB或LABVIEW下用EMD分析。同时请说明大小齿轮齿数，哪个是故障齿轮及其转速，采样频率和采样点数，谢谢！
200431145@163.com

看了Hilbert谱的一些文章，很奇怪Hilbert谱表现形式有很多，二维三维的都有，法国人的那个程序试了下，是二维的点的图像，弱弱的问句，怎样得到二维直线的Hilbert谱，三维的呢，谢谢

输入为y＝3*sin(2*pi*20*t)+2*sin(2*pi*40*t)+0.5*sin(2*pi*80*t)的信号，要求设计：
1能滤除20hz40hz的低通滤波器
2能滤除80hz 的高通滤波器
3能滤除20hz、80hz的带通滤波器
4能滤除40hz的带阻滤波器。
小弟才疏学浅，望各位大侠帮助，滤波器的设计参数不限，只要能达到这些功能就行，在下在此先谢谢了！！
联系方式：邮箱：qchyde◎163.com  QQ:462752320

引用:

原帖由 小光仔 于 2008-4-29 17:04 发表


关于fir滤波器的设计的matlab程序输入为y＝3*sin(2*pi*20*t)+2*sin(2*pi*40*t)+0.5*sin(2*pi*80*t)的信号，要求设计：
1能滤除20hz40hz的低通滤波器
2能滤除80hz 的高通滤波器
3能滤除20hz、80hz的带通滤波器
4能滤除40hz的带阻滤波器。
小弟才疏学浅，望各位大侠帮助，滤波器的设计参数不限，只要能达到这些功能就行，在下在此先谢谢了！！

实际设计这种数字滤波器并不难，这在这里只给出（3）的设计方法，其它几种可以仿照设计。
fs=300;
N=256;
n=1:N;
t=(n-1)/fs;
x1=3*sin(2*pi*20*t);
x2=2*sin(2*pi*40*t);
x3=0.5*sin(2*pi*80*t);
x=x1+x2+x3;
subplot 211; plot(t,x1,'r',t,x2,'g',t,x3,'b',t,x,'k')
legend('x1','x2','x3','x'); axis tight
b3=fir1(51,[30,70]/150);
y3=filtfilt(b3,[1],x);
%freqz(b3,1);
%break
subplot 212; plot(t,x2,'r','linewidth',2); hold on
plot(t,y3); hold off;
legend('x2','y3'); axis tight
原始信号和滤波结果如图，因为按（3）的要求滤波后只留下40hz的信号，把40hz信号（红线）和滤波器输出（蓝线）比较，可看出完全重合。

非常感谢你的答复，这里有两个小问题：
1.你给出的是带通滤波器，我修改你程序上哪几个参数可以变成要求的低通，高通，带阻滤波器？
2.我是否简单地修改你设定地频率范围的话就能对别的范围的信号也能进行滤波，我还要做个对一个有10hz、80hz、与160hz的信号作分析。
看见本贴的话望大侠帮助！！谢谢！

引用:

原帖由 小光仔 于 2008-5-2 22:23 发表
1.你给出的是带通滤波器，我修改你程序上哪几个参数可以变成要求的低通，高通，带阻滤波器？
2.我是否简单地修改你设定地频率范围的话就能对别的范围的信号也能进行滤波，我还要做个对一个有10hz、80hz、与160hz的信号作分析。

楼主只要看一下fir1函数调用的参数，以上两个问题都能得到解决的。

不好意思，希望楼主作出最后帮助。前面的fir1的频率设置我已经理解了，现在有个最后的问题。当我把fs调整至300以上时，滤波信号的幅值好像就会衰减，我要怎样修改它？我做的设计FS参数一般是256hz、和512hz，望楼主最后帮助，麻烦了，谢谢！！

fir1的矢用方法为：b=fir1(n,Wn,'type');
Wn是截止频率的参数，不同的采样频率，该参数也是变化的。例如在上例中fs=300时，带通的参数是
[30 70]/150=[30 70]/(300/2)=[fc1 fc2]/(fs/2)，
当fs=256时该参数变为[30 70]/(256/2)。这样设计出的滤波器在通带内不会衰减。

[ 本帖最后由 songzy41 于 2008-5-3 20:10 编辑 ]