用matlab画线段图,基于matlab的直接序列扩频通信系统课程设计报告

基于matlab的直接序列扩频通信系统课程设计报告

1《扩频通信原理》课程设计报告题目：直接扩频系统仿真班级：0110910 和 0110911姓名： 詹晓丹(2009210432)姜 微(2009210503)张建华(2009210336)指导老师：李兆玉21. 课程设计目的(1) 了解、掌握直接扩频通信系统的组成、工作原理；(2) 了解、熟悉扩频调制、解调、解扩方法，并分析其性能；(3) 学习、掌握 Matlab 相关编程知识并用其实现仿真的直接扩频通信系统；2. 课程设计实验原理直接扩频通信系统工作原理：直接序列扩频，就是直接用高码率的扩频码序列在发端去扩展信号的频谱，在收端用相同的扩频码去解扩，把展宽的扩频信号还原成原始的基带信号。 在发端输入的信息与扩频码发生器产生的伪随机码序列(这里使用的是 m 序列)进行波形相乘，得到复合信号，实现信号频谱的展宽，展宽后的信号再调制射频载波发送出去。由于采用平衡调制可以提高系统抗侦波的能力，所以直接序列扩频调制一般都采用二相平衡调制方式。一般扩频调制时一个信息码包含一个周期的伪码，用扩频后的复合信号对载波进行二相相移监控(BPSK)调制，当 gt 从“0”变成“1”或从“1”变到“0”时，载波相位发生 180 度相移。 接收端的本振信号与发射端射频载波相差一个中频，接收端收到的宽带射频信号与本振信号混频、低频滤波后得到中频信号，然后与本地产生的与发端相同并且同步的扩频码序列进行波形相乘，实现相关解扩，再经信息解调，恢复出原始信号。3. 建立模型描述(1) 直接扩频通信系统组成框图：(2) 直接扩频通信系统波形图：34. 模块功能分析(1) 直扩系统的调制功能模块：(都包含模块框图和不同调制、解调方式介绍、分析)(a) 扩频调制模块用扩频码发生器产生一个伪随机码 pn(这里用的是 m 序列) ，与信源信息码序列 xt 相乘，实现频谱的展宽 (b) BPSK 调制模块4调制的方式可以有二相相移监控 BPSK、四相相移键控 QPSK、偏移四相相移监控OQPSK、最小频移监控 MSK。QPSK 调制的目的是节省频谱，但在扩频系统中有时候带宽的利用并不是最重要的；OQPSK 的优点就是调制信号的相位改变没有倒 π 现象；MSK 调制信号时可以避免相位突变，由于以上调制方式实现比较复杂，所以我们选用扩频系统中最常用的 BPSK 调制方式。(2) 直扩系统的解调功能模块：(a) BPSK 解调模块在常规数字通信中，解调可以用锁相环解调器、平方环解调器、科斯塔斯环解调器。在直扩系统中，一般扩频调制方式是用抑制载波双平衡调制来产生 BPSK 信号的，而对于BPSK 信号，不管是绝对相移还是差分相移。其载波分量都被抑制了几十分贝，并且直扩信号谱密度都很低，而大气噪声及接收机内部噪声又很大，有用信号常淹没在噪声中，所以用一般的锁相环难以提取载波。而平方环虽然便于载波提取，但环路工作在二倍频后的频率上，工作频率高，环路稳定性较差。我们选用的是科斯塔斯环，因为它的突出优点是能够解调移相键控信号和抑制了载波的信号，且环路的工作频率与输入信号载波频率完全相同。(b) 扩频解调模块解扩方式有相关解扩、直接式相关器解扩、外差式相关器解扩、序列匹配滤波器解扩。5直接式相关器的优点是结构简单，缺点是对干扰信号有直通和码速率泄露现象外差式相关器的抗干扰能力较低；由于相关解扩在性能上很好，在接收端产生的本地 pn’码，可以用科斯塔斯环实现与发端的 pn 码精确的同步。5. 模块源代码及调试过程(1)直扩系统的调制模块(a)信息码生成模块code_length=20; %信息码元个数N=1:code_length;rand( seed ,0);x=sign(rand(1,code_length)-0.5); %信息码从 0、1 序列变成-1、1 序列x1=rectpulse(x,800);%每个码元内采样 800 个点plot(x1);axis([0 16000 -1.5 1.5]);title( 信源信息码序列 );grid on;生成信息码的波形图(b)伪随机码生成模块function y = mgen(g,state,N)%输入 g:m 序列生成多项式(10 进制输入)%state:寄存器初始状态(10 进制输入)%N:输出序列长度figure(1)g=19;state=8;N=2000;gen = dec2bin(g)-48;M = length(gen);curState = dec2bin(state,M-1) - 48;for k =1:Ny(k) = curState(M-1);a = rem(sum( gen(2:end).*curState),2);curState = [a curState(1:M-2)];6endx_code=sign(y-0.5); pn=rectpulse(x_code,8);%每个伪码元内采样 8 个点plot(pn);axis([0 600 -1.5 1.5]);title( 伪随机码序列 );grid on;生成伪随机码的波形图(c)扩频调制模块gt=x1.*pn;plot(gt);axis([0 1000 -1.5 1.5]);title( 复合码序列 );grid on;生成的复合码波形图7(d)BPSK 调制模块%用 BPSK 调制fs=20e6; f0=30e6;for i=1:2000AI=2;dt=fs/f0;n=0:dt/7:dt; %一个载波周期内采样 8 个点cI=AI*cos(2*pi*f0*n/fs);x_bpsk((1+(i-1)*8):i*8)=gt((1+(i-1)*8):i*8).*cIendplot(x_bpsk);axis([0 200 -2.5 2.5]);title( BPSK 调制后的波形 );grid on;生成 BPSK 调制后的波形图(2)加噪模块sigma=0.1;nt=sigma*randn(1,20);nt1=rectpulse(nt,800);gt1=gt+nt1;fs=20e6; f0=30e6;for i=1:2000AI=2;dt=fs/f0;n=0:dt/7:dt; %一个载波周期内采样 8 个点cI=AI*cos(2*pi*f0*n/fs);x_bpsk1((1+(i-1)*8):i*8)=gt1((1+(i-1)*8):i*8).*cIendplot(x_bpsk1);axis([0 200 -2.5 2.5]);8title( 加噪后已调波的波形 );grid on;生成加噪后已调波的波形图(3)直扩系统的解调模块(a)BPSK 解调模块AI=1;dt=fs/f0;n=0:dt/7:dt; %一个载波周期内采样八个点cI=AI*cos(2*pi*f0*n/fs);for i=1:2000s((1+(i-1)*8):i*8)= x_bpsk1((1+(i-1)*8):i*8).*cI;endplot(s);axis([0 200 -2.5 2.5]);title( 解调后的波形 );grid on;生成 BPSK 解调后的波形图9(b)解扩模块%相关解扩jk_code=s.*pn;%低通滤波wn=5