求助随机减量法程序

tumugong521

各位信号处理高手：在下急需一基于matlab的随机减量法程序，希望有这个程序的同仁可以提供。邮箱：whut_jh@163.com。谢谢！

yangkeer10

tumugong521兄有没有搞到RDT的程序啊？有的话能不能给我发一份啊？邮箱seudengyang@yahoo.com.cn

freeplus

%随机减量技术提取自由振动响应
% y 为随机响应信号（实测信号）
% A 为割点值
% L 为所需的自由响应信号长度
% z 为提取出来的自由响应信号
% freeplus, 2010,1 调试通过
function z=random_reduce(y,A,L)
n=length(y);
t1=[];
for k=1:(n-L)  %求割线点集
    if ( ((y(k)<=A)&&(y(k+1)>=A)) || ((y(k)>=A)&&(y(k)<=A)) )
        t1=[t1; k];
    end;
end;
m=length(t1);
if (floor(m/2)*2 == (m-1)), t1=t1(1:(m-1)); m=length(t1); end;  %保证偶数个割点

zz=zeros(L,m);
for k=1:m, zz(:,k)=y(t1(k):(t1(k)+L-1)); end;
z=sum(zz,2)/m;