《课程设计（论文）基于Matlab的FM信号分析.doc》由会员分享，可在线阅读，更多相关《课程设计（论文）基于Matlab的FM信号分析.doc（25页珍藏版）》请在三一办公上搜索。
1、目 录摘 要IAbstractII1. FM调制系统的介绍及设计11.1 通信系统的简介11.2 FM调制模型的建立21.3 FM调制仿真结果42. FM解调系统的介绍及设计72.1 FM解调模型的建立72.2 解调过程分析72.3 FM解调的部分程序及波形83. 高斯白噪声信道特性114. 调频系统的抗噪声性能分析145. MATLAB源代码196心得体会217参考文献22摘 要调制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的影响，调制方式往往决定
2、了一个通信系统的性能。FM信号的调制属于频谱的非线性搬移，它的解调也有相干和非相干解调两种方式。本课程设计使用的仿真软件为MATLAB 7.0，利用MATLAB集成环境下的M文件，编写程序来实现FM信号的仿真分析，并分别绘制出基带信号、载波信号、已调信号的时域波形和频域波形；再进一步分别绘制出对已调信号叠加噪声后信号、同步解调前信号和解调后基带信号的时域波形；最后绘出FM基带信号调制和解调系统后的输入输出信噪比的关系，并通过与理论结果波形对比来分析该仿真调制与解调过程的正确性及噪声对FM信号解调的影响。在课程设计中，系统开发平台为Windows Xp，在该平台运行MATLAB程序完成了对FM调
3、制和解调以及对叠加噪声后解调结果的观察和分析。关键词 FM；调制；解调；MATLAB 7.0；噪声AbstractModulation in a communication system has extremely important role. Through the modulation, not only can frequency shift, the modulation signal spectrum is moved to want, thus will position modulation signal is converted into suitable for trans
4、mission of the signal, and it has attune to the transmission effectiveness and reliability of transferring having very big effect, modulation method often determines a communication system performance. FM signals modulation of nonlinear shift, belong to the spectrum of its demodulation also have coh
5、erent and incoherent demodulation in two ways. The course is designed to use simulation software for MATLAB 7.0, use of MATLAB integration environment of M files, write a program to realize FM signals of the simulation analysis, and draw the baseband signal respectively, carrier signal, already adju
6、stable signal and the time domain waveform frequency domain waveform; Further respectively to map out already adjustable signal after adding noise signal, synchronous demodulation signal demodulation before and after baseband signal temporal profile, Finally draw FM baseband signal modulation and de
7、modulation system after the input/output SNR, and the relationship with the theoretical results wave contrast to analyze the simulation modulation and demodulation process accuracy and noise on FM signal demodulation of the influence. In the curriculum design, system developing platform for Windows
8、Xp, on this platform operation MATLAB finished to FM modulation and demodulation of adding noise and after the observation and analysis results demodulation. Key words ：FM, Modulation, Demodulation, MATLAB 7.0，noise1. FM调制系统的介绍及设计1.1 通信系统的简介通信的目的是传输信息。通信系统的作用就是将信息从信息源发送到一个或多个目的地。通信系统对信号进行两种基本变换： 第一、
9、要把发送的消息要变换成原始电信号。 第二、将原始电信号调制到频率较高的载频上，使其频带适合信道的传输。 调制前和解调后的信号称为基带信号，已调信号也称为频带信号。 对于任何一个通信系统，均可视为由发送端、信道和接收端三大部分组成（如图2-1所示）。信息源发送设备信 道接受设备信息源噪声源发送端接收端信道图1-1 通信系统一般模型信息源（简称信源）的作用是把各种信息转换成原始信号。根据消息的种类不同信源分为模拟信源和数字信源。发送设备的作用产生适合传输的信号，即使发送信号的特性和信道特性相匹配，具有抗噪声的能力，并且具有足够的功率满足原距离传输的需求。信息源和发送设备统称为发送端。发送端将信息直
10、接转换得到的较低频率的原始电信号称为基带信号。通常基带信号不宜直接在信道中传输。因此，在通信系统的发送端需将基带信号的频谱搬移（调制）到适合信道传输的频率范围内进行传输。这就是调制的过程。信号通过信道传输后，具有将信号放大和反变换功能的接收端将已调制的信号搬移（解调）到原来的频率范围，这就是解调的过程。信号在信道中传输的过程总会受到噪声的干扰，通信系统中没有传输信号时也有噪声，噪声永远存在于通信系统中。由于这样的噪声是叠加在信号上的，所以有时将其称为加性噪声。噪声对于信号的传输是有害的，它能使模拟信号失真。在本仿真的过程中我们假设信道为高斯白噪声信道。模拟系统框图如下图1-2所示： 图1-2
11、模拟系统框图调制在通信系统中具有十分重要的作用。一方面，通过调制可以把基带信号的频谱搬移到所希望的位置上去，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号。另一方面，通过调制可以提高信号通过信道传输时的抗干扰能力，同时，它还和传输效率有关。不同的调制方式产生的已调信号的带宽不同，因此调制影响传输带宽的利用率。可见，调制方式往往决定一个通信系统的性能。在本仿真的过程中我们选择用调频调制方法进行调制。在本仿真的过程中我们选择用同步解调方法进行解调。1.2 FM调制模型的建立为基带调制信号，设调制信号为 m(t)=cos(2*pi*fm*t) 公式(2-1) 设正弦载波为 c(t)=
12、cos(2*pi*fc*t) 公式(2-2)信号传输信道为高斯白噪声信道，其功率为。在调制时，调制信号的频率去控制载波的频率的变化，载波的瞬时频偏随调制信号成正比例变化，即: 公式(2-3) 式中，为调频灵敏度。这时相位偏移为 公式(2-4) 公式(2-5)则可得到调频信号为 公式(2-6) FM 的频谱的计算： 公式(2-7) 公式(2-8) 公式(2-9)可以看出FM 的频谱与的值有关。其信号带宽为 公式(2-10)FM 的频谱理论值无穷大，但可根据调频指数分为宽带调频和窄带调频。在MATLAB程序中利用自建的积分函数实现对调制信号的间接调频，部分程序如下所示：%*FM调制*t0=0; t
13、1=0.01; %调制信号的时域范围 fs=100000; %抽样频率 ts=1/fs; %采样率t=t0:ts:t1;fc=30000; %设定载波频率x=cos(2*pi*fc*t); %生成载波fm=300; %设定调制信号频率mf=0.5; %设定调频指数kf=2*pi*fm*mf;f=fm*mf; %频偏最大值%*正弦波调制信号* m=cos(2*pi*fm*t); %生成调制信号% %*三角波调制信号*%m=sawtooth(2*pi*fm*t,0.5);u=GetFMSignal(m,t0,t1,ts,fc,kf); %已调制信号function u=GetFMSignal(m,
14、t0,t1,ts,fc,kf)%由调制信号和载波，产生已调制的FM信号%因为是调频，载波幅度设置为1，即为cos(2*pi*fc*t)t=t0:ts:t1;int_m=m_INT(m,t0,t1,ts);u=cos(2*pi*fc*t+kf*int_m);1.3 FM调制仿真结果在调频指数不同的情况下的调制波形的时域和频域图如下列图所示： 图1-3 时的各波形图 图1-4 时的各波形图图1-5 时的各波形图三角波调制时的波形与三角函数一样，下面只给出了在调频指数为3时的各波形图。图1-6 时的各波形图2. FM解调系统的介绍及设计2.1 FM解调模型的建立调制信号的解调分为相干解调和非相干解调
15、两种。相干解调即同步解调，仅仅适用于窄带调频信号，且需同步信号，故应用范围受限；而非相干解调不需同步信号，且对于NBFM信号和WBFM信号均适用，因此是FM系统的主要解调方式。但在本仿真的过程中我们对窄带信号进行调制与解调，选择用同步解调方法进行解调。图2-1 FM同步解调解调模型限幅器输入为已调频信号和噪声，限幅器是为了消除接收信号在幅度上可能出现的畸变；带通滤波器的作用是用来限制带外噪声，使调频信号顺利通过。鉴频器中的乘法器把调频信号与相干载波相乘，然后由低通滤波器和微分器取出调制信号。2.2 解调过程分析由上述公式(2-6)知道输入调频信号为 设相干载波为c(t)=cos(2*pi*fc
16、*t) 公式(2-1)乘法器的作用是把调频信号变成有多种频率的波的混合，乘法器输出为 公式(2-2)经低通滤波器后取出器低频分量为 公式(2-3)在经过微分器，即得出解调出的基带信号： 公式(2-4)相干解调可以恢复出原来的基带信号，而且要求本地载波与调制载波同步，否则会使解调信号失真。2.3 FM解调的部分程序及波形 在这次的解调设计中，采用MATLAB自带的解调函数fmdemod对已调信号进行解调。该函数使用形式如下所示：Z = FMDEMOD(Y,Fc,Fs,FREQDEV)其中，Y为输入的调频波函数，Fc、Fs分别为载波频率和抽样频率，而且要满足，FREQDEV为最大频偏。各种情况下的
17、解调波形如下所示：图2-2 时的各波形图图2-3 时的各波形图图2-4 时的各波形图图2-5 时的各波形图图2-6 时的各波形图3. 高斯白噪声信道特性设正弦波通过加性高斯白噪声信道后的信号为 公式（3-1） 其中，白噪声的取值的概率分布服从高斯分布。MATLAB本身自带了标准高斯分布的内部函数awgn。awgn函数产生的随机序列服从均值为，方差的高斯分布。正弦波通过加性高斯白噪声信道后的信号为 公式（3-2）故其有用信号功率为 公式（3-3） 噪声功率为 公式（3-4） 公式（3-5） 则可得到公式 公式（3-6）我们可以通过这个公式方便的设置高斯白噪声的方差。在本仿真过程中，我们选择了5d
18、b、10db、15db、20db、25db五种不同信噪比以示区别，并绘制输入输出信噪比关系曲线。图3-1 三角函数的输入输出信噪比曲线图3-2 三角波的输入输出信噪比曲线4. 调频系统的抗噪声性能分析从前面的分析可知，调频信号的解调有相干解调和非相干解调两种。相干解调仅适用于窄带调频信号，且需同步信号；而非相干解调适用于窄带和宽带调频信号，而且不需同步信号，因而是FM系统的主要解调方式，所以这里仅仅讨论非相干解调系统的抗噪声性能，其分析模型如图9所示。 图4-1 调频系统抗噪声性能分析模型图中带通滤波器的作用是抑制信号带宽以外的噪声。是均值为零，单边功率谱密度为的高斯白噪声，经过带通滤波器后变
19、为窄带高斯噪声 。限幅器是为了消除接收信号在幅度上可能出现的畸变。设调频信号为 故其输入功率为 公式（4-1）输入噪声功率为 公式（4-2）因此输入信噪比为 公式（4-3）在大信噪比条件下，信号和噪声的相互作用可以忽略，这时可以把信号和噪声分开来算，这里，我们可以得到解调器的输出信噪比 公式（4-4）上式中，为载波的振幅，为调频器灵敏度，为调制信号的最高频率，为噪声单边功率谱密度。我们如若考虑为单一频率余弦波时的情况，可得到解调器的制度增益为 公式（4-5）考虑在宽带调频时，信号带宽为 公式（4-6）则可以得到 公式（4-7）可以看出，大信噪比时宽带调频系统的信噪比增益是很高的，它与调频指数的
20、立方成正比。可见，加大调频指数，可使调频系统的抗噪声性能迅速改善。5. MATLAB源代码echo offclose allclear allclc%*FM调制*t0=0; t1=0.01; %调制信号的时域范围 fs=100000; %抽样频率 ts=1/fs;t=t0:ts:t1;fc=30000; %设定载波频率x=cos(2*pi*fc*t); %生成载波fm=300; %设定调制信号频率% mf=0.5; %设定调频指数% mf=1;mf=3;kf=2*pi*fm*mf;f=fm*mf; %频偏最大值%*正弦波调制信号* m=cos(2*pi*fm*t); %生成调制信号% %*三角
21、波调制信号*% m=sawtooth(2*pi*fm*t,0.5);function u=GetFMSignal(m,t0,t1,ts,fc,kf)%由调制信号和载波，产生已调制的FM信号%因为是调频，载波幅度设置为1，即为cos(2*pi*fc*t)t=t0:ts:t1; int_m=m_INT(m,t0,t1,ts);u=cos(2*pi*fc*t+kf*int_m);function M,m_new,df_new=GetFFT(m,t0,t1,ts,df)%对信号求FFT，先得出FFT的点数t=t0:ts:t1;fs=1/ts;n1=length(t); if nargin=4 n2=0
22、; else n2=fs/df; endn=2(max(nextpow2(n1),nextpow2(n2);M=fft(m,n); %m的n点FFTm_new=m,zeros(1,n-length(t); %新的m，有可能对m进行扩展df_new=fs/n; %新的频率分辨率u=GetFMSignal(m,t0,t1,ts,fc,kf); %已调制信号%*FM解调*yn0=awgn(u,5); %加入高斯白噪声yn1=awgn(u,10);yn2=awgn(u,15); yn3=awgn(u,20); yn4=awgn(u,25); m_dem_new=fmdemod(u,fc,fs,f);
23、yyn0=fmdemod(yn0,fc,fs,f); yyn1=fmdemod(yn1,fc,fs,f); yyn2=fmdemod(yn2,fc,fs,f); yyn3=fmdemod(yn3,fc,fs,f); yyn4=fmdemod(yn4,fc,fs,f);dy10=yn0-u; %高斯白噪声snr10=var(u)/var(dy10); %输入信噪比dy20=yyn0-m_dem_new; %解调后噪声snr20=var(m_dem_new)/var(dy20); %输出信噪比 dy11=yn1-u; snr11=var(u)/var(dy11); dy21=yyn1-m_dem_
24、new; snr21=var(m_dem_new)/var(dy21); dy12=yn2-u; snr12=var(u)/var(dy12); dy22=yyn2-m_dem_new; snr22=var(m_dem_new)/var(dy22);dy13=yn3-u; snr13=var(u)/var(dy13); dy23=yyn3-m_dem_new; snr23=var(m_dem_new)/var(dy23);dy14=yn4-u; snr14=var(u)/var(dy14); dy24=yyn4-m_dem_new; snr24=var(m_dem_new)/var(dy24)
25、; in=snr10,snr11,snr12,snr13,snr14; out=snr20,snr21,snr22,snr23,snr24; %*求傅里叶变换* df=0.25; %FFT时的频率分辨率f=0:df:df*(length(t)-1)-fs/2; %与t对应的频率向量M,m_new,df_new_m=GetFFT(m,t0,t1,ts,df);U,u_new,df_new_u=GetFFT(u,t0,t1,ts,df);f1=0:df_new_m:df_new_m*(length(m_new)-1)-fs/2;f2=0:df_new_u:df_new_u*(length(u_ne
26、w)-1)-fs/2;%画调制信号和已调制信号的波形disp(按任意键可以看到原调制信号和已调信号在时域和频域内的图形)pausefigure(1);subplot(2,2,1);plot(t,m);title(基带信号);subplot(2,2,2);plot(t,u);title(已调制信号);axis(0 0.001 -1 1)subplot(2,2,3);plot(f1,abs(fftshift(M);xlabel(基带信号频谱);% axis(0 5000 0 5000)subplot(2,2,4);plot(f2,abs(fftshift(U);xlabel(已调制信号频谱);%
27、axis(0 5000 0 10000)%绘制加噪前后的解调信号的波形disp(按任意键可以看到加噪前后的原基带信号和解调信号在时域和频域内的图形)pausefigure(2);subplot(3,2,1);plot(t,m);xlabel(原基带信号);subplot(3,2,2);plot(t,m_dem_new);xlabel(解调后基带信号);axis(0 0.01 -1 1)subplot(3,2,3);plot(t,yn0,r);xlabel(加噪后基带信号);axis(0 0.001 -1 1)subplot(3,2,4);plot(t,yyn0,r);xlabel(加噪解调后基
28、带信号);m_dem_new,m_dem_new_new,df_new_m_dem_new=GetFFT(m_dem_new,t0,t1,ts,df);f3=0:df_new_m_dem_new:df_new_m_dem_new*(length(m_dem_new_new)-1)-fs/2;subplot(3,2,5);plot(f1,abs(fftshift(M);xlabel(基带信号频谱);subplot(3,2,6);plot(f1,abs(fftshift(M);xlabel(无噪声解调后基带信号频谱);%绘制输入输出信噪比曲线disp(按任意键可以看到输入输出信噪比关系曲线的图形)
29、pausefigure(3);plot(in,out,*)hold onplot(in,out)xlabel(输入信噪比)ylabel(输出信噪比) function int_m=m_INT(m,t0,t1,ts);% 计算离散信号的积分t=t0:ts:t1;int_m(1)=0;for j=1:length(t)-1 int_m(j+1)=int_m(j)+m(j)*ts;end6心得体会我这次所做的课程设计是基于MATLAB的FM调制实现，虽然在制作此课程设计中遇到许多的困难，但主要还是自己对所学的知识并没有完全的消化，以至对FM原理都是似懂非懂，浪费了较多的时间，另外通过此次的课程设计，对MATLAB软件的使用也有了更深一步的了解与掌握，特别是如何对基带信号如何进行调制，时域到频域变换，及加入噪声方面有了深入的了解。通过此次课程设计不仅对原来通信原理知识得到了巩固，而且还学会了许多原来不会的东西，所以这次的课程设计对我的帮助很大。7.参考文献1通信原理（第六版） 樊昌信等。 北京：国防工业出版社。2MATLAB7.0在数字信号处理中的应用 罗军辉等。北京：机械工业出版社。3MATLAB程序设计教程 刘卫国等。 北京：中国水利水电出版社。4 MATLAB及在电子信息类课程中的应用唐向宏等编著。北京：电子工业出版社
