数字频带通信系统仿真

1课程设计说明书设计及说明主要结果引言数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。通常使用键控法来实现数字调制，比如对载波的振幅、频率和相位进行键控。其次，还有4PSK、16QAM等调制方式。1、任务和要求1.1、课程设计目的☞通过课程设计加深对《通信原理》课程基本内容的理解☞学习与掌握通信系统仿真的基本方法1.2、课程设计要求1）完成2ASK、2FSK、2PSK、4PSK、16QAM五种调制系统的仿真，要求有调制、信道、解调三个部分；2）信道为高斯信道；3）系统参数可自行设置；4）仿真结果要求给出已调信号时域波形和频谱，解调后信号的时域波形与频谱。21.3、MATLAB仿真平台简介MATLAB是矩阵实验室（MatrixLaboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。32、2ASK调制解调原理与系统仿真2.1、2ASK调制原理与仿真数字幅度调制又称幅度键控（ASK），二进制幅度键控记作2ASK。2ASK是利用代表数字信息“0”或“1”的基带矩形脉冲去键控一个连续的载波，使载波时断时续地输出。有载波输出时表示发送“1”，无载波输出时表示发送“0”。2ASK信号可表示为twtsteccos)()(0（2-1）式中，cw为载波角频率，s(t)为单极性NRZ矩形脉冲序列)()(bnnnTtgats（2-2）其中，g(t)是持续时间bT、高度为1的矩形脉冲，常称为门函数；na为二进制数字PPan101，出现概率为，出现概率为（2-3）2ASK/OOK信号的产生方法通常有两种：模拟调制（相乘器法）和键控法。本课程设计运用模拟幅度调制的方法，用乘法器实现。相应的调制如图2-3：图2-3模拟相乘法乘法器)(2teASK)(tstccos4调制代码：clear;%清除空间变量m=[10110110];%数字信号序列Lm=length(m);%序列的长度F=150;%数字信号的带宽f=600;%正弦载波信号的频率A=1;%载波的幅度Q=f/F;%频率比，即一个码元宽度中的正弦周期个数,为适配下面滤波器参数%%选%%取，推荐将Q设为Q=1/3M=500;%一个正弦周期内的采样点数t=(0:M-1)/M/f;%一个正弦信号周期内的时间carry1=repmat(A*sin(2*pi*f*t),1,Q);%一个码元宽度内的正弦载波信号Lcarry1=length(carry1);%一个码元宽度内的信号长度carry2=kron(ones(size(m)),carry1);%载波信号ask=kron(m,carry1);%调制后的信号V=0.002;ask1=imnoise(ask,'gaussian',0,V);%加高斯白噪声后的已调信号Fask=fft(ask);%对信号进行傅里叶变换Fask=abs(Fask(1:length(Fask)/2+1));%已调信号频谱f1=(0:length(Fask)-1)*f/length(Fask)/2;%已调信号频谱带宽N=length(ask);%长度tau=(0:N-1)/(M-1)/f;%时间Tmin=min(tau);%最小时刻5Tmax=max(tau);%最大时刻T=ones(size(carry1));%一个数字信号‘1’dsig=kron(m,T);%数字信号波形subplot(2,2,1)%子图分割plot(tau,dsig),title('数字信号波形')%画出数字信号的波形图gridonaxis([TminTmax-0.21.2])%设置坐标范围subplot(2,2,2)%子图分割plot(tau,carry2),title('载波波形')%画出载波波形gridonaxis([TminTmax-1.2*A1.2*A])%设置坐标范围subplot(2,2,3),plot(tau,ask)%画出调制后的波形title('调制信号')%子图分割gridonaxis([TminTmax-1.2*A1.2*A])%设置坐标范围subplot(224);plot(f1,Fask),title('已调信号频谱');figure(2)subplot(111),plot(tau,ask1)%画出调制加噪波形title('调制后加噪声波形');gridonaxis([TminTmax-1.2*A1.2*A])%设置坐标范围运行结果：6信号经过信道后，会受到噪声的干扰，假设信道噪声为高斯白噪声。信号加噪后波形：2.2、2ASK解调原理与仿真2ASK有两种基本解调方法：相干解调法（同步检测法）和非相干解调法（包络检波法）。相干解调需要将载频位置的已调信号频谱重新搬回原始基带位置，因此用相乘器与载波相乘来实现。为确fc7带通整流低通抽样判决判决定时脉冲S(t)eo(t)abcd保无失真还原信号，必须在接收端提供一个与调制载波严格同步的本地载波，这是整个解调过程能否顺利完好进行的关键。解调过程如图5所示。图5相干解调包络检波器通常由整流器和低通滤波器组成，可以直接从已调波中提取原始基带信号，结构简单，如图6所示。经过各个模块后波形变化如图7所示。图6非相干解调图7非相干解调过程的时间波形本实验采用非相干解调方式。解调代码：S(t)定时脉冲带通相乘低通抽样判决eo(t)cosωct2fs11100000101abcd8clear;%清除空间变量m=[10110110];%数字信号序列Lm=length(m);%序列的长度F=150;%数字信号的带宽f=600;%正弦载波信号的频率A=1;%载波的幅度Q=f/F;%频率比即一个码元宽度中的正弦周期个数,为适配下面滤波器参数选%取M=500;%一个正弦周期内的采样点数t=(0:M-1)/M/f;%一个正弦信号周期内的时间carry1=repmat(A*sin(2*pi*f*t),1,Q);%一个码元宽度内的正弦载波信号Lcarry1=length(carry1);%一个码元宽度内的信号长度carry2=kron(ones(size(m)),carry1);%载波信号ask=kron(m,carry1);%调制后的信号V=0.002;ask1=imnoise(ask,'gaussian',0,V);%加高斯白噪声后的已调信号sig_mul=ask.*carry2;%已调信号与载波信号相乘N=length(ask);%长度tau=(0:N-1)/(M-1)/f;%时间Tmin=min(tau);%最小时刻Tmax=max(tau);%最大时刻T=ones(size(carry1));%一个数字信号‘1’dsig=kron(m,T);%数字信号波形Fask=fft(ask);%对信号进行傅里叶变换Fask=abs(Fask(1:length(Fask)/2+1));%已调信号频谱f1=(0:length(Fask)-1)*f/length(Fask)/2;%已调信号频谱带宽、%Butterworth模拟低通滤波器设计%获得Butterworth模拟低通原型滤波器的阶数及3dB截止频率[Ord,omega_c]=buttord(4*pi*f*0.6,4*pi*f*0.8,2,30,'s');9%由原型滤波器向实际滤波器转换，获得滤波器的分子，分母系数[num,den]=butter(Ord,omega_c,'s');h=tf(num,den);%获得滤波器传递函数%滤波x=lsim(h,sig_mul,tau);%运用模拟滤波器对信号进行滤波Fx=fft(x);%对信号进行傅里叶变换Fx=abs(Fx(1:length(Fx)/2+1));%解调信号频谱f1=(0:length(Fx)-1)*f/length(Fx)/2;%解调信号频谱带宽subplot(3,2,2)%子图分割plot(tau,x),title('解调输出波形')%画出滤波后的波形gridon%添加网格axis([TminTmax-0.30.8])%设置坐标范围subplot(323);plot(f1,Fx),title('解调信号频谱');%解调信号频谱gridon%添加网格axis([TminTmax-0.30.8])%设置坐标范围th=0.25;%抽样判决的阈值设置t_judge=(0:Lm-1)*Lcarry1+Lcarry1/2;%抽样判决点的选取y=(x(t_judge))';%抽样判决时刻时的信号值y_judge=1*(y=th)+0*(y=th);%抽样判决信号值的0阶保持y_value=kron(y_judge,ones(size(carry1)));%抽样判决后的数字信号波形n_tau=tau+0.5/F;%抽样判决后的信号对应的时间subplot(3,2,4)%子图分割plot(n_tau,y_value),title('抽样判决后的数字信号波形')%画出抽样判决%后%%的数字%%信号波形axis([min(n_tau)max(n_tau)-0.21.2])%设置坐标范围gridonsubplot(3,2,5)%子图分割plot(tau,dsig),title('原始信号波形与解调后的信号作对比')%画出原始%%信号波形与%解调后的信号作对比10axis([TminTmax-0.21.2])%设置坐标范围gridon%添加网格subplot(321),plot(tau,sig_mul)title('已调信号与载波信号相乘波形')%已调信号与%载波信号相乘波形gridonaxis([TminTmax-1.2*A1.2*A])%设置坐标范围运行结果：分析结果：解调信号的频谱没出来，检查代码终究还是没找到问题所在。解调抽样后恢复原二进制序列。3、2FSK调制解调原理与系统仿真3.1、2FSK调制原理112FSK是利用数字基带信号控制在波的频率来传送信息。例如，1码用频率f1来传输，0码用频率f2来传输，而其振幅和初始相位不变。故其表示式为：)cos()cos(21122)(tAtAFSKt时发送时发送10式中，假设码元的初始相位分别为1和2；112fπ和222fπ为两个不同的码元的角频率；幅度为A为一常数，表示码元的包络为矩形脉冲。2FSK信号的产生方法有两种：（1）模拟法，即用数字基带信号作为调制信号进行调频。如图1-1（a）所示。（2）键控法，用数字基带信号)(tg及其反)(tg相分别控制两个开关门电路，以此对两个载波发生器进行选通。如图1-