伪随机序列仿真实验报告

伪随机序列仿真仿真实验报告

报告时间：2013/5/13

姓名：刘敬亚

学号：2010101026

电子邮件：403555418@https://www.360docs.net/doc/9310551089.html,

一、实验目的

1、熟悉MATLAB 仿真的应用；

2、掌握伪随机序列的原理、软件产生，并仿真分析其相关特性。

二、实验内容

1、MATLAB中产生127位gold序列，并分析其自相关和互相关特性。

2、采用MATLAB的伪随机函数，生成127位二进制伪随机序列，并分析其自

相关和互相关特性。

3、收集学习技术资料，列举五种以上通信系统中采用的伪随机序列。

三、实验过程记录

根据搜集的资料在matlab中编写代码，然后运行，得到相应的特性图。

1、产生127位的gold序列

2、gold序列的自相关性

close all;

clear all;

clc;

n=7;%移位寄存器级数

P=2^n-1%m列的周期

%产生m1序列%

c1=[1 0 1 0 0 1 1];%反馈系数==211

reg1=[zeros(1,n-1) 1];%初始化移位寄存器

m1(1)=reg1(n);

for i=2:P

new_reg1(1)=mod(sum(c1.*reg1),2);

for j=2:n

new_reg1(j)=reg1(j-1);

end

reg1=new_reg1;

m1(i)=reg1(n);

end

subplot(311);

stairs(m1);

xlabel('k');

title('m1序列');

%产生m2序列%

c2=[1 0 1 1 0 0 1];%反馈系数==217

reg2=[zeros(1,n-1) 1];%初始化移位寄存器m2(1)=reg2(n);

for i=2:P

new_reg2(1)=mod(sum(c2.*reg2),2); for j=2:n

new_reg2(j)=reg2(j-1);

end

reg2=new_reg2;

m2(i)=reg2(n);

end

subplot(312);

stairs(m2);

xlabel('k');

title('m2序列');

%产生gold序列%

for i=0:P-1

shift_m2=[m2(i+1:P) m2(1:i)];%m2

g(i+1,:)=mod(m1+shift_m2,2);%模2加end

subplot(313)

stairs(g(1,:));

xlabel('k');

title('gold序列的第1行');

%gold序列的自相关特性和互相关特性%

g1=g(1,:);

g2=g(2,:);%取gold序列的第2行

% for i=-P+1:-1

% shift=[g1(i+P+1:P) g1(1:i+P)];

% r1(P+i)=sum(g1(P+i).*shift(P+i)); % end

% for i=0:P-1

% shift=[g1(i+1:P) g1(1:i)];

% r1(P+i)=sum(g1(i+1).*shift(i+1)); % end

r1=conv(g1,g1);

figure(2)

subplot(211)

stem(r1/P);

xlabel('k');

title('g1的自相关特性');

% for i=-P+1:-1

% shift=[g1(i+P+1:P) g1(1:i+P)]; % r2(P+i)=g2(P+i).*shift(P+i); % end

% for i=0:P-1

% shift=[g1(i+1:P) g1(1:i)];

% r2(P+i)=g2(i+1).*shift(i+1); % end

r2=conv(g1,g2);

subplot(212);

stem(r2/P);

xlabel('k');

title('g1与g2的互相关特性');

产生的图形：

四、分析及结论

从实验结果可以得出，Gold序列间不仅有良好的互相关性能，而且个数远比m序列多。但是Gold序列的周期自相关函数在非零点也是三值的，这一点远不如m序列。