m序列产生及自相关和互相关函数曲线.docxVIP

1、 m 序列产生及自相关和互相关函数曲线 function PN=makem(x) %m 序列产生函数 ss1=num2str(x); ss2=dec2bin(base2dec(ss1,8)); %先把八进制转换为十进制， 再把十进制转换为 二进制 G=2^(length(ss2)-1)-1; %最大周期 sd=[]; for j=1:(length(ss2)-2) sd=[sd 0]; end sd=[sd 1]; %寄存器初始状态 0 0 0...0 1 PN=[]; for j=1:G PN=[PN sd(length(sd))]; %m 序列输出的第一位 onenum=[]; for jj=1:length(ss2) if str2num(ss2(jj))==1 onenum=[onenum jj-1]; %存储二进制反馈系数里面“ 1”的位置 -1，即进行异或的位置 end end temp=sd(onenum(2)); for jj=3:length(onenum) %根据“ 1”的位置进行异或运算 temp=xor(temp,sd(onenum(jj))); end for jj=length(ss2)-1:-1:2 %移位（序列后一位值等于前一位值） sd(jj)=sd(jj-1); end sd(1)=temp; %序列第一位等于反馈出来的值 end function mandzi(ss) %m 序列曲线及自相关函数曲线绘图函数 ss1=num2str(ss); ss2=dec2bin(base2dec(ss1,8)); %转换为二进制 G=2^(length(ss2)-1)-1; %最大周期 PN=makem(ss); %调用函数计算 m 序列 pp=(-2).*PN+1; %0→1 1→-1 pp2=[]; for tao=-(G-1):G-1 pp1=circshift(pp,[0,tao]); pp2=[pp2 sum(pp.*pp1)/G]; %计算自相关函数 end subplot(2,1,1) stem(PN) grid on; title([ 使用生成多项式 (,num2str(ss),)8=(,ss2,)2产生的 m 序列 ]) subplot(2,1,2) tao=-(G-1):G-1; plot(tao,pp2) grid on; title( 自相关函数曲线 ) function huxg(x,y) x1=num2str(x); x2=dec2bin(base2dec(x1,8)); G1=2^(length(x2)-1)-1; y1=num2str(y); y2=dec2bin(base2dec(y1,8)); G2=2^(length(y2)-1)-1; if G1~=G2  %m 序列互相关绘图函数 %转换为二进制 %最大周期 %转换为二进制 %最大周期 error(周期不同，无法计算  ) return end pn1=makem(x);  %分别调用函数计算出  m 序列 pn2=makem(y); pp=[]; for tao=-(G1-1):G1-1 pn1tao=circshift(pn1,[0,tao]); %pp=[pp sum(pn2.*pn1tao)/G1]; pp=[pp sum(pn2.*pn1tao)]; end  %计算互相关函数 tao=-(G1-1):G1-1; plot(tao,pp) grid on; title([ 反馈系数 (,num2str(x),)8 和(,num2str(y),)8 的互相关函数曲线 ]) 2、 Rake 接收机仿真 clear all; clc Numusers=1; Nc=16; %扩频因子 ISI_Length=1; %每径延时为 ISI_Lengh/2 EbN0db=[0:1:30]; %信噪比，单位 dB Tlen=8000; %数据长度 Bit_Error_Number1=0; %误比特率初始值 Bit_Error_Number2=0; Bit_Error_Number3=0; power_unitary_factor1=sqrt(6/9); %每径功率因子 power_unitary_factor2=sqrt(2/9); power_unitary_factor3=sqrt(1/9); s_initial=randsrc(1,Tlen); %数据源 wal2=[1 1;1 -1]; %产生 walsh 矩阵 wal4=[wal2 wal2;wal2 wal2*(-1)]; wal8=[wal4 wal4;wal4 wal4*(-1)]; wal16=[wal8 wal8;wal8 wal