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

1、 . 电子产品辅助设计与仿真 课程考核报告基于MATLAB频分复用系统的研究与仿真设计系 部：电子与信息工程系专业班级：通信工程10秋2 班姓 名：学 号：指导教师：完 成 日 期 2013年 12月目 录 TOC o 1-3 h z u HYPERLINK l _Toc375937800 1引言 PAGEREF _Toc375937800 h 1 HYPERLINK l _Toc375937801 2课程设计要求 PAGEREF _Toc375937801 h 2 HYPERLINK l _Toc375937802 2.1课程设计题目 PAGEREF _Toc375937802 h 2 HY

2、PERLINK l _Toc375937803 2.2课程设计目的 PAGEREF _Toc375937803 h 2 HYPERLINK l _Toc375937804 2.3设计要求 PAGEREF _Toc375937804 h 2 HYPERLINK l _Toc375937805 3设计过程及原理 PAGEREF _Toc375937805 h 3 HYPERLINK l _Toc375937806 3.1频分复用通信系统模型建立 PAGEREF _Toc375937806 h 3 HYPERLINK l _Toc375937807 3.2频分复用通信系统理论原理 PAGEREF _

3、Toc375937807 h 4 HYPERLINK l _Toc375937808 4 MATLAB仿真 PAGEREF _Toc375937808 h 5 HYPERLINK l _Toc375937809 4.1 语音信号的时域和频域仿真 PAGEREF _Toc375937809 h 5 HYPERLINK l _Toc375937810 4.2 复用信号的频谱仿真 PAGEREF _Toc375937810 h 6 HYPERLINK l _Toc375937811 4.3传输信号的仿真 PAGEREF _Toc375937811 h 6 HYPERLINK l _Toc375937

4、812 4.4 解调信号的频谱仿真 PAGEREF _Toc375937812 h 7 HYPERLINK l _Toc375937813 4.5恢复信号的时域与频域仿真 PAGEREF _Toc375937813 h 8 HYPERLINK l _Toc375937814 5体会与收获 PAGEREF _Toc375937814 h 9 HYPERLINK l _Toc375937815 参考文献 PAGEREF _Toc375937815 h 10 HYPERLINK l _Toc375937816 附录 PAGEREF _Toc375937816 h 111引言MATLAB是matrix

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

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

7、，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。2课程设计要求2.1课程设计题目 基于MATLAB频分复用系统的研究与仿真设计2.2课程设计目的通过本次综合设计，运用已学的课程知识，根据题目要求，使用Matlab软件进行通信系统的设计和仿真,对电子产品辅助设计与仿真课程中涉及的知识原理等方面有一定的感性认识和实践操作能力，从而加深对本课程知识点的理解，使学生应用知识能力、设计能力、调试能力以及报告撰写能力等方面有显著提高。2.3设计要求用Matlab设计一个通信系统，调用M文件编程界面完成设计，M文件编程中可以调用Matlab通信模块内置的功能模块，也可以自

8、己编写功能模块。设计完成后，点击M文件界面上的运行按钮，可进行必要的交互式操作，最终通过绘图能显示运行结果。（1）使用Matlab软件画出采样后语音信号的时域波形和频谱图。（2）选择合适的高频载波，对采样信号进行调制。（3）使用Matlab软件画出复用信号的频谱图。（4）设计合适的带通滤波器，并画出带通滤波器的频率响应。（5）对滤波后的信号进行解调，画出解调后各路信号的频谱图。（6）设计低通滤波器，画出低通滤波器的频率响应。恢复信号的时域波形和频谱图。3设计过程及原理3.1频分复用通信系统模型建立频分多址（FDMA）是使用最早、目前使用较多的一种多址接入方式，广泛应用于卫星通信、移动通信、一点

9、多址微波通信系统中。FDMA通信系统核心的思想是频分复用（FDM），复用是一种将若干个彼此独立的信号合并为一个可在同一个信道上传送的复合信号的方法。例如，在电话通信系统中，语音信号频谱在3003400Hz内，而一条干线的通信资源往往远大于传送一路语音信号所需的带宽。这时，如果用一条干线只传一路语音信号会使资源大大的浪费，所以常用的方法是“复用”，使一条干线上同时传输几路电话信号，提高资源利用率。频分复用（FDM）是信道复用按频率区分信号，即将信号资源划分为多个子频带，每个子频带占用不同的频率，如图3.1所示。然后把需要在同一信道上同时传输的多个信号的频谱调制到不同的频带上，合并在一起不会相互影

10、响，并且能再接收端彼此分离开。混频过程的时域表示式为： （3-1）图 3.1 频分复用的子频带划分3.2频分复用通信系统理论原理抽样的理论基础是抽烟定理，它说明在什么条件下能从抽样输出信号f0(t)中恢复输入信号f(t)。根据频谱分析理论，只有抽样信号的频率不发生重叠现象时，抽样的频谱才能与信号频谱相一致。因此，抽样定理可表述为：为了使抽样信号f0(t)能完全恢复连续信号f(t)，抽样信号重复频率 必须大于等于2倍的 ， 为包含任何干扰在内的信号f(t)的最高有效频率。由于实际滤波器特性的不理想，抽样频率 通常都有高于 ，一般取3到5倍 。语音信号频谱在3003400Hz内，语音采样频率 必须

11、大于6.8KHz。在MATLAB数据采集箱中提供语音采集wavrecord命令，wavrecord命令利用Windows 音频输入设备记录声音,其调用形式为:wavrecord (n ,fs ,ch)。利用Windows音频输入设备记录n个音频采样, 频率为fs Hz ,通道数为ch。采样值返回到一个大小为n*ch 的矩阵中。缺省时,fs = 11025 ,ch = 1。其中MATLAB提供的标准音频采样频率有:8000、11025、22050 和44100Hz。为了保证语音的质量，本次设计中取语音信号的采用频率为44100Hz，该采样频率为语音信号CD音质。语音信号采集后，可以用MATLAB

12、数据采集箱中wavwrite命令保存采集的语音信号。如图3.2所示图 3.2 抽样过程波形4 MATLAB仿真4.1 语音信号的时域和频域仿真 （1） 信号的时域仿真使用MATLAB软件可以对采集的语音信号进行时域和频域分析。可以使用subplot（m,n,p）或者subplot（m n p）将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，也就是整个figure中有n个图是排成一行的，一共m行，p则是指要把曲线画到figure中哪个图上。MATLAB中绘图命令plot(x,y)，其含义是以x为横坐标，y为纵坐标，绘制图形。可得到如图4.1所示的时域分析图图4.1 声音

13、样本的时域分析（2）信号频域仿真频域分析主要是将3个声音样本信号sd1、sd2和sd3用MATLAB软件进行快速傅里叶变换后，再画出3个信号的频谱图。其中快速傅里叶变换可以直接用MATLAB中的fft命令，然后通过abs得到经过快速傅里叶变换后信号的振幅。最后用MATLAB中stem命令对于得到的离散序列实现其频谱图的绘制。4.2 复用信号的频谱仿真在MATLAB软件中将采样的3路语音信号经过混频处理得到3路已调信号x1、x2和x3，再通过加法器将3路信号变为一路复用信号s，通过MATLAB软件中stem(t,abs(fft(s),.)命令对复用信号s进行了频谱分析，其频谱分析如图4.2所示。

14、图4.2 声音样本的频谱分析4.3传输信号的仿真我们都知道FDMA通信系统的复用信号传输是通过空气介质传输的，复用信号在空气传输中会有很多的噪声，其中主要是以高斯白噪声为主，所以在信号传输的设计仿真中，主要对复用信号加入高斯白噪声。在MATLAB中可以通过awgn函数在某一信号中加入高斯白噪声，其调用方式为：y = awgn(x,SNR)，其意义是在信号x中加入高斯白噪声；信噪比SNR以dB为单位，x的强度假定为0dBW。如果x是复数，就加入复噪声。通过前面的调制和信号复用设计后，得到了复用信号s，使用MATLAB中的awgn函数加入高斯白噪声后复用信号变为ys。为了使后面能够较好的恢复语音信

15、号，所以在这里加入白噪声时，信噪比不能设置的太小。仿真发现大于20dB时失真比较小。图4.3为加入高斯白噪声后，复用信号ys的频谱图。图4.3 加入高斯白噪声后复用信号的频谱分析4.4 解调信号的频谱仿真 信号解调前，首先通过3个带通滤波器对复用信号s进行滤波，得到3路调制的语音信息y1、y2和y3，然后在对这三路信号进行解调，解调过程与调制的过程相同，使用与原来调制载波相同的信号分别与滤波后的3路信号相乘。得到3路解调信号y01、y02和y03。然后对各路信号使用MATLAB软件中的快速傅里叶变换函数fft进行变换，并通过MATLAB软件，得到的3路解调信号的频谱如图4.4所示。图4.4 解

16、调后信号的频谱图4.5恢复信号的时域与频域仿真语音信号的恢复就是将前面解调所得到的3路信号y01、y02和y03再通过低通滤波器使用filter函数滤波后，分别得到3路恢复的语音信号。然后调用MATLAB中的plot(t,yy1)函数和subplot函数对恢复的3路语音信号进行时域分析，其时域分析波形如图4.5所示。图4.5 恢复信号的时域波形5体会与收获经这次课程设计，我不仅复习巩固了课堂所学的理论知识，还提高了对所学知识的综合应用。同时，在以前课本学习中没有弄懂的问题，通过这次课程设计，我都有了更深入的理解。比如通信原理中的时域采样定理、滤波器参数设计等。在设计经过不断的修改调试，在MAT

17、LAB上仿真频分多址通信技术取得了较好的效果。录音的声音再经过调试和解调后的信号与原来相比较为接近。我觉得仿真的成功关键在于载波频率的选择以及带通和低通滤波器的参数设计。另外在低通滤波阶段，得到的恢复信号与原始信号基本一致，但是在t=0附近有所失真，这是由于频谱混叠所致，各信号频谱混叠部分均为高频部分，恢复信号在附近的波峰变换最快。即频率最高的区域，引起高频部分失真，这是因为录音期间引入频率高于语音信号的噪声，所以如果在完全无噪音的环境中进行录音，可得无失真的恢复信号。仿真结果分析表明，信号在频分复用时还存在着频间干扰的问题，对此，采用了适当加大采样频率的方法，在较大程度上使该问题得以解决至于

18、完全消除频谱间的干扰，还有待进一步研究与完善。参考文献1 张志涌 精通MATLAB 6.5版教程.北京：北京航天航空大学出版社, 2003(56) 2 胡鹏,徐会燕.基于matlab的图像去噪算法的研究与实现福建电脑，2009(12)3 李彦军,苏红旗等.改进的中值滤波图像去噪方法研究计算机工程与设计，2009(12)4 孙宏琦，施维颖,巨永峰.利用中值滤波进行图像处理长安大学学报（自然科学版） ，2003(2)5 周建兴,MATLAB从入门到精通. 人民邮电出版社 2008(60)附录放入程序源代码附录一：MATLAB仿真程序%（1）获取录音文件 pause fs=44100; %声音的采样

19、频率为44.1Khz duration=3; %录音时间为3s fprintf(按任意键开始录音1：n); pause fprintf(录音中n); sd1=wavrecord(duration*fs,fs); %duration*fs每次获得总的采样数为132300，保存声音文件名为sd1 fprintf(放音中n); wavplay(sd1,fs); fprintf(录音1播放完毕。n); wavwrite(sd1,fs,sound1.wav); %将录音文件保存为WAV格式的声音文件 fprintf(按任意键开始录音2：n); pause fprintf(录音中n); sd2=wavre

20、cord(duration*fs,fs); fprintf(放音中n); wavplay(sd2,fs); fprintf(录音2播放完毕。n); wavwrite(sd2,fs,sound2.wav); fprintf(按任意键开始录音3：n); pause fprintf(录音中n); sd3=wavrecord(duration*fs,fs); fprintf(放音中n); wavplay(sd3,fs); fprintf(录音3播放完毕。n); wavwrite(sd3,fs,sound3.wav);%（2）声音样本的时域和频域分析 fprintf(按任意键开始声音样本的时域分析：n)

21、; pause fs=44100; %声音的采样频率为44.1Khz duration=3; t=0:duration*fs-1; %总的采样数 sd1,fs=wavread(sound1.wav); %打开保存的录音文件 sd2,fs=wavread(sound2.wav); sd3,fs=wavread(sound3.wav); figure(1) %图一为三个声音样本的时域波形 subplot(311) plot(t,sd1);xlabel(单位：s);ylabel(幅度); title(三个声音样本的时域波形); subplot(312) plot(t,sd2);xlabel(单位：s

22、);ylabel(幅度); subplot(313) plot(t,sd3);xlabel(单位：s);ylabel(幅度); fprintf(按任意键开始声音样本的频域分析：n); pause figure(2) %图二为三个声音样本的频谱分析 subplot(311) stem(t,abs(fft(sd1),.); %fft对声音信号进行快速傅里叶变换xlabel(单位：Hz);ylabel(幅度); title(三个声音样本的频谱分析); subplot(312) stem(t,abs(fft(sd2),.);xlabel(单位：Hz);ylabel(幅度); subplot(313)

23、stem(t,abs(fft(sd3),.);xlabel(单位：Hz);ylabel(幅度); %（3）调制，将三个声音信号用高频载波进行调制 fprintf(按任意键开始信号的调制和复用信号频域分析：n); pause x1=4*sd1.*cos(2*pi*4000*t/fs); x2=4*sd2.*cos(2*pi*11000*t/fs); x3=4*sd3.*cos(2*pi*18000*t/fs); s=x1+x2+x3; figure(3) stem(t,abs(fft(s),.);xlabel(单位：Hz);ylabel(幅度); title(复用信号的频谱分析); %（4）信号

24、传输仿真设计 fprintf(按任意键开始信道仿真设计：n); %加入高斯白噪声 pause ys=awgn(s,20); snr=10*log10(s*s)/(s-ys)*(s-ys); snr %计算信噪比 figure(4) stem(t,abs(fft(ys),.);xlabel(单位：Hz);ylabel(幅度); title(加入高斯白噪声后复用信号的频谱分析); %（5）带通滤波器的设计 fprintf(按任意键开始带通滤波器的设计：n); pause Rp=0.5; Rs=40; Wp1=4000 8000/22050; Ws1=3800 8500/22050; n1,Wn1=

25、cheb2ord(Wp1,Ws1,Rp,Rs); b1,a1=cheby2(n1,Rs,Wn1); h1,w1=freqz(b1,a1); mag1=abs(h1); db1=20*log10(mag1+eps)/max(mag1); Wp2=9000 13000/22050; Ws2=8000 14000/22050; n2,Wn2=cheb2ord(Wp2,Ws2,Rp,Rs); b2,a2=cheby2(n2,Rs,Wn2); h2,w2=freqz(b2,a2); mag2=abs(h2); db2=20*log10(mag2+eps)/max(mag2); Wp3=14500 185

26、00/22050; Ws3=14000 19000/22050; n3,Wn3=cheb2ord(Wp3,Ws3,Rp,Rs); b3,a3=cheby2(n3,Rs,Wn3); h3,w3=freqz(b3,a3); mag3=abs(h3); db3=20*log10(mag3+eps)/max(mag3); figure(5); subplot(3,1,1); plot(w1/pi,db1);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|); title(用切比雪夫2型设计三个带通滤波器); subplot(3,1,2); plot(w

27、2/pi,db2);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|); subplot(3,1,3); plot(w3/pi,db3);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|); y1=filter(b1,a1,ys); y2=filter(b2,a2,ys); y3=filter(b3,a3,ys); %(6)解调 fprintf(按任意键开始信号的解调和3路信号频域分析：n); pause fs=44100; y01=y1.*cos(2*pi*4000*t/fs); y02=y2

28、.*cos(2*pi*11000*t/fs); y03=y3.*cos(2*pi*18000*t/fs); figure(6) subplot(311) stem(t,abs(fft(y01),.);xlabel(单位：Hz);ylabel(幅度); title(解调后的3路信号各自的频谱图); subplot(312) stem(t,abs(fft(y02),.);xlabel(单位：Hz);ylabel(幅度); subplot(313) stem(t,abs(fft(y03),.);xlabel(单位：Hz);ylabel(幅度); %（7）低通滤波 fprintf(按任意键开始低通滤波

29、器的设计：n); pause Rp=0.5; Rs=40; Wp1=3400/22050; Ws1=4000/22050; n1,Wn1=cheb2ord(Wp1,Ws1,Rp,Rs); b1,a1=cheby2(n1,Rs,Wn1); h1,w1=freqz(b1,a1); mag1=abs(h1); db1=20*log10(mag1+eps)/max(mag1); figure(7); plot(w1/pi,db1); axis(0 1 -50 20); xlabel(w/pi); ylabel(20lg|H(ejw)|); title(低通滤波器的频率响应); %（8）恢复信号的时域波

30、形和频谱分析 fprintf(按任意键开始恢复信号以及3路信号时域分析：n); pause yy1=filter(b1,a1,y01); yy2=filter(b1,a1,y02); yy3=filter(b1,a1,y03); figure(8) %图七恢复信号的时域波形 subplot(311) plot(t,yy1);xlabel(单位：s);ylabel(幅度);title(恢复信号的时域波形); subplot(312) plot(t,yy2);xlabel(单位：s);ylabel(幅度); subplot(313) plot(t,yy3);xlabel(单位：s);ylabel(

31、幅度); fprintf(按任意键开始3路信号频域分析：n); pause figure(9) %图八恢复信号的频谱分析 subplot(311) stem(t,abs(fft(yy1);xlabel(单位：Hz);ylabel(幅度); title(恢复信号的频谱分析); subplot(312) stem(t,abs(fft(yy2);xlabel(单位：Hz);ylabel(幅度); subplot(313) stem(t,abs(fft(yy3);xlabel(单位：Hz);ylabel(幅度); fprintf(按任意键使声音再现：n); pause fprintf(放音中n); w

32、avplay(yy1,fs); wavplay(yy2,fs); wavplay(yy3,fs); fprintf(放音完毕，仿真结束！n);附录资料：MATLAB的30个方法1 内部常数pi 圆周率 exp(1)自然对数的底数ei 或j 虚数单位Inf或 inf 无穷大 2 数学运算符a+b 加法a-b减法a*b矩阵乘法a.*b数组乘法a/b矩阵右除ab矩阵左除a./b数组右除a.b数组左除ab 矩阵乘方a.b数组乘方-a负号 共轭转置.一般转置3 关系运算符=等于大于=大于或等于=不等于4 常用内部数学函数 指数函数exp(x)以e为底数对数函数log(x)自然对数，即以e为底数的对数lo

33、g10(x)常用对数，即以10为底数的对数log2(x)以2为底数的x的对数开方函数sqrt(x)表示x的算术平方根绝对值函数abs(x)表示实数的绝对值以及复数的模三角函数（自变量的单位为弧度）sin(x)正弦函数cos(x)余弦函数tan(x)正切函数cot(x)余切函数sec(x)正割函数csc(x)余割函数反三角函数 asin(x)反正弦函数acos(x)反余弦函数atan(x)反正切函数acot(x)反余切函数asec(x)反正割函数acsc(x)反余割函数双曲函数 sinh(x)双曲正弦函数cosh(x)双曲余弦函数tanh(x)双曲正切函数coth(x)双曲余切函数sech(x)

34、双曲正割函数csch(x)双曲余割函数反双曲函数 asinh(x)反双曲正弦函数acosh(x)反双曲余弦函数atanh(x)反双曲正切函数acoth(x)反双曲余切函数asech(x)反双曲正割函数acsch(x)反双曲余割函数求角度函数atan2(y,x)以坐标原点为顶点，x轴正半轴为始边，从原点到点（x，y）的射线为终边的角，其单位为弧度，范围为（ ， 数论函数gcd(a,b)两个整数的最大公约数lcm(a,b)两个整数的最小公倍数排列组合函数factorial(n)阶乘函数，表示n的阶乘 复数函数 real(z)实部函数imag(z)虚部函数abs(z)求复数z的模angle(z)求复

35、数z的辐角，其范围是（ ， conj(z)求复数z的共轭复数求整函数与截尾函数ceil(x)表示大于或等于实数x的最小整数floor(x)表示小于或等于实数x的最大整数round(x)最接近x的整数最大、最小函数max(a，b，c，)求最大数min(a，b，c，)求最小数符号函数 sign(x)5 自定义函数-调用时：“返回值列=M文件名（参数列）”function 返回变量=函数名（输入变量） 注释说明语句段（此部分可有可无）函数体语句 6进行函数的复合运算compose(f,g) 返回值为f(g(y)compose(f,g,z) 返回值为f(g(z)compose(f,g,x,.z) 返回

36、值为f(g(z)compose(f,g,x,y,z) 返回值为f(g(z)7 因式分解syms 表达式中包含的变量 factor(表达式) 8 代数式展开syms 表达式中包含的变量 expand(表达式)9 合并同类项syms 表达式中包含的变量 collect(表达式，指定的变量)10 进行数学式化简syms 表达式中包含的变量 simplify(表达式)11 进行变量替换syms 表达式和代换式中包含的所有变量 subs(表达式，要替换的变量或式子，代换式)12 进行数学式的转换调用Maple中数学式的转换命令，调用格式如下：maple(Maple的数学式转换命令) 即：maple(co

37、nvert(表达式，form)将表达式转换成form的表示方式 maple(convert(表达式，form, x) 指定变量为x，将依赖于变量x的函数转换成form的表示方式（此指令仅对form为exp与sincos的转换式有用） 13 解方程solve(方程，变元) 注：方程的等号用普通的等号： = 14 解不等式调用maple中解不等式的命令即可，调用形式如下： maple(maple中解不等式的命令)具体说，包括以下五种：maple( solve（不等式）) maple( solve（不等式，变元） ) maple( solve（不等式，变元） ) maple( solve（不等式，变

38、元） ) maple( solve（不等式，变元） )15 解不等式组调用maple中解不等式组的命令即可，调用形式如下： maple(maple中解不等式组的命令) 即：maple( solve（不等式组，变元组） )16 画图方法：先产生横坐标的取值和相应的纵坐标的取值，然后执行命令： plot(x,y) 方法2：fplot(f(x)，xmin,xmax) fplot(f(x)，xmin,xmax,ymin,ymax) 方法3：ezplot(f(x) ezplot(f(x) ,xmin,xmax) ezplot(f(x) ,xmin,xmax,ymin,ymax) 17 求极限（1）极限：syms x limit(f(x), x, a) （2）单侧极限：左极限：syms x limit(f(x), x, a,left)右极限：syms x limit(f(x), x, a,right) 18 求导数diff(f(x) diff(f(x),x) 或者：syms x diff(f(x) syms x diff(f(x), x) 19 求高阶导数 diff(f(x)，n) diff(f(x),x，n)或者：