求助以下程序

sushujing

那位知道去希尔伯特包络趋势的算法在那可以找到？下面程序中pd为零，怎么循环？fs_kHz=8
% Function to remove the trend in the hilbert envelope
%--------------------------------------------------------------------------

function [Ret] = rem_trend(h,fs_kHz);

pd = floor(fs_kHz*2/1000);

for i = 1:pd
    temp = h(1:i+pd);
    Ret(i) = h(i)/mean(temp);
    Ret(i) = Ret(i)*h(i);
end;

for i = pd+1:length(h)-pd
    temp1 = h(i-pd:i+pd);
    Ret(i) = h(i)/mean(temp1);
    Ret(i) = Ret(i)*h(i);
end;

for i = (length(h)-pd+1):length(h)
    temp2 = h(i:length(h));
    Ret(i) = h(i)/mean(temp2);
    Ret(i) = Ret(i)*h(i);
end;

zhlong

下面程序中pd为零，怎么循环？fs_kHz=8

是不是你给的fs_kHz=8太小了呢？

你程序中pd的作用似乎就是把h分为几块，然后分段求均值，再处理

[ 本帖最后由 zhlong 于 2007-5-9 11:41 编辑 ]