rdt.m

来自「随机减量方法」· M 代码 · 共 20 行

function [t,y]=rdt(x)
sf=20;         %采样频率
np=601;               %输出数据长度
t=(0:1/sf:(np-1)/sf)';    %建立离散输出时间向量
nt=length(x);
s=1.5*std(x);
%为获取输入信号的子样本函数
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;