M序列产生及其特性仿真实验报告.doc

m序列产生及其特性仿真实验报告一、三种扩频码序列简介1.1 m序列 它是由多级移位寄存器或其他延迟元件通过线性反馈产生的最长的码序列。 m序列的特性1、最长周期序列：n=2n-12、功率平衡性：1的个数比0的个数多13、0、1随机分布：近似高斯噪声4、相移不变性：任意循环移位仍是m序列，仅初相不同5、离散自相关函数：0-+1，1-11.2 gold序列 gold序列是两个等长m序列模二加的复合序列 两个m序列应是“优选对”特点： 1、包括两个优选对m序列，一个gold序列族中共有2n+1个gold序列2、 gold序列族中任一个序列的自相关旁瓣及任意两个序列的互相关峰值均不超过两个m序列优选对的互相关峰值 1.3ovsf序列 又叫正交可变扩频因子，系统根据扩频因子的大小给用户分配资源，数值越大，提供的带宽越小，是一个实现码分多址(cdma)信号传输的代码，它由walsh函数生成，ovsf码互相关为零，相互完全正交。ovsf序列的特点1、序列之间完全正交2、极适合用于同步码分多址系统3、序列长度可变，不影响正交性，是可变速率码分系统的首选多址扩频码4、自相关性很差，需与伪随机扰码组合使用2、 三种扩频码序列产生仿真一、m序列的产生代码：x1=1;x2=0;x3=1;x4=0; %移位寄存器输入xi初t态（0101）， yi为移位寄存器各级输出m=60; %置m序列总长度for i=1:m %1# y4=x4; y3=x3; y2=x2; y1=x1; x4=y3; x3=y2; x2=y1; x1=xor(y3,y4); %异或运算 if y4=0 u(i)=-1; else u(i)=y4;endendm=u%绘图i1=ik=1:1:i1;plot(k,u,k,u,rx)xlabel(k)ylabel(m序列)title(移位寄存器产生的m序列)用阶梯图产生表示：x1=1;x2=0;x3=1;x4=0; %移位寄存器输入xi初t态（0101）， yi为移位寄存器各级输出m=60; %置m序列总长度for i=1:m %1# y4=x4; y3=x3; y2=x2; y1=x1; x4=y3; x3=y2; x2=y1; x1=xor(y3,y4); %异或运算 if y4=0 u(i)=-1; else u(i)=y4;endendm=u%绘图stairs(m);二、gold序列的产生： m序列a的生成：x1=1;x2=0;x3=1;x4=0; %移位寄存器输入xi初t态（1010）， yi为移位寄存器各级输出m=60; %置m序列总长度for i=1:m %1# y4=x4; y3=x3; y2=x2; y1=x1; x4=y3; x3=y2; x2=y1; x1=xor(y3,y4); %异或运算 if y4=0 a(i)=0; else a(i)=y4;endendm=a%绘图i1=ik=1:1:i1;plot(k,a,k,a,rx)xlabel(k)ylabel(m序列)title(移位寄存器产生的m序列) m序列b的生成： x1=0;x2=1;x3=0;x4=1; %移位寄存器输入xi初t态（0101）， yi为移位寄 存器各级输出m=60; %置m序列总长度for i=1:m %1# y4=x4; y3=x3; y2=x2; y1=x1; x4=y3; x3=y2; x2=y1; x1=xor(y3,y4); %异或运算 if y4=0 b(i)=0; else b(i)=y4;endendn=b%绘图i1=ik=1:1:i1;plot(k,b,k,b,rx)xlabel(k)ylabel(m序列)title(移位寄存器产生的m序列)生成gold序列： c=xor(a,b); stairs(c); 三、ovsf序列的产生： %function ovsf_codes=ovsf_generator(spread_fator,code_number)%code_number=-1 表示生成所有扩频因子=spread_factor的ovsf码code_number=-1;spread_fator=8;ovsf_codes=1;if spread_fator=1 return;endfor i=1:1:log2(spread_fator) temp=ovsf_codes; for j=1:1:size(ovsf_codes,1) if j=1 ovsf_codes=temp(j,:),temp(j,:) temp(j,:),(-1)*temp(j,:); else ovsf_codes=ovsf_codes temp(j,:),temp(j,:) temp(j,:),(-1)*temp(j,:); end endend%if code_number-1 % ovsf_codes=ovsf_codes(code_number+1),:);%endfigure(3) b4,t4=stairs(1:length(ovsf_codes),ovsf_codes); plot(b4,t4); axis(0 130 -1.1 1.1);title(ovsf序列)三、三种扩频码序列特性仿真 （一）m序列自相关函数x1=1;x2=0;x3=1;x4=0; %移位寄存器输入xi初t态（0101）， yi为移位寄存器各级输出m=28-1; %置m序列总长度for i=1:m %1# y4=x4; y3=x3; y2=x2; y1=x1; x4=y3; x3=y2; x2=y1; x1=xor(y3,y4); %异或运算 if y4=0 u(i)=-1; else u(i)=y4;y = xcorr(u);stairs(y);end互相关函数：输入两个m序列clcclear allclose allm1 = 0 1 0 0 1 1 0 0 1 0 1 0 1 0 0 1 1 0 1 1 0 0 0 1m2 = 1 0 1 0 1 1 1 0 0 1 1 1 0 0 1 0 1 0 0 1 1 1 0 1y = xcorr(m1,m2,unbiased);stairs(y)（二）gold码的自相关函数x2=(2*c)-1;%将运行结果gold序列c从单极性序列变为双极性序列y1=xcorr(x2,unbiased);%求自相关性stairs(y1);gridxlabel(t)ylabel(相关性)title(移位寄存器产生的gold序列的相关性)互相关性gold序列和m序列的互相关性y1=xcorr(c,m1,unbiased);stairs(y1);（3） ovsf码的互相关和自相关 a=1 -1 1 1 -1 1 -1 -1;b=1 -1 -1 1 1 -1 -1 1;p=length(a);%求序列a的自相关函数ra(1)=sum(a.*a);for k=1:p-1 ra(k+1)=sum(a.*circshift(a,0,k);end%求序列b的自相关函数rb(1)=sum(b.*b);for k=1:p-1 rb(k+1)=sum(b.*circshift(b,0,k);end%求序列a和b的互相关函数rab(1)=sum(a.*b);for k=1:p-1 rab(k+1)=sum(a.*circshift(b,0,k);endx=0:p-1;figure(9)subplot(3,1,1);stem(x,rab);ylabel(a和b的互相关函数);axis(0 p-1 -10 12);grid;xlabel(偏移量);subplot(3,1,2);stem(x,ra);ylabel(a自相关函数);xlabel(偏移量);%axis(0 p-1 -5 30);subplot(3,1,3);ste