通信系统课程设计之基于MATLAB的FM通信系统

西南科技大学课程设计报告课程名称： _______ 通信系统课程设计_________设计名称： ______ 基于MATLAB勺FM通信系统设计___________ 姓名：容晓庆____________________________学号：20096025 ________________________班级： ____________ 通信0901班___________________________ 指导教师： __________ 侯宝林_________________________起止日期：2012617-2012625 ____________西南科技大学信息工程学院课程设计任务书学生班级：通信0901班学生姓名：容晓庆学号：设计名称：基于MATLAB的FM通信系统仿真起止日期：2012.06.10-2012.06.25 指导教师：侯宝林课程设计评语表指导教师评语:成绩: 指导教师:基于MATLAB的FM通信系统仿真「、设计目的和意义(1) 熟悉MATLAB文件中M文件的使用方法，包括函数、原理和方法的应用。

(2) 加深对FM信号调制原理的理解。

(3) 增强在通信原理仿真方面的动手能力与自学能力。

(4) 在做完FM调制仿真之后，在今后遇到类似的问题，学会对所面对的问题进行系统的分析，并能从多个层面进行比较、设计原理图1模拟通信系统模型调制器：使信号与信道相匹配，便于频分复用等。

发滤波器：滤除调制器输出的无用信号。

收滤波器：滤除信号频带以外的噪声，一般设N(t)为高斯白噪声，则Ni(t)为窄带白噪声。

在通信系统中一般需要将信号进行相应调制，以利于信号在信道上的传输，调制是将用原始信号去控制高频振荡信号的某一参数，使之随原始信号的变化而成规律变化。

调制可分为线性调制和非线性调制。

线性调制有AM、DSB等，非线性调制有FM、PM等，这里主要讨论FM调制通信系统1. FM调制原理角调制不是线性调制，角调制中已调信号和调制信号频谱之间不是线性关系而是产生出新的与频谱搬移不同的新的频率分量，呈现非线性特性，故又成为非线性调制。

FM调制中瞬时角频率是关于调制信号的线性函数，瞬时角频率偏移量w k FM f (t),则，瞬时角频率为w w c k FM f (t)k FM为频偏指数则，调频信号为S FM (t)cos(w c t k FM f(t)dt)当调制信号是单频余弦时，调制信号为S FM(t) Acos(w c t k FM A m COSW m tdt) Acos(w c t FM sin w m t)，调制信号的信号的产生图2m(t) | ⑺dt|[PM S FM (t)图3图2为调频信号的直接产生，图3为调频信号的间接产生，在间接产生方法中，因为实 际调相的相位调节器的范围在(-,)，故而间接调频只能用于相位偏移和频率偏移 不大的窄调制情况，而直接调频则常用于宽带调制情况。

窄带角频率调制时，最大瞬时相位偏移远小于 30°即 k FM f (t)dt max /6调频信号的时域表达式S FM(t) COS (W c t k FMf(t)dt)，当为窄带调制时，S Z M (t) Acosw c t [Ak FM f (t)dt]sinw c tF(W),而且假设飞f (t )的平均值为零，有傅氏变换可得2. 解调原理 1)非相干解调由于调频信号的瞬时频率正比于调制信号的幅度，因而调制信号的调节器必须产生正比 于输入频率的输出电压，即输入为调制信号为S NBFM (w)A[ (w W c )此时调频信号的带宽为调制信号的两倍 其相应图形如下:-叱+ %Aht .% -馮『叱些亠％：调频信号频谱假设调制线号的频谱为 (W W c )]W c )]信号频谱构成鉴频器。

2)相干解调由于窄带调频信号可分解为同相分量和正交分量，因而可以用线性调制中的相干解调 来进行解调。

原理图如下：n i(t )S NBFM (t)------- » 带通滤波Si(t). 相乘器Sp(t) 微分器低通滤波So(t)-------►-------►------------- .带通滤波器的作用是抑制信号带宽以外的噪声，低通滤波器输出 S o (t) Ak FM f(t)/2 ；上图所示的相干解调法只适用与窄带调频3. 调频系统的抗噪声性能主要讨论非相干解调的抗噪声性能，已知输入信号为 S FM (t) cos(w c t k FM f (t)dt)，输况下，解调器的输出端的信噪比为S o 3A 2K ：m 2(t) N o 8 2n °f ；单频时，带宽B FM 2( f f m ),所以增益可化解为噪比增益和调频指数的三次方成正比。

加大调频指数 FM,可使调频系统的抗噪声性能 迅速改善。

二. 详细设计步骤1.产生调制信号和载波信号，画出相应的时域图和频域图入功率 S i A 2，输入噪声N in o B FM ，信噪比为Si/NiA* A/2n o B FM ,非相干解调的解调器的输入端是调频信号与噪声的叠加，为S(t) n i (t)S FM (t) n i (t)，在大信噪比情考虑单频余弦信号调制，故可得大信噪比情况下的信噪比增益为G FM50 N o 51N i2m fA 2 2n o fmG FM 3 FM ( FM 1).由此可看出，性原理图如下:包络检波输出为So(t)= k FM k d f(t) , k d 为鉴频器的灵敏度，微分器和包络检波最简单的解调器是鉴频器片：二七FH 寸砂凹1 1 11. . . . J 111■L. . . .30002500调制信号的频谱图00O 2 00 5 00 00 6 A)luruoadso60o40002oo4—■oo -20030-500cen u载波信号时域囲2.利用modulate函数产生调制信号说频信号的频谱圏0 _______________________________ ^00^400 ^300 -200 -100 0 100 200 300 400 500Frequence (Hz)3.对调制信号进行解调无噪声条件下解调信号的时域圏1 --------------------------------- 1 --------------------------------- 1 --------------0 8 I0 0.5 1 1.5250200150100504.对信道中加入加性高斯白噪声，绘制调频信号的时域图和频域图O 05□ 1F r uqucji ic}Its)O 1 5(rw m s ^s d5•绘制加噪声的调频信号的解调时域图和频域图吋冋t设计结果及分析1. 对已调制未加噪声的调频信号而言，在时域图中出现疏密交错的变化规律，有理论和实际都可以看到，在瞬时最大正频偏处，波形最密；在瞬时最低处，波形最疏..2. 由图可知，输入的调制信号通过调制之后，波形发生了明显的变化，原本规则的正弦信号变成了不规则的上下起伏波动的图形，而且调制后的图形也没有原本正弦信号般圆滑，出现了十分尖锐的突起。

说明正弦信号通过FM调制之后波形发生了明显的改变。

通过频谱图的对照比较我们可以看出FM调制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形，调制后的波形与调制前的完全不同，这证明FM 调制并不是线性的，而是非线性的.3. 改变抽样频率，可以看到当抽样频率低于调制信号的时，无法得出信号波形，由奈奎斯特定义可以知道，当抽样频率低于调制信号时，信号会混叠失真•4. 加入噪声后时域波形与原来的时域波形相比，波形明显失真，波形不仅不如原本波形般规则，而且曲线之间还出现了未链接在一起的断裂，但随着信噪比的增大，与原有的波形的相似度也增大了，说信噪比越大，噪声对信号的影响也变小了。

从加入噪声的图形与未加入噪声的对比中我们还可以看出噪声对时域图的变化明显比频域图的变化更为突出，白噪声在整个频谱内每个频点的能量为常数，且基本恒定，所以他对于时域的影响更大。

三、体会此次课程设计让自己对FM 调制及MATLAB 软件有了较深了解，设计时遇到许多难题，但通过自己努力，以及与同学一起探讨，以及在老师的帮助下，逐步对自己的课题有了更加深入的了解，刚开始对自己所要做的FM 调制原理没有搞清楚，书籍知识也只理解到皮毛而且对MATLAB 软件的使用和其语言的使用也不是很会，虽然在数字信号处理过程中学过，但仅仅停留在表面上，只会简单的调试，对其未进行深入学习，这也增加了此次设计的对自己的难度，但通过相关书籍的阅读和与老师的交流，在老师，同学的帮助下，完成了课程设计的基本要求。

在课程设计的这段时间内主要遇到了以下问题：1、对MATLAB 中的函数的理解不够，在信号的调制时，一开始并不知道要用到modulate 函数来对未调制信号进行调制，还很辛苦的在在网上找了一大堆程序来完成调制过程，后来经过同学的提醒，简化了程序，简单的完成了对信号的调制。

2、在做时域到频域的变化步骤时，由于要利用傅立叶函数进行变化，刚开始完全不知如何下手，很多书籍上有DFT，IFFT 函数，不能明白其意思，后来通过老师的指点，使用FFT可以使时域直接到频域进行变化，从而完成时域到频域的转换。

3、加噪声是对系统抗噪声性能的一种测试，用awgn 这个函数时，不知是什么原因，系统总是报错，后来经过对软件的重新启动又可以执行相关功能了，软件的性能也很重要呀。

4、在此次课程设计中遇到的棘手问题只要在MATLAB 软件语言上的使用，虽然对FM 调制有了一定理解，但却不能将其用相应的语言来表示出FM 调制过程，在程序方面主要还是通过网络方面收集，并对其进行整理和分析，综合成相关的程序要求，简单的完成了这次的课程设计。

四、参考文献[1] 曹志刚钱亚生《现代通信原理》. 清华大学出版社[2] 程培青《数字信号处理》清华大学出版社[3] 聂翔飞、王海宝《MATLAB 程序设计及其在信号处理中的应用》大学出西南交通版社[4] 沈伟慈《通信电路》西安电子科技大学出版社[5] 百度文库等网上资源.title （' 调制信号时域图 '）;fc=50;% ct=cos(2*pi*fc*t);% figure(2) plot(t,ct) ; % xlabel('t')title （' 载波信号时域图 '）; fs=1000; % sfm=modulate(mt,fc,fs,'FM');ts=0.001;% fs=1/ts;% df=0.25; %设定载波频率生成载波 设定调频指数调制，产生已调信号 抽样间隔 抽样频率所需的频率分辨率， 用在求傅里叶变换%寸，它表示FFT 的最小频率间隔%***** 对调制信号 m(t) 求傅里叶变换 *****mt=am*cos(2*pi*fm*t); % fs=1/ts; if nargin==2 n1=0; elsen1=fs/df; endn2=length(mt);n=2A (max( nextpow2( n1), nextpow2( n2))); M=fft(mt,n);mt=[mt,zeros(1,n-n2)];df1=fs/n;%里变换M=M/fs;%f=[0:df1:df1*(length(mt)-1)]-fs/2; %figure(11)以上程序是对调制后的信号 u 求傅 缩放，便于在频铺图上整体观察时间向量对应的频率向量........ % 调制信号的产生 .......dt=0.001; %t=0:dt:1.5; %am=5; % fm=5; % echo offclose all clear all clcmt=am*cos(2*pi*fm*t); % 生成调制信号 figure(1) plot(t,mt); xlabel('t'); 设定时间步长 产生时间向量 设定调制信号幅度 设定调制信号频率原调信号plot(f,abs(fftshift(M))) %fftshift: 将FFT 中的DC分量移到频谱中心xlabel(' 频率f') title(' 调制信号的频谱图');figure(3)plot(t,sfm);xlabel('t')title(' 调频信号的时域图'); axis([0 0.15 -1 1]);b=fft(sfm,1024); f=(0:length(b)-1)*fs/length(b) -fs/2;figure(4)plot(f,abs(fftshift(b))); %FM 信号频谱图xlabel('Frequence(Hz)'); % 频率分量ylabel('Power Spectrum(y)'); % 频率谱能量title(' 调频信号的频谱图');..... % 调频信号的解调.......nsfm=sfm;for i=1:length(t)-1 % 接受信号通过微分器处理diff_nsfm(i)=(nsfm(i+1)-nsfm(i))./dt;enddiff_nsfmn = abs(hilbert(diff_nsfm)); %hilbert 变换，求绝对值得到瞬时幅度(包络检波)zero=(max(diff_nsfmn)-min(diff_nsfmn))/2; diff_nsfmn1=diff_nsfmn-zero;figure(5)..... % 绘制无噪声条件下解调信号的时域图plot((1:length(diff_nsfmn1))./1000,diff_nsfmn1./400,'r');xlabel(' 时间t');title(' 无噪声条件下解调信号的时域图'); %对已调信号进行傅氏变换x=fft(diff_nsfmn1,1024); f=(0:length(x)-1)*fs/length(x) -fs/2;x=x/1000; figure(6) plot(f,abs(fftshift(x)));xlabel(' 频率f')title(' 已解调信号的频谱图') %有噪声情况下的调频信号fc = 50; % 载波频率fs = 1000; % 采样频率t = (0 :0.001:0.15); % 时间区域am=5;fm=5;mt=am*cos(2*pi*fm*t) %输入信号 sfm = modulate(mt,fc,fs,'FM'); %调制信号 y1 = sfm + awgn(sfm,10,0); % FM信号加入噪声 figure(7)plot(t,y1); %FM信号时域图 xlabel('t(s)');ylabel('sfm');title(' 加噪声之后的调频信号时域图 , 性噪比为 10') axis([0 0.15 -1 1]); b=fft(y1,1024);f=(0:length(b)-1)*fs/length(b) -fs/2; figure(8)plot(f,abs(fftshift(b))); %FM 信号频谱图 xlabel('Frequence(Hz)'); % 频率分量 ylabel('PowerSpectrum(y)'); % 频率谱能量title(' 加噪声后的调频信号频谱图 ');%对加性噪声调频信号进行解调 ny1=y1;for i=1:length(t)-1 % diff_ny1(i)=(ny1(i+1)-ny1(i))./dt; enddiff_ny1n = abs(hilbert(diff_ny1)); %hilbert 络检波) zero=(max(diff_ny1n)-min(diff_ny1n))/2; diff_ny1n1=diff_ny1n-zero; figure(9).... % 绘制无噪声条件下解调信号的时域图plot((1:length(diff_ny1n1))./1000,diff_ny1n1./400,'r'); xlabel(' 时间 t');title(' 加性噪声条件下解调信号的时域图 '); %对已调信号进行傅氏变换 x=fft(diff_ny1n1,1024); f=(0:length(x)-1)*fs/length(x) -fs/2; x=x/1000; figure(10) plot(f,abs(fftshift(x))); xlabel(' 频率 f');title(' 加噪声的已解调信号的频谱图 ');接受信号通过微分器处理变换，求绝对值得到瞬时幅度（包。