shawnzahk的博客

(考试必备模板1) fs=1000; t=0:1/fs:1; x=2sin(2pi20t)+2sin(2pi100t); y=x+1*randn(size(t)); plot(y(1:50)); Y=fft(y,256); Pyy=Y.conj(Y)/256; f=fs/256(0:127); figure(2); plot(f,Pyy(1:128)); xlabel(‘Freguency(HZ)’) [b,a]=ellip(4,0.1,40,[15 25]*2/fs); sf=filter(b,a,y);