MATLAB仿真m序列和Gold序列自相关与互相关2

来自恶搞哥哥的MATLAB仿真

MATLAB仿真m序列和Gold序列自相关与互相关

MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可

视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括

MATLAB和Simulink两大部分。MATLAB可以进行矩阵运算、绘制函数和数据、

实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、

控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

本章正是运用MATLAB来仿真m序列和Gold序列的相关特性，以及OCDMA

系统的误码率同用户数的关系曲线。

1仿真过程

在理论分析的基础上，下面使用附录上两段程序，通过MATLAB仿真得出m

序列和Gold序列的自匕相关性。

这段m序列产生程序采用了8个移位寄存器，将最后两个移位寄存器的值进

行异或处理反馈给第一个移位寄存器，然后向前移位，输出最后一个移位寄存

器的值，Gold序列的产生只是将两个m序列中的一个进行延时移位，再进行异或,

产生的主要原理和m序列并无较大差异。

在这两段程序个前半部分m序列和Gold序列生成的基础上，只要将

yl=xcorr(xl)改为yl=xcorr(xl,x2)即可求出它们的互相关仿真。

2仿真结果

在系统中采用上述序列仿真得到自相关和互相关特性曲线如图1,2及图3

所示。

—

200

150

阳100

50

图1m序列的自相关曲线

移位寄存器产生的Gold序列的相关性

250

200

豺150

*

5n

100

50

0

102030405060708090100110120

图4-2m序列与Gold序列的自相关曲线

m序列和Gold序列的互相关性

o

乃

o

210

—5o

阳

o

10

50

0

102030405060708090100110120

t

S4-3m序列与Gold序列的互相关曲线

附录

程序1

Xl=l;X2=0;X3=l;X4=0;%移位寄存器输入乂1初丁态(0101),Yi为移位寄存器

各级输出

m=120;%置乂序列总长度

fori=l:m%1#

Y8=X8;Y7=X7;Y6=X6;Y5=X5;Y4=X4;Y3=X3;Y2=X2;Y1=X1;

X8=Y7;X7=Y6;X6=Y5;X5=Y4;X4=Y3;X3=Y2;X2=Y1;

Xl=xor(Y7,Y8);为异或运算

ifY8==0

U(i)=-1:

else

U(i)=Y8;

end

end

M—U

%绘图

il=i

xl=将运行结果m序列M从单极性序列变为双极性序列

yl=xcorr(xl)；%求自相关性

t=l:1:i1；

plot(t,yl(l:il))口*15([1,120,-12,288])%绘出信号的相关图

grid

xlabelCt')

ylabelC相关性')

title('移位寄存器产生的M序列的相关性')

程序2

functionc=gold()

n=7;

a=[11111111];

co=[];

forv=l:2^n-l

co=[co,a(l)];

a(8)=mod(a(5)+a(l),2);

a⑴二@出；

a(2)=a(3);

a(3)=a(4);

a(4)=a(5);

a(5)=a(6);

a(6)=a(7);

a(7)=a(8);

end

ml二co;

b=[l0100001];

co=[];

forv=l:2^n-l

co=[co,b(l)];

m=mod(b(5)+b(l),2);

p=mod(b(6)+m,2);

b(8)=mod(b(5)+b(l),2);

b⑴二b⑵；

b(2)=b(3);

b(3)=b(4);

b(4)=b(5);

b(5)=b(6);

b(6)=b(7);

b(7)=b(8);

end

m2=co;

c=xor(ml,m2);

x2=[(2*c)T]';%将运行结果Gold序列c从单极性序列变为双极性序列

yl=xcorr(x2)；