声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 11267|回复: 39

[HHT] 关于二维EMD的matlab程序中FastRBF的问题

  [复制链接]
发表于 2009-11-1 22:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
在网上下了Nunes J.C的BEMD程序,同时也下载了FastRBF™的工具箱,地址是http://www.fastrbf.com/,但在程序执行过程中出现了下述问题,看到前面有好几位也出现了同样的现象,但是始终没有给出具体的答案,就想拿出来请各位高人再指点一下,谢谢了~
在MATLAB中执行BEMD之后出现了下述错误:
..................................................................................................................................................
FastRBF - data interpolation toolbox V 1.4.2 (build: Oct  4 2006)
Copyright �1997-2001, FarField Technology Ltd
http://www.farfieldtechnology.com/

type 'help fastrbf' for more information

-----------------------------License------------------------------
  ERROR: file is invalid
----------------------------------------------------------------------

??? License Invalid

Error in ==> toolbox\private\FastRBF_MEX at 27
    [out{:}] = FastRBF_MEXr14(varargin{:});

Error in ==> fastrbf_fit at 103
Sol = FastRBF_MEX( 'Fit', varargin{:} );

Error in ==> Env_Rbf at 64
    rbfmin=fastrbf_fit(datamin,0.00000001);

Error in ==> BEMD_RBF at 47
        [Ienv_I,Ienv_S,nb_extr]=Env_Rbf(I,ext_lpe);

Error in ==> BEMD at 39
[If,Ires,Imodes]=BEMD_RBF(I,nb_modes,SDmax,ext_lpe,precision);
.............................................................................................................................................................
希望各位大师多提宝贵意见,小弟静候佳音了!~

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2010-9-22 08:58 | 显示全部楼层
关于二维EMD的matlab程序中FastRBF的问题
(1)首先matlab的版本应低于matlab2006A,或者是matlab7.1;
(2)其次给support@farfieldtechnology.com 发邮件所要FastRBF License;

评分

1

查看全部评分

发表于 2010-9-27 19:42 | 显示全部楼层
各位高手 我的MATLAB是7.0版本 怎么办啊?没有 提示错误 就是一堆英文 翻译还翻译不出来 我都崩溃了
发表于 2010-9-27 22:02 | 显示全部楼层

matlab再安装一个要求的版本
发表于 2011-1-28 11:07 | 显示全部楼层
zhuxiaofei 发表于 2010-9-22 08:58
关于二维EMD的matlab程序中FastRBF的问题
(1)首先matlab的版本应低于matlab2006A,或者是matlab7.1;
( ...

他们会给么,我已经写信了~
发表于 2011-1-31 07:54 | 显示全部楼层
zhuxiaofei 发表于 2010-9-22 08:58
关于二维EMD的matlab程序中FastRBF的问题
(1)首先matlab的版本应低于matlab2006A,或者是matlab7.1;
( ...

给他们官方写信了,今天早上刚回的我,原来真会给LIC文件~
发表于 2011-4-16 17:28 | 显示全部楼层
问下,LIC文件是不是只能用一次??
发表于 2011-4-17 16:02 | 显示全部楼层
回复 7 # cindy1127 的帖子

一般他们发你邮箱里的LICENSE只能用一个月,处理的点数也有限。
可以下个月换个邮箱继续申请一个~
发表于 2011-4-18 23:24 | 显示全部楼层
回复 8 # chenlu1986 的帖子

了了,谢谢!
发表于 2011-4-18 23:41 | 显示全部楼层
回复 8 # chenlu1986 的帖子

请问下下面这个程序是不是实现BEMD分解的局部取极值点,如果是的话,请问用的是什么方法??麻烦解释下代码好吗,谢谢!!
%   BEMD_RBF
%   Interpolation by radial basis functions from extrema or crest lines
%
%   J.C. Nunes, Y. Bouaoune, E. Del閏helle, N. Oumar, and Ph. Bunel.
%   "Image analysis by bidimensional empirical mode decomposition".
%   Image and Vision Computing Journal (IVC), (to appear), 2003.

function [I,Ires,If,nbr_extrema]=BEMD_RBF(Iori,nb_modes,SDmaxori,ext_lpe,precision)

SDmax=SDmaxori;             %   SDmax de d閜art     Crit鑢e d'arret d'obtention d'un mode
SDmax_prec=SDmaxori*2;      %   SDmax pr閏閐ent
kzeros=1;                   %   on ne tient pas compte des surfaces inf閞ieures � ce scale
difference=0;

if ext_lpe==1
    disp('Interpolation par les extrema')
else
    disp('Interpolation par les lignes cr鑤es')
end
t=size(Iori);
If=[];
alpha=1;                    % coefficient multiplicatif du SDmax entre deux modes

nbextrema=0;

for j= 1:nb_modes           % nombre de modes
    I=Iori;
    chaine='mode non trouv�';
    for i= 1:300            % nombre d'it閞ations par modes
        SDmax=SDmaxori;     %   standard deviation
        if i==1
            difference=0;
        end
        if j==1
            if i==1
                SDmaxori=SDmaxori;
            end
            SDmax=SDmaxori;
        end
        if j>1
            if i==1
                SDmaxori=SDmaxori*alpha;
            end
            SDmax=SDmaxori;
        end
            
        [Ienv_I,Ienv_S,nb_extr]=Env_Rbf(I,ext_lpe);

        [If(:,:,j),difference,SDmax_prec]=moy_env(I,Ienv_S,Ienv_I,SDmax,SDmax_prec,precision);
        chaine='mode non trouv�';
%        disp(['SDmax_prec= ' num2str(SDmax_prec)])
        
        if difference==1
            chaine='mode trouv�';
            difference=0;
            nbr_extrema(j)=nb_extr;
            SDmax=1;
            disp('!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!   Obtention d''un Mode')
            break;
        end
        I=If(:,:,j);
    end
    difference=0;
    nbr_extrema(j)=nb_extr;
   
    Iori=Iori-If(:,:,j);

    Ires=Iori;
%    figure
%    subplot(121),imagesc(If(:,:,j)), colormap(gray),  title(['Mode n�' num2str(j) chaine]),
%    subplot(122),imagesc(Ires), colormap(gray), title('R閟idu'),truesize, drawnow
   
end

return
发表于 2011-4-19 09:22 | 显示全部楼层
回复 10 # cindy1127 的帖子

是的 插值方法在注释的第一行有写:Interpolation by radial basis functions from extrema or crest lines
这是二维EMD分解的程序,你找副图片做下试验就会深入理解点的。
发表于 2011-4-19 19:42 | 显示全部楼层
回复 11 # chenlu1986 的帖子

谢谢!!
发表于 2011-5-9 22:40 | 显示全部楼层
回复 8 # chenlu1986 的帖子

不好意思!!我忘记怎么重新申请了?我换了个邮箱,发了邮件,可是没有回复,是不是漏了什么步骤??麻烦告知一下,谢谢
发表于 2011-5-11 16:38 | 显示全部楼层
发表于 2011-5-14 19:10 | 显示全部楼层
回复 13 # cindy1127 的帖子

support@farfieldtechnology.com
一般马上就有回复,在附件里有一个月试用的LINCESE~
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-17 07:03 , Processed in 0.084499 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表