library(waveslim)
fs=128#采样频率
#模拟波形
t=seq(1/fs,2,1/fs)#样本数要是2^J
S=1.5*sin(2*pi*6*t)+2*sin(2*pi*50*t)+3*sin(2*pi*100*t)
d<-dwt(S,, wf="d8",n.levels=3)
d$d2<-rep(0,length(d$d2))
d$d1<-rep(0,length(d$d1))
d$d3<-rep(0,length(d$d3))
id<-idwt(d)
op<-par(mfrow=c(3,1))
plot(S~t,t='l',ylab='原始信号')
plot(id~t,t='l',ylab='除高频')
plot(t,1.5*sin(2*pi*6*t),t='l')
par(op)
离散小波变化