调幅通信系统数字仿真课程设计.ppt

1、数字通信计算机仿真课程设计,主讲人：胡正平,燕山大学电子与通信工程系,课程设计的价值,利用所学知识解决一些问题！计算机编程进行仿真！数字信号处理！通信原理基本知识！通过这个课程设计能让你学习不少东西，学会思考，学会动手，学会查资料，学会相互配合。,一、基础,什么是调幅？使载波振幅按照调制信号改变的调制方式叫调幅。经过调幅的电波叫调幅波。它保持着高频载波的频率特性，但包络线的形状则和信号波形相似。调幅波的振幅大小，由调制信号的强度决定。,为什么要调制？信号调制可以将信号的频谱搬移到任意位置，从而有利于信号的传送，并且使频谱资源得到充分利用。有效辐射；频分复用。,二、目的,通过本次课程设计使学生深

2、入理解和掌握调幅通信系统的各个关键环节，包括调制、解调、滤波、传输、噪声对通信质量的影响等。在数字信号处理实验课的基础上更加深入地掌握数字滤波器的设计原理及实现方法。使学生对系统各关键点的信号波形及频谱有深刻的认识。,三、意义,本次设计是对学生综合能力的检验，它涉及三门主干课程，包括通信原理、数字信号处理、C/C+语言程序设计。通过本次设计对学生的综合运用专业基础知识及软件设计能力也会有较大提高。,四、设计内容,本次设计的主要内容是用软件模拟调幅（AM）通信系统。原理如下所示：,作用？,作用？,作用？,作用？,作用？,输入的信号可采用学号的后4位（例如3199）进行编码作为基带信号。基带信号经

3、过低通滤波后再与载波相乘进行调制，经过调制后的信号送入信道传输，在传输过程中，有用信号会受到各种信道噪声影响，这里用高斯白噪声模拟信道噪声。在接收端先经过带通滤波器提取信号，再采用相干解调恢复基带信号，进而进行抽样判决，解调出传输的学号信息。,五、系统各测试点信号,输入信号,低通滤波器,低通滤波后的信号,调制后的信号,白噪声,加入白噪声,带通滤波器,带通滤波后的信号,相干解调,低通滤波后的信号,抽样判决后的信号,六、设计要求及关键问题,基带输入信号基带输入信号采用各位同学的学号的后四位，采用双极性码表示，例如：3199D=0011000110011001B-1-1+1+1-1-1-1+1+1-

4、1-1+1+1-1-1+1共16个码元。基带码元宽度为1/200s。基带信号带宽为200Hz。理想冲击信号的带宽是无穷大。对于我们能够实现的脉冲信号，如下图所示：,但是对于本设计中的数字信号来讲，想恢复出“-1”、“+1”信息只需判决其是否过零即可。因此正确判决并不需要全部信号精确重建，允许存在失真。这样我们可以舍掉部分谐波成分，也就是使用低通滤波器，降低信号的带宽，避免了采样的混叠失真（当信号最高频率超过折叠频率就会发生混叠），这样有利于传输，并且不影响重构数字信号。,采样频率fs的确定,由于本次课程设计可以选用了TurboC2.0的开发环境，画图功能挺面，所以你选择采样频率时就要受到限制。

5、在不修改寄存器的前提下，TurboC2.0的开发环境支持最大640480的分辨率。建议利用visualC+进行设计。,如果fs=8000Hz，那么基带信号的数据点数为：168000200=640点，正好达到最大分辨率。一个载波周期8个点。如果fs=6400Hz，那么基带信号的数据点数为：512点，一个载波周期6.4个点。如果fs=4000Hz，那么基带信号的数据点数为：320点，一个载波周期4个点。,建议大家使用VC+或VC+.net，编写托管（如.net）或非托管（如MFC）程序都可以，那样就没有图形界面上的限制了。,滤波器设计要求采用线性相位FIR滤波器，可采用两种方法设计：频率采样法和窗

6、函数法。,频率采样法已知频率响应求时域脉冲响应h(n)。低通和带通滤波器频率响应如下图所示：,频率对应关系如下：频率间隔=fs/N，如果fs=8000Hz，那么4000Hz对应频率(最高频率！)。,频率采样法设计FIR滤波器示意图,(N-1)/2,(N-1)/2,(N-1)/2,N为奇数时：利用图形对应关系，求出Hg(k)，然后利用式，求Hd(k)。,窗函数法确定理想频率响应,带通滤波器：,加窗,载波产生,卷积,voidconvolution(floatxM,floath127,floatyM)intl;floatb11000;l=M+N;for(i=0;iM)b1i=b1i+0;elseb1

7、i+=hj*xi-j;for(i=0;iM;i+)yi=b1i+(N-1)/2;,噪声,VisualC+版本,voidNoise_Gen(double*noise,doubleinten,intn)inten=sqrt(inten*fs/2000)/8;for(inti=0;in;i+)noisei=0;for(intj=0;j12;j+)doubler=rand();doubleri=(2*r-32767)/32767;noisei+=ri;noisei*=inten;,C语言版本,floatGenWhiteNoise(floata)floatb,q,randi=0.0;q=RAND_MAX;for(i=0;i12;i+)randi+=rand()/q;b=a*(randi-6);returnb;,C语言绘图,画图操作intgdriver,gmode;gdriver=DETECT;initgraph(,DFT特点：速度慢但是灵活。,voidDft(floatxM,intN1,floatyM)floatarM,aiM;for(k=0;kN1;k+)ark=0.0;aik=0.0;for(n=0;nN1;n+)ark=ar