lmfchirp2.m

来自「这是有关锁相环的一些必用知识.希望能对大家有所帮助哈」· M 代码 · 共 37 行

m=300;
for i=-m:m
    mix_cum(:,i+m+1)=cum4est(mix,m,1600,60,'unbiased',i,i);
end

for i=1:2*m+1
    bb(i)=mix_cum(i,i);
end
bb(m+1)=bb(m+1)/10;

bb_spect=abs(fft(bb,length(mix)));
figure
subplot(2,1,1)
plot(bb)
subplot(2,1,2)
plot(bb_spect)



%---------------------------------------

mix=[zeros(1,(n-length(mix))/2),mix,zeros(1,(n-length(mix))/2)];
x=[zeros(1,n/2),mix,zeros(1,n/2)];
r=0;
for tao=0:n/2
    r1=zeros(1,2*n);
    r1(n/2+1:n/2+n)=x(n/2+1+tao:n/2+n+tao).*x(n/2+1-tao:n/2+n-tao);
    r=r+r1;
end
r=r(n/2+1+250:n/2+n-250);
r_spect=abs(fft(r));

figure
subplot(2,1,1)
plot(r)
subplot(2,1,2)
plot(r_spect)