基于MATLAB的直接序列扩频通信系统性能仿真分析研究

直接序列扩频通信系统仿真

直接序列扩频通信系统仿真

直接序列扩频通信系统仿真 一、实验的背景及内容 1、直接扩频通信背景 扩频通信，即扩展频谱通信(Spread Spectrum Communication)，它与光纤通信、卫星通信，一同被誉为进入信息时代的三大高技术通信传输方式。 有关扩频通信技术的观点是在1941年由好莱坞女演员Hedy Lamarr和钢琴家George Antheil提出的。解决了短距离数据收发信机、如：卫星定位系统(GPS)、移动通信系统、WLAN(IEEE802.11a, IEEE802.11b, IEE802.11g)和蓝牙技术等应用的关键问题。扩频技术也为提高无线电频率的利用率(无线电频谱是有限的因此也是一种昂贵的资源)提供帮助。 扩频通信技术自50年代中期美国军方便开始研究，一直为军事通信所独占，广泛应用于军事通信、电子对抗以及导航、测量等各个领域。直到80年代初才被应用于民用通信领域。为了满足日益增长的民用通信容量的需求和有效地利用频谱资源，各国都纷纷提出在数字峰窝移动通信、卫星移动通信和未来的个人通信中采用扩频技术，扩频技术现已广泛应用于蜂窝电话、无绳电话、微波通信、无线数据通信、遥测、监控、报警等等的系统中。 2、实验的内容及意义 本次实验主要研究了直接序列扩频系统，建立了直接序列扩频系统的matlab仿真模型，在信道中存在高斯白噪声和干扰的情况下，对系统误码率性能进行了仿真及分析。 近年来，随着超大规模集成电路技术、微处理器技术的飞速发展，以及一些新型元器件的应用，扩频通信在技术上已迈上了一个新的台阶，不仅在军事通信中占有重要地位，而且正迅速地渗透到了个人通信和计算机通信等民用领域，成为新世纪最有潜力的通信技术之一因此研究扩频通信具有很深远的意义。本人通过此次实验，进行深入地研究学习扩频通信技术及对它进行仿真应用，将所学的知识进行归纳与总结，从而巩固通信专业基础知识，为以后的个人学习和工作打下基础。

MATLAB中产生高斯白噪声

MATLAB中产生高斯白噪声，涉及到awgn和wgn函数 MATLAB中产生高斯白噪声非常方便，可以直接应用两个函数，一个是WGN，另一个是AWGN。WGN用于产生高斯白噪声，AWGN则用于在某一信号中加入高斯白噪声。 1. WGN：产生高斯白噪声 y = wgn(m,n,p) 产生一个m行n列的高斯白噪声的矩阵，p以dBW为单位指定输出噪声的强度。 y = wgn(m,n,p,imp) 以欧姆(Ohm)为单位指定负载阻抗。 y = wgn(m,n,p,imp,state) 重置RANDN的状态。 在数值变量后还可附加一些标志性参数： y = wgn(…,POWERTYPE) 指定p的单位。POWERTYPE可以是'dBW', 'dBm'或 'linear'。线性强度(linear power)以瓦特(Watt)为单位。 y = wgn(…,OUTPUTTYPE) 指定输出类型。OUTPUTTYPE可以是'real'或 'complex'。 2. AWGN：在某一信号中加入高斯白噪声 y = awgn(x,SNR) 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数，就加入复噪声。 y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值，则其代表以dBW为单位的信号强度；如果SIGPOWER为'measured'，则函数将在加入噪声之前测定信号强度。y = awgn(x,SNR,SIGPOWER,STATE) 重置RANDN的状态。 y = awgn(…,POWERTYPE)指定SNR和SIGPOWER的单位。POWERTYPE可以是'dB'或'linear'。如果POWERTYPE是'dB'，那么SNR以dB为单位，而SIGPOWER以dBW为单位。如果POWERTYPE是'linear'，那么SNR作为比值来度量，而SIGPOWER 以瓦特为单位。 注释 1. 分贝(decibel,dB)：分贝（dB）是表示相对功率或幅度电平的标准单位，换句话说，就是我们用来表示两个能量之间的差别的一种表示单位，它不是一个绝对单位。例如，电子系统中将电压、电流、功率等物理量的强弱通称为电平，电平的单位通常就以分贝表示，即事先取一个电压或电流作为参考值（0dB），用待表示的量与参考值之比取对数，再乘以20作为电平的分贝数（功率的电平值改乘10）。 2. 分贝瓦(dBW, dB Watt)：指以1W的输出功率为基准时，用分贝来测量的功率放大器的功率值。 3. dBm (dB-milliWatt)：即与1milliWatt（毫瓦）作比较得出的数字。 0 dBm = 1 mW 10 dBm = 10 mW 20 dBm = 100 mW 也可直接用randn函数产生高斯分布序列，例如： 程序代码 y=randn(1,2500); y=y/std(y);

matlab实现：常见的离散时间信号

1． 单位抽样序列，或称为离散时间冲激，单位冲激： ? ??=01)(n δ 00≠=n n 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： ???=-01)(k n δ 0≠=n k n 2．单位阶跃序列 ? ??01)(n u 00<≥n n 在MATLAB 中可以利用ones( )函数实现。 );,1(N ones x = 3．正弦序列 )(cos )(0φω+=n A n x 这里， ,,0ωA 和φ都是实数，它们分别称为本正弦信号)(n x 的振幅，角频率和初始相位。 πω200=f 为频率。 4．复正弦序列 n j e n x ω=)( 5．实指数序列 n A n x α=)( 6. 随机序列 长度为N 的随机序列 基本数学函数参考教材P69页以及随后的使用说明。 注意使用行向量，特别是冒号运算符。 举例，长度为N 的实指数序列在MATLAB 中实现： n a x N n .^1 :0=-= 1. 单位采样 长度为N 的单位采样序列u(n)可以通过下面的MATLAB 命令获得：

u=[1 )1,1(-N zeros ]； 延迟M 个采样点的长度为N 的单位采样序列ud(n)(M

直接序列扩频通信

ＭＡＴＬＡＢ仿真直接序列扩频通信 1.摘要 直接序列扩频通信系统（DS-CDMA）因其抗干扰性强、隐蔽性好、易于实现码分多址（CDMA）、抗多径干扰、直扩通信速率高等众多优点，而被广泛应用于许多领域中。针对频通信广泛的应用，本文用MATLAB工具箱中的SIMULINK通信仿真模块和MATLAB函数对直接序列扩频通信系统进行了分析和仿真，使其更加形象和具体。 关键字：扩频通信m序列gold正交序列matlab仿真 2.引言 直接序列扩频(DSSS— Direct Sequence Spread Spectrum)技术是当今人们所熟知的扩频技术之一。这种技术是将要发送的信息用伪随机码（PN码）扩展到一个很宽的频带上去，在接收端，用与发端扩展用的相同的伪随机码对接收到的扩频信号进行相关处理，恢复出发送的信息。 它是二战期间开发的，最初的用途是为军事通信提供安全保障, 是美军重要的无线保密通信技术。这种技术使敌人很难探测到信号。即便探测到信号，如果不知道正确的编码，也不可能将噪声信号重新汇编成原始的信号。有关扩频通信技术的观点是在1941年由好莱坞女演员Hedy Lamarr 和钢琴家George Antheil提出的。基于对鱼雷控制的安全无线通信的思路，他们申请了美国专利#2.292.387。不幸的是，当时该技术并没有引起美国军方的重视，直到十九世纪八十年代才引起关注，将它用于敌对环境中的无线通信系统。 直序扩频解决了短距离数据收发信机、如：卫星定位系统(GPS)、3G移动通信系统、WLAN (IEEE802.11a, IEEE802.11b, IEE802.11g)和蓝牙技术等应用的关键问题。扩频技术也为提高无线电频率的利用率(无线电频谱是有限的因此也是一种昂贵的资源)提供帮助。 3.直接序列扩频DS-SS是直接用具有高码率的扩频码序列在发送端去扩展信 号的频谱。而在收端，用相同的扩频码序列去进行解扩，把展宽的扩频信号还原成原始的信息。

时间序列ARIMA模型的SAS程序编写

goptions vsize=7cm hsize=10cm; data b; format time monyy5.; input monyy7. asr; dif=dif(asr) ; keep time asr dif; cards; Jan1999 50 Feb1999 54.5 Mar1999 51 Apr1999 49 May1999 50 Jun1999 52 Jul1999 49 Aug1999 49 Sep1999 55 Oct1999 58 Nov1999 60 Dec1999 67.6 Jan2000 62 Feb2000 58.4 Mar2000 55 Apr2000 52.7 May2000 54.4 Jun2000 55.9 Jul2000 53.6 Aug2000 53.4 Sep2000 58.7 Oct2000 62.8 Nov2000 64.2 Dec2000 73.9 Jan2001 66.9 Feb2001 61.7 Mar2001 58.5 Apr2001 56.3 May2001 60.1 Jun2001 60.3 Jul2001 58 Aug2001 58.5 Sep2001 64.3 Oct2001 68.5 Nov2001 70.6 Dec2001 79.2 Jan2002 72.4

Feb2002 67.3 Mar2002 62.9 Apr2002 60.7 May2002 65.9 Jun2002 65.8 Jul2002 62.9 Aug2002 63.6 Sep2002 70.5 Oct2002 76 Nov2002 79 Dec2002 85.1 Jan2003 79.9 Feb2003 73.5 Mar2003 69.5 Apr2003 64.8 May2003 67.6 Jun2003 73.4 Jul2003 70.2 Aug2003 71.6 Sep2003 79.3 Oct2003 85.5 Nov2003 88.5 Dec2003 98.4 Jan2004 90.8 Feb2004 81.8 Mar2004 78.8 Apr2004 75 May2004 81 Jun2004 83.9 Jul2004 80.1 Aug2004 81.1 Sep2004 89.7 Oct2004 98.7 Nov2004 101.7 Dec2004 116.3 Jan2005 103.7 Feb2005 94.2 Mar2005 89.1 Apr2005 86.2 May2005 91.9 Jun2005 98.6 Jul2005 92.2 Aug2005 96.1 Sep2005 103.5

基于matlab的直接序列扩频通信系统仿真

基于MATLAB的直接序列扩频通信系统仿真 1.实验原理：直接序列扩频（DSSS）是直接利用具有高码率的扩频码系列采用各种调 制方式在发端与扩展信号的频谱，而在收端，用相同的扩频码序去进行解扩，把扩展宽的扩频信号还原成原始的信息。它是一种数字调制方法，具体说，就是将信源与一定的PN码（伪噪声码）进行摸二加。例如说在发射端将"1"用11000100110，而将"0"用00110010110去代替，这个过程就实现了扩频，而在接收机处只要把收到的序列是11000100110就恢复成"1"是00110010110就恢复成"0"，这就是解扩。这样信源速率就被提高了11倍，同时也使处理增益达到10DB以上，从而有效地提高了整机倍噪比。 1.1 直扩系统模型 直接序列扩频系统是将要发送的信息用伪随机码(PN码)扩展到一个很宽的频带上去，在接收端用与发送端相同的伪随机码对接收到的扩频信号进行相关处理，恢复出发送的信号。对干扰信号而言，与伪随机码不相关，在接收端被扩展，使落入信号通频带的干扰信号功率大大降低，从而提高了相关的输出信噪比，达到了抗干扰的目的。直扩系统一般采用频率调制或相位调制的方式来进行数据调制，在码分多址通信中，其调制多采用BPSK、DPSK、QPSK、MPSK等方式，本实验中采取BPSK方式。 直扩系统的组成如图1所示，与信源输出的信号a(t)是码元持续时间为Ta的信息流，伪随机码产生器产生伪随机码c(t)，每个伪随机码的码元宽度为Tc (Tc<

常用信号的MATLAB表示

5 常用信号的MATLAB表示5.1单位冲激函数、单位冲激序列 示例7： t = -5:0.01:5; y = (t==0); subplot(121); plot(t, y, 'r'); n = -5:5; x = (n==0); subplot(122); stem(n, x); 图5 运行结果如图5所示。

程序说明： (1)由n = -5:5得到一个1×11数组n；而在x = (n==0)中，n==0是一个向量运算，即向量n中的每一个元素与0比较是否相等，其比较结果0或1放在x中。这样得到的向量x也是1×11数组，且正好就是单位冲激序列。 (2)在MATLAB中，任何向量x的下标是从1开始的，不能取零或负值，而x(n)中的时间变量n则不此受限制。因此向量x的下标与时间变量n是两个概念，如本例中向量x(n)的下标是从1到11，而时间变量n是从-5到5。所以必须用一个与向量x等长的定位时间变量n，以及向量x，才能完整地表示序列x(n)。在信号的表示和运算中，这一点请务必注意；只有当序列x(n)的时间变量正好是从1开始时，才能省去时间变量n，因为此时向量的下标与时间变量相同。 (3)单位冲激函数的实现方法实际上与单位冲激序列是完全相同的，都是用序列表示。只不过表示连续时间信号的序列中两相邻元素所对应的时间间隔更小，如本例中t的间隔为0.01，而表示离散时间信号的序列中两相邻元素所对应的时间间隔一般为1。 由于单位冲激序列在信号与系统中经常使用，我们专门编制一个函数文件delta.m，在后面的实验部分直接调用该函数即可产生需要的波形。 % delta.m function [x, n] = delta(n1,n2,k) % 产生冲激序列δ(n-k),其中n1<=n<=n 2, n1<=k<=n2

基于MATLAB的太阳黑子时间序列与仿真

2012.3 26 基于MATLAB 的太阳黑子时间序列 分析与仿真 周园 肖洪祥 董俊飞 桂林理工大学信息科学与工程学院 广西 541004 摘要：本文研究了时间序列的分析方法，具体分析了基于最大Lyapunov 指数的方法在太阳黑子时间序列分析中的应用。介绍利用MATLAB 对太阳黑子时间序列进行分析与仿真的方法，并给出相关的流程、程序和相应的仿真结果。最终证明太阳黑子时间序列是一个混沌时间序列。 关键词：混沌时间序列；最大Lyapunov 指数；太阳黑子数；仿真 0 引言 在非线性系统中，初始条件的微小变化，往往会导致结果以指数级的大小发生分离，这时我们称这个系统存在混沌。时间序列是非线性动力系统的一种模型。如果时间序列对初始条件敏感，采用传统线性时间序列分析方法将很难予以分析，因此传统时间序列预测模型对混沌时间序列的拟合和预测准确度都很差。经过混沌学的发展，可以使用序列本身的规律对其进行预测。Lyapunov 指数法即是其中之一。通过最大Lyapunov 指数的数值，可以判断一个时间序列是否是混沌时间序列，亦即该非线性系统中是否存在着混沌。本文对太阳黑子序列进行分析，证明其是一个混沌时间序列。 1 基于Lyapunov 指数的时间序列分析方法 对时间序列进行分析，首先必须进行相空间重构。根据有限的数据重构吸引子以研究系统动力行为的方法即是相空间重构。主要思想为：系统中每个分量的演化皆是由与之联系的其他分量所决定的，相关分量的信息隐含在任意其他分量的变化过程中，即是运用系统的任何一个观察量可以重构出整个系统的模型。 设时间序列为{}t x ，其中1,2,...,t N =。重构相空间m R 的元素组为： (1)(,,)(,,...,),T=1,2,3,...,T T T T m X m N X X X p τττ++-= (1) 其中，N 为重构相空间维数；τ为延迟时间间隔数，且为正整数；(1)p N m τ=--为时间序列嵌入相空间的向量数，N 为时间序列的数据点数。 由Tokens 定理，在理论条件下可任选τ。但在现实条件下时间序列都是有限长且有噪声的。因而在重构相空间时， τ的选取至关重要。目前所采用的方法大多是通过经验来选 择τ， 从而使得T X 和T X τ+相互独立并不完全相关。 Lyapunov 指数是描述奇异吸引子性质的数据量。在m 维离散系统中存在m 个Lyapunov 指数，即Lyapunov 指数族。正的Lyapunov 指数意为在此维度方向，系统以指数级速度分离。1983年，G.Grebogi 证明了若最大Lyapunov 指数 max 0λ>，则系统一定存在着混沌。因此要判断一个时间序 列是否为混沌时间序列，必须求出其最大Lyapunov 指数。为了保证领域点沿着不同的轨道运动，最近邻域点间必须有分离间隔。此处取分离间隔为/w T t =?，其中T 为用FFT 计算出的序列平均周期；t ?为序列的采样周期。 2 计算机仿真步骤 仿真步骤如图1所示。输入太阳黑子年平均序列，通过 FFT 算法计算得到其平均周期T 。计算分离间隔作为时间窗 ωτ。由公式得到嵌入维数m 。运用所得的参数使用Wolf 法 算出最大Lyapunov 指数。进而判断该序列是否是混沌序列。

白噪声的测试MATLAB程序

白噪声的测试MATLAB程序 学术篇 2009-11-13 22:18:03 阅读232 评论0 字号：大中小订阅 clear; clc; %生成各种分布的随机数 x1=unifrnd(-1,1,1,1024);%生成长度为1024的均匀分布 x2=normrnd(0,1,1,1024);%生成长度为1024的正态分布 x3=exprnd(1,1,1024);%生成长度为1024的指数分布均值为零 x4=raylrnd(1,1,1024);%生成长度为1024的瑞利分布 x5=chi2rnd(1,1,1024);%生成长度为1024的kaifang分布%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %求各种分布的均值 m1=mean(x1),m2=mean(x2),m3=mean(x3),m4=mean(x4),m5=mean(x5) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %求各种分布的方差 v1=var(x1),v2=var(x2),v3=var(x3),v4=var(x4),v5=var(x5) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %求各种分布的自相关函数 figure(1);title('自相关函数图'); cor1=xcorr(x1);cor2=xcorr(x2);cor3=xcorr(x3);cor4=xcorr(x4);cor5=xcorr(x5); subplot(3,2,1),plot(1:2047,cor1);title('均匀分布自相关函数图'); subplot(3,2,2),plot(1:2047,cor2);title('正态分布'); subplot(3,2,3),plot(1:2047,cor3);title('指数分布'); subplot(3,2,4),plot(1:2047,cor4);title('瑞利分布'); subplot(3,2,5),plot(1:2047,cor5);title('K方分布'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %求各种分布的概率密度函数 y1=unifpdf(x1,-1,1); y2=normpdf(x2,0,1); y3=exppdf(x3,1); y4=raylpdf(x4,1); y5=chi2pdf(x5,1); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %各种分布的频数直方图 figure(2); subplot(3,2,1),hist(x1);title('均匀分布频数直方图'); subplot(3,2,2),hist(x2,[-4:0.1:4]);title('正态分布'); subplot(3,2,3),hist(x3,[0:.1:20]);title('指数分布'); subplot(3,2,4),hist(x4,[0:0.1:4]);title('瑞利分布'); subplot(3,2,5),hist(x5,[0:0.1:10]);title('K方分布'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %各种分布的概率密度估计 figure(3);

实验1 常见离散信号的MATLAB产生和图形显示

实验1 常见离散信号的MATLAB 产生和图形显示 一、实验目的：加深对常用离散信号的理解 二、实验原理： 1.单位抽样序列：???=01)(n δ 00 ≠=n n 在MATLAB 中可以利用zeros()函数实现。 ; 1)1();,1(==x N zeros x 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n ?δ即：? ??=?01)(k n δ 0≠=n k n 2．单位阶越序列：???0 1 )(n u 00<≥n n 在MATLAB 中可以利用ones()函数实现。(1,)x ones N = 3．正弦序列：)/2sin()(?π+=Fs fn A n x 在MATLAB 中 ) /***2sin(*1 :0fai Fs n f pi A x N n +=?= 4．复正弦序列：n j e n x ?=)( 在MATLAB 中 ) **exp(1 :0n w j x N n =?= 5．指数序列：n a n x =)( 在MATLAB 中 n a x N n .^1:0=?= 三、实验内容： 1、编制程序产生上述5种信号（长度可输入确定），并绘出其图形。 2、讨论正弦序列、复指数序列的性质： （1）绘出信号()zn x n e =，当1126z j π=?+、1126z j π=+、112 z =、62πj z +=、6π j z =时 信号的实部和虚部图；当6 π j z =时信号的周期为多少？

（2）绘出信号() 1.5sin(2*0.1)x n n π=的频率是多少？周期是多少？产生一个数字频率为0.9的正弦序列，并显示该信号，说明其周期。 3、使用帮助功能学习square(方波),sawtooth(锯齿波)和sinc 函数，并绘图。 四、实验要求： 1、预先阅读MATLAB 基础； 2、讨论复指数序列的性质。

直接序列扩频通信系统的误码率仿真培训讲学

直接序列扩频通信系统的误码率仿真

直接序列扩频通信系统的误码率仿真 1．引言 扩展频谱通信系统是将基带信号的频谱扩展至很宽的频带上，然后再进行 传输的一种通信系统,即将待传送的信息数据用伪随机编码调制，实现频谱扩展后再传输，接收端则采用同样的编码进行解调及相关处理，恢复原始信息数据。 扩频通信的基础理论根据信息论中的shannon 公式 ） （N S B C /1log 2+= 式中，C 是系统的信道容量，B 是系统信道带宽,N 是噪声功率,S 为信号的功率,S/N 即为信噪比。 Shannon 公式表明了一个系统信道无误差的传输信息的能力与存在于信道中的信噪比以及用于传输信息的系统信道带宽之间的关系。该公式说明了两个极为重要的概念：一是在一定的信道容量条件下，可以用减少发送信号功率、增加带宽的方法来达到信道容量的要求；另一个是可以采用减少带宽而增加信号功率的方法来达到信道容量的要求。这也就说明了信道容量可以通过带宽与信噪比的互换来保持不变。在实际的工程应用中，改变信号的功率并不容易，相比较而言，扩展信号的带宽更容易操作，所以，要提高信道容量，采用增加信号的带宽比提高信号功率的方法要有效的多。 由于扩频通信系统可以在信号功率远低于噪声功率的环境中工作，因此扩 频通信系统具有抗干扰能力强，保密性强等优点，在现在通信领域内的应用越 来越广泛。 2．系统概述 本次仿真实验是以MATLAB 为仿真平台，信号是8位双极性二进制信号,由 1和-1组成。随后对产生的双极性信号进行时域抽样，得到基带信号s ，是一组1024位的信息码。伪随机序列由mgen 函数产生，共有1024个码元。对已得到的基带信号进行扩频调制，直接把基带信号S 与产生的伪随机序列相乘，得到扩频信号。然后对已作扩频处理的信号作BPSK 载波调制，得到发射信号。发射信号通过存在高斯白噪声的信道，到达接到端，接收端首先对信号进

利用SPSS和Matlab进行时间序列预测

§7.利用SPSS 和Matlab 进行时间序列预测 1.移动平均和滑动平均计算 例1：下表给出了某地区1990～2004年粮食产量数据（表1）。试分别用Matlab 和SPSS 软件，对该地区的粮食产量进行移动平均和和滑动平均计算。 表1 某地区1990～2004年粮食产量及其平滑结果 移动平均 滑动平均 年份 自然序号 粮食产量y （单位：104 t ）三点移动 五点移动 三点滑动 五点滑动 1990 1 3149.44 1991 2 3303.66 3154.47 1992 3 3010.30 3010.30 3141.19 3242.44 1993 4 3109.61 3154.47 3253.04 3263.32 1994 5 3639.21 3141.19 3334.21 3295.88 1995 6 3253.80 3253.04 3242.44 3453.17 3461.80 1996 7 3466.50 3334.21 3263.32 3520.07 3618.81 1997 8 3839.90 3453.17 3295.88 3733.69 3692.89 1998 9 3894.66 3520.07 3461.80 3914.72 3892.78 1999 10 4009.61 3733.69 3618.81 4052.51 4019.78 2000 11 4253.25 3914.72 3692.89 4121.45 4075.78 2001 12 4101.50 4052.51 3892.78 4158.21 4148.58 2002 13 4119.88 4121.45 4019.78 4160.01 4227.01 2003 14 4258.65 4158.21 4075.78 4260.11 2004 15 4401.79 4160.01 4148.58 利用spss 进行移动平均计算主要有以下步骤： （1） 在菜单中依次选择transform->Create time series…,在弹出的对话框中 的单击function 下面的下拉条，选择Prior moving average 方法，span 框中输入数值为3（表示进行三点滑动平均）。

直接序列扩频系统的Simulink仿真

直接序列扩频系统的Matlab/Simulink仿真 摘要：本文利用Matlab/Simulink对直接序列扩频系统进行了仿真，对其原理进 行了相关的说明。读者可以通过对本文的阅读对直接序列扩频的相关原理有一定的了解。 关键字：扩频通信直接序列扩频 一、仿真的意义 随着信息技术的发展，通信技术变得越来越复杂，技术更新的周期也越来越短。对于大部分学者，特别是我们学生来说，在学习通信技术时，若对每一个系统都要实体研究是不现实的。此时通信系统仿真对我们来说可以说是必不可少的。通过建立相应的通信系统的模型，对其进行仿真，可以使我们把琐碎的知识联系在一起，形成一个个通信系统的概念，可以让我们对各个知识点的原理有更加深刻的理解和掌握。 二、直接序列扩频的原理 扩频通信，即扩展频谱通信(Spread Spectrum Communication)是将待传送的信息数据用伪随机编码(扩频序列：Spread Sequence)调制，实现频谱扩展后再传输而接收端则采用相同的编码进行解调及相关处理，恢复原始信息数据。扩频通信具有抗干扰能力强、抗噪声、保密性强、功率谱密度低，具有隐蔽性和较低的截获概率、可多址复用和任意选址、高精度测量等优点。 根据扩展频谱方式的不同，可以将扩频通信系统分为直接序列扩频(Direct Sequence Spread Spectrum)工作方式，简称直扩(DS)方式；跳变频率(Frequency Hopping)工作方式，简称跳频(FH)方式；跳变时间(Time Hopping)工作方式，简称跳时(TH)方式；宽带线性调频(Chirp Modulation)工作方式，简称Chirp方式和各种混合方式。 直接序列(DS-Direct Scquency)扩频，就是直接用具有高码率的扩频码序列在发端去扩展信号的频谱，而在收端，用相同的扩频码序列去进行解扩，把展宽的扩频信号还原成原始的信息。直接序列扩频是扩频通信系统最基本的工作方式。 假设信源序列对应的双极性波形为a(t)，其电平取值为±1 ，码元速率为Rabps，码元宽度为Ta=1/Ra/秒。扩频所使用的伪随机序列c(t)也是电平取值为±1 的双极性波形，伪随机序列（PN序列）的码元也称为码片（chip），码片速率设为Rcchip/s，对应的码片宽度就是Tc=1/Rc/秒。对于双极性波形而言，扩频过程等价于数据流a(t)与伪随机序列c(t)相乘的过程，扩频输出序列设为d(t)，也是取值为±1 的双极性波形，其速率等于码片速率。扩频序列经过调制后得到调制输出信号s(t)送入信道。对于BPSK调制，发送的信号就相当于是数据流与伪随机序列相乘后再乘于一个高频的余弦信号。在接收端，接收到的信号中有包含了有用信号s(t)及各种干扰J(t)和噪声n(t)。由于接收端采用相关解扩，即将s(t)J(t)n(t)和本地PN序列c(t)相乘，只有有用信号的频谱能够被还原为窄带信号，其他的噪声和干扰的频谱只会被展宽，当信号通过窄带滤波器后只有一小部分被展宽了的频谱会混进有用信号中，由此大大增强了其抗干扰的能力。 三、仿真的系统与结果 此处是对直接序列扩频通信系统的仿真。假设该系统以BPSK方式调制，数

白噪声产生程序

第二章的白噪声产生程序 例2.2 用乘同余法产生(见光盘FLch2bzsheg2.m) ①编程如下： A=6; x0=1; M=255; f=2; N=100；%初始化； x0=1; M=255; for k=1: N %乘同余法递推100次； x2=A*x0; %分别用x2和x0表示x i+1和x i-1； x1=mod (x2,M); %取x2存储器的数除以M的余数放x1（x i）中； v1=x1/256; %将x1存储器中的数除以256得到小于1的随机数放v1中； ）减去0.5再乘以存储器f中的系数，存放在v(:,k)=(v1-0.5 )*f; %将v1中的数（ i 矩阵存储器v的第k列中，v(:,k)表示行不变、列随递推循环 次数变化； x0=x1; % x i-1= x i； v0=v1; end %递推100次结束； v2=v %该语句后无‘；’，实现矩阵存储器v中随机数放在v2中，且 可直接显示在MATLAB的window中； k1=k; %grapher %以下是绘图程序； k=1:k1; plot(k,v,k,v,'r'); xlabel('k'), ylabel('v');tktle(' (-1,+1)均匀分布的白噪声') ②程序运行结果如图2.6所示。 图2.6 采用MA TLAB产生的(-1,+1)均匀分布的白噪声序列 ③产生的(-1，1)均匀分布的白噪声序列 在程序运行结束后，产生的(-1，1)均匀分布的白噪声序列，直接从MATLAB的window 界面中copy出来如下（v2中每行存6个随机数）：

v2 = -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 *另外，书中图2.3白噪声的产生如下： 显然，只要在例2.2程序的初始化部分中给N=300，f=6，运行程序就可以得到如图2.3所示的（-3，3）的白噪声过程. ①编程如下： A=6; x0=1; M=255; f=6; N=300；%初始化； x0=1; M=255; for k=1: N %乘同余法递推100次； x2=A*x0; %分别用x2和x0表示x i+1和x i-1； x1=mod (x2,M); %取x2存储器的数除以M的余数放x1（x i）中； v1=x1/256; %将x1存储器中的数除以256得到小于1的随机数放v1中； ）减去0.5再乘以存储器f中的系数，存放在v(:,k)=(v1-0.5 )*f; %将v1中的数（ i 矩阵存储器v的第k列中，v(:,k)表示行不变、列随递推循环 次数变化； x0=x1; % x i-1= x i； v0=v1; end %递推100次结束； v2=v %该语句后无‘；’，实现矩阵存储器v中随机数放在v2中，且 可直接显示在MATLAB的window中； k1=k; %grapher %以下是绘图程序； k=1:k1; plot(k,v,k,v,'r'); xlabel('k'), ylabel('v');tktle(' (-1,+1)均匀分布的白噪声')

直接序列扩频通信系统开题报告

哈尔滨工业大学华德应用技术学院毕业设计（论文）开题报告 题目：直接序列扩频通信系统的设计与仿真实现 系（部）应用电子与通信技术 专业通信工程 学生薛光宇 学号24 班号0992222 指导教师周凯 开题报告日期2012.10,22 哈工大华德学院

说明 一、开题报告应包括下列主要内容： 1．通过学生对文献论述和方案论证，判断是否已充分理解毕业设计（论文）的内容和要求 2．进度计划是否切实可行； 3．是否具备毕业设计所要求的基础条件。 4．预计研究过程中可能遇到的困难和问题，以及解决的措施； 5．主要参考文献。 二、如学生首次开题报告未通过，需在一周内再进行一次。 三、开题报告由指导教师填写意见、签字后，统一交所在系（部）保存，以备检查。指导教师评语： 指导教师签字：检查日期：

一、课题题目和课题研究现状 课题题目：直接序列扩频通信系统的设计与仿真实现。 研究现状：目前扩频技术中研究最多的对象是CDMA技术，其中又以码捕获技术和多用户检测（MUD）技术代表了目前扩频技术研究的现状。 1．码捕获 同步的实现是直扩系统中一个关键问题。只有在接收机将本地产生的伪码和接收信号中调制信息的伪码实现同步以后，才有可能实现直序扩频通信的各种优点。同步过程分为两步来实现:首先是捕获阶段，实现对接收信号中伪码的粗跟踪；然后是跟踪阶段，实现对伪码的精确跟踪。目前的研究主要集中在码捕获过程。 2, 多用户检测 CDMA系统容量受到来自其他用户的多址干扰的限制，多用户检测能够利用这些多址干扰来改善接收机的性能，因此是一种提高系统容量的有效方法。传统的CDMA 接收机是由一系列单用户检测器组成，每个检测器都是与特定扩频码对应的相关器，它并没有考虑多址干扰的结构，而是把来自其它用户的干扰当成加性噪声，因此当用户数量增加时，其性能急剧下降。通过对所有用户的联合译码可以极大地改善CDMA系统的性能。但是最优的多用户接收机，其复杂度随用户数量成指数增长，因此在实际通信系统中几乎不可能实现。这样寻找在性能和复杂度之间折中的次最优多用户检测器成为研究的热点 二、目的及意义 通过对该课题的研究，了解科研学术论文的撰写流程，并且将自己所学的理论知识运用到论文中，全面多角度的分析该领域的发展现状，同时提高自己的思维能力，对搜集的数据进行恰当处理和准确分析，对大学本科四年学习成果进行有效的检验，并且进一步提高自学能力和自主进行科学研究的水平。 三、课题的基本内容 所谓直接序列扩频（DS），就是直接用具有高速率的扩频码序列在发送端去扩展信号的频谱。而接收端，用相同的扩频码序列进行解扩，把展宽的扩频信号还原成原始信息。

时间序列MATLAB程序

时间序列MATLAB程序

时间序列 移动平均法 clc,clear y=[533.8 574.6 606.9 649.8 705.1 772.0 816.4 892.7 963.9 1015.1 1102.7]; m=length(y); n=[4,5];%n为移动平均的项数 for i=1:length(n) %由于n的取值不同，下面使用了细胞数组 for j=1:m-n(i)+1

yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i); end y12(i)=yhat{i}(end);%提出第12月份的预测值 s(i)=sqrt(mean((y(n(i)+1:end)-yhat{i}(1:end-1)).^2));%求预测的标准误差end y12, s %分别显示两种方法的预测值和预测的标准误差 指数平滑 一次指数平滑 程序： clc,clear yt=load('dianqi.txt'); %读取dianqi.txt的数据， n=length(yt); %求yt的长度 alpha=[0.2 0.5 0.8]; %输入a的值 m=length(alpha) yhat(1,[1:m])=(yt(1)+yt(2))/2; %求第一个预测值索引

for i=2:n yhat(i,:)=alpha*yt(i-1)+(1-alpha).*yhat(i-1,:); end yhat %求预测值 err=sqrt(mean((repmat(yt,1,m)-yhat).^2)) %求预测的标准误差 xlswrite('dianqi.xls',yhat) %把预测数据写到Excel文件，准备在word表格中使用yhat1988=alpha*yt(n)+(1-alpha).*yhat(n,:) %求1988的预测值 二次指数平滑

matlab 正弦波 高斯白噪声 均匀白噪声 功率谱密度 自相关函数

现代通信原理作业一 姓名：张英伟学号：8036 班级：13级理工部3班 利用matlab完成： ●产生正弦波信号、均匀白噪声以及高斯白噪声并分别将两种噪声叠加到正弦 波信号上，绘出波形。 ●分别求取均匀白噪声序列和高斯白噪声序列的自相关及功率谱密度，绘出波 形。 一、白噪声区别及产生方法 1、定义： 均匀白噪声：噪声的幅度分布服从均匀分布，功率谱密度在整个频域内均匀分布的噪声。 高斯白噪声：噪声的幅度分布服从正态分布，功率谱密度在整个频域内均匀分布的噪声。 2、matlab仿真函数： rand函数默认产生是区间在[0,1]的随机数，这里需要利用公式： z2=a+(b-(a))*rand(m,n)............（公式1） randn函数默认产生均值是0、方差是1的随机序列，所以可以用其来产生均值为0、方差为1的正态分布白噪声，即N(0，12)。利用公式： z1=a+b*randn(1,n).................（公式2） 可以产生均值为a，方差为b2 高斯白噪声，即N(a，b2)。 二、自相关函数与功率谱密度之间的关系 1、功率谱密度：每单位频率波携带的功率，这被称为信号的功率谱密度。 2、自相关函数：描述随机信号X(t)在任意两个不同时刻t1，t2的取值之间的相关程度。 3、维纳-辛钦定理： 由于平均值不为零的信号不是平方可积的，所以在这种情况下就没有傅里叶变换。幸运的是维纳-辛钦定理提供了一个简单的替换方法，如果信号可以看作是平稳随机过程，那么功率谱密度就是信号自相关函数的傅里叶变换。 4、平稳随机过程：是在固定时间和位置的概率分布与所有时间和位置的概率分布相同的随机过程。（就是指得仅一个随机过程，中途没有变成另外一个统计特性的随机过程）