gold序列实验资料报告材料
移动通信 实验 解扩实验
实验十二解扩实验一．实验目的：1、通过本实验掌握载波已调信号m序列解扩原理及方法，掌握解扩前后信号在时域及频域上的变化。
2、通过本实验掌握载波已调信号GOLD序列解扩原理及方法，掌握解扩前后信号在时域及频域上的变化。
二．实验内容：1、观察解扩时本地扩频码与扩频时扩频码的同步情况。
2、观察已调信号在解扩前后的频域变化。
三．基本原理：m序列解扩的是在接收到的RF信号上进行的，其实解扩的原理很简单，即用一个与发送端完全相同的m序列与接收到的信号直接相乘就可以完成信号的解扩，两个m序列的相位必须一致，即接收端产生的m序列必须进行捕获和跟踪，以使其速率和相位与发送端m序列保持一致。
四．实验原理：1、实验模块简介（1）CDMA发送模块：本模块的主要功能：产生PN31伪随机序列，将伪随机序列或外部输入的其它数字序列扩频，扩频增益为32，扩频后输出码速率为512kbps，可输出两条不同扩频码信号。
（2）CDMA接收模块：本模块的主要功能：完成10.7MHz射频信号的选频放大，当本地扩频码设置为与发送端扩频码相同时，可完成扩频码的捕获及跟踪，进而完成射频信号的解扩。
（3）IQ调制解调模块：本模块的主要功能：产生调制及解调用的正交载波；完成射频正交调制及小功率线性放大；完成射频信号正交解调。
2、扩频后的PSK已调信号分为三路送入CDMA接收模块中，分别与结婚搜模块中产生的m序列的超前、同相、滞后序列相乘。
在扩频码没有捕获到时，同相支路的捕获输出为低电平，扣码电路工作，每周期扣掉1/4个码元，使发送端和接收端的两个PN序列产生相对滑动，当滑动到两个序列的相位差小于一个码元时，电平，扣码电路停止工作，系统进入跟踪状态。
此时超前-滞后支路产生的复合相关特性出现，经低通滤波后控制VCO，使收发端PN序列完全同步，此后跟踪过程一直存在，维持PN序列的同步。
PN码同相支路的相乘信号经带通滤波后即为解扩后的信号。
该信号时一个基带信元的PSK调制信号，扩频码调制部分已经被去除。
南昌大学完美的GOLD序列实验程序参考
%产生GOLD序列clc;clear all;close all;r=6;N=2^r-1;s1(1:6)=[1 0 0 0 0 1];%初始化S1s2(1:6)=[1 0 0 0 0 0];%初始化S2f1=[1 0 0 0 0 1 1];f2=[1 1 0 0 1 1 1];for n=r+1:N;s1(n)=mod(sum(s1(n-r:n-1).*f1(1:r)),2); end;for n=r+1:N;s2(n)=mod(sum(s2(n-r:n-1).*f2(1:r)),2); end;for n=r+1:N;s=mod(s1+s2,2);gold_sequence_1=s;figure(1);stem(gold_sequence_1);end;图为产生的GOLD序列%自相关分析clc;clear all;close all;r=6;N=2^r-1;s1(1:6)=[1 0 0 0 0 0];%初始化S1s2(1:6)=[1 0 0 0 0 0];%初始化S2f1=[1 0 0 0 0 1 1];f2=[1 1 0 0 1 1 1];for n=r+1:N;s1(n)=mod(sum(s1(n-r:n-1).*s1(1:r)),2);end;for n=r+1:N;s2(n)=mod(sum(s2(n-r:n-1).*s2(1:r)),2);end;s=mod(s1+s2,2);gold_sequence_1=s;figure(1);stem(gold_sequence_1);grid on;gold_sequence_1=2*gold_sequence_1-1;for j=0:N-1s3(j+1)=sum(gold_sequence_1.*[gold_sequence_1(1+j:N),gold_sequence_1( 1:j)])/N;end;j=-N+1:N-1;rho=[fliplr(s3(2:N)),s3];figure(2)plot(j,rho);axis([-30 30 -1 1.2])title('第一个gold序列的自相关函数')图为产生的GOLD序列图为GOLD序列对应的自相关函数图%互相关分析clc;clear all;close all;r=6;N=2^r-1;s1(1:6)=[1 0 0 0 0 0];%初始化S1s2(1:6)=[1 0 0 0 0 0];%初始化S2f1=[1 0 0 0 0 1 1];f2=[1 1 0 0 1 1 1];for n=r+1:N;s1(n)=mod(sum(s1(n-r:n-1).*f1(1:r)),2);end;for n=r+1:N;s2(n)=mod(sum(s2(n-r:n-1).*f2(1:r)),2);end;s=mod(s1+s2,2);gold_sequence_1=s;figure(1);stem(gold_sequence_1);grid on;r=6;M=2^r-1;s1(1:6)=[1 0 0 0 1 0];%初始化S1s2(1:6)=[1 0 0 0 1 0];%初始化S2f1=[1 0 0 0 0 1 1];f2=[1 1 0 0 1 1 1];for n=r+1:M;s1(n)=mod(sum(s1(n-r:n-1).*f1(1:r)),2);end;for n=r+1:N;s2(n)=mod(sum(s2(n-r:n-1).*f2(1:r)),2);end;s=mod(s1+s2,2);gold_sequence_2=s;figure(2);stem(gold_sequence_2);grid on;gold_sequence_1=2*gold_sequence_1-1;gold_sequence_2=2*gold_sequence_2-1;for j=0:M-1s3(j+1)=sum(gold_sequence_1.*[gold_sequence_2(1+j:N),gold_sequence_2( 1:j)])/N;end;j=-M+1:M-1;rho=[fliplr(s3(2:N)),s3];figure(3)plot(j,rho);axis([-30 30 -1 1.2])title('第一个gold序列的互相关函数')图为GOLD序列1图为GOLD序列2图为GOLD序列的互相关函数。
基于Gold序列的直接扩频通信系统仿真
频 技术 中应 用最 为广 泛 的就是 直接 扩频 技术 （ｉ ｃ Ｓｑｅｃ ｐｅｄＳ ｅｔ ｍ， Ｄｒ ｔ ｅｕｎｅ ｒａ ｐｃ ｕ 简称 Ｄ Ｓ ） ｅ Ｓ ｒ ＳＳ。
１ Ｇｏｄ序 列 ｌ
ｍ序列是一种应用较为广泛的伪随机序列 ， 是 带 由线 性 反馈 的移 位 寄 存 器 产 生 的 周期 最 长 的一种二进制序列 ， 然而其优选特性较差 ， 能 不
第 ２ 卷 第 ２期 ６
２ １ 年 ４月 ０２
商洛 学院学报
Ｊ ｕｎ ｌｏ ｈ ｎ ｈｏＵｎｖ ｒｉ ｏ ｒ ａ ｆ ａ ｇ ｉｅｓ ｙ Ｓ ｔ
Ｖ０．６ Ｎｏ２ １ ． ２
Ａｐ ． ２ ２ ｒ ０１
基 于 Ｇｏｄ序 列的直接扩频通信 系统仿真 ｌ
李 亚 文
（ 商洛学院物理与电子信息工程系， 陕西商洛 ７６０） ２００
定的， 除了军事 中应用外， 民用 中扩频技术也 在 正在显示 出它的强大生命力。 扩频技术正广泛用
７６ ０） ２０ ０
Ａｂ ｔ ａ ｔ Ｔ ｅ ｗｏ ｋ ｎ ｒｎ ｉ ｌ ｆ ｄ ｒ ｃ ｐ ｅ ｄ ｓ ｅ ｔ ｍ ｏ ｓｒ ｃ ： ｈ ｒ ｉ ｇ ｐ ｉ ｃｐ ｅ ｏ ｉｅ ｔ ｓ ｒ ａ ｐ ｃｒ ｕ ｃ ｍｍｕ ｉａ ｉ ｎ ｓ ｓｅ ａ ａ ａｙ ｅ ｉ ｈ ｉ ｎ ｃ ｔ ｙ ｔ ｍｉ ｎ ｌ ｓ ｄ ｗｈ ｃ ｓ ｏ ｓ ｒ ａ ｅ ｐ ｃｒ ｍ ｙ Ｇ ｌ ｏ ｅ・ ｑ ｅ ｃ ， ｎ ｉ ｌ ｔ ｇ ｔ ｉ ｃ ｍｍｕ ｉ ａ ｏ ｔ ｔ ｂ Ｓ ｍｕ ｉ ｋ ｙ ｐｅｄ ｄ ｓｅ ｔ ｕ ｂ ｏ ｄ ｃ ｄ ｓ ｕ ｎ ｅ ａ ｄ ｓｍｕ ａ ｉ ｈ ｓ ｏ ｅ ｎ ｎ ｃ ｔｎ ｗｉ Ｍａ ｌ ／ ｉ ｌ ．Ｂ ｈ ａ ｎ
无线通信原理与应用-实验二 Gold序列及截短的Gold序列相关特性
实验三、Gold序列及截短的Gold序列相关特性一、实验目的了解常用正交序列--Gold序列及截短Gold序列的自相关及互相关特性。
测量实验系统在异步CDMA工作方式下作为基站地址码的中的截短Gold序列。
二、实验内容1. 用示波器测量常用正交序列--Gold序列及截短Gold序列的波形及其相关运算后的自相关函数及互相关函数，了解其相关特性。
2. 用示波器测量实验系统在异步CDMA工作方式下作为基站地址码的中的截短Gold 序列(长32位)。
三、基本原理见实验一的”三、基本原理”。
下面是本实验待测量的Gold序列及截短Gold序列。
1. Gold序列（1）5阶Gold序列表3-3-1 5阶Gold序列的自相关特性测量（序列长25-1=31位）PN i(t) 0000,0000,1001,0100,1001,1110,1010,110.用实验一表3-1-2相位的二个5阶m序列优选对模二加产生PN j(t) 同上同上表3-3-2 5阶Gold序列的互相关特性测量（序列长25-1=31位）PN i(t) 0000,0000,1001,0100,1001,1110,1010,110.用实验一表3-1-2相位的二个5阶m序列优选对模二加产生PN j(t) 0110,1010,1010,1111,0111,1010,0110,111.用实验一表3-1-2的第一个序列与延时27位(即超前4位)的第二个序列模二加产生这就是本实验系统异步CDMA方式的二个基站地址码，只是相位不同(见式(2-2))。
54（2）7阶Gold序列表3-3-3 7阶Gold序列的自相关特性测量（序列长27-1=127位）PN i(t) 0000,0000,0011,1111,0000,1100,1110,1111,用实验一表3-1-4相位的二个7阶m序列优选对模二加产生0100,1100,0000,0010,0001,1010,1010,1100,0111,1100,1001,0011,0101,1101,0111,0101,0000,1100,0000,1000,0111,1000,1011,010.PN j(t) 同上同上表3-3-4 7阶Gold序列的互相关特性测量（序列长27-1=127位）PN i(t) 0000,0000,0011,1111,0000,1100,1110,1111,用实验一表3-1-4相位的二个7阶m序列优选对模二加产生0100,1100,0000,0010,0001,1010,1010,1100,0111,1100,1001,0011,0101,1101,0111,0101,0000,1100,0000,1000,0111,1000,1011,010.PN j(t) 0001,1110,0111,1010,1001,0001,1010,0000,用实验一表3-1-4第一个序列与延时123位(即超前4位)的第二个序列模二加产生1110,1101,1100,0110,1001,0001,1111,0111,1010,0100,0100,0001,1011,0011,0001,0000,0101,0011,1100,1000,1111,1011,1011,111.2. 截短的Gold序列（1）截短的Gold序列一：Gc1序列表3-3-5 Gc1序列自相关特性测量（序列长32位）PN i(t) 0010,1101,1110,0111,0010,1011,0011,0000.从实验一表3-1-4 PN i的第40位码片开始截取32位PN j(t) 同上同上表3-3-6 Gc1序列互相关特性测量（序列长32位）PN i(t) 0010,1101,1110,0111,0010,1011,0011,0000.从实验一表3-1-4 PN i的第40位码片开始截取32位PN j(t) 0010,1100,1110,1010,0111,1101,0000,1110.从实验一表3-1-4 PN j的第40位码片开始截取32位55（2）截短的Gold序列二：Gc2序列表3-3-7 Gc2序列自相关特性测量（序列长64位）PN i(t) 0010,1101,1110,0111,0010,1011,0111,0000, 从实验一表3-1-4 PN i的第40位码片开始截取64位0110,1101,0111,0100,0110,0100,0100,0000.PN j(t) 同上同上表3-3-8 Gc2序列互相关特性测量（序列长64位）PN i(t) 0010,1101,1110,0111,0010,1011,0111,0000, 从表3-1-4 PNi的第40位码片开始截取64位0110,1101,0111,0100,0110,0100,0100,0000.PN j(t) 0010,1100,1110,1010,0111,1101,0000,1110, 从表3-1-4 PNj的第40位码片开始截取64位0010,0100,1101,1010,1101,1110,1100,0110.四、实验步骤1. 实验箱不要插天线，打开电源。
Gold序列产生及特性分析实验
2、在测试点TP201测试输出的时钟，在测试点TP202、TP203、TP204测试用于产生Gold序列的周期为31的m序列优选对。
实验二、Gold序列产生及特性分析实验
1、实验目的
1、了解Gold序列的性质和特点。
2、熟悉Gold序列的产生方法。
2、实验内容
1、熟悉Gold序列的产生方法。
2、测试Gold序列的波形。
3、实验原理
m序列虽然性能优良，但同样长度的m序列个数不多，且m序列之间的互相关函数并不理想（为多值函数）。1967年，R.Gold提出和讨论了一种新的序列，即Gold序列。这种序列有较为优良的自相关和互相关特性，构造简单，产生的序列数多，因而得到广泛的应用。
TP202测试点输出的m序列为：1 1 0 0 1 1 0 1 1 0 1 1.......
TP203测试点输出的m序列为：1 0 1 0 0 0 1 1 0 1 1 1.......
经验证符合实验结果。
1.m序列优选对
m序列优选对是指在m序列集中，其互相关函数最大值的绝对值满足下式的两条n介m序列：
2.Gold序列的产生方法
Gold序列是m序列的组合序列，由同步时钟控制的两个码元不同的m序列优选对逐位模2加得到。这两个序列发生器的周期相同，速率相同，因而两者保持一定的相位关系，这样产生的组合序列与这两个自序列的周期也相同。当改变两个序列的相对位移，会得到一个新的Gold序列。Gold序列具有以下性质：
（1）两个m序列优选对经不同移位相加产生的新序列都是Gold序列，两个n级移位寄存器可以产生 个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:Nnewregister(1)= mod(sum(fbconnection.*register),2); %寄存器与反馈的模2和for j = 2:n,newregister(j)= register(j-1);end;register = newregister; %移位后的寄存器mseq(i) = register(n); %新的寄存器输出endclear 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;endclear all;close all;clc;connection1=[0 0 0 0 1 1];connection2=[1 0 0 1 1 1];goldseq = g_sequence(connection1,connection2);stem(goldseq);axis([0 65 -0.2 1.2]);grid on;THANKS !!!致力为企业和个人提供合同协议，策划案计划书，学习课件等等打造全网一站式需求欢迎您的下载，资料仅供参考。
m序列和Gold序列特性研究要点上课讲义
m序列和G o l d序列特性研究要点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 =++++∈L (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 ----=++++∈L (2) 特征多项式(1)与递归多项式(2)是r 级线性移位寄存器反馈逻辑的两种不同种表示法，因其应用的场合不同而采用不同的表示方法。
m序列和Gold序列特性研究
扩频通信实验报告Harbin Institute of Technology扩频通信实验报告课程名称：扩频通信实验题目：Gold码特性研究院系：电信学院班级：通信一班姓名：学号：指导教师：迟永钢时间： 2012年5月8日哈尔滨工业大学- I-第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 级线性移位寄存器反馈逻辑的两种不同种表示法，因其应用的场合不同而采用不同的表示方法。
Gold码
通信08-1 艾盼盼0850283101设计Gold序列发生器姓名：艾盼盼学号：0850283101 班级：通信08-1摘要：m序列，尤其是m序列优选对，是特性很好的伪随机序列。
但是，它们能彼此构成优选对的数目很少，不便于在码分多址系统中应用。
R.Gold于1967年提出了一种基于m 序列优选对的码序列，称为Gold序列。
它是m序列的组合码，由优选对的两个m序列逐位模2加得到，当改变其中一个m序列的相位（向后移位）时，可得到一新的Gold序列。
Gold 序列虽然是由m序列模2加得到的，但它已不是m序列，不过它具有与m序列优选对类似的自相关和互相关特性，而且构造简单，产生的序列数多，因而获得广泛的应用。
【关键词】：m序列优选对，Gold序列，模2加，自相关1. Gold码的概述1.1 gold码定义R.Gold于1967年提出了一种基于m序列优选对的码序列，称为Gold序列。
它是m序列的组合码，由优选对的两个m序列逐位模2加得到，当改变其中一个m序列的相位（向后移位）时，可得到一新的Gold序列。
Gold序列虽然是由m序列模2加得到的，但它已不是m 序列，不过它具有与m序列优选对类似的自相关和互相关特性，而且构造简单，产生的序列数多，因而获得广泛的应用。
1.2 gold码基本功能单元Gold码发生器的基本功能单元为线性反馈移位寄存器LFSR（Linear Fdddback Bhift Register）。
2.Gold序列的设计2.1 m序列优选对寻找方法产生gold序列的必要条件是m序列优选对，设A是对应于n级本原多项式f（x）所产生的m序列，B是对应于n级本原多项式g（x）所产生的m序列，当它们的互相关函数|Ra.b（k）|满足：则f（x）和g（x）所产生的m序列A和B构成一对优选对。
寻找m序列优选对的方法还有硬件计算法，分圆陪集法，逐步移位模2加法，三值判别法。
2.2gold序列设计的理论证明证明，若F1（x），F2（x）为两个不同的本原多项式，令F1（x）产生的序列为G（F1），F2（x）产生的序列为G（F2），F1（x）. F2（x）所产生的序列为G（F1，F2），则有上式表明两本原多项式乘积所产生的序列等于两个本原多项式分别产生的模2和序列。
移动通信原理 实验报告
南京邮电大学实验报告实验名称__CDMA扩频与解扩_ 呼叫实验_____课程名称现代移动通信 _ _班级学号姓名开课时间 2011 /2012 学年，第二学期实验一 CDMA扩频与解扩一、实验目的1. 了解扩频调制的基本概念；2.掌握PN码的概念以及m序列的生成方法；3.掌握扩频调制过程中信号频谱的变化规律。
4. 了解CDMA解扩的基本概念；5. 掌握解扩的基本方法；6. 掌握解扩过程中信号频谱的变化规律。
二、实验设备1. 移动通信实验机箱一台2. 微型计算机一台三、实验原理1. 扩频实验原理m序列是最长线性反馈移位寄存器序列的简称，它是由带线性反馈的移位器产生的周期最长的一种序列。
如果把两个m序列发生器产生的优选对序列模二相加，则产生一个新的码序列，即Gold码序列。
实验中三种可选的扩频序列分别是长度为15的m序列、长度为31的m序列以及长度为31的Gold序列。
1.长度为15的m序列由4级移存器产生，反馈器如图所示。
初始状态 1 0 0 01 1 0 01 1 1 01 1 1 10 1 1 11 0 1 10 1 0 11 0 1 01 1 0 10 1 1 00 0 1 11 0 0 10 1 0 00 0 1 00 0 0 1……………………………….1 0 0 02.长度为31的m 序列由5级移存器产生，反馈器如图所示。
a4a3a2a1+a03. 长度为31的gold 序列:Gold 码是Gold 于1967年提出的，它是用一对优选的周期和速率均相同的m 序列模二加后得到的。
其构成原理如图2.1.3所示。
两个m 序列发生器的级数相同，即n n n ==21。
如果两个m 序列相对相移不同，所得到的是不同的Gold 码序列。
对n 级m 序列，共有12-n 个不同相位，所以通过模二加后可得到12-n 个Gold 码序列，这些码序列的周期均为12-n ，如图2.1.4所示。
两组数据为： 1 0 0 0 0 1 0 0 0 00 1 0 0 0 0 1 0 0 00 0 1 0 0 0 0 1 0 01 0 0 1 0 0 0 0 1 0m 序列发生器 n 级 m 序列发生器n 级 初态设置 时钟 Gold 码 21m m ⊕ 1m 2m0 1 0 0 1 0 0 0 0 11 0 1 0 0 1 1 1 0 11 1 0 1 0 1 0 0 1 10 1 1 0 1 1 0 1 0 00 0 1 1 0 0 1 0 1 01 0 0 1 1 0 0 1 0 11 1 0 0 1 1 1 1 1 11 1 1 0 0 1 0 0 1 01 1 1 1 0 0 1 0 0 11 1 1 1 1 1 1 0 0 10 1 1 1 1 1 0 0 0 10 0 1 1 1 1 0 1 0 10 0 0 1 1 1 0 1 1 11 0 0 0 1 1 0 1 1 01 1 0 0 0 0 1 0 1 10 1 1 0 0 1 1 0 0 01 0 1 1 0 0 1 1 0 01 1 0 1 1 0 0 1 1 01 1 1 0 1 0 0 0 1 10 1 1 1 0 1 1 1 0 01 0 1 1 1 0 1 1 1 00 1 0 1 10 0 1 1 11 0 1 0 1 1 1 1 1 00 1 0 1 0 0 1 1 1 10 0 1 0 1 1 1 0 1 00 0 0 1 0 0 1 1 0 10 0 0 0 1 1 1 0 1 1……………………………………………………………………………….所以生成长度为31的Gold序列为：{0,0,0,0,0,1,1,1,0,0,0,0,1,0,0,0,0,1,1,0,0,1,0,0,1,0,1,1,1,1,0} 在硬件上，扩频调制是通过单片机和学生平台软件联合实现的。
基于Gold序列的直接扩频通信系统仿真
基于Gold序列的直接扩频通信系统仿真李亚文【摘要】The working principle of direct spread spectrum communication systemia analysed which is spreaded spectrum by Gold code sequence,and simulating this communicaton with Matlab/Simulink.By the experiment,it follows the complete process of spread spectrum communication.The result of simulation indicates that the Direct spread spectrum communication system base on Gold sequence has advantages of good performance of anti-interference and simple circuit,it can implementate users need of communication transmission and have the better reliability.%分析了直接序列扩频通信系统的工作原理,选用Gold序列作为扩频序列码,利用MATLAB/Simulink对直接扩频通信系统进行仿真,通过仿真实验,实时跟踪了扩频通信的完整过程。
分析实验结果可知,基于Gold序列的直接扩频通信系统抗干扰性较好,电路实现简单,能实现用户通信传输的需要,并且可靠性较好。
【期刊名称】《商洛学院学报》【年(卷),期】2012(026)002【总页数】4页(P51-54)【关键词】Gold序列;直接扩频通信系统;simulink;可靠性【作者】李亚文【作者单位】商洛学院物理与电子信息工程系,陕西商洛726000【正文语种】中文【中图分类】O453通信已经成为人们传递信息的主要途径，为了提高通信过程中信息传输的速率，改进通信系统传输的有效性，因而引入扩频通信技术。
黄金的实验报告
一、实验目的1. 了解黄金的物理性质，包括颜色、密度、延展性、导电性等。
2. 掌握黄金的化学性质，包括耐腐蚀性、与酸碱反应等。
3. 熟悉黄金的提取、提纯方法。
二、实验原理黄金（Au）是一种化学性质非常稳定的金属，具有良好的导电性、导热性和延展性。
黄金的密度较大，约为19.32g/cm³。
在常温下，黄金不与空气中的氧气、氮气、二氧化碳等气体发生反应，也不与大多数酸碱反应，具有良好的耐腐蚀性。
三、实验仪器与试剂1. 仪器：天平、量筒、试管、酒精灯、蒸馏烧瓶、冷凝管、烧杯、玻璃棒、试管架、铁架台、加热套、试管夹等。
2. 试剂：黄金、盐酸、氢氧化钠、稀硫酸、硝酸、氯水、硫酸铜溶液、硫酸锌溶液等。
四、实验步骤1. 黄金的物理性质测定（1）颜色观察：将黄金放在白纸上，观察其颜色，记录下来。
（2）密度测定：使用天平称取一定质量的黄金，用量筒量取一定体积的水，将黄金放入量筒中，记录水的体积变化，根据密度公式计算黄金的密度。
（3）延展性测定：将黄金用锤子锤成薄片，观察其变化，记录下来。
（4）导电性测定：将黄金制成导线，接入电路中，观察电流表指针的偏转情况，记录下来。
2. 黄金的化学性质测定（1）耐腐蚀性测定：将黄金放入盐酸、氢氧化钠、稀硫酸、硝酸等溶液中，观察黄金的变化，记录下来。
（2）与酸碱反应：将黄金放入氯水、硫酸铜溶液、硫酸锌溶液等溶液中，观察黄金的变化，记录下来。
3. 黄金的提取与提纯（1）黄金的提取：将含有黄金的矿石破碎、磨细，加入硝酸溶解，过滤、洗涤、干燥，得到含有黄金的沉淀。
（2）黄金的提纯：将含有黄金的沉淀用氯水氧化，过滤、洗涤、干燥，得到纯黄金。
五、实验结果与分析1. 黄金的物理性质（1）颜色：黄金呈金黄色。
（2）密度：实验测得黄金的密度为19.32g/cm³。
（3）延展性：黄金可锤成薄片。
（4）导电性：黄金具有良好的导电性。
2. 黄金的化学性质（1）耐腐蚀性：黄金在盐酸、氢氧化钠、稀硫酸、硝酸等溶液中均无变化。
m序列对及平衡Gold序列的产生与搜索
平衡Gold序列的产生与搜索方法
平衡Gold序列是一类具有优良伪随机性和自相关性的二元序列。产生平衡 Gold序列的方法主要有两种：一是通过平衡Gold多项式产生平衡Gold序列； 二是通过选择适当的m序列和反馈函数产生平衡Gold序列。搜索方法主要有基 于差分编码的搜索算法和基于自相关函数的搜索算法。
3、信息隐藏：平衡Gold序列可以用于信息隐藏中的隐写术，以隐藏敏感信息 在公开信息中。由于平衡Gold序列具有优良的自相关性，它可以用于构造高效 的隐写算法，提高信息隐藏的安全性和鲁棒性。
未来展望
随着技术的发展，平衡Gold序列的研究将面临新的挑战和机遇。以下是一些未 来可能的研究方向：
1、高效产生方法：目前平衡Gold序列的产生方法还比较耗时，因此研究更高 效的产生方法以提高搜索速度和降低计算成本是未来的一个研究方向。
1、反馈函数设计：反馈函数是m序列发生器的核心部分，其设计的好坏直接影 响到m序列的性能。一个好的反馈函数应该具有较低的误码率、较高的线性复 杂度、易于实现等优点。常用的反馈函数包括异或、同或、与等运算。
2、LFSR级数和反馈位数选择：m序列发生器的性能与LFSR的级数和反馈位数 密切相关。一般来说，增加LFSR的级数和反馈位数可以提高m序列的性能，但 同时也会增加硬件实现的复杂度和功耗。因此，在设计中需要根据实际需求选 择合适的LFSR级数和反馈位数。
在软件实现中，可以使用各种编程语言如C、C++、Python等编写m序列发生器。 软件实现具有灵活度高、易于调试和修改等优点，但也存在运行速度较慢和需 要运行额外的处理器等缺点。
m序列发生器的应用
m序列发生器在信号处理、通信技术等领域具有广泛的应用。以下是m序列发生 器的一些典型应用：
gold序列实验资料报告材料
Gold 序列实验报告一、实验过程及分析1 m 序列&G old 序列的产生（1）选取本原多项式367()1f x x x x x =++++，产生m 序列。
fbconnection_m=[1 0 1 0 0 1 1]; mseq=m_sequence(fbconnection_m);（2）由本原多项式37237()1&()1f x x x f x x x x x =++=++++生成的m 序列为m 序列优选对，以此优选对来产生Gold 序列。
fbconnection_op_1=[0 0 1 0 0 0 1]; fbconnection_op_2=[1 1 1 0 0 0 1];goldseq=gold_seq(fbconnection_op_1,fbconnection_op_2);2 m 序列的相关特性-150-100-50050100150kR (k )autocorrelation function of m sequence图1 m 序列自相关函数kR (k )图 2 m 序列互相关函数3 G old 序列的相关特性kR (k )autocorrelation function of gold sequence图3 G old 序列自相关函数-20-15-10-5051015kR (k )cross-correlation function of gold sequence图4 G old 序列互相关函数图2-5 m 序列&G old 序列最小，最大相关值4 m 序列&G old 序列相关特性比较-200-1000100200-20020406080100120140kR (k )autocorrelation function of m sequence -200-1000100200-20020406080100120140kR (k )autocorrelation function of gold sequence图6 m 序列&G old 序列自相关函数比较-200-1000100200-50-40-30-20-100102030kR (k )cross-correlation function of m sequence -200-1000100200-20-15-10-551015kR (k )cross-correlation function of gold sequence图7 m 序列&G old 序列互相关函数比较gold 序列在时延不为0是三值的，这一点不如m 序列，在时延为0处具有与m 序列相同的峰值特性。
实验一 GOLD序列特性实验
实验 GOLD 序列特性实验一、实验目的1、掌握GOLD 序列的特点。
2、了解GOLD 序列在直接扩频通信中所起的作用二、实验器材1、移动通信原理实验箱一台 2、20M 双踪示波器一台 3、频谱分析仪或带FFT 功能的数字示波器（选配） 一台三、实验内容1、观察GOLD 序列的波形（频谱）。
2、观察GOLD 序列的自相关和互相关特性。
四、实验原理1、伪随机序列工程上常用二元{0，1}序列来产生伪噪声码。
它具有如下特点：（1） 每一周期内“0”和“1”出现的次数近似相等。
（2） 每一周期内，长度为n 比特的游程出现的次数比长度为n+1比特的游程出现的次数多一倍。
（游程是指相同码元的码元串）（3） 序列具有双值自相关函数，即：⎪⎩⎪⎨⎧-≤≤=11k 01)(p p R τττ当－＝当（4.1-1）在（4.1-1）式中，p 为二元序列周期，又称码长，k 为小于p 的整数，τ为码元延时。
2、m 序列二元m 序列是一种基本的伪随机序列，有优良的自相关函数，易于产生和复制，在扩频技术中得到了广泛的应用。
长度为2n -1位的m 序列可以用n 级线性移位寄存器来产生。
如图4.1-1所示：图4.1-1 线性移位寄存器m 序列的特性如下（1） 在每一周期p= 2n -1内，“0”出现2n －1-1次，“1”出现2n －1次，“1”比“0”多出现一次。
（2） 在每一周期内共有2n －1个元属游程，其中“0”的游程和“1”的游程数目各占模二加法器一半。
并且，对n>2，当1≤k ≤n-1时，长为k 的游程占游程总数的1/ 2 k ，其中“0”的游程和“1”的游程各占一半。
长为n –1的游程只有一个，为“0”的游程；长为n 的游程也只有一个，为“1”的游程。
（3） m 序列（a k ）与其位移序列（τ-k a ）的模二和仍然是m 序列的另一位移序列（τ'-k a ），即：{}{}{}ττ'--=+k k k a a a（4） m 序列的自相关函数为：⎪⎩⎪⎨⎧≠-=p p p R mod 01mod 01)(τττ当＝当 （4.1-2）3、Gold 序列虽然m 序列有优良的自相关特性，但是使用m 序列作CDMA （码分多址）通信的地址码时，其主要问题是由m 序列组成的互相关特性好的互为优选的序列集很少，对于多址应用来说，可用的地址数太少了。
扩频通信实验报告
中南大学扩频通信实验报告实验一：扩频与解扩观测实验时间：4月9号一、实验目的1、了解直接序列扩频的原理。
2、了解扩频前后信号在时域及频域上的变化。
二、实验器材⒈主控&信号源模块、2号、14号、11号模块各一块⒉双踪示波器一台⒊连接线若干三、实验原理1、实验原理框图数字信号源扩频解扩DoutMUXBSOUT2# 模块14# 模块11# 模块NRZ1NRZ-CLK1扩频序列1序列1（TP8）扩频序列2序列2（TP8）CDMA1AD 输入1AD 输入2CDMA2Dout实验框图2、实验框图说明本实验选择【扩频与解扩观测实验】菜单。
如框图所示，我们用2号模块作为信号源，DoutMUX 输出32K 数字信号，送入至14号模块的NRZ1。
14号模块此时完成扩频功能，扩频序列由14号模块内部产生，将开关S1设置为0000，开关S2设置为0111，即可设置该路扩频序列1的码型（测试点为TP8序列1）。
扩频信号由端口CDMA1输出。
同时，当14号模块的开关S3设置为0111、开关S4设置为0000且端口NRZ2和NRZ-CLK2无信号输入时，端口CDMA2输出的伪随机序列与14号模块的扩频序列1相同，本实验中将该序列“CDMA2”可作为后续的解扩序列。
此时的11号模块完成解扩功能，其中扩频信号从端口“AD 输入1”输入，解扩序列从“AD 输入2”输入，解扩信号从11号模块的“Dout ”输出。
该实验【扩频与解扩观测实验】中扩频序列的长度可通过PN 序列长度设置开关S6进行选择15位或16位。
当开关S6拨至“127位”时，表示该实验的扩频为15位；当开关S6拨至“128位”时，表示该实验的扩频为16位。
注：为配合示波器调节，为了较好的对比观测扩频前和扩频后的码元，建议选择16位。
四、实验步骤1、按框图所示连线。
源端口目标端口连线说明 模块2：DoutMUX模块14：TH3(NRZ1)数据送入扩频单元模块2：BSOUT 模块14：TH1(NRZ-CLK1) 时钟送入扩频单元模块14：TH4(CDMA1) 模块11：TH2(AD输入1) 送入解扩单元模块14：TH5(CDMA2) 模块11：TH3(AD输入2) 提供解扩序列2、选择主菜单【移动通信】→【扩频与解扩观测实验】，此时2号模块DoutMUX输出速率为32K。
GOLD序列特性实验
实验一GOLD序列特性实验②用示波器观察测试点“GOLD1”处的波形。
改变拨位开关“扩频码速率”的设置，按“发射机复位”键，再观察“GOLD1”处的波形。
“GOLD1置位”设置为10000000时100kbit/s的GOLD1波形“GOLD1置位”设置为10000000时200kbit/s的GOLD1波形观察Gold序列的自相关和互相关特性用示波器测“TX3”处波形，该波形即为Gold序列的GOLD序列的自相关特性用示波器测“TX3”处波形，该波形即为Gold序列的互相关特性。
GOLD序列的互相关特性1、顺时针将“跟踪”电位器旋到底，用示波器测“VCO-C”处波形，该波形即为延迟锁相环的鉴相特性曲线。
2、用示波器双踪分别观察“G1-BS”和“G3-BS”处的波形，调节“跟踪”旋钮，直到二个波形完全一致，没有相差为止。
此时表明接收机的Gold序列和发射机的Gold序列在相位与码速率上都一致。
3、用示波器双踪分别观察“GOLD1”和“GD-TX”处的波形，二者的波形应完全一致。
说明：由于本系统的Gold序列频率较高，且周期很长，模拟双踪示波器应在“断续（CHOP）”模式下比较“GOLD1”和“GD-TX”处的波形，如果在“交替（ALT）”模式下即使两者输出波形一致，观察结果也可能不一致。
数字示波器则不存在该问题。
实验三扩频与解扩实验①将“SIGN1置位”设置成不为全0或全1的码字，设置“GOLD1置位”。
用示波器分别观察“SIGN1”和“S1-KP”的波形，并做对比。
信码速率为1kbit/s、扩频码速率为100kbit/s时“SIGN1”和“S1-KP”处波形②（选做）用带FFT功能的数字示波器分别观察“SIGN1”和“S1-KP”的频谱，并做对比。
③分别改变发射机的信码速率和扩频码速率，重复上一步骤。
5、（选做）观察扩频前后PSK调制频谱的实验①码字设置不变，将“扩频”开关拨下，用频谱仪观察“PSK1”的频谱。
GOLD码产生与特性分析实验
实验八GOLD 码特性实验一、实验目的1、掌握GOLD 码的编解码原理。
2、掌握GOLD 码的软件仿真方法。
3、掌握GOLD 码的硬件仿真方法。
4、掌握GOLD 码的硬件设计方法。
二、预习要求1、掌握GOLD 码的编解码原理和方法。
2、熟悉matlab 的应用和仿真方法。
3、熟悉Quatus 的应用和FPGA 的开发方法。
三、实验原理1、GOLD 序列简介GOLD 序列是由m 序列的“优选对”构成的。
所谓优选对是指m 序列中互相关值为[－1，－t(n),t(n)-2]的一对序列。
其中(1)/21(2)/212,2,(){n n n n t n ++++=为奇数，为偶数下表为部分m 序列的部分优选对 表1 部分m 序列的部分优选对上表中的m 序列采用8进制（可参见PN 码实验）。
2、GOLD 序列由m 序列中的优选对{xi}和{yi}本身加上它们的相对移位模二相加构成的2n －1个序列组成，序列总数为2n ＋1。
任一队序列之间的互相关函数都是三值的，即(1)/2(1)/21(1)/2(1)/211{(21)(),(21)()2R {1{(21)(),(21)()2n n cn n t n n t n t n n t n ++++++--+=--=-=--+=--=-为奇数为偶数，但不被4整除即，GOLD 序列的最大互相关值为|R |()c m t n下表为GOLD 序列的t(n)值及其与自相关峰值Rs （0）的比值，同时给出GOLD 序列族中的序列数。
表为 部分GOLD 序列的t(n)值、Rs （0）、序列数表1、建立GOLD 的仿真文件(GOLD.MDL)GOLD1…GOLD7的Sample Time 均设置为SampleTime ；Preferred polynomial （1）设置为[1 0 1 1]；Initial states （1）设置为[0 0 1]; Preferred polynomial （2）设置为[1 1 0 1]；Initial states （2）设置为[0 0 1]。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Gold 序列实验报告一、实验过程及分析1 m 序列&G old 序列的产生（1）选取本原多项式367()1f x x x x x =++++，产生m 序列。
fbconnection_m=[1 0 1 0 0 1 1]; mseq=m_sequence(fbconnection_m);（2）由本原多项式37237()1&()1f x x x f x x x x x =++=++++生成的m 序列为m 序列优选对，以此优选对来产生Gold 序列。
fbconnection_op_1=[0 0 1 0 0 0 1]; fbconnection_op_2=[1 1 1 0 0 0 1];goldseq=gold_seq(fbconnection_op_1,fbconnection_op_2);2 m 序列的相关特性-150-100-50050100150kR (k )autocorrelation function of m sequence图1 m 序列自相关函数kR (k )图 2 m 序列互相关函数3 G old 序列的相关特性kR (k )autocorrelation function of gold sequence图3 G old 序列自相关函数-20-15-10-5051015kR (k )cross-correlation function of gold sequence图4 G old 序列互相关函数图2-5 m 序列&G old 序列最小，最大相关值4 m 序列&G old 序列相关特性比较-200-1000100200-20020406080100120140kR (k )autocorrelation function of m sequence -200-1000100200-20020406080100120140kR (k )autocorrelation function of gold sequence图6 m 序列&G old 序列自相关函数比较-200-1000100200-50-40-30-20-100102030kR (k )cross-correlation function of m sequence -200-1000100200-20-15-10-551015kR (k )cross-correlation function of gold sequence图7 m 序列&G old 序列互相关函数比较gold 序列在时延不为0是三值的，这一点不如m 序列，在时延为0处具有与m 序列相同的峰值特性。
显然，m 序列的自相关特性比Gold 序列的自相关特性要好。
Gold 序列具有更小的互相关峰值，显然Gold 序列的互相关特性比m 序列的互相关特性性能要好。
5 m 序列优选对m 序列优选对也是gold 序列族中的gold 序列，所以其互相关函数与gold 序列一致。
kR (k )cross-correlation function of m sequence optimum pairs图8 m 序列优选对的互相关函数6 任选6个G old 序列，观察其自相关函数kR (k )kR (k )kR (k )kR (k )kR (k )kR (k )图9 任选6个G old 序列的自相关函数7 任选6对G old 序列，观察其互相关函数kR (k )=17kR (k )=17kR (k )=17kR (k )=17kR (k )=17kR (k )=17图 10 任选6对G old 序列的互相关函数任意两个Gold 序列互相关函数幅度小于或等于17，符合Gold 理论。
源码(1) PN_Correlation.m% ======================================= % PN_Correlation.m% Correlation analysis of PN sequence %% Author: jzh (04011646 seu) % Date: 2014.12.29% =======================================% cleanclose all;clear all;clc;% optimum pairsfbconnection_op_1=[0 0 1 0 0 0 1];fbconnection_op_2=[1 1 1 0 0 0 1];% any m sequence for testfbconnection_m=[1 0 1 0 0 1 1];% gold sequencemseq_op_1=m_sequence(fbconnection_op_1);mseq_op_2=m_sequence(fbconnection_op_2);goldseq=gold_seq(fbconnection_op_1,fbconnection_op_2);% m sequencemseq=m_sequence(fbconnection_m);% mappingmseq_op_1((mseq_op_1==0))=-1;mseq_op_2((mseq_op_2==0))=-1;goldseq((goldseq==0))=-1;mseq((mseq==0))=-1;m_auto=ccorr(mseq);% any two m sequences except optimum pairsm_cross=ccorr(mseq,mseq_op_1);gold_auto=ccorr(goldseq(1,:));gold_cross=ccorr(goldseq(1,:),goldseq(2,:));m_op_cross=ccorr(mseq_op_1,mseq_op_2);m_auto_max = max(m_auto)m_auto_min = min(m_auto)gold_auto_max = max(gold_auto)m_cross_amplitude_max = max(abs(m_cross)) gold_cross_amplitude_max = max(abs(gold_cross))N=length(mseq); % sequence lengthabscissa=-N+1:N-1; % abscissa axisfigure(1);plot(abscissa,m_auto);xlabel('k'); ylabel('R(k)');legend('autocorrelation function of m sequence');title('autocorrelation function of m sequence');figure(2);plot(abscissa,m_cross);xlabel('k'); ylabel('R(k)');legend('cross-correlation function of m sequence'); title('cross-correlation function of m sequence');figure(3);plot(abscissa,gold_auto);xlabel('k'); ylabel('R(k)');legend('autocorrelation function of gold sequence'); title('autocorrelation function of gold sequence');figure(4);plot(abscissa,gold_cross);xlabel('k'); ylabel('R(k)');legend('cross-correlation function of gold sequence'); title('cross-correlation function of gold sequence');figure(5);subplot(1,2,1);plot(abscissa,m_auto);xlabel('k'); ylabel('R(k)');title('autocorrelation function of m sequence');subplot(1,2,2);plot(abscissa,gold_auto);xlabel('k'); ylabel('R(k)');title('autocorrelation function of gold sequence');figure(6);subplot(1,2,1);plot(abscissa,m_cross);xlabel('k'); ylabel('R(k)');title('cross-correlation function of m sequence');subplot(1,2,2);plot(abscissa,gold_cross);xlabel('k'); ylabel('R(k)');title('cross-correlation function of gold sequence');figure(7);plot(abscissa,m_op_cross);xlabel('k'); ylabel('R(k)');legend('cross-correlation function of m sequence optimum pairs'); title('cross-correlation function of m sequence optimum pairs');% test autocorrelation of any gold sequencefor test_time = 1:6ii = ceil(rand()*(N+2));figure();temp_gold_auto = ccorr(goldseq(ii,:));plot(abscissa,temp_gold_auto);xlabel('k'); ylabel('R(k)');legend(strcat( 'seq number = ',num2str(ii)));title('autocorrelation function of a random gold sequence');end% test cross-correlation of any 2 gold sequencefor test_time = 1:6ii = ceil(rand()*(N+2));jj=ceil(rand()*(N+2));while jj == iijj=ceil(rand()*(N+2));endfigure();temp_gold_cross = ccorr(goldseq(ii,:),goldseq(jj,:));plot(abscissa,temp_gold_cross);xlabel('k'); ylabel('R(k)');legend(strcat( 'seq1 number=',num2str(ii),', seq2 number=',num2str(jj)));title(strcat('cross-correlation function of 2 random gold sequence:',...' cross correlation amplitude MAX = ',...num2str(max(abs(temp_gold_cross)))));end(2) m_sequence.m% ==============================% m_sequence.m% generating m sequence%% Author: jzh (04011646 seu)% Date: 2014.12.29% ==============================function mseq = m_sequence(fbconnection)% ============================== m_sequence ============================== % generating m sequence%% Input:% fbconnection: feedback connection responding to primitive polynomial% Output:% msequ: generation of m sequence, 1-by-N% ========================================================================n = length(fbconnection);N = 2^n-1; % length of m sequenceregister = [zeros(1,n-1) 1]; % initial state of shift register% preallocatingnewregister = zeros(1,n);mseq = zeros(1,N);mseq(1)=register(n);for ii=2:N% feedbacknewregister(1) = mod(sum(fbconnection.*register),2);% shiftfor jj=2:nnewregister(jj)=register(jj-1);endregister=newregister; % updatemseq(ii)=register(n); % output m sequenceend(3) gold_seq.m% ==================================% gold_seq.m% generating gold sequence family%% Author: jzh (04011646 seu)% Date: 2014.12.29% ==================================function goldseq = gold_seq(fbconnection1,fbconnection2)% =============================== gold_seq =============================== % generating gold sequence family%% Input:% fbconnection1: feedback connection 1% fbconnection2: feedback connection 2% Note: fbconnection1 & fbconnection2 are optimum pairs% Output:% goldseq: generation of gold sequence family, (N+2)-by-N% ========================================================================mseq1 = m_sequence(fbconnection1);mseq2 = m_sequence(fbconnection2);N=2^length(fbconnection1) - 1;goldseq = zeros(N+2,N); % preallocating% different shift amount of two m sequencefor shift_amount = 0:N-1shift_mseq2 = [mseq2(shift_amount+1:N) mseq2(1:shift_amount)];goldseq(shift_amount+1,:) = mod(mseq1+shift_mseq2,2);endgoldseq(N+1,:)=mseq1;goldseq(N+2,:)=mseq2;(4) ccorr.m% =======================================% ccorr.m% calculate cross-correlation function%% Author: jzh (04011646 seu)% Date: 2014.12.29% =======================================function r = ccorr(seq1,seq2)% ================================ ccorr ================================= % calculate cross-correlation function%% Input:% seq1: sequence 1% seq2: sequence 2% Note: (a) length of two sequence must be the same% (b) two input sequence: calculate cross-correlation function% (c) one input sequence: calculate autocorrelation function% Output:% r: correlation function, 1-by-2*N-1% ========================================================================% one input sequence: calculate autocorrelation if nargin==1seq2 = seq1;endN = length(seq1);r = zeros(1,2*N-1); % preallocatingfor k=-N+1:-1seq2_shift = [seq2(k+N+1:N) seq2(1:k+N)];r(N+k) = seq1*seq2_shift';endfor k=0:N-1seq2_shift = [seq2(k+1:N) seq2(1:k)];r(N+k) = seq1*seq2_shift';end。