gonglvpu.m

来自「对噪声及信号的功率谱估计,生物医学信号处理,matlab编程实现」· M 代码 · 共 17 行

clear;              %功率谱估计
Fs=1000;            %采样频率    
L=2;                %数据长度
n=0:1/Fs:L;         
N=length(n);        
x1=sin(2*pi*40*n);
e=randn(1,N);       %白噪声
x=x1+e;
z=zeros(1,N);
X=[x z];            %补零
Rx=xcorr(X,X);      %自相关 
M=abs(fft(Rx,2*N)); %做FFT并取模
plot((1:2*N)*Fs/(2*N),M);  %把横坐标转换成频率，并使数据一致