clf;clear;
Ts=0.005;?
fs=1/Ts;?
NLen=1024;?
n=0:NLen-1;
s=sin（2*pi*4.3*n*Ts）


x=sin（2*pi*4.3*n*Ts）+0.5*randn（11024）%sin（2*pi*42.5*n*Ts）+sin（2*pi*85*n*Ts）+sin（2*pi*20*n*Ts）+

subplot（311）;
plot（s）;
subplot（312）;
plot（x）;

s1=sgolayfilt（x341）

subplot（313）;
plot（s1）;
xfft=fft（xNLen）;??
%mag=abs（xfft）
xfft=xfft.*conj（xfft）/NLen;?%求幅值
y1=ff