clear close all
Fs=256; %%采样率 N=128; n=N/2:3/2N-1;%fc=100;w=2pi*fc; %载波是fc t=n/Fs; %%时长 f1=61;
x=exp(1i*(2pif1t+(2rand(1,1)-1)2pi)); %x信号 x=awgn(x,5); h=Blackman(N); h1=Blackman(2*N); % for m=1:N % x(m)=x(m)h(m); % end nn=0:2N-1; tt=nn/Fs;
f2=20; x_delay=exp(1i*(2pif2tt+(2rand(1,1)-1)2pi)); %x延迟信号 x_delay=awgn(x_delay,0); % for m=1:2*N % x_delay(m)=x_delay(m)*h1(m); % end
xor_len=length(x_delay)-length(x); xor_value=zeros(xor_len,N);
x_cor=0; for ii=1:length(x)
x_cor=x_cor+x(ii)*x(ii);
end x_cor=x_cor/length(x);
x_cor_conj=0; for jj=1:length(conj(x))
x_cor_conj