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

1、基于MATLAB的直接序列扩频通信系统仿真1.实验原理：直接序列扩频（DSSS）是直接利用具有高码率的扩频码系列采用各种调制方式在发端与扩展信号的频谱，而在收端，用相同的扩频码序去进行解扩，把扩展宽的扩频信号还原成原始的信息。它是一种数字调制方法，具体说，就是将信源与一定的PN码（伪噪声码）进行摸二加。例如说在发射端将1用11000100110，而将0用00110010110去代替，这个过程就实现了扩频，而在接收机处只要把收到的序列是11000100110就恢复成1是00110010110就恢复成0，这就是解扩。这样信源速率就被提高了11倍，同时也使处理增益达到10DB以上，从而有效地提高了整

2、机倍噪比。1.1 直扩系统模型直接序列扩频系统是将要发送的信息用伪随机码(PN码)扩展到一个很宽的频带上去，在接收端用与发送端相同的伪随机码对接收到的扩频信号进行相关处理，恢复出发送的信号。对干扰信号而言，与伪随机码不相关，在接收端被扩展，使落入信号通频带内的干扰信号功率大大降低，从而提高了相关的输出信噪比，达到了抗干扰的目的。直扩系统一般采用频率调制或相位调制的方式来进行数据调制，在码分多址通信中，其调制多采用BPSK、DPSK、QPSK、MPSK等方式，本实验中采取BPSK方式。直扩系统的组成如图1所示，与信源输出的信号a(t)是码元持续时间为Ta的信息流，伪随机码产生器产生伪随机码c(t

3、)，每个伪随机码的码元宽度为Tc (Tc=0 xo(i)=1; else xo(i)=-1; endend figure(8)plot(t0,xo)xlabel(t(ms);ylabel(xo(t)(V);title(解调后输出的信号);axis(0,400,-1.5,1.5);grid on for a=1:15 r= awgn(s,SNR(a),measured); rs=r.*cos(2*fc*pi*td); jk_code=rs.*k_code; f,rf=T2F(td,jk_code); tl,l=lpf(f,rf,fc/10000); liremod=cos(2*fs*pi*tq)

4、*sqrt(2).*l; fi,lif=T2F(tq,liremod);tji,iremod=lpf(fi,lif,fs/10000); iremodo=sign(iremod); lqremod=(-1)*sin(2*fs*pi*tq)*sqrt(2).*l; fq,lqf=T2F(tq,lqremod); tjq,qremod=lpf(fq,lqf,fs/10000); qremodo=sign(qremod); zl=qremodo,iremodo; z2=0;xo=zeros(size(x);for i=1:100 z2=0 for j=1+(i-1)*400:i*400 z2=z2+

5、zl(j); end if z2=0 xo(i)=1; else xo(i)=-1; endend z=0;for i=1:100 if(x(i)=xo(i) z=z+1; endendz1(a)=z/100endfigure(8)plot(SNR,z1)xlabel(SNR);ylabel(z1);title(误码率);axis(0,15,10(-4),1);grid on 程序运行结果：3 仿真参数：信源速率2kbps，扩频因子：64，调制方式：BPSK，信道：高斯信道4总结与感想：扩频通信以其较强的抗干扰、抗衰落、抗多径性能而成为第三代通信的核心技术，本文阐述了扩频通信的理论基础和实现方

6、法，利用MATLAB，详细讲述了各模块的设计，并给出了仿真建模中需注意的问题。在给定仿真条件下，运行了仿真系统，验证了所建仿真模型的正确性。通过仿真研究了扩频增益和输出端信噪比的关系，结果表明，在相同误码率下，增大扩频增益，可以提高系统输出端的信噪比，从而提高系统的抗干扰能力。参考文献1 邵保华. 软件无线电在扩频通信发射中的应用研究D. 哈尔滨工程大学，2004年.2 张蕾，郑实勤. 基于MATLAB的直接序列扩频通信系统性能仿真分析研究J. 电气传动自动化，2007年第3期.3 查光明，熊贤祚. 扩频通信M. 西安电子科技大学出版社，1990年.附录资料：不需要的可以自行删除Pascal/

7、C/C+语句对比（补充版）一、Hello world 先看三种语言的样例：Pascalbegin writeln(Hello world);end.C#include int main() printf(Hello world!n); return 0;C+#include using namespace std;int main()cout Hello world! endl; return 0; 从这三个程序可以看到一些最基本的东西。在Pascal中的begin和end，在C/C+里就是；Pascal主程序没有返回值，而C/C+返回0（好像在C中可以为NULL）。在C/C+中，main函数

8、以前的是头文件，样例中C为stdio.h，C+除了iostream还有第二行的using namespace std，这个是打开命名空间的，NOIP不会考这个，可以不管，只要知道就行了。 此外说明 注释单行用/，段落的话Pascal为，C/C+为/* */。* 常用头文件（模板）#include #include #include #include #include #include using namespace std;int main() system(“pause”);return 0;二、数据类型及定义 这里只列出常用的类型。1、整型PascalC/C+范围shortint-128

9、127integershort-32768 32767longintInt -2147483648 2147483647int64long long-9223372036854775808 9223372036854775807byte-0 255wordunsigned short0 65535longwordunsigned int0 4294967295qwordunsigned long long0 18446744073709551615 * 当对long long 变量赋值时，后要加LLLong long x=6327844632743269843LL* 如果位移 x2LL* Li

10、nux: printf(“%lldn”,x);* Windows: printf(“%I64dn”,x);2、实型PascalC/C+范围realfloat2.9E-39 1.7E38single-1.5E-45 3.4E38doubledouble5.0E-324 1.7E3083、字符即字符串 字符在三种语言中都为char，C里没有字符串，只有用字符数组来代替字符串，Pascal和C+均为string。Pascal中字符串长度有限制，为255，C+则没有。 字符串和字符在Pascal中均用单引号注明，在C/C+中字符用单引号，字符串用双引号。4、布尔类型 Pascal 中为 boolean

11、，C/C+ 为 bool。值均为True 或 False。C/C+中除0外bool都为真。5、定义 常量的定义均为 const，只是在C/C+中必须要注明常量的类型。在C/C+中还可以用宏来定义常量，此时不注明类型。PascalC/C+const a = 60; b = -a + 30; d = ;const int a = 60;const int b = - a + 30;const string d = “”;define MAXN 501 /这个是宏 * 宏定义其实就是直接在程序相应的位置替换： #define randomize srand(unsigned time(NULL) #

12、define wait for(int w=0;w a;cout a;cout a endl;特别说明C+中cin一个字符的话会自动跳过空格和回车，Pascal和C则会读入空格和回车。在Pascal中writeln(a:n:m) 表示在n个字符宽的输出域上输出a保留m位小数。例如：pascal write(a:6) c/c+ printf(“%6d”,a) Pascal write(a:6:2) c/c+ printf(“%6.2f”,a) C+ 如果用 cout ? (繁琐！) 需要加头文件 #inlude cout setprecision(2)a; /作用永久 cout setw(6)a

13、; /作用临时 以下三个进制设定都是永久作用： cout deca; 相当 printf(“%d”,a); /十进制 cout hexa; 相当 printf(“%X”,a); /十六进制 cout octa; 相当 printf(“%o”,a); /八进制例如：cout 12hex12oct1212endl;输出：12c1414 C 的输入输出里面的字符串中%表示变量，%后面的字目表示变量类型。下面是类型表：%hd1个short型整数%d1个int型整数%u1个unsigned int型整数%I64d1个long long型整数%c1个字符%s1个C字符串%f1个float型实数%lf1个d

14、ouble型实数%10.4f输出1个总宽度为10,保留4位小数的实数 文件输入输出：Pascalassign(input, test.in);assign(output, test.out);reset(input);rewrite(output);read(a, b);writeln(a, b);close(input);close(output);CFILE *fin = fopen(“test.in”, “r”);FILE *fout = fopen(“test.out”, “w”);fscanf(fin, “%d%d”, &a, &b);fprintf(fout, “%d%d”, a,

15、 b);fclose(fin); fclose(fout);C+#include using namespace std;ifstream fin(“test.in”);ofstream fout(“test.out”);fin a b;fout a b endl;fin.close(); fout.close();因为C+的读入较慢，个人建议C+的话使用C的输入方式。当然也有人用C的读入，C+的输出的，这种方式我们称之为城乡结合。*中国计算机学会竞赛须知发布的C读写程序：（C+ 也能用，cin,cout,scanf,printf 可混用）#include int main() int a,b

16、; freopen(“sum.in”,”r”,stdin);freopen(“sum.out”,”w”,stdout); scanf(“%d%d”,&a,&b); printf(“%dn”,a+b); return 0; 或者：freopen(“sum.in”,”r”,stdin);freopen(“sum.out”,”w”,stdout);ios:sync_with_stdio(false); 取消同步，cin,cout的速度就不慢了！ cinab;couta+bendl; return 0;以下扩充c/c+混用是可行的：#include #include using namespace s

17、td;int main() int a,b,c,d; freopen(sum.in,r,stdin); freopen(sum.out,w,stdout); scanf(%d%d,&a,&b); cincd; printf(%dn,a+b); couta+b+c+dsn).Cwhile(scanf(%s%d,s,&n)!=EOF).四、赋值语句及运算符号 一一对应的关系PascalC/C+赋值运算赋值:=基本运算加+减-乘*除（实数）/ (double)除法取整div(int) / (int)取余mod%比较等于=不等于!=大于大于等于=小于小于等于=逻辑且and&或or|非not!位运算左移

18、(*2)shl且and&或or|非not异或xor其他增一inc(x)x+减一dec(x)x- 在C/C+中对某个变量自身进行运算可以简写为 变量名 运算符号= 改变量 如 x += 8 就表示 x = x + 8， 即 inc(x, 8)。 在 C/C+里还存在一种三目运算 变量名 = 条件 ? 值A : 值B 如 x = x 0 ? x : -x; /表示若x 0 则取 x， 否则取 x， 同 if x 0 then x := x else x := -x;五、条件语句1、if C/C+中if 语句的条件必须要用括号括起来，后面不使用then。PascalC/C+if a b then f

19、lag := true else flag := false;if (a b) flag = true;else flag = false;2、多种分支 C/C+中为switch，Pascal为case：PascalC/C+case x of 1: inc(x); 2: dec(x); else x := x * x;end;switch (x) case 1: x +; break; case 2: x -; break; default: x *= x; 切记C/C+中一定要写break，后果你可以去掉break，运行看看就知道了。六、循环语句1、forPascalC/C+for 变量名

20、:= 初始值 to(downto) 终止值 dofor (变量名=初始值;条件;改变方式)for i := 5 to 10 do dec(a);/终止值大于初始值用 tofor i := 5 downto 1 do dec(a);/终止值小于于初始值用 downtofor (i = 5; i = 1; i-) a-;/*只要i 满足条件就会一直循环。C/C+中i是实数、指针都可以*/C/C+中for的特殊用法：/变量为实数for (double i = 1; i 符号为间接引用，后面会提到。for (type1 *p = head - next; p; p = p - next) printf

21、(“%d”, p - k);2、whilePascalC/C+while 条件 dowhile (条件)while i 0 do dec(i);while (i != 0) i-;/也可写作 while (i) i-;/在C/C+中非0即为真。3、repeat-until & do-whilePascalC/C+repeat 语句 until 结束条件;do while (运行条件)repeat int(i) until i 100;do i+; while (i = 100);七、数组 Pascal中数组的下标可以随意定义，而C/C+下标始终为从0开始到(数组大小1)。PascalC/C+定

22、义a : array 1.100 of integer;b :array 1.10,1.10 of int64;int a100;int b1010;含义a 为大小为100的integer数组，合法下标为1到100b 为大小为10*10的int64数组，合法下标为1,1到10,10a 为大小为100的int数组，合法下标为0到99b 为大小为10*10的int数组，合法下标为0,0到9,9;使用inc(a21);b2,2:=b1,1+b1,2+b2,1;a21+;b11=b01+b00+b10; 数组清零PascalC/C+Fillchar(a, sizeof(a), 0);memset(a,

23、 0, sizeof(a);/头文件包含 string.h*如果要填最大： memset(a,127,sizeof(a) （但达不到 INT_MAX） 如果要填最小： memset(a,128,sizeof(a) (但达不到 INT_MIN） 如果填0： memset(a,0,sizeof(a) 如果填-1： memset(a,-1,sizeof(a)八、字符串 C风格的字符串就是字符数组。 C+和Pascal的字符串使用基本相同，只是C+中字符串下标以0开始，Pascal以1开始。字符串处理很多这里不一一列举，只写最常用的几个。PascalC (包含)定义用：char sC+（包含）定义用：

24、string s输入输出Readln(s);Writeln(s);Scanf(“%s”,s);Printf(“%sn”,s);注：不能输入输出c+的字符串Cins;Couts = s 的区别： getline(cin,s)cins一次性整行读入，直至行末尾。只读入一个“单词”，遇空格和行末停止。例如输入；How are you?s=” How are you?”读入整串含空格例如输入；How are you?s=”How”如果三个都读：cins1s2s3*C+ 数字与数值之间的转换：#include #include #include /必须加入using namespace std;int

25、main() string text = 152; int number; stringstream ss;ss number; /string - int coutnumber+100endl; ss string string str = ss.str(); return 0;九、过程和函数1、过程 在C/C+中没有过程，但可以把返回值为“空”的函数理解为过程。PascalC/C+无参过程procedure 过程名;说明部分begin 语句部分 end;/说明部分、begin、end语句部分统称为过程体void 函数名(); 主体部分; return ;带参过程procedure 过程名(

26、形参表)过程体void 函数名(形参表)过程体 值传和址传：当一个参数是值传时，形参在子过程中相当于一个局部变量，对它的改变不影响实在的参数值。址传则会影响。下例中a为值传，b为址传。初始a = 5,b = 5，运行后a = 5,b = 10；PascalC/C+var a, b:integer;procedure doit(a:integer; var b:integer);begin b := a + b; a := a + b;end;begina := 5;b := 5;doit(a, b);writeln(a, , b);end.void doit(int a, int &b) HY

27、PERLINK a a认为值参，b认为变量传参 b += a; a += b; return ;int main()int a = 5, b = 5;doit(a, b);cout a b;return 0;* 用若干地址传参可以给调用者传回若干值 Void tryit(int &x,int &y,int &z) 调用时: tryit(a,b,c) ,可以传回 a,b,c的值。* 用数组名（也是地址）传参可以传回整组的数据 Void tryit( int a) 调用时: tryit(x),可以传回整个数组。例如：void tryit(int a) for(int i=0;i=10;i+) ai

28、=i*2; return; int main() int x10; tryit(x); for(int i=0;i=10;i+) coutxiendl; system(pause); return 0;*用指向函数的指针作为参数，可以执行指定的函数。(略)STL 的两个应用：* C+ 快排函数#include Bool com(int a,int b) Return ab;Int main() Int a10=5,7,3,2,6,8,4,3,5,7;Sort(a,a+10,com); /如果升序可以省略com.For(int i=0;i10;i+) Coutai” “;* 优先队列（以堆排为例

29、）#include #include using namespace std;priority_queue Q;int main() int n,a; cinn; while (n-) cina; Q.push(a); while (!Q.empty() cout Q.top() ; Q.pop(); return 0;* 队列、栈、优先队列 三种数据结构汇总：#include #include priority_queue Q;queue Q1;stack S;int main()Q.push(5); x = Q.top(); Q.pop();Q.empty()；Q.size()Q1.pus

30、h(5); Q1.size(); Q1.front(); Q1.empty(); Q1.pop();S.push(5); S.size(); S.top(); S.emtpy(); S.pop();2、函数 . PascalC/C+funtion 函数名(形参表):返回值类型;函数体返回值类型 函数名(形参表) 主体; return 返回值; 注意在Pascal中： 在函数体的语句部分中，必须有对应的函数标识符赋值的语句，并且这些语句在函数被引用时至少要有一句被执行，最后赋给函数标识符的值就是函数值。或者exit(返回值);下面为求n阶乘的样例：PascalC/C+function fac(n:integer):longint;var sum, i:longint;beginsum := 1;for i := 2 to n do sum := sum * i;fac := sum; /exit(sum);end;int fac(int n)int sum = 1, i; for (i = 2; i = n; i+) sum *= i; return sum;*关于随机函数 （包含） Rand() 产生范围在1,INT_MAX之间的一个整数。