format?long
y=load（‘1.txt‘）;
n=length（y）;
var1=std（y）;
subplot（411）
plot（y）
grid
for?m=1:floor（n/1000）;
x（m）=mean（y（（m-1）*1000+1:m*1000））;
end
var2=std（x）;
%y=y（500:1700）;
%y=S1;
sampling_rate=1000;
y0=mean（y）
fy=fft（y-y0）;
fre=（0:length（fy）-1）‘*sampling_rate/length（fy）;

%Pyy?=?fy.*?conj（fy）?;
%Pyy=（Pyy/sum（Pyy））*2;??
%plot（frePyy）;
%Pyy?=?fy.*?conj（fy）;?
%plot（frePyy）;%功率；
%grid
fw=2*fy./n;
subplot（412）
plot（freabs（fw））
grid
i=0;
fw=abs（