跪求小波多脊提取matlab程序

xuxuesan

求一个小波多脊提取程序，matlab的，参考一下，xuesan.2007@163.com

xuxuesan

怎么没有人回复呢？

航仔

dsp2008

跪在哪里？

xuxuesan

呵呵，先拿出段程序来再说啊

xuxuesan

dsp2008 发表于 2012-11-4 18:41

                               
登录/注册后可看大图

跪在哪里？

四楼有没有啊？大神帮个忙吧

寂寞的部落

不知楼主解决问题了没，我这右端程序，也不是我自己写的，别人的，你可以参考下
%多脊提取
clear
SampFreq = 10;
x=1/SampFreq:1/SampFreq:20;%原理想信号
%sig=exp(-0.045*x).*cos(2*pi*1.5*x+pi)+exp(-0.098*x).*cos(2*pi*1.0*x+pi/2);
sig=exp(-0.045*x).*cos(2*pi*2*x+pi)+exp(-0.098*x).*cos(3*pi*1.0*x+pi/2);
% sig=sig';
% noise=randn(1,200);
% noise=noise';
% dB=20;
% sig=sigmerge(sig,noise,dB);
% sig=awgn(sig,20)
% sig=sig';

fmax = 0.5;                                    % 最高分析频率（归一化频率）
fmin = 0.005;                                  % 最低分析频率（归一化频率）
fb = 4 ;                                       % 取cmor4-2小波进行实验，带宽参数为4
fc = 2;                                        % 中心频率2Hz                                      
totalscal =512 ;                               % 所取尺度的数目                                         
FreqBins = linspace(fmin,fmax,totalscal);      % 将频率轴在分析范围内等间隔划分
Scales = fc./ FreqBins;                        % 计算相应的尺度参数     
RealFreqBins = FreqBins * SampFreq;            % 尺度所对应的实际频率
swd=cwt(sig,Scales,'cmor4-2');
swd=abs(swd);
figure(1)
imagesc(x,Scales,swd);
[a1,b1]=max(swd(2:120,:));
[a2,b2]=max(swd(100:250,:));
figure(2)
mesh(x,RealFreqBins,swd);
colormap jet;
shading interp;
colorbar;
axis([min(x) max(x) min(RealFreqBins) max(RealFreqBins)]);
ylabel('Frequency / Hz');
xlabel('Time / sec');
figure(3)
m=0.495/512*10*(b1+2);
plot(m)
axis([1 200 0.05 5])
hold on
n=0.495/512*10*(b2+100);
plot(n);
axis([1 200 0.05 5])

yghit08

寂寞的部落 发表于 2013-1-20 09:54

                               
登录/注册后可看大图

不知楼主解决问题了没，我这右端程序，也不是我自己写的，别人的，你可以参考下
%多脊提取
clear

比较简单的一段程序：分区域提取小波嵴（如果不知道每个嵴的大概范围呢？），优化算法这边是利用的Matlab自带的函数，这个不知道性能怎么样。再者，对于信号自身来说是不是还有负嵴（如果没做swd=abs(swd)的话）。当然，如果没做swd=abs(swd)会更困难些，因为优化函数不好选（有正有负），那么相应的优化算法就要考虑能提取局部区域内最值点的算法了。
大家接着分享自己知道的程序啊

寂寞的部落

你好，楼上，不知你是参数识别用什么方法？向你请教问题

yghit08

寂寞的部落 发表于 2013-1-20 10:31

                               
登录/注册后可看大图

你好，楼上，不知你是参数识别用什么方法？向你请教问题

目前，用过盲信号分离，子空间类方法，和ＥＭＤ，其他没做过。（ＡＲＭＡ模型类用过，不了解具体内容）。我提到的这三个内容，只子空间类方法我没有在本版发过（因为论文还没有被接收），不过你可以问我，如果你也做的话！

寂寞的部落

对你来说我应该是初学者，我用小波，目前遇到好多问题

yghit08

寂寞的部落 发表于 2013-1-20 11:01

                               
登录/注册后可看大图

对你来说我应该是初学者，我用小波，目前遇到好多问题

你可以把你的问题在这个帖子里列一列（先搜索论坛再做）

yesterday

学习下