RDT随机减量程序的编制

wht0211

下面是EMD后得到的本征模态进行随机减量求自由相应的程序：但看了关于随机减量的具体方法，还是感觉课本与程序对应不起来，哪位高手给讲解一下，谢谢啦
sf=255;
np=1024;                               %%%输出数据长度
t=0:1/sf:(np-1)/sf;                     %建立离散输出时间向量
nt=length(x);                          %取输入数据长度
s=1.5*std(x);                          %设置截取振幅为输入信号标准差的1.5倍
m=0;                                   %获取输入信号的子样本函数进行叠加
y=zeros(1,np);
for k=2:nt-np
a=abs(x(k-1)-s);
b=abs(x(k)-s);
c=abs(x(k+1)-s);
if b<a&b<c
y(1:np)=y(1:np)+x(k:k+np-1);
m=m+1;
end
end
y=y./m;                                % 对叠加结果做平均
figure(3);
plot(t,y);                             % 绘制自由振动时程曲线图

[ 本帖最后由 wht0211 于 2010-8-17 15:43 编辑 ]

wht0211

RDT方法.doc (84.5 KB, 下载次数: 37)

2010-8-17 15:45 上传

点击文件名下载附件