gold序列的生成与相关特性仿真

Gold序列生成与相关性仿真

1.1 references

[1] 基于Matlab的Gold码序列的仿真与实现.

[2] Code Selection for CDMA Systems.

1.2 m序列的生成原理

1.2.1生成本原多项式

利用Matlab编程环境求解本原多项式，其运行结果如表1所示．选择n=7，采用7级移位寄存器，产生的序列周期是127，其程序如下所示．

N=7; %以7级寄存器为例，并组其中的一组优选对：211，,217

connections=gfprimfd(N,'all');

表（1）n=7 本原多项式

上面的多项式中，仅有9个是独立的．因为第一行和第十行，第二行和四行，第三行和第十六行，第五行和第八行，第六行和第十四行，第七行和第十三行，第九行和第十八行，第十一行和第十二行，第十五行和第十七行是两两对称的．用八进制数表示时，所选择的本原多项式为211、217、235、367、277、325、203、313和345共9条．在这9条本原多项式中，选择一个基准本原多项式，再按要求选择另一本原多项式与之配对，构成m序列优选对，对7级m序列优选对如下表：

表（2）n=7 m序列所以优选对

1.2.2构成移位寄存器

根据产生Gold码序列的方法，从上述本原多项式中选择一对m序列优选对，以211作为基准本原多项式，217作为配对本原多项式，通过并联结构形式来产生Gold序列，生成gold 序列的结构如图（6）所示：

图（6）Gold序列生成结构

1.3 自相关函数

仿真参数及初始值设定如下：

N=7; %以7级寄存器为例，并组其中的一组优选对：211，,217

connections=gfprimfd(N,'all');

f1=connections(4,:); %取一组本原多项式序列，211

f2=connections(16,:); %取另一组本原多项式序列,217

registers1=[1 0 0 0 0 0 0];%给定寄存器的初始状态

registers2=[1 0 0 0 0 0 0];%取相同的初始状态

生成的gold 序列自相关函数如图（7）、（8）所示

图（7） Gold 序列周期自相关函数

结论：自相关函数取值集合{127，15，-1，-17}

图（8）Gold 序列非周期自相关函数

020406080100120140

gold 序列周期自相关函数

020406080100120140

-40

-20

20

40

60

80

100

120

140

gold 序列非周期自相关函数

1.4 互相关函数

仿真时改变m序列寄存器初始状态，从而生成两个gold序列，求得互相关函数如图（9）（10）所示。

N=7; %以7级寄存器为例，并组其中的一组优选对：211，,217

connections=gfprimfd(N,'all');

f1=connections(4,:); %取一组本原多项式序列，211

f2=connections(16,:); %取另一组本原多项式序列,217

registers1=[1 0 0 0 0 0 0];%给定寄存器的初始状态

registers2=[1 0 0 0 0 0 0];%取相同的初始状态

registers3=[0 0 0 0 0 1 0];%给定寄存器的初始状态

registers4=[0 1 0 0 0 0 0];%取不相同的初始状态

gold序列周期互相关函数

图（9）gold序列周期互相关函数

gold序列非周期互相关函数

图（10）gold序列非周期互相关函数

结论：仿真得到gold序列具有三值互相关函数{-1，-17, 15}，从而验证了上述定理。

遍历第2个寄存器，得到127个gold序列，和第一个做互相关，记录最大值。

for m=1:L

registers4=dec2bin(m,7);%遍历第2个寄存器，得到127个gold序列，和第一个做互相关，记录最大值。

…

maxR=max(R);

minR=min(R);%记录最大最小值

end

仿真结果：对于每个循环，第二个寄存器初始值不同，求得自相关函数最值：

maxR=15

minR=-17

1.5 matlab仿真代码

clear all;

clc;

N=7; %以7级寄存器为例，并组其中的一组优选对：211，,217

connections=gfprimfd(N,'all');

f1=connections(4,:); %取一组本原多项式序列，211

f2=connections(16,:); %取另一组本原多项式序列,217

registers1=[1 0 0 0 0 0 0];%给定寄存器的初始状态

registers2=[1 0 0 0 0 0 0];%取相同的初始状态

L=2^N-1; %周期长度

sum2=0;

sum1=0;

for k=1:L

seq1(k)=registers1(N); %第一组m序列

seq2(k)=registers2(N); %第二组m序列

for j=1:N %进行模2加

sum1=sum1+f1(j+1)*registers1(j); %各级寄存器送参与模2加的值sum1=mod(sum1,2);

sum2=sum2+f2(j+1)*registers2(j); %各级寄存器送参与模2加的值sum2=mod(sum2,2);

end

for t=N:-1:2 %寄存器移位

registers1(t)=registers1(t-1);

registers2(t)=registers2(t-1);

end

registers1(1)=sum1;

registers2(1)=sum2;

sum2=0;

sum1=0;

end

z=seq1+seq2; %m序列的相加

gold=mod(z,2); %模2运算

%-------周期自相关函数------%%

a=1-2*gold;

b=a;

for i=1:L

R(i)=sum(a.*b)

b=[b(127) b(1:126)]

end

plot(R);

title('gold序列周期自相关函数');

% % %-------非周期自相关函数--------%%

% a=1-2*gold;

% b=a;

% for i=1:L

% R(i)=sum(a.*b)

% b=[0 b(1:126)]

% end

% plot(R)

% title('gold序列非周期自相关函数');

生成m序列与gold序列

一、生成m序列 function [mseq] = m_sequence(fbconnection); n = length(fbconnection); N = 2^n-1; %m序列的长度 register = [zeros(1,n - 1) 1]; %定义移位寄存器的初始状态 mseq(1)= register(n); %m序列的第一个输出码元 for i = 2:N newregister(1)= mod(sum(fbconnection.*register),2); %寄存器与反馈的模2和 for j = 2:n, newregister(j)= register(j-1); end; register = newregister; %移位后的寄存器 mseq(i) = register(n); %新的寄存器输出 end clear all; close all; clc; fbconnection=[0 0 1 0 1]; %输入本原多项式系数，从C1开始 m_sequence=m_sequence(fbconnection); stem(m_sequence); %对m序列绘图 axis([0 35 -0.2 1.2]); grid on;

二、生成gold序列 function goldseq = g_sequence(connection1,connection2); msequence1 = m_sequence(connection1); %生成第一个m序列 msequence2 = m_sequence(connection2); %生成第二个m序列 N=2^length(connection1)-1; %gold序列长度 for i = 1:N; s = mod(msequence1+msequence2,2); %两个m序列模二加产生gold序列 goldseq = s; end clear all; close all; clc; connection1=[0 0 0 0 1 1]; connection2=[1 0 0 1 1 1]; goldseq = g_sequence(connection1,connection2);

时间序列分析实验报告(3)

《时间序列分析》课程实验报告

一、上机练习（P124） 1.拟合线性趋势 12.79 14.02 12.92 18.27 21.22 18.81 25.73 26.27 26.75 28.73 31.71 33.95 程序： data xiti1; input x@@; t=_n_; cards; 12.79 14.02 12.92 18.27 21.22 18.81 25.73 26.27 26.75 28.73 31.71 33.95 ; proc gplot data=xiti1; plot x*t; symbol c=red v=star i=join; run; proc autoreg data=xiti1; model x=t; output predicted=xhat out=out; run; proc gplot data=out; plot x*t=1 xhat*t=2/overlay; symbol2c=green v=star i=join; run; 运行结果：

分析：上图为该序列的时序图，可以看出其具有明显的线性递增趋势，故使用线性模型进行拟合：x t=a+bt+I t,t=1,2,3,…,12 分析：上图为拟合模型的参数估计值，其中a=9.7086,b=1.9829，它们的检验P值均小于0.0001，即小于显著性水平0.05，拒绝原假设，故其参数均显著。从而所拟合模型为：x t=9.7086+1.9829t.

分析：上图中绿色的线段为线性趋势拟合线，可以看出其与原数据基本吻合。 2.拟合非线性趋势 1.85 7.48 14.29 23.02 37.42 74.27 140.72 265.81 528.23 1040.27 2064.25 4113.73 8212.21 16405.95 程序： data xiti2; input x@@; t=_n_; cards; 1.85 7.48 14.29 23.02 37.42 74.27 140.72 265.81 528.23 1040.27 2064.25 4113.73 8212.21 16405.95 ; proc gplot data=xiti2; plot x*t; symbol c=red v=star i=none; run; proc nlin method=gauss; model x=a*b**t; parameters a=0.1 b=1.1; der.a=b**t; der.b=a*t*b**(t-1); output predicted=xh out=out; run; proc gplot data=out; plot x*t=1 xh*t=2/overlay;

扩频编码M序列和gold序列

M序列 由n级移位寄存器所能产生的周期最长的序列。这种序列必须由非线性移位寄存器产生,并且周期为2n（n 为移位寄存器的级数）。例如，考察图中a的非线性反馈移位寄存器,其状态转移关系如表：

状态(a k-3，a k-2，a k-1)的接续状态是(a k-2,a k-1，a k)，其中a k＝a k-3嘰a k-1嘰1嘰a k-2a k-1是一种非线性逻辑。从任一状态出发，例如从(000)出发,其接续状态恰好构成一个完全循环(图b)，由此产生一个周期为23＝8的3级序列。M序列最早是用抽象的数学方法构造的。它出现于组合数学的一些数学游戏中，例如L.欧拉关于哥尼斯堡的七桥问题等。后来发现这种序列具有某些良好的伪随机特性。例如,M序列在一个周期中，0与1的个数各占一半。同时,同样长度的0游程与1游程也各占一半。所有这些性质在数据通信、自动控制、光学技术和密码学诸领域中均有重要应用。 隐蔽通信内容的通信方式。为了使非法的截收者不能理解通信内容的含义，信息在传输前必须先进行各种形式的变化，成为加密信息，在收信端进行相应的逆变化以恢复原信息。电报通信、电话通信、图像通信和数据通信,都有相应的保密技术问题。另一方面,为了从保密通信中获得军事、政治、经济、技术等机密信息，破译技术也在发展。保密技术和破译技术是在相互对立中发展起来的。 1881年世界上出现了第一个电话保密专利。电话保密开始是采用模拟保密或置乱的方法，即把话音的频谱或时间分段打乱。置乱后的信号仍保持连续变化的性质。在第二次世界大战期间，频域和时域的置乱器在技术上已基本成熟。70年代以来，由于采用集成电路，电话保密通信得到进一步完善。但置乱器仍是有线载波和短波单边带电话保密通信的主要手段。模拟保密还可以采用加噪声掩盖、人工混响或逆向混响等方法，但因恢复后话音的质量大幅度下降或保密效果差，这些方法没有得到推广应用。数字保密是由文字密码发展起来的。数字信号（包括由模拟信号转换成的数字信号），由相同速率的密码序列加密，成为数字保密信号；保密信号传输到收信端后由同一密码序列去密，恢复原数字信号。随着集成电路的发展，数字保密通信已成为保密通信的主要发展方向。话音、图像等模拟信号都可以用数字保密方式。一般来说，数字破译要比模拟破译困难得多。数字保密的主要限制是传输数字信号所需带宽要比传输模拟信号的带宽大好多倍。 模拟保密通信话音信号置乱后的带宽基本保持不变，这是模拟保密通信的一个特点。但是，置乱后恢复的话音质量有所下降。置乱的过程越复杂，则话音质量下降的程度越大。 倒频用倒频器（图1）把话音频谱颠倒过来，使高频变为低频,低频变为高频,这是最简单的一种频域置乱方法。频域置乱器的基本电路是平衡调制器和带通滤波器。平衡调制器可以搬移和倒置频谱，而滤波器可以滤取所需要的频谱成分。输入的话音信号经过平衡调制器后输出上、下两个边带。适当地选择

基于5级m序列的探究及仿真

通信系统综合设计与实践 题目基于5级m序列 的反馈系数的探究 院（系）名称信息工程学院 专业名称通信工程 学生姓名金宇、张艳丽、赵春阳 学生学号090110079、090110085、090110026 指导教师赵春雨 2012年05月21日

目录 1背景及原理 (1) 1.1探究n级m序列的反馈系数背景及目的 (1) 1.2 生成m序列的原理及方法 (1) 2 确定反馈系数的方法 (3) 2.1 判断本原多项式的方法 (3) 2.2 基于5级循环序列发生器特征方程组中满足本原多项式的反馈系数的分 析 (3) 2.3 基于5级循环序列发器反馈系数的程序 (5) 3 m序列的相关性质 (8) 3.1 m序列的性质 (8) 3.2 不同反馈系数对应m序列的性质的分析 (8) 3.2.1、平衡性 (8) 3.2.2、自相关性 (9) 3.2.3 互相关性分析 (12) 4 不同反馈系数对应的m序列对扩频通信系统抗干扰性能影响 (16) 4.1 m序列在直接序列扩频通信系统应用的simulink的仿真观察 (16) 4.2 不同/相同长度的不同m序列对扩频通信系统性能影响的matlab的仿 真 (19) 总结 (22) 参考文献 (23) 附录 (24) 附录A (24) 附录B (24) 附录C (25) 附录D (26) 指导教师评语： (31)

基于5级m序列的反馈系数的探究 摘要 m序列易于从干扰信号中被识别和分离出来,又可以方便地产生和重复，有随机噪声的优点，易于实现相关接受或匹配接受, 因此伪随机序列在相关辩识、伪码测距、扩频通信、多址通信、分离多径、误码测试、数据加扰、信号同步等方面均有广泛的应用。n级循环序列生成器生成m序列和自身的反馈系数密切相关，本文我们提供了n级循环序列发生器能产生m序列的反馈系数的判断方法，及分析了不同的反馈系数对扩频通信系统性能的影响，并在matlab环境下运行了模拟仿真。 首先，我们利用本原多项式生成算法，确定一个伪随机序列的特征方程中对应的m序列，进而确定相关m序列的反馈系数。又对m序列的性质进行的相关分析，我们得出了m序列平衡性为1，m序列越尖锐自相关性越明显，且仅在k=0时出现峰值。我们又进一步利用matlab仿真并分析不同m序列在直接扩频通信系统中的仿真，我们可知在实际应用中选择自相关性大、互相关小的m序列作为扩频序列，另外m序列越长可以提高系统的抗干扰能力、降低系统的误码率及增加系统的容量，由此选择最佳的反馈系数生成最优的m序列对于提高扩频通信系统性能指标具有非常重要的意义。 关键词: 伪随机序,m序列,反馈系数,仿真

m序列和Gold序列特性研究

扩频通信实验报告 - I- Harbin Institute of Technology 扩频通信实验报告 课程名称： 扩频通信 实验题目： Gold 码特性研究 院 系： 电信学院 班 级： 通信一班 姓 名： 学 号： 指导教师： 迟永钢 时 间： 2012年5月8日 哈尔滨工业大学

第1章实验要求 1.以r=5 1 45E为基础，抽取出其他的m序列，请详细说明抽取过程； 2.画出r=5的全部m序列移位寄存器结构，并明确哪些序列彼此是互反多项式； 3.在生成的m序列集中，寻找出m序列优选对，请确定优选对的数量，并画 出它们的自相关和互相关函数图形； 4.依据所选取的m序列优选对生成所有Gold序列族，确定产生Gold序列族的 数量，标出每个Gold序列族中的所有序列，并实例验证族内序列彼此的自相关和互相关特性； 5.在生成的每个Gold序列族内，明确标出平衡序列和非平衡序列，并验证其 分布关系。 6.完整的作业提交包括：纸质打印版和电子版两部分，要求两部分内容统一， 且在作业后面附上源程序，并加必要注释。 7.要求统一采用Matlab软件中的M文件实现。

第2章 实验原理 2.1 m 序列 二元m 序列是一种伪随机序列，有优良的自相关函数，是狭义伪随机序列。m 序列易于产生于复制，在扩频技术中得到了广泛应用。 2.1.1 m 序列的定义 r 级非退化的移位寄存器的组成如图1所示，移位时钟源的频率为c R 。r 级线性移位寄存器的反馈逻辑可用二元域GF(2)上的r 次多项式表示 2012() {0,1}r r i f x c c x c x c x c =++++∈ (1) 图 2-1 r 级线性移位寄存器 式(1)称为线性移位寄存器的特征多项式，其给出的表示反馈网络的而逻辑关系式是现行的。因此成为线性移位寄存器。否则称为，非线性移位寄存器。 对于动态线性移位寄存器，其反馈逻辑也可以用线性移位寄存器的递归关系式来表示 112233 {0,1}i i i i r i r i a c a c a c a c a c ----=++++∈ (2) 特征多项式(1)与递归多项式(2)是r 级线性移位寄存器反馈逻辑的两种不同种表示法，因其应用的场合不同而采用不同的表示方法。以式(1)为特征多项式的r 级线性反馈移位寄存器所产生的序列，其周期21r N ≤-。假设以GF(2)域上r 次多项式(1)为特征多项式的r 级线性移位寄存器所产生的非零序列{}i a 的周期为21r N =-，称序列为{}i a 是最大周期的r 级线性移位寄存器序列，简称m 序列。

m序列Simulink仿真实现

m序列Simulink仿真Verilog实现1. 4阶m序列生成器

Simulink模型如下： 其中，可以在Unit Delay属性中设置初始值为1000，由于Unit Delay输出为double，所以要将其转为Boolean以便进行模二加运算，使用XOR实现。

下面分别是最后一级和所有级的输出波形，可以看出，与上面的是一致的。

Verilog实现 `timescale 1ns / 1ps ////////////////////////////////////////////////////////////////////////////////// // Company: // Engineer: // // Create Date: 11:02:17 05/01/2012 // Design Name: // Module Name: PNcode // Project Name: ////////////////////////////////////////////////////////////////////////////////// module PNcode( clk, rst, PNstate, PNout

); input clk; input rst; output PNstate; output PNout; // PN code n = 4, f(x) = 1 + x + x^4 parameter order = 4; reg PNout = 0; reg [order-1 : 0] PNstate = 0; always @ (posedge clk) if(rst == 1) begin PNout <= 0; PNstate <= 4'b1000; // PN seed = b1000 end else begin PNout <= PNstate[0]; PNstate <= {PNstate[3]+PNstate[0], PNstate[3:1]}; end endmodule 测试文件： `timescale 1ns / 1ps //////////////////////////////////////////////////////////////////////////////// // Company: // Engineer: // // Create Date: 14:37:43 05/01/2012 // Design Name: PNcode // Module Name: E:/me/CAST/Project/FPGAcomm/PNcode/PNcode_tb.v

Gold序列与m序列仿真应用

1. 绪论 m 序列具有优良的双值自相关特性，但互相关特性不是很好。作为CDMA 通信地址码时，由于互相关特性不理想，使得系统内多址干扰影响增大，且可用地址码数量较少。在某些应用场合，利用狭义伪随机序列复合而成复合序列更为有利。这是因为通过适当方法构造的复合序列具有某些特殊性质。Gold 序列就是一种复合序列，而且具有良好的自相关与互相关特性，地址码数量远大于m 序列，且易于实现、结构简单，在工程上得到广泛应用。 表1是m 序列和Gold 序列的主要性能比较，表中max ?为m 序列的自相关峰值，(0)s ?为自相关主峰；()t n 为Gold 序列的互相关峰值，(0)g ?为其自相关主峰。从表1中可以看出：当级数n 一定时，Gold 序列中可用序列个数明显多于m 序列数，且Gold 序列的互相关峰值和主瓣与旁瓣之比都比m 序列小得多，这一特性在实现码分多址时非常有用。 表1. m 序列和Gold 序列性能比较 在引入Gold 序列概念之前先介绍一下m 序列优选对。m 序列优选对，是指在m 序列集中，其互相关函数绝对值的最大值(称为峰值互相关函数)max ()R τ最接近或达到互相关值下限(最小值)的一对m 序列。 设{a i }是对应于r 次本原多项式F 1(x )所产生的m 序列， {b i } 是另一r 次本原多项式F 2(x )产生的m 序列，峰值互相关函数满足 12 max 2 221()214r ab r r R τr ++?+?≤??+? 为奇数 为偶数但不是的整倍数 （1） 则m 序列{a i }与{b i }构成m 序列优选对。 例如：6r =的本原多项式61()1F x x x =++与6522()1F x x x x x =++++所产生的m 序列{}i a 与{}i b ，其峰值互相关函数2622 2 max ()172 12117r ab R τ++=≤+=+=。满足式（1） ，故{}i a 与{}i b 构成m 序列优选对。而本原多项式65323()1F x x x x x =++++所产生的m 序列 {}i c ，与m 序列{}i a 的峰值互相关函数max ()2317ac R τ=>，不满足上式，故{}i a 与{}i c 不 是m 序列优选对。 2. Gold 序列 1967年，R·Gold 指出：“给定移位寄存器级数r 时，总可找到一对互相关函数值是最小的码序列，采用移位相加方法构成新码组，其互相关旁瓣都很小，且自相关函数和互相关函数均有界”。这样生成的序列称为Gold 码（Gold 序列）。 Gold 序列是m 序列的复合序列，由两个码长相等、码时钟速率相同的m 序列优选对的模2

基于m序列的直接序列扩频

扩频通信实验 实验名称：基于m序列的直接序列扩频 专业班级：通信111501班 学生姓名：穆琦沈傲立孙琳王瑞学熊晓倩

学号：201115040111 13 16 20 27 指导教师：郑秀萍 时间：2014.10.29 1 需求分析 在通信发射端将载波信号展宽到较宽的频段上;在接收端,用同样的扩频码序列进行解扩和解调,把展宽的信号还原成原始信息.通过扩展频谱的相关处理,大大降低了频谱的平均能量密度,可在负信噪比条件下工作,获得了高处理增益,从而降低了被截获和检测的概率,避免了干扰影响.通过仿真模型结果分析抗噪声性能结果。 2 概要设计 扩频通信系统分为直接序列扩频系统、跳频扩频系统、跳时扩频系统和混合式扩频系统。直接序列扩频系统，又称“平均”系统或伪噪声系统，就是采用高码率的扩频码序列PN 码(伪随机码),在发送端与编码数据信号进行模2 加,产生一扩频序列,这一码序列由于码元很窄,占用了很宽的频带,达到扩频的目的,然后用扩频序列去调制载波并予以传输。在接收端接收到的扩频信号经高频放大混频之后,用与发端相同且同步的伪随机码对扩频信号进行相关解扩,由于收发端伪随机码的相关系数为1,故可以完全恢复所传的信息,而干扰和噪声由于与接收机伪

随机码不相关,在相关解调时大大降低进入信号通频带内的干扰。它是目前应用较广泛的一种扩展频谱系统。在国外已获得成功的空间探测器“喷气推进实验室（JPL）测距技术”就是一种直接序列调制，TATS-1 军用卫星中的扩展频谱多址(SSMA)系统等都使用DSSS。 直接序列扩频系统的接收一般采用相关接收，并分成两步，即解扩和解调。在接收端，接收信号经过数控振荡器放大混频后，用与发射端相同且同步的由M 序列发生器产生的伪随机码对中频信号进行相关解扩，把扩频信号恢复成窄带信号，然后再由基带滤波器进行解调，最后恢复出原始信息序列。扩频与解扩过程中,利用PN序列生成器模块( PN Sequence Generator ) ,产生6级、传输速率500b/s的PN伪随机序列来达到扩频和多址接入效果,这里扩频增益为50倍.扩频的运算是信息流与PN码相乘或模二加的过程.解扩的过程与扩频过程完全相同,即将接收的信号用PN码进行第二次扩频处理.要求使用的PN码与发送端扩频用PN码不仅码字相同,而且相位相同.否则会使有用信号自身相互抵消.解扩处理将信号压缩到信号频带内,由宽带信号恢复为窄带信号.同时将干扰信号扩展,降低干扰信号的谱密度,使之进入到信息频带内的功率下降,从而使系统获得处理增益,提高系统的抗干扰能力.调制与解调使用二相相移键控PSK方式. 为了方便分析, 我们可对系统作如下假设: 系统各用户同步;系统各用户功率相同;仅考虑系统MAI和白噪声干扰引起的误码, 忽略信号传输、调制解调过程中的误码。 3 开发工具和编程语言 开发工具：

基于MATLAB的m序列产生

第一章设计内容及要求 基于MATLAB产生m序列 要求： 1．通过matlab编程产生m序列的产生原理及其产生方法。 2．对特定长度的m序列，分析其性质，及其用来构造其它序列的方法。

第二章m序列设计方案的选择 2.1 方案一 MATLAB编程非常简单，无需进行变量声明，可以很方便的实现m序列。 2.2 方案二 图2.1 Simulink实现m序列 Simulink是MATLAB最重要的组件之一，它提供了一个动态系统建模，仿真和综合分析的集成环境。在此环境中无需大量书写程序，而只需通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应性广，结构及流程清晰及仿真精细等优点，基于以上优点，Simulink已被广泛的运用到控制理论和数字信号处理的复杂仿真和设计。 通过比较方案一和方案二，发现方案一的有点具有通用性而方案二利用MATLAB的Simulink直接搭建模块，在移位寄存器较少的情况下利用此方法比较简单，可是当移位寄存器的个数增多时，要搭建那么多的模块就显的很繁琐了，缺乏通用性，因此本次实验选择方案一。

第三章m序列的产生及性质 3.1 m序列的产生原理、结构及产生 m序列是最长线性反馈移位寄存器序列的简称，m序列是由带线性反馈的移位寄存器产生的。 由n级串联的移位寄存器和反馈逻辑线路可组成动态移位寄存器，如果反馈逻辑线路只由模2和构成，则称为线性反馈移位寄存器。 带线性反馈逻辑的移位寄存器设定初始状态后，在时钟触发下，每次移位后各级寄存器会发生变化，其中任何一级寄存器的输出，随着时钟节拍的推移都会产生一个序列，该序列称为移位寄存器序列。 n级线性移位寄存器的如图3.1所示： ◇A 图3.1 n级线性移位寄存器 图中C i表示反馈线的两种可能连接方式，C i=1表示连线接通，第n-i 级输出加入反馈中；C i=0表示连线断开，第n-i级输出未参加反馈。 因此，一般形式的线性反馈逻辑表达式为 ------表达式3.1将等式左边的a n移至右边，并将a n=C0a n(C0=1)带入上式，则上式可以 写成 -------表达式3.2 定义一个与上式相对应的多项式 --------表达式3.3 其中x的幂次表示元素的相应位置。该式为线性反馈移位寄存器的特征

spss时间序列模型

《统计软件实验报告》SPSS软件的上机实践应用 时间序列分析

数学与统计学学院 一、实验内容： 时间序列是指一个依时间顺序做成的观察资料的集合。时间序列分析过程中最常用的方法是：指数平滑、自回归、综合移动平均及季节分解。 本次实验研究就业理论中的就业人口总量问题。但人口经济的理论和实践表明，就业总量往往受到许多因素的制约，这些因素之间有着错综复杂的联系，因此，运用结构性的因果模型分析和预测就业总量往往是比较困难的。时间序列分析中的自回归求积分移动平均法（ARIMA）则是一个较好的选择。对于时间序列的短期预测来说，随机时序ARIMA是一种精度较高的模型。 我们已辽宁省历年（1969-2005）从业人员人数为数据基础建立一个就业总量的预测时间序列模型,通过spss建立模型并用此模型来预测就业总量的未来发展趋势。 二、实验目的： 1.准确理解时间序列分析的方法原理 2.学会实用SPSS建立时间序列变量 3.学会使用SPSS绘制时间序列图以反应时间序列的直观特征。

4.掌握时间序列模型的平稳化方法。 5.掌握时间序列模型的定阶方法。 6.学会使用SPSS建立时间序列模型与短期预测。 7.培养运用时间序列分析方法解决身边实际问题的能力。 三、实验分析： 总体分析： 先对数据进行必要的预处理和观察，直到它变成稳态后再用SPSS对数据进行分析。 数据的预处理阶段，将它分为三个步骤：首先，对有缺失值的数据进行修补，其次将数据资料定义为相应的时间序列，最后对时间序列数据的平稳性进行计算观察。 数据分析和建模阶段：根据时间序列的特征和分析的要求，选择恰当的模型进行数据建模和分析。 四、实验步骤： SPSS的数据准备包括数据文件的建立、时间定义和数据期间的指定。 SPSS的时间定义功能用来将数据编辑窗口中的一个或多个变量指定为时间序列变量，并给它们赋予相应的时间标志，具体操作步骤是： 1.选择菜单：Date→Define Dates，出现窗口：

Gold序列产生及其特性实验

湖南科技大学 移动通信实验报告 姓名：吴文建 学号：1208030104 专业班级：应用电子技术教育一班 实验名称：GOLD序列产生及其特性实验 实验目的：1）掌握Gold序列的特性、产生方法及应用。 2) 掌握Gold序列与m序列的区别。 实验仪器：1、pc机一台2、 实验原理： m序列虽然性能优良，但同样长度的m序列个数不多，且m序列之间的互相关函数并不理想（为多值函数）。 1.m序列优选对 m序列优选对是指在m序列集中，其互相关函数最大值的绝对值满足下式的两条n介m序列： 2.Gold序列的产生方法 Gold序列是m序列的组合序列，由同步时钟控制的两个码元不同的m序列优选对逐位模2加得到。这两个序列发生器的周期相同，速率相同，因而两者保持一定的相位关系，这样产生的组合序列与这两个自序列的周期也相同。当改变两个序列的相对位移，会得到一个新的Gold序列。Gold序列具有以下性质： （1）两个m序列优选对经不同移位相加产生的新序列都是Gold序列，两个n级移位寄存器可以产生2n+1个Gold序列，周期均为2n?1。 （2）Gold序列的周期性自相关函数是一个三值函数，与m序列相比，具有良好的互相关特性。 Gold序列的产生有两种形式：并联形式和串联形式 实验步骤： 1.预习Gold序列的产生原理及性质及独立设计Glod序列产生方法。 2.画出Gold序列仿真流程图。

3.编写MATLAB程序并上机调试。 4.比较m序列与Glod序列的异同。 5.撰写实验报告。 实验数据、结果表达及误差分析： 实验仿真图形如图所示 实验编写程序（此程序在实验五编写程序之上方可运行）：function c=gold() n=7; a=[1 1 1 1 1 1 1 1]; co=[]; for v=1:2^n-1 co=[co,a(1)]; a(8)=mod(a(5)+a(1),2); a(1)=a(2); 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 m1=co; b=[1 0 1 0 0 0 0 1];

m序列产生及其特性

一、实验目的 通过本实验掌握m 序列的特性、产生方法及应用。 二、实验内容 1、观察m 序列，识别其特征。 2、观察m 序列的自相关特性。 三、基本原理 m 序列是有n 级线性移位寄存器产生的周期为21n -的码序列，是最长线性移位寄存器序列的简称。码分多址系统主要采用两种长度的m 序列：一种是周期为1521-的m 序列，又称短PN 序列；另一种是周期为 4221-的m 序列，又称为长PN 码序列。m 序列主要有两个功能：①扩展调制信号的带宽到更大的传输带宽， 即所谓的扩展频谱；②区分通过多址接入方式使用同一传输频带的不同用户的信号。 3、m 序列的互相关函数 两个码序列的互相关函数是两个不同码序列一致程度（相似性）的度量，它也是位移量的函数。当使 用码序列来区分地址时，必须选择码序列互相关函数值很小的码，以避免用户之间互相干扰。 研究表明，两个长度周期相同，由不同反馈系数产生的m 序列，其互相关函数（或互相关系数）与自 相关函数相比，没有尖锐的二值特性，是多值的。作为地址码而言，希望选择的互相关函数越小越好，这 样便于区分不同用户，或者说，抗干扰能力强。 在二进制情况下，假设码序列周期为P 的两个m 序列，其互相关函数R xy (τ)为 ()xy R A D τ=- （9-9） 式中，A 为两序列对应位相同的个数，即两序列模2加后“0”的个数；D 为两序列对应位不同的个数， 即两序列模2加后“1”的个数。 为了理解上述指出的互相关函数问题，在此以5n =时由不同的反馈系数产生的两个m 序列为例计算它 们的互相关系数，以进一步讲述m 序列的互相关特性。将反馈系数为8(45)和8(75)时产生的两个5级m 序 列分别记做：1m ：1000010010110011111000110111010和2m ：111110111000101011010000110100，序列1m 和 2m 的互相关函数如表9-3所示。 表9-3序列1m 和2m 的互相关函数表

GOLD 序列码产生及特性分析实验

实验二 GOLD 序列码产生及特性分析实验 一、实验目的 1. 了解Gold 码的性质和特点； 2. 熟悉Gold 码的产生方法； 二、实验内容 1. 熟悉Gold 码的的产生方法； 2. 测试Gold 码的的波形； 三、实验原理 m 序列虽然性能优良，但同样长度的m 序列个数不多，且m 序列之间的互相关函数值并不理想（为多值函数）。1967年，R ．Gold 提出和讨论了一种新的序列，即Gold 码序列。这种序列有较为优良的自相关和互相关特性，构造简单，产生的序列数多，因而得到广泛的应用。 a) m 序列优选对 m 序列优选对是指在m 序列集中，其互相关函数最大值的绝对值满足下式的两条n 阶m 序列： 表2－1给出了部分m 序列优选对。 表2－1 部分优选对码表 级数 基准本原多项式 配对本原多项式 7 211 217,235,277,325,203,357,301,323 9 1021 1131,1333 10 2415 2011,3515,3177 11 4445 4005,5205,5337,5263 2．Gold 码的产生方法 Gold 码是m 序列的组合码，由同步时钟控制的两个码字不同的m 序列优选对逐位模2加得到，其原理如图2-1所示。这两个码发生器的周期相同，速率也相同，因而两者保持一整除为偶数，但不能被位奇数41212)(2/)2(2/)1(n n R n n xy ???++≤++τ

定的相位关系，这样产生的组合码与这两个子码序列的周期也相同。当改变两个m 序列的相对位移时，会得到一个新的Gold 码。Gold 码虽然是m 序列模2加得到的，但它已不再是m 序列，不过仍具有与m 序列近似的优良特性，各个码组之间的互相关特性与原来两个m 序列之间的互相关特性一样，最大的互相关值不会超过原来两个m 序列间最大互相关值。Gold 码最大的优点是具有比m 序列多得多的独立码组。 图2-1 Gold 码序列发生器 Gold 码序列具有以下性质： （1）两个m 序列优选对经不同移位相加产生的新序列都是Gold 序列，两个n 级移位寄存器可以产生2n +1个Gold 序列，周期均为2n -1。 （2）Gold 码序列的周期性自相关函数是一个三值函数，与m 序列相比，具有良好的互相关特性。 Gold 码的产生有两种形式：并联形式和串联形式。例如m 序列本原多项式为：61)(x x x f ++=和6521)(x x x x x f ++++=，构成的并联和串联形式的Gold 码发生器如2－2图所示。（a ）为并联形式，（b ）为串联形式。 （a ）并联结构 （b ）串联结构 图2－2Gold 码发生器 （a ） 并联形式（b ）串联形式 为了观测方便，本实验用两个周期为31的m 序列优选对采用并联结构产生一个Gold

m序列产生要点

设计内容及要求 基于MATLAB产生m序列 要求： 1．通过matlab编程产生m序列的产生原理及其产生方法。 2．对特定长度的m序列，分析其性质，及其用来构造其它序列的方法。 第二章m序列设计方案的选择 2.1 方案一 MATLAB编程非常简单，无需进行变量声明，可以很方便的实现m序列。 2.2 方案二 图2.1 Simulink实现m序列 Simulink是MATLAB最重要的组件之一，它提供了一个动态系统建模，仿真和综合分析的集成环境。在此环境中无需大量书写程序，而只需通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应性广，结构及流程清晰及仿真精细等优点，基于以上优点，Simulink已被广泛的运用到控制理论和数字信号处理的复杂仿真和设计。

通过比较方案一和方案二，发现方案一的有点具有通用性而方案二利用MATLAB的Simulink直接搭建模块，在移位寄存器较少的情况下利用此方法比较简单，可是当移位寄存器的个数增多时，要搭建那么多的模块就显的很繁琐了，缺乏通用性，因此本次实验选择方案一。 第三章m序列的产生及性质 3.1 m序列的产生原理、结构及产生 m序列是最长线性反馈移位寄存器序列的简称，m序列是由带线性反馈的移位寄存器产生的。 由n级串联的移位寄存器和反馈逻辑线路可组成动态移位寄存器，如果反馈逻辑线路只由模2和构成，则称为线性反馈移位寄存器。 带线性反馈逻辑的移位寄存器设定初始状态后，在时钟触发下，每次移位后各级寄存器会发生变化，其中任何一级寄存器的输出，随着时钟节拍的推移都会产生一个序列，该序列称为移位寄存器序列。 n级线性移位寄存器的如图3.1所示： ◇A 图3.1 n级线性移位寄存器 图中C i表示反馈线的两种可能连接方式，C i=1表示连线接通，第n-i 级输出加入反馈中；C i=0表示连线断开，第n-i级输出未参加反馈。 因此，一般形式的线性反馈逻辑表达式为 ------表达式3.1将等式左边的a n移至右边，并将a n=C0a n(C0=1)带入上式，则上式可以 写成

应用时间序列实验报告

河南工程学院课程设计 《时间序列分析课程设计》学生姓名学号： 学院：理学院 专业班级： 专业课程：时间序列分析课程设计指导教师： 2017年 6 月 2 日

目录 1. 实验一澳大利亚常住人口变动分析..... 错误!未定义书签。 实验目的............................................... 错误!未定义书签。 实验原理............................................... 错误!未定义书签。 实验内容............................................... 错误!未定义书签。 实验过程............................................... 错误!未定义书签。 2. 实验二我国铁路货运量分析........... 错误!未定义书签。 实验目的............................................... 错误!未定义书签。 实验原理............................................... 错误!未定义书签。 实验内容............................................... 错误!未定义书签。 实验过程............................................... 错误!未定义书签。 3. 实验三美国月度事故死亡数据分析...... 错误!未定义书签。 实验目的............................................... 错误!未定义书签。 实验原理............................................... 错误!未定义书签。 实验内容............................................... 错误!未定义书签。 实验过程............................................... 错误!未定义书签。课程设计体会 ............................ 错误!未定义书签。

gold序列的生成与相关特性仿真

Gold序列生成与相关性仿真 1.1 references [1] 基于Matlab的Gold码序列的仿真与实现. [2] Code Selection for CDMA Systems. 1.2 m序列的生成原理 1.2.1生成本原多项式 利用Matlab编程环境求解本原多项式，其运行结果如表1所示．选择n=7，采用7级移位寄存器，产生的序列周期是127，其程序如下所示． N=7; %以7级寄存器为例，并组其中的一组优选对：211，,217 connections=gfprimfd(N,'all'); 表（1）n=7 本原多项式 上面的多项式中，仅有9个是独立的．因为第一行和第十行，第二行和四行，第三行和第十六行，第五行和第八行，第六行和第十四行，第七行和第十三行，第九行和第十八行，第十一行和第十二行，第十五行和第十七行是两两对称的．用八进制数表示时，所选择的本原多项式为211、217、235、367、277、325、203、313和345共9条．在这9条本原多项式中，选择一个基准本原多项式，再按要求选择另一本原多项式与之配对，构成m序列优选对，对7级m序列优选对如下表：

表（2）n=7 m序列所以优选对 1.2.2构成移位寄存器 根据产生Gold码序列的方法，从上述本原多项式中选择一对m序列优选对，以211作为基准本原多项式，217作为配对本原多项式，通过并联结构形式来产生Gold序列，生成gold 序列的结构如图（6）所示： 图（6）Gold序列生成结构 1.3 自相关函数 仿真参数及初始值设定如下：

N=7; %以7级寄存器为例，并组其中的一组优选对：211，,217 connections=gfprimfd(N,'all'); f1=connections(4,:); %取一组本原多项式序列，211 f2=connections(16,:); %取另一组本原多项式序列,217 registers1=[1 0 0 0 0 0 0];%给定寄存器的初始状态 registers2=[1 0 0 0 0 0 0];%取相同的初始状态 生成的gold 序列自相关函数如图（7）、（8）所示 图（7） Gold 序列周期自相关函数 结论：自相关函数取值集合{127，15，-1，-17} 图（8）Gold 序列非周期自相关函数 020406080100120140 gold 序列周期自相关函数 020406080100120140 -40 -20 20 40 60 80 100 120 140 gold 序列非周期自相关函数

M序列的产生和性能分析

M序列的产生和性能分析 摘要 在扩频函数中，伪随机信号不但要求具有尖锐的互相关函数，互相关函数应接近于零，而且具有足够长的码周期，以确保抗侦破、抗干扰的要求；由足够多的独立地址数，以实现码分多址的要求。M序列是伪随机序列的一种，可由m序列添加全0状态而得到。m序列与M序列对比得出在同级移位寄存器下M序列的数量远远大于m序列数量，其可供选择序列数多，在作跳频和加密码具有极强的抗侦破能力。 本文在matlab中的Simulink下用移位寄存器建立了4级、5级、6级M 序列的仿真模型，进行了仿真，画出其时域图、频谱图、互相关性图。通过时域图和频域图可看出，经过扩频后的信号频带明显的被扩展；由M序列互相关性图，得出M序列有较小的互相关性，较强的自相关性，但相关性略差于m序列。最后，本文又将M序列应用于CDMA扩频通信仿真系统中，得到下列结论：当使用与扩频时相同的M序列做解扩操作与用其他序列做解扩的输出有巨大的差别。使用相同的序列进行解扩时系统输出值很大，而使用其他序列解扩时输出值在零附近变化。这就是扩频通信的基础。 关键词：伪随机编码, 扩频通信自相关函数，互相关函数

M SEQUENCE GENERATION AND PERFORMANCE ANALYSIS ABSTRACT In spread-spectrum communication, pseudo-random sequence must have high autocorrelation value, low cross correlation, long code period and lots of dependent address to satisfy code division multiple access(CDMA). M sequence is one kind of the pseudo-random sequences. It can be may obtained through adding entire 0 states to m sequence. The number of M sequence is greater than the m-sequence under the same level shift register. It may supply the more choice. The M-sequence is often applied to the frequency hopping and adds the password to have greatly strengthened anti- solves the ability. At first, M sequences which has n=4、5、7 levels of shift registers are produced under Simulink of Matlab. The time domain chart, the spectrograph, the mutual correlation chart are plotted. Through the time domain chart and the spectrograph, we could see how the bandwidth of the information signal is expanded. The pseudo-random symbol speed rate higher noise signal frequency spectrum is proliferated widely, the output power spectrum scope is lower. This can explain the spread-spectrum communication system principle from the frequency range. Through the M sequence’s auto correlation chart we can see that the auto correlation of M-sequence is quite good but is inferior to the m sequence. Finally, the M sequence is applied to the code division multiple access (CDMA) communication system. This is the spread-spectrum communication foundation. KEY WORDS：Pseudo-random code, auto-correlation, cross-correlation