RDT随机减量程序的编制
下面是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 编辑 ]
RDT方法叙述
页:
[1]