直扩信号伪码序列估计要点讲解.docVIP

直扩信号的伪码估计 一、短码直扩信号的伪码估计 clear all; close all; clc; %---------初始化设置------- randn(state,sum(100*clock)); rand(state,sum(100*clock)); Ns=500; Nc=31; Sa=1; p=Nc*Sa; SNR=-15; delay=1; %------------------------------------- %--------产生PN序列----------- PN=2*(rand(1,Nc)0.5)-1; n=0; for i=1:Nc for j=1:Sa n=n+1; pn(n)=PN(i); end end ppn=zeros(1,p) %------------------------------------------------------- %---------产生随机信号，并对该信号通过PN码序列进行扩频------ signal=zeros(1,Nc*Sa*(Ns+2)); symbol=sign(rand(1,Ns+5)-0.50); ppn=[pn(1:Nc)]; %对ＰＮ码进行分用户，分为Ｇ组 s1=ppn.*symbol; %相应的用户上的PN码对信息位进行扩频 s2=s1(:).; z=size(s2); signal=s2(delay:delay+Nc*Sa*(Ns+2)-1); y=signal; %------------------------------------------------------- %-------对扩频后的信号加噪-------------------------------- L=length(y); snr=10.^(SNR/10); Eb=sum(y.^2)/L; N0=Eb/snr; sgma=sqrt(N0); noise=sgma*randn(1,L); yt=y+noise; %----------------------------------------------------- %-------------以一个伪码周期为窗口长度对500组数据求自相关------- for i=1:Ns x(:,i)=yt((i-1)*p+1:i*p).; end R=(x*x.)/Ns; %---------------------------------------------------------------------------------------------- %--------------对自相关函数进行奇异值分解，得到的U中的第一列即为所求伪码序列------------------------ [U,S,V]=svd(R); Us=U(:,1); figure(1) subplot(212) stem(Us); subplot(211) stem(pn); 图1 上图为伪码序列的真实值，下图为 信噪比为-15DB时，采样率为1/pin ，信息码长为500时的伪码序列的估计值，从图中可以看出估计值和真实值完全相反。 二、长码直扩信号的伪码估计 clear all; close all; clc; %------------初始化设置---------- randn(state,sum(100*clock)); rand(state,sum(100*clock)); Ns=500; Nc=31; Sa=5;G=5; p=Nc*Sa; m=p/G; SNR=-10; %------产生伪码序列-------------- delay=1; PN=2*(rand(1,Nc)0.5)-1; % % % % % % % % % % % % % % % % % % % % % % % b=[0 0 0 1 1]; %生成m序列 % % % % % % % % % % % % % % % % % % % % % % % PN=m_sequence(b); % % % % % % % % % % % % % % % % % % % % % % % for i=1:2^G-1 % % % % % % % % % % % % % % % % % % % % % % % if PN(i)==0 % % % % % % % % % % % % % % % % % % % % % % % PN(i)=1; % % %