faultnoise.m

此为基于二代提升小波的阈值压缩算法源码

function fault=faultnoise()
fs=100000;
sym t;
w=314;
t1=0:1/fs:0.1-1/fs;
k1=0*t1;
t2=0.1:1/fs:0.3;
k2=ones(size(t2));
t3=0.3+1/fs:1/fs:0.4;
k3=0*t3;
t=[t1 t2 t3];
k=[k1 k2 k3];
y=k.*(0.2*sin(3*w*t)+0.2*sin(7*w*t))+sin(w*t)+0.01*randn(size(t));
fault=y;
%plot(t,y);