FM的调制解调

1、目目 录录 目目 录录.2 摘摘 要要.3 ABSTRACT.4 1.引言引言.5 2. FM 调制系统的介绍及设计调制系统的介绍及设计.6 2.1 通信系统的简介.6 2.2 FM 调制模型的建立.8 2.3 FM 调制仿真结果.11 3. FM 解调系统的介绍及设计解调系统的介绍及设计.13 3.1 FM 解调模型的建立.13 3.2 解调过程分析.14 3.3 FM 解调的部分程序及波形.15 4. 高斯白噪声信道特性高斯白噪声信道特性.18 5. 调频系统的抗噪声性能分析调频系统的抗噪声性能分析.20 6. MATLAB 源代码源代码 .22 7心得体会心得体会.27 8参考文献参考文

2、献.28 摘摘 要要 调制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬 移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合 于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的 影响，调制方式往往决定了一个通信系统的性能。 FM 信号的调制属于频谱的 非线性搬移，它的解调也有相干和非相干解调两种方式。 本课程设计使用的仿真软件为 MATLAB 7.0，利用 MATLAB 集成环境下的 M 文件，编写程序来实现 FM 信号的仿真分析，并分别绘制出基带信号、载波信号、 已调信号的时域波形和频域波形；再进一步分别绘制出对已调信号叠加噪声后 信号、同步解调前

3、信号和解调后基带信号的时域波形；最后绘出 FM 基带信号调 制和解调系统后的输入输出信噪比的关系，并通过与理论结果波形对比来分析 该仿真调制与解调过程的正确性及噪声对 FM 信号解调的影响。在课程设计中， 系统开发平台为 Windows Xp，在该平台运行 MATLAB 程序完成了对 FM 调制和解 调以及对叠加噪声后解调结果的观察和分析。 关键词 FM；调制；解调；MATLAB 7.0；噪声 AbstractAbstract Modulation in a communication system has extremely important role. Through the modul

4、ation, not only can frequency shift, the modulation signal spectrum is moved to want, thus will position modulation signal is converted into suitable for transmission of the signal, and it has attune to the transmission effectiveness and reliability of transferring having very big effect, modulation

5、 method often determines a communication system performance. FM signals modulation of nonlinear shift, belong to the spectrum of its demodulation also have coherent and incoherent demodulation in two ways. The course is designed to use simulation software for MATLAB 7.0, use of MATLAB integration en

6、vironment of M files, write a program to realize FM signals of the simulation analysis, and draw the baseband signal respectively, carrier signal, already adjustable signal and the time domain waveform frequency domain waveform; Further respectively to map out already adjustable signal after adding

7、noise signal, synchronous demodulation signal demodulation before and after baseband signal temporal profile, Finally draw FM baseband signal modulation and demodulation system after the input/output SNR, and the relationship with the theoretical results wave contrast to analyze the simulation modul

8、ation and demodulation process accuracy and noise on FM signal demodulation of the influence. In the curriculum design, system developing platform for Windows Xp, on this platform operation MATLAB finished to FM modulation and demodulation of adding noise and after the observation and analysis resul

9、ts demodulation. Key words ：FM, Modulation, Demodulation, MATLAB 7.0，noise 1.1.引言引言 在通信系统中，调制过程是一个频谱搬移的过程，它是将低频信号的频谱 搬移到载频位置。解调是调制的逆过程，即是将已调制的信号还原成原始基带 信号的过程。信号的接收端就是通过解调来还原已调制信号从而读取发送端发 送的信息。因此信号的解调对系统的传输有效性和传输可靠性有着很大的影响。 角度及解调电路不同于频谱线性搬移电路。它是用低频信号去调制高频振 荡的相角，或是从已调波中解出调制信号所进行的频谱变换，这种变换不是线 性变换，而是非线性

10、变换。因此，我们把角度调制及调角波的解调电路称为频 谱非线性变换电路。FM 在通信系统中的使用非常广泛。FM 广泛应用于高保真音 乐广播、电视伴音信号的传输、卫星通信和蜂窝电话系统等。本次课设中使用 功能强大的 MATLAB 仿真软件对 FM 信号进行仿真分析。 MATLAB 的名称源自 Matrix Laboratory，它的首创者是在数值线性代数领 域颇有影响的 Cleve Moler 博士，他也是生产经营 MATLAB 产品的美国 Mathworks 公司的创始人之一。 MATLAB 是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB 将高性 能的数值计算和可视化集成在一起，并提

11、供了大量的内置函数，从而使其被广 泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作中。 MATLAB 与其它计算机高级语言如 C，C+等相比，MATLAB 语言编程要简洁得多， 编程语句更加接近数学描述，可读性好，其强大的图形功能和可视化数据处理 能力也是其它高级语言望尘莫及的。对于具有任何一门高级语言基础的学生来 说，学习 MATLAB 十分容易。 MATLAB 使得人们摆脱了常规计算机编程的繁琐，让人们能够将大部分精力 投入到研究问题的数学建模上。可以说，应用 MATLAB 这一数学计算和系统仿真 的强大工具，可以使科学研究的效率得以成百倍的提高。 2.2. FMFM 调

12、制系统的介绍及设计调制系统的介绍及设计 2.1 通信系统的简介通信系统的简介 通信的目的是传输信息。通信系统的作用就是将信息从信息源发送到一个 或多个目的地。通信系统对信号进行两种基本变换： 第一、要把发送的消息要变换成原始电信号。 第二、将原始电信号调制到频率较高的载频上，使其频带适合信道的传 输。 调制前和解调后的信号称为基带信号，已调信号也称为频带信号。 对于任何一个通信系统，均可视为由发送端、信道和接收端三大部分组成 （如图 2-1 所示） 。 信息源发送设备信 道接受设备信息源 噪声源 发送端发送端接收端接收端信道信道 图 2-1 通信系统一般模型 信息源（简称信源）的作用是把各种信

13、息转换成原始信号。根据消息的种 类不同信源分为模拟信源和数字信源。发送设备的作用产生适合传输的信号， 即使发送信号的特性和信道特性相匹配，具有抗噪声的能力，并且具有足够的 功率满足原距离传输的需求。 信息源和发送设备统称为发送端。 发送端将信息直接转换得到的较低频率的原始电信号称为基带信号。通常 基带信号不宜直接在信道中传输。因此，在通信系统的发送端需将基带信号的 频谱搬移（调制）到适合信道传输的频率范围内进行传输。这就是调制的过程。 信号通过信道传输后，具有将信号放大和反变换功能的接收端将已调制的 信号搬移（解调）到原来的频率范围，这就是解调的过程。 信号在信道中传输的过程总会受到噪声的干扰

14、，通信系统中没有传输信号 时也有噪声，噪声永远存在于通信系统中。由于这样的噪声是叠加在信号上的， 所以有时将其称为加性噪声。噪声对于信号的传输是有害的，它能使模拟信号 失真。在本仿真的过程中我们假设信道为高斯白噪声信道。 模拟系统框图如下图 2-2 所示： 图 2-2 模拟系统框图 调制在通信系统中具有十分重要的作用。一方面，通过调制可以把基带信 号的频谱搬移到所希望的位置上去，从而将调制信号转换成适合于信道传输或 便于信道多路复用的已调信号。另一方面，通过调制可以提高信号通过信道传 输时的抗干扰能力，同时，它还和传输效率有关。不同的调制方式产生的已调 信号的带宽不同，因此调制影响传输带宽的利

15、用率。可见，调制方式往往决定 一个通信系统的性能。在本仿真的过程中我们选择用调频调制方法进行调制。 角度调制及调角波的解调电路是频谱非线性变换电路，所以要采用与调幅电路 不同的方法来对基带信号进行调制和解调。在本仿真的过程中我们选择用同步 解调方法进行解调。 2.2 FM 调制模型的建立 从频率调制的相位与频率关系可以看出，调频信号可通过直接调频和间接 调频两种方法得到，所谓间接调频就是先对调制信号积分再调相而得到。同样， 6 max )(K dttm f 调相信号也可以通过直接调相和间接调相两种方法得到，间接调相就是先对调 制信号进行微分再进行频率调制。 根据调制后已调信号的瞬时相位偏移的大

16、小，可将角度调制分为宽带调制 (宽带调频和宽带调相)和窄带调制(窄带调频和窄带调相)。如果调频信号或调 相信号的最大瞬时相位偏移保持在很小的范围内，一般小于 30即满足条件： 公式(2-1) 时，则称为窄带调频或窄带调相。当上述条件不满足时，就称为宽带调频或宽 带调相。 图 2-3 FM 间接调制模型 其中，为基带调制信号，设调制信号为 ( )m t m(t)=cos(2*pi*fm*t) 公式(2-2) 设正弦载波为 c(t)=cos(2*pi*fc*t) 公式(2-3) 信号传输信道为高斯白噪声信道，其功率为。 2 在调制时，调制信号的频率去控制载波的频率的变化，载波的瞬时频偏随调制 信号

17、成正比例变化，即: ( )m t 公式(2-4) )()( 0 tukt f 式中，为调频灵敏度。 fK 这时相位偏移为 公式(2-5) dttuktt t f )()( 0 0 公式(2-6) 则可得到调频信号为 公式(2-6) )(cos()(u 0 0FM dttuktUt t fc FM 的频谱的计算： 公式(2-7) 公式(2-8) 公式(2-9) 可以看出 FM 的频谱与 )( fn mJ 的值有关。 其信号带宽为 公式(2-9) ) 1(2 f mB FM 的频谱理论值无穷大，但可根据调频指数分为宽带调频和窄带调频。 在 MATLAB 程序中利用自建的积分函数实现对调制信号的间接

18、调频，部分程 序如下所示： %*FM 调制* t0=0; t1=0.01; %调制信号的时域范围 fs=; %抽样频率 ts=1/fs; %采样率 t=t0:ts:t1; fc=30000; %设定载波频率 m m f m 1 20 420 2cos)(2)( 4cos)(22cos)(2)()sincos( n fnf ffff tnmJmJ tmJtmJmJtm 0 )12( 531 ) 12sin()(2 5sin)(23sin)(2sin)(2)sinsin( n fn ffff tnmJ tmJtmJtmJtm tnmJA tVk tAt c n fn f cFM )(cos)( s

19、in cos)(U _ 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); %生成调制信号 % %*三角波调制信号* % m=sawtooth(2*pi*fm*t,0.5); u=GetFMSignal(m,t0,t1,ts,fc,kf); %已调制信号 function u=GetFMSignal(m,t0,t1,ts,fc,kf) %由调制信号和载波，产生已调制的 FM 信号 %因

20、为是调频，载波幅度设置为 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); 2.3 FM 调制仿真结果 在调频指数不同的情况下的调制波形的时域和频域图如下列图所示： 00.0050.01 -1 -0.5 0 0.5 1 信 信 信 信 00.51 x 10 -3 -1 -0.5 0 0.5 1 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 50 100 150 200 250 信 信 信 信 信 信

21、信 图 2-4 时的各波形图3 f m 00.0050.01 -1 -0.5 0 0.5 1 信 信 信 信 00.51 x 10 -3 -1 -0.5 0 0.5 1 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 100 200 300 400 信 信 信 信 信 信 信 图 2-5 时的各波形图1 f m 00.0050.01 -1 -0.5 0 0.5 1 信 信 信 信 00.51 x 10 -3 -1 -0.5 0 0.5 1 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信

22、信 信 信 信 -505 x 10 4 0 100 200 300 400 500 信 信 信 信 信 信 信 图 2-6 时的各波形图5 . 0 f m 三角波调制时的波形与三角函数一样，下面只给出了在调频指数为 3 时的 各波形图。 00.0050.01 -1 -0.5 0 0.5 1 信 信 信 信 00.51 x 10 -3 -1 -0.5 0 0.5 1 信 信 信 信 信 -505 x 10 4 0 100 200 300 400 500 信 信 信 信 信 信 -505 x 10 4 0 100 200 300 信 信 信 信 信 信 信 图 2-7 时的各波形图3 f m 3.

23、3. FMFM 解调系统的介绍及设计解调系统的介绍及设计 3.1 FM 解调模型的建立 调制信号的解调分为相干解调和非相干解调两种。相干解调即同步解调，仅仅 适用于窄带调频信号，且需同步信号，故应用范围受限；而非相干解调不需同 步信号，且对于 NBFM 信号和 WBFM 信号均适用，因此是 FM 系统的主要解 调方式。但在本仿真的过程中我们对窄带信号进行调制与解调，选择用同步解 调方法进行解调。 图 3-1 FM 同步解调解调模型 限幅器输入为已调频信号和噪声，限幅器是为了消除接收信号在幅度上可能出 现的畸变；带通滤波器的作用是用来限制带外噪声，使调频信号顺利通过。鉴 频器中的乘法器把调频信号

24、与相干载波相乘，然后由低通滤波器和微分器取出 调制信号。 3.2 解调过程分析 由上述公式(2-6)知道输入调频信号为 )(cos()(u 0 0FM dttuktUt t fc 设相干载波为 c(t)=cos(2*pi*fc*t) 公式(3-1) 乘法器的作用是把调频信号变成有多种频率的波的混合，乘法器输出为 公式(3-2)2cos1 ()( 2 1 2sin 2 1 )(tdttmKtts cfcp 经低通滤波器后取出器低频分量为 公式(3-3)dttmKts fd )( 2 1 )( 在经过微分器，即得出解调出的基带信号： 公式(3-4)( 2 1 )( 0 tmKtm f 相干解调可以

25、恢复出原来的基带信号，而且要求本地载波与调制载波同步，否 则会使解调信号失真。 3.3 FM 解调的部分程序及波形 在这次的解调设计中，采用 MATLAB 自带的解调函数 fmdemod 对已调信 号进行解调。该函数使用形式如下所示： Z = FMDEMOD(Y,Fc,Fs,FREQDEV) 其中，Y 为输入的调频波函数，Fc、Fs 分别为载波频率和抽样频率，而且要满 足，FREQDEV 为最大频偏。 Fc2Fs 各种情况下的解调波形如下所示： 00.0050.01 -1 0 1 信 信 信 信 信 00.0050.01 -1 0 1 信 信 信 信 信 信 信 00.51 x 10 -3 -

26、1 0 1 信 信 信 信 信 信 信 00.0050.01 -100 0 100 信 信 信 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 信 信 信 信 信 信 图 3-2 时的各波形图3 f m 00.0050.01 -1 0 1 信 信 信 信 信 00.0050.01 -1 0 1 信 信 信 信 信 信 信 00.51 x 10 -3 -1 0 1 信 信 信 信 信 信 信 00.0050.01 -200 0 200 信 信 信 信 信 信 信 信 信

27、 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 信 信 信 信 信 信 图 3-3 时的各波形图1 f m 00.0050.01 -1 0 1 信 信 信 信 信 00.0050.01 -1 0 1 信 信 信 信 信 信 信 00.51 x 10 -3 -1 0 1 信 信 信 信 信 信 信 00.0050.01 -500 0 500 信 信 信 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 200 40

28、0 600 信 信 信 信 信 信 信 信 信 信 信 信 图 3-4 时的各波形图5 . 0 f m 00.0050.01 -1 0 1 信 信 信 信 信 00.0050.01 -1 0 1 信 信 信 信 信 信 信 00.51 x 10 -3 -1 0 1 信 信 信 信 信 信 信 00.0050.01 -500 0 500 信 信 信 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 信 信 信 信 信 信 图 3-5 时的各波形图5 . 0 f m 00.

29、0050.01 -1 0 1 信 信 信 信 信 00.0050.01 -1 0 1 信 信 信 信 信 信 信 00.51 x 10 -3 -1 0 1 信 信 信 信 信 信 信 00.0050.01 -100 0 100 信 信 信 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 -505 x 10 4 0 200 400 600 信 信 信 信 信 信 信 信 信 信 信 信 图 3-6 时的各波形图3 f m 4.4. 高斯白噪声信道特性高斯白噪声信道特性 设正弦波通过加性高斯白噪声信道后的信号为 ( )cos()( )cr tAt

30、n t 其中，白噪声的取值的概率分布服从高斯分布。 ( )n t MATLAB 本身自带了标准高斯分布的内部函数 awgn。awgn 函数产生的随机 序列服从均值为，方差的高斯分布。 0m 2 1 正弦波通过加性高斯白噪声信道后的信号为 ( )cos()( )cr tAtn t 故其有用信号功率为 2 2 A S 噪声功率为 2 N 信噪比满足公式 S N 1010log () S B N 则可得到公式 2 2 10 2 10 B A 我们可以通过这个公式方便的设置高斯白噪声的方差。 在本仿真过程中，我们选择了 5db、10db、15db、20db、25db 五种不同信 噪比以示区别，并绘制输

31、入输出信噪比关系曲线。 020406080100120140160180 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 信 信 信 信 信 信 信 信 信 信 图 4-1 三角函数的输入输出信噪比曲线 020406080100120140160180 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 信 信 信 信 信 信 信 信 信 信 图 4-2 三角波的输入输出信噪比曲线 . 5.5. 调频系统的抗噪声性能分析调频系统的抗噪声性能分析 从前面的分析可知，调频信号的解调有相干解调和非相干解调两种。相干 解调仅适用于窄带调频信号，且需同步信

32、号；而非相干解调适用于窄带和宽带 调频信号，而且不需同步信号，因而是 FM 系统的主要解调方式，所以这里仅仅 讨论非相干解调系统的抗噪声性能，其分析模型如图 9 所示。 图 9 调频系统抗噪声性能分析模型 图中带通滤波器的作用是抑制信号带宽以外的噪声。是均值为零，单 ( )n t 边功率谱密度为的高斯白噪声，经过带通滤波器后变为窄带高斯噪声 0n ( )in t 。限幅器是为了消除接收信号在幅度上可能出现的畸变。 设调频信号为 ( )cos( ) t FMcfStAtKmd 故其输入功率为 2 2 i A S 输入噪声功率为 ioFMNn B 因此输入信噪比为 2 2 i iFM SA NB

33、在大信噪比条件下，信号和噪声的相互作用可以忽略，这时可以把信号和 噪声分开来算，这里，我们可以得到解调器的输出信噪比 222 23 3( ) 8 o f oo m A K m t S Nn f 上式中，为载波的振幅，为调频器灵敏度，为调制信号的最 A fK mf( )m t 高频率，为噪声单边功率谱密度。 on 我们如若考虑为单一频率余弦波时的情况，可得到解调器的制度增益 ( )m t 为 2 2 2 3 2 A o o FM f i o m i S N Gm S n f N 考虑在宽带调频时，信号带宽为 2(1)2()FMfmmBmfff 则可以得到 2 3(1)FMf f Gmm 可以看出

34、，大信噪比时宽带调频系统的信噪比增益是很高的，它与调频指 数的立方成正比。可见，加大调频指数，可使调频系统的抗噪声性能迅速改 fm 善。 6.6. MATLABMATLAB 源代码源代码 echo off close all clear all clc %*FM 调制* t0=0; t1=0.01; %调制信号的时域范围 fs=; %抽样频率 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

35、*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); %已调制信号 %*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); yyn0=fmdemod(yn0,fc,fs,f); yyn1=f

36、mdemod(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_new; snr21=var(m_dem_new)/var

37、(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); in=snr10,snr11,snr12,snr1

38、3,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_new)-1)-fs/2; %画调制信号和已

39、调制信号的波形 disp(按任意键可以看到原调制信号和已调信号在时域和频域内的图形) pause figure(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(已调制信号频谱); % ax

40、is(0 5000 0 10000) %绘制加噪前后的解调信号的波形 disp(按任意键可以看到加噪前后的原基带信号和解调信号在时域和频域 内的图形) pause figure(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)

41、; xlabel(加噪解调后基带信号); m_dem_new,m_dem_new_new,df_new_m_dem_new=GetFFT(m_dem_new,t0,t1,t s,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(无噪声解调后基带信号频谱); %绘制输入输出信噪比曲线 dis

42、p(按任意键可以看到输入输出信噪比关系曲线的图形) pause figure(3); plot(in,out,*) hold on plot(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; end 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=t