1.相关函数和增强的谱线.m

来自「自适应谱线增强问题」· M 代码 · 共 28 行

clear;
clc;

Fi=2*pi*rand(1,1);                  % 产生随机相位Fi
X1=[];
for n=1:20000                       % 产生信号序列X1(n)
    x=sin(0.05*pi*n+Fi);
    X1=[X1;x];
end
XX1=xcorr(X1,'unbiased');

V=randn(30002,1);                   % 产生白噪声序列V(n)
X2=[];
for n=1:30000                       % 产生信号序列X2(n)
    x=V(n+2)+2*V(n+1)+V(n);
    X2=[X2;x];
end
X2=X2(10001:30000);
XX2=xcorr(X2,'unbiased');

figure(1)
plot(XX1(19901:20100))
figure(2)
plot(XX2(19901:20100))
figure(3)
hold on;
plot(XX1(19901:20100),'r')
plot(XX2(19901:20100))