基于MATLAB频分复用系统的研究与仿真设计课程报告

《基于MATLAB频分复用系统的研究与仿真设计课程报告》由会员分享，可在线阅读，更多相关《基于MATLAB频分复用系统的研究与仿真设计课程报告（28页珍藏版）》请在装配图网上搜索。

1、. 电子产品辅助设计与仿真 课程考核报告基于MATLAB频分复用系统的研究与仿真设计完 成 日 期 20XX 12月25 / 28目录1引言12课程设计要求22.1课程设计题目22.2课程设计目的22.3设计要求23设计过程及原理33.1频分复用通信系统模型建立33.2频分复用通信系统理论原理44 MATLAB仿真54.1 语音信号的时域和频域仿真54.2 复用信号的频谱仿真64.3传输信号的仿真64.4 解调信号的频谱仿真74.5恢复信号的时域与频域仿真85体会与收获9参考文献10附录111引言MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂矩阵实验室MATLAB工

2、作界面。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言如C、Fortran的编辑模式,代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面

3、、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN,C+,JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以

4、用。2课程设计要求2.1课程设计题目基于MATLAB频分复用系统的研究与仿真设计2.2课程设计目的通过本次综合设计,运用已学的课程知识,根据题目要求,使用Matlab软件进行通信系统的设计和仿真,对电子产品辅助设计与仿真课程中涉及的知识原理等方面有一定的感性认识和实践操作能力,从而加深对本课程知识点的理解,使学生应用知识能力、设计能力、调试能力以及报告撰写能力等方面有显著提高。2.3设计要求用Matlab设计一个通信系统,调用M文件编程界面完成设计,M文件编程中可以调用Matlab通信模块置的功能模块,也可以自己编写功能模块。设计完成后,点击M文件界面上的运行按钮,可进行必要的交互式操作,最终

5、通过绘图能显示运行结果。1使用Matlab软件画出采样后语音信号的时域波形和频谱图。2选择合适的高频载波,对采样信号进行调制。3使用Matlab软件画出复用信号的频谱图。4设计合适的带通滤波器,并画出带通滤波器的频率响应。5对滤波后的信号进行解调,画出解调后各路信号的频谱图。6设计低通滤波器,画出低通滤波器的频率响应。恢复信号的时域波形和频谱图。3设计过程及原理3.1频分复用通信系统模型建立频分多址FDMA是使用最早、目前使用较多的一种多址接入方式,广泛应用于卫星通信、移动通信、一点多址微波通信系统中。FDMA通信系统核心的思想是频分复用FDM,复用是一种将若干个彼此独立的信号合并为一个可在同

6、一个信道上传送的复合信号的方法。例如,在通信系统中,语音信号频谱在3003400Hz,而一条干线的通信资源往往远大于传送一路语音信号所需的带宽。这时,如果用一条干线只传一路语音信号会使资源大大的浪费,所以常用的方法是复用,使一条干线上同时传输几路信号,提高资源利用率。频分复用FDM是信道复用按频率区分信号,即将信号资源划分为多个子频带,每个子频带占用不同的频率,如图3.1所示。然后把需要在同一信道上同时传输的多个信号的频谱调制到不同的频带上,合并在一起不会相互影响,并且能再接收端彼此分离开。混频过程的时域表示式为： 3-1图 3.1 频分复用的子频带划分3.2频分复用通信系统理论原理抽样的理论

7、基础是抽烟定理,它说明在什么条件下能从抽样输出信号f0中恢复输入信号f。根据频谱分析理论,只有抽样信号的频率不发生重叠现象时,抽样的频谱才能与信号频谱相一致。因此,抽样定理可表述为：为了使抽样信号f0能完全恢复连续信号f,抽样信号重复频率 必须大于等于2倍的 , 为包含任何干扰在的信号f的最高有效频率。由于实际滤波器特性的不理想,抽样频率 通常都有高于 ,一般取3到5倍 。语音信号频谱在3003400Hz,语音采样频率 必须大于6.8KHz。在MATLAB数据采集箱中提供语音采集wavrecord命令,wavrecord命令利用Windows 音频输入设备记录声音,其调用形式为:wavreco

8、rd 。利用Windows音频输入设备记录n个音频采样, 频率为fs Hz ,通道数为ch。采样值返回到一个大小为n*ch 的矩阵中。缺省时,fs = 11025 ,ch = 1。其中MATLAB提供的标准音频采样频率有:8000、11025、22050 和44100Hz。为了保证语音的质量,本次设计中取语音信号的采用频率为44100Hz,该采样频率为语音信号CD音质。语音信号采集后,可以用MATLAB数据采集箱中wavwrite命令保存采集的语音信号。如图3.2所示图 3.2 抽样过程波形4MATLAB仿真4.1 语音信号的时域和频域仿真1 信号的时域仿真使用MATLAB软件可以对采集的语音

9、信号进行时域和频域分析。可以使用subplotm,n,p或者subplotm n p将多个图画到一个平面上的工具。其中,m表示是图排成m行,n表示图排成n列,也就是整个figure中有n个图是排成一行的,一共m行,p则是指要把曲线画到figure中哪个图上。MATLAB中绘图命令plot,其含义是以x为横坐标,y为纵坐标,绘制图形。可得到如图4.1所示的时域分析图图4.1 声音样本的时域分析2信号频域仿真频域分析主要是将3个声音样本信号sd1、sd2和sd3用MATLAB软件进行快速傅里叶变换后,再画出3个信号的频谱图。其中快速傅里叶变换可以直接用MATLAB中的fft命令,然后通过abs得到

10、经过快速傅里叶变换后信号的振幅。最后用MATLAB中stem命令对于得到的离散序列实现其频谱图的绘制。4.2 复用信号的频谱仿真在MATLAB软件中将采样的3路语音信号经过混频处理得到3路已调信号x1、x2和x3,再通过加法器将3路信号变为一路复用信号s,通过MATLAB软件中stemt,absfft,.命令对复用信号s进行了频谱分析,其频谱分析如图4.2所示。图4.2 声音样本的频谱分析4.3传输信号的仿真我们都知道FDMA通信系统的复用信号传输是通过空气介质传输的,复用信号在空气传输中会有很多的噪声,其中主要是以高斯白噪声为主,所以在信号传输的设计仿真中,主要对复用信号加入高斯白噪声。在M

11、ATLAB中可以通过awgn函数在某一信号中加入高斯白噪声,其调用方式为：y = awgn,其意义是在信号x中加入高斯白噪声；信噪比SNR以dB为单位,x的强度假定为0dBW。如果x是复数,就加入复噪声。通过前面的调制和信号复用设计后,得到了复用信号s,使用MATLAB中的awgn函数加入高斯白噪声后复用信号变为ys。为了使后面能够较好的恢复语音信号,所以在这里加入白噪声时,信噪比不能设置的太小。仿真发现大于20dB时失真比较小。图4.3为加入高斯白噪声后,复用信号ys的频谱图。图4.3 加入高斯白噪声后复用信号的频谱分析4.4 解调信号的频谱仿真信号解调前,首先通过3个带通滤波器对复用信号s

12、进行滤波,得到3路调制的语音信息y1、y2和y3,然后在对这三路信号进行解调,解调过程与调制的过程相同,使用与原来调制载波相同的信号分别与滤波后的3路信号相乘。得到3路解调信号y01、y02和y03。然后对各路信号使用MATLAB软件中的快速傅里叶变换函数fft进行变换,并通过MATLAB软件,得到的3路解调信号的频谱如图4.4所示。图4.4 解调后信号的频谱图4.5恢复信号的时域与频域仿真语音信号的恢复就是将前面解调所得到的3路信号y01、y02和y03再通过低通滤波器使用filter函数滤波后,分别得到3路恢复的语音信号。然后调用MATLAB中的plot函数和subplot函数对恢复的3路

13、语音信号进行时域分析,其时域分析波形如图4.5所示。图4.5 恢复信号的时域波形5体会与收获经这次课程设计,我不仅复习巩固了课堂所学的理论知识,还提高了对所学知识的综合应用。同时,在以前课本学习中没有弄懂的问题,通过这次课程设计,我都有了更深入的理解。比如通信原理中的时域采样定理、滤波器参数设计等。在设计经过不断的修改调试,在MATLAB上仿真频分多址通信技术取得了较好的效果。录音的声音再经过调试和解调后的信号与原来相比较为接近。我觉得仿真的成功关键在于载波频率的选择以及带通和低通滤波器的参数设计。另外在低通滤波阶段,得到的恢复信号与原始信号基本一致,但是在t=0附近有所失真,这是由于频谱混叠

14、所致,各信号频谱混叠部分均为高频部分,恢复信号在附近的波峰变换最快。即频率最高的区域,引起高频部分失真,这是因为录音期间引入频率高于语音信号的噪声,所以如果在完全无噪音的环境中进行录音,可得无失真的恢复信号。仿真结果分析表明,信号在频分复用时还存在着频间干扰的问题,对此,采用了适当加大采样频率的方法,在较大程度上使该问题得以解决至于完全消除频谱间的干扰,还有待进一步研究与完善。参考文献1 志涌 精通MATLAB 6.5版教程.：航天航空大学, 2003 2 胡鹏,徐会燕.基于matlab的图像去噪算法的研究与实现电脑,20093 彦军,红旗等.改进的中值滤波图像去噪方法研究计算机工程与设计,2

15、0094 宏琦,施维颖,巨永峰.利用中值滤波进行图像处理长安大学学报自然科学版 ,20035 周建兴,MATLAB从入门到精通. 人民邮电 2008附录放入程序源代码附录一：MATLAB仿真程序%1获取录音文件 pause fs=44100; %声音的采样频率为44.1Khz duration=3; %录音时间为3s fprintf; pause fprintf; sd1=wavrecord; %duration*fs每次获得总的采样数为132300,保存声音文件名为sd1 fprintf; wavplay; fprintf; wavwrite; %将录音文件保存为WAV格式的声音文件 fpr

16、intf; pause fprintf; sd2=wavrecord; fprintf; wavplay; fprintf; wavwrite; fprintf; pause fprintf; sd3=wavrecord; fprintf; wavplay; fprintf; wavwrite;%2声音样本的时域和频域分析 fprintf; pause fs=44100; %声音的采样频率为44.1Khz duration=3; t=0:duration*fs-1; %总的采样数 sd1,fs=wavread; %打开保存的录音文件 sd2,fs=wavread; sd3,fs=wavread

17、; figure %图一为三个声音样本的时域波形 subplot plot;xlabel;ylabel; title; subplot plot;xlabel;ylabel; subplot plot;xlabel;ylabel; fprintf; pause figure %图二为三个声音样本的频谱分析 subplot stemt,absfft,.; %fft对声音信号进行快速傅里叶变换xlabel;ylabel; title; subplot stemt,absfft,.;xlabel;ylabel; subplot stemt,absfft,.;xlabel;ylabel; %3调制,将

18、三个声音信号用高频载波进行调制 fprintf; pause x1=4*sd1.*cos; x2=4*sd2.*cos; x3=4*sd3.*cos; s=x1+x2+x3; figure stemt,absfft,.;xlabel;ylabel; title; %4信号传输仿真设计 fprintf; %加入高斯白噪声 pause ys=awgn; snr=10*log10/*; snr %计算信噪比 figure stemt,absfft,.;xlabel;ylabel; title; %5带通滤波器的设计 fprintf; pause Rp=0.5; Rs=40; Wp1=4000 800

19、0/22050; Ws1=3800 8500/22050; n1,Wn1=cheb2ord; b1,a1=cheby2; h1,w1=freqz; mag1=abs; db1=20*log10/max; Wp2=9000 13000/22050; Ws2=8000 14000/22050; n2,Wn2=cheb2ord; b2,a2=cheby2; h2,w2=freqz; mag2=abs; db2=20*log10/max; Wp3=14500 18500/22050; Ws3=14000 19000/22050; n3,Wn3=cheb2ord; b3,a3=cheby2; h3,w3

20、=freqz; mag3=abs; db3=20*log10/max; figure; subplot; plot;axis;xlabel;ylabel20lg|H|; title; subplot; plot;axis;xlabel;ylabel20lg|H|; subplot; plot;axis;xlabel;ylabel20lg|H|; y1=filter; y2=filter; y3=filter; %解调 fprintf; pause fs=44100; y01=y1.*cos; y02=y2.*cos; y03=y3.*cos; figure subplot stemt,absf

21、ft,.;xlabel;ylabel; title; subplot stemt,absfft,.;xlabel;ylabel; subplot stemt,absfft,.;xlabel;ylabel; %7低通滤波 fprintf; pause Rp=0.5; Rs=40; Wp1=3400/22050; Ws1=4000/22050; n1,Wn1=cheb2ord; b1,a1=cheby2; h1,w1=freqz; mag1=abs; db1=20*log10/max; figure; plot; axis; xlabel; ylabel20lg|H|;title; %8恢复信号的

22、时域波形和频谱分析 fprintf; pause yy1=filter; yy2=filter; yy3=filter; figure %图七恢复信号的时域波形 subplot plot;xlabel;ylabel;title; subplot plot;xlabel;ylabel; subplot plot;xlabel;ylabel; fprintf; pause figure %图八恢复信号的频谱分析 subplot stemt,absfft;xlabel;ylabel; title; subplot stemt,absfft;xlabel;ylabel; subplot stemt,absfft;xlabel;ylabel; fprintf; pause fprintf; wavplay; wavplay; wavplay; fprintf;