[资料共享]costas序列的模糊函数

% costas序列的模糊函数 %dd B$(  
clear all; 6m(+X MS  
close all; %,8 "cM`D  
seq1=[1 8 3 7 6 0 2 5];%costas序列 #K-O<:s=y  
T=1;%脉冲宽度 DM)Re~*  
Tr=5*T;%脉冲重复周期 )ARV>(  
fs=1/Tr;%采样频率 4e[ 0.2?  
i=0; oO &%&;[/A  
for tau=-1*T:.01:1*T %时延 1D"EF  
    i=i+1; Y,4?>:39J  
    j=0; $Vs5d=
B  
    for fd=-5/T:.05:5/T %多普勒频移 S}/ZHo  
        j = j + 1; !_>/ r  
        ambig(j,i)=a2(seq1,T,Tr,fs,tau,fd); ?$\sMkn  
    end ZI :wJU:f  
end v~dUH0P<>e  
tau=-1*T:.01:1*T; S[/D._5QD%  
fd=-5/T:.05:5/T; W-1Ub
|8C  
figure SrWmV@"y  
mesh(tau,fd,abs(ambig)/max(max(abs(ambig)))); }p9#Bzc  
xlabel('tau'); }vF=XA  
ylabel('fd'); "Q.C1#W}.  
子函数： ?ZKIs9E[m  
%跳频编码信号的模糊函数 {oo(HD;5  
function ambig=amb_hop(seq,T,Tr,fs,tau,fd) ZQ\O| n8  
%seq-跳频序列,T-脉冲宽度，Tr-脉冲重复周期 ll"6KI'X  
%fs-采样率,tau-时延，fd-频移 ,0,& L  
L=length(seq);%序列长度 >I/~)B`jhE  
eps=0.0000001; J rYL8 1  
l=sqrt(-1);%虚数l #zn`)n  
s=0; XS0NjZW  
    for n=1:1:L y>J6)F =  
        for m=1:1:L L6=5]
?B=  
            val1=exp(-l*2*pi*seq(m)*fs*tau); d\ 7OtM  
            val2=exp(l*2*pi*(fd-(seq(m)-seq(n))*fs)*n*Tr); ` gor  
            val3=exp(l*pi*(fd-(seq(m)-seq(n))*fs)*(T-tau+(m-n)*Tr)); w*Ze5j4@ \  
            val4=pi*(fd-(seq(m)-seq(n))*fs)*(T-abs(tau-(m-n)*Tr)); cn_KHz=  
            val5=(T-abs(tau-(m-n)*Tr))/T; RBeQT=B8~  
            amb(n,m)=val1*val2*val3*val5*sin(val4)/(val4+eps); {4R;C~E8  
 &n .. #@M'*X_%}K  
Jj1lAg0  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册