声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1235|回复: 1

[HHT] extr函数

[复制链接]
发表于 2010-5-7 09:58 | 显示全部楼层 |阅读模式

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

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

x
请问各位高手,在这个函数中,带红色的部分
是什么意思呢?通过x1=x(1:m-1);
x2=x(2:m);
indzer = find(x1.*x2<0);
这几个语句,不是可以求出来极值点的位置吗?

function [indmin, indmax, indzer] = extr(x,t);
% [indmin, indmax, indzer] = EXTR(x,t) finds extrema and zero-crossings
%
% inputs : - x : analyzed signal
%          - t (optional) : sampling times, default 1:length(x)
%
% outputs : - indmin = indices of minima
%           - indmax = indices of maxima
%           - indzer = indices of zero-crossings
if(nargin==1)
  t=1:length(x);
end
m = length(x);
x1=x(1:m-1);
x2=x(2:m);
indzer = find(x1.*x2<0);
if any(x == 0)
  iz = find( x==0 );
  indz = [];
  if any(diff(iz)==1)
    zer = x == 0;
    dz = diff([0 zer 0]);
    debz = find(dz == 1);
    finz = find(dz == -1)-1;
    indz = round((debz+finz)/2);
  else
    indz = iz;
  end
  indzer = sort([indzer indz]);
end
  
d = diff(x);
n = length(d);
d1 = d(1:n-1);
d2 = d(2:n);
indmin = find(d1.*d2<0 & d1<0)+1;
indmax = find(d1.*d2<0 & d1>0)+1;
if any(d==0)
  
  imax = [];
  imin = [];
   
  bad = (d==0);
  dd = diff([0 bad 0]);
  debs = find(dd == 1);
  fins = find(dd == -1);
  if debs(1) == 1
    if length(debs) > 1
      debs = debs(2:end);
      fins = fins(2:end);
    else
      debs = [];
      fins = [];
    end
  end
  if length(debs) > 0
    if fins(end) == m
      if length(debs) > 1
        debs = debs(1:(end-1));
        fins = fins(1:(end-1));
      else
        debs = [];
        fins = [];
      end      
    end
  end
  lc = length(debs);
  if lc > 0
    for k = 1:lc
      if d(debs(k)-1) > 0
        if d(fins(k)) < 0
          imax = [imax round((fins(k)+debs(k))/2)];
        end
      else
        if d(fins(k)) > 0
          imin = [imin round((fins(k)+debs(k))/2)];
        end
      end
    end
  end
  
  if length(imax) > 0
    indmax = sort([indmax imax]);
  end
  if length(imin) > 0
    indmin = sort([indmin imin]);
  end
  
end
回复
分享到:

使用道具 举报

发表于 2010-5-8 15:45 | 显示全部楼层
if any(diff(iz)==1) % find nonconsecutive zeros
    zer = x == 0;
    dz = diff([0 zer 0]);
    debz = find(dz == 1); % find peak-side of a zero
    finz = find(dz == -1)-1; % find trough-side of a zero
    indz = round((debz+finz)/2); % avg peaks and troughs to get zero crossings
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 18:31 , Processed in 0.057695 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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