信号瞬时频率估计方法的研究： 在信号处理中，信号本身有很多重要的属性，频率特性有：带宽、各频率分量的相对幅值、频率分量间的相对相位关系等；时域特性有信号时宽等。在很多时候，对信号的处理都涉及到需要对平稳或者非平稳信号的频率特性进行估计。平稳信号的频率特性是时不变的，而非平稳信号的频率特性往往是时变的，因此，瞬时频率的定义主要是针对非平稳信号而提出的。Ville 给出了一种统一的瞬时频率的定义： 1()[arg ()] 2i d f t z t dt π= 其中，z(t)是实信号()cos(())s t A t φ=的解析信号。 瞬时频率估计的方法可以分为时频分析和时域分析两类。 就平稳信号而言，由于其功率谱密度函数是不随时间变化的，因此可以直接用参数化或者非参数化谱估计的方法来得到其功率谱，将功率谱中峰值所对应的频率值作为组成该平稳信号的各频率分量的频率的估计值。但是，对于非平稳信号而言，由于其功率谱密度函数是时变的，因此如果要在频域估计其瞬时频率，最简单的方法就是先将其视为短时平稳的信号，每次都用足够短的时间内的数据来构建其功率谱密度函数，将估计得到的结果作为该短时间内的信号瞬时频率，这也就是时频分析中的短时傅立叶变换方法。当然，时频分析还有诸如小波变换等其他的性能更好的变换方法这里不再展开叙述。 下图是用短时傅立叶变换得到的一个非线性调频信号的时频分布图：
时域处理方法则主要是根据信号瞬时频率的定义，先将实信号变换为复信号，再通过对复信号的相位进行求导（模拟）或者差分（数字）的方法来求得瞬时频率。时频分析处理的好处是对于有多个频率分量的信号可以根据功率谱密度函数的各个峰值点估计出对应分量的瞬时频率。而基于相位求导或者差分的时域处理方法却是无法对多频率分量的信号进行瞬时频率估计的。针对这一问题，HUANG. N. E 提出了局域波分解方法，首先将复杂的信号分解成有限个基本模式分量，再对这些基本模式进行相位求导或者差分以估计各分量的瞬时频率。通过局域波分解的方法可以很好的解决相位求导或差分方法的缺陷。时域处理的好处是计算量远小于时频分析处理。 这里主要讨论时域的处理。而要进行时域处理，则通常要首先将物理上的实信号变换为复信号以便取其相位。现有的两种的方法分别是正交变换和hilbert变换。
2DPSK调制与解调系统的仿真 1、 2DPSK基本原理 2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图所示。 图 2DPSK信号 在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。 定义为本码元初相与前一码元初相之差，假设： →数字信息“0”； →数字信息“1”。 则数字信息序列与2DPSK信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK信号相位： 或：0 2DPSK信号的调制原理
一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。 图1.2.1 模拟调制法 2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。 图1.2.2 键控法调制原理图 2DPSK信号的解调原理 2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。 2DPSK信号解调的极性比较法 它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判
基于S 变换的信号瞬时频率特征提取 摘要： S 变换是一种优越的时频分析方法，能够清晰表达信号瞬时频率的变化特征。与传统时频分析方法相对比，S 变换的抗噪性较强，无交叉项干扰。本文提出了采用S 变换来提取调制信号的瞬时频率。仿真实验结果表明，S 变换时频谱能够清晰表示出不同信号的瞬时频率特征。 关键词：时频分析；S 变换；时频图；调制信号；瞬时频率 1 引言 信号的瞬时频率特征可以反映信号在不同时刻的频率变化规律。与传统的时频分析方法相比较，S 变换的时频分析方法具有频率分辨率高、抗噪性强、无交叉项干扰等优点，这使得S 变换能够准确提取信号的瞬时频率。 2S 变换的基本原理 2.1S 变换的提出 S 变换由短时傅里叶变换发展而来，借鉴了短时傅里叶变换加窗的思想。将短时傅里叶变换中的高斯窗函数进行相关伸缩和平移，从而使信号的频率分辨率具备随频率的适应性。这个特点使得S 变换在信号的时频分析中具有明显的优势。 S 变换[1]是由地球物理学家Stockwell 于1996年首次提出的。它可由短时傅里叶变换推导而来，对于连续信号()h t 的短时傅里叶变换为： 2(,)()()j ft STFT f x t w t e dt π+∞ --∞τ=-τ?(1) 其中， 22()t t -δω= (2) 若窗函数为归一化的高斯函数，且对窗函数进行依赖频率的伸缩和平移，那么 22()2(,)t f t f τ τ--ω-= (3) 这样就得到了连续信号()h t 的S 变换定义式： 22()22(,)(f t i ft ST f h t e dt πτ-+∞---∞τ=? (4) 其中，τ为时移因子。 利用S 变换与傅里叶变换之间的紧密联系，可实现信号从S 变换中的无损恢复。S 变换的逆变换形式如式(5)所示： {} 2()(,)j ft h t S f d e df πττ+∞ +∞-∞-∞=?? (5) S 变换还可以看成是信号的小波变换与相位因子的乘积。它采用平移、伸缩的局部高斯窗函数作为母小波，具有频率分辨率高、抗噪性强的优点，且不需满足小波变换的容许性条件。因此，S 变换并不是严格意义上的小波变换，但可以看成是小波变换的一种扩展。 2.2S 变换的瞬时频率表达 由于S 变换为复数，包含实部和虚部，所以S 变换可以表示为： (,)(,)(,)j f S f A f e τττΦ= (6) 其中(,)A f τ为振幅谱，(,)f τΦ为相位谱： (,)f τA =[][]Im (,)(,)arctan Re (,)S f f S f τττ????Φ=?????? (8)
实验题目——PSK（DPSK）调制与解调 一、实验目的 1、掌握绝对码、相对码的概念以及它们之间的变换关系和变换方法。 2、掌握产生PSK(DPSK)信号的方法。 3、掌握PSK(DPSK)信号的频谱特性。 二、实验内容 1、观察绝对码和相对码的波形。 2、观察PSK(DPSK)信号波形。 3、观察PSK(DPSK)信号频谱。 4、观察PSK(DPSK)相干解调器各点波形。 三、实验仪器 1、信号源模块 2、数字调制模块 3、数字解调模块 4、20M双踪示波器 5、导线若干 四、实验原理 1、2PSK(2DPSK)调制原理 2PSK信号是用载波相位的变化表征被传输信息状态的，通常规定0相位载波和π相位载波分别代表传1和传0，其时域波形示意图如图所示。 2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在这种绝对移相的方式中，由于发送端是以某一个相位作为基准的，因而在接收系统也必须有这样一个固定基准相位作参考。如果这个参考相位发生变化，则恢复的数字信息就会与发送的数字信息完全相反，从而造成错误的恢复。这种现象常称为2PSK的“倒π”现象，因此，实际中一般
不采用2PSK 方式，而采用差分移相（2DPSK ）方式。 2DPSK 方式即是利用前后相邻码元的相对载波相位值去表示数字信息的一种方式。如图为对同一组二进制信号调制后的2PSK 与2DPSK 波形。 0 0 0 1 0 1 1 1 0 0 0 1 1 1 0 0 1数字信息（绝对码）PSK 波形 DPSK 波形 相对码 从图中可以看出，2DPSK 信号波形与2PSK 的不同。2DPSK 波形的同一相位并不对应相同的数字信息符号，而前后码元相对相位的差才唯一决定信息符号。这说明，解调2DPSK 信号时并不依赖于某一固定的载波相位参考值。只要前后码元的相对相位关系不破坏，则鉴别这个关系就可以正确恢复数字信息，这就避免了2PSK 方式中的“倒π”现象发生。同时我们也可以看到，单纯从波形上看，2PSK 与2DPSK 信号是无法分辨的。这说明，一方面，只有已知移相键控方式是绝对的还是相对的，才能正确判定原信息；另一方面，相对移相信号可以看成是把数字信息序列（绝对码）变换成相对码，然后再根据相对码进行绝对移相而形成。 2DPSK 的调制原理与2FSK 的调制原理类似，也是用二进制基带信号作为模拟开关的控制信号轮流选通不同相位的载波，完成2DPSK 调制，其调制的基带信号和载波信号分别从“PSK 基带输入”和“PSK 载波输入”输入，差分变换的时钟信号从“PSK-BS 输入”点输入，其原理框图如图所示： 2DPSK 调制原理框图 2、2PSK （2DPSK ）解调原理
中南民族大学 软件课程设计报告 电信学院级通信工程专业 题目2PSK数字信号的调制与解调学生学号 42 指导教师 2012年4月21日
基于MATLAB数字信号2PSK的调制与解调 摘要：为了使数字信号在信道中有效地传播，必须使用数字基带信号的调制与解调，以使得信号与信道的特性相匹配。基于matlab实验平台实现对数字信号的2psk的调制与解调的模拟。本文详细的介绍了PSK波形的产生和仿真过程加深了我们对数字信号调制与解调的认知程度。 关键字：2PSK;调制与解调；MATLAB 引言 当今社会已经步入信息时代，在各种信息技术中，信息的传输及通信起着支撑作用。而对于信息的传输，数字通信已经成为重要的手段。因此，数字信号的调制就显得非常重要。 调制分为基带调制和带通调制。不过一般狭义的理解调制为带通调制。带通调制通常需要一个正弦波作为载波，把基带信号调制到这个载波上，使这个载波的一个或者几个参量上载有基带数字信号的信息，并且还要使已调信号的频谱倒置适合在给定的带通信道中传输。特别是在无线电通信中，调制是必不可少的，因为要使信号能以电磁波的方式发送出去，信号所占用的频带位置必须足够高，并且信号所占用的频带宽度不能超过天线的的通频带，所以基带信号的频谱必须用一个频率很高的载波调制，使期带信号搬移到足够高的频率上，才能够通过天线发送出去。 主要通过对它们的三个参数进行调制，振幅，角频率，和相位。使这三个参量都按时间变化。所以基带的数字信号调制主要有三种方式：FSK，PSK，ASK。在这三种调制的基础上为了得到更高的效果也出现了很多其它的调制方式，如：DPSK，MASK，MFSK，MPSK，APK。它们其中有的一些是将基本的调制方式用在多进制上或者引入了一些新的方式来解决基本调制的一些问题如相位模糊和无法提取位定时信号，另外一些由是组合多种基本的调制方式来达到更好的效果。 基带信号的调制主要分为线性调制和非线性调制，线性调制是指已调信号的频谱结构与原基带信号的频谱结构基本相同，只是占用的频率位置搬移了。而非线性调制则是指它们的结构完全不同不仅仅是频谱搬移，在接收方会出现很多新的频谱分量。在三种基本的调制中，ASK 属于线性调制，而FSK和PSK属于非线性调制。已调信号会在接收方通过各种方式通过解调得到，但是由于噪声和码间串扰，总会有一定的失真。所以人们总是在寻找不同的接收方式来降低误码率，其中的接收方式主要有相干接收和非相干接收。在接收方通过载波的相位信号去检测信号的方法称为相干检测，反之若不利用就称为非相干检测，而对于一些特别的调制有特别的解调方式，如过零检测法。 系统的性能好坏取决于传输信号的误码率，而误码率不仅仅与信道、接收方法有关还和发送端采用的调制方式有很大的关系。我们研究的ASK，FSK，PSK等就主要是发送方的调制方式。
第32卷第1期2009年3月 辽宁师范大学学报(自然科学版)Journal of Liaoning Normal University (Natural Science Edition ) Vol.32 No.1Mar. 2009 文章编号:100021735(2009)0120051207 基于EMD 的信号瞬时频率估计 刘小丹, 孙晓奇, 沈 滨 (辽宁师范大学计算机与信息技术学院,辽宁大连 116029) 收稿日期:2008209224基金项目:辽宁省教育厅科学技术研究项目(20060466) 作者简介:刘小丹(19572),男,吉林蛟河人,辽宁师范大学教授,硕士.E 2mail :xdliu @https://www.doczj.com/doc/9712718275.html, 摘 要:分析了信号瞬时频率的定义及其两种主要的获得信号相位的方法:解析信号法和正交模型法.提出了一种基 于经验模式分解的新的瞬时频率估计方法———正交包络法.该方法计算简单,克服了正交模型法无法由一个时间函数 确定两个时间函数的困难.与Hilbert 变换方法相比,正交包络法使边界问题得到了明显改善.实验证明这是一种有效 的瞬时频率估计方法. 关键词:瞬时频率;正交包络法;EMD ;Hilbert 变换 中图分类号:TP202.4 文献标识码:A 根据Fo urier 分析理论,任何一个平稳信号都可以表示为多个谐波的加权和,对于谐波的某一特定频率,其幅值和相位是常数.而对于非平稳信号,由于其谱特性是随时间变化的,因此不能简单地用Fourier 变换作为非平稳信号的分析工具[1],平稳信号的频率概念也就无法准确解释非平稳信号的时变特性,于是就需要引入一个随时间变化的频率的概念,即瞬时频率. 瞬时频率的一个重要特性是作为时间的函数,用它可以确定信号谱峰的位置.基于这一特性,瞬时频率的概念有着极其重要的应用,因此瞬时频率的估计也就成为许多实际的信号处理应用中一项很有意义的工作.一些信息探测系统只要系统与目标之间有相对运动,多普勒效应就会使频率改变,传播媒质的扰动也会使频率变化,雷达、声呐、移动通信、医疗设备和天文观测都存在这一问题.以雷达信号处理为例,其主要目的是对目标实行检测、跟踪和成像,而像军用飞机一类的目标为了逃避被跟踪,其径向速度是随时间改变的,这使得雷达的多普勒频率具有非平稳的谱.因此,跟踪这类目标需要用到瞬时频率估计技术.瞬时频率估计技术也应用于生物医学.例如,血流的多普勒变化直接关系到心脑血管疾病的诊断.同时,在地震信号处理中,可以利用瞬时频率来确定不同的地质构造.在语音处理等其他诸多领域都有瞬时频率估计技术的应用,详见文献[223]. 从物理学的角度,信号可以分为单分量信号和多分量信号.单分量信号在任意时刻都只有一个频率,该频率称为信号的瞬时频率,而多分量信号则在某些时刻具有多个不同的瞬时频率. 瞬时频率的定义最早是由Carson 和Fry 在研究调频信号时分别提出的,在Gabor 提出了解析信号的概念之后,Ville 将二者结合起来,提出了现在普遍接受的实信号的瞬时频率的定义[4],即:实信号的瞬时频率就是该信号所对应的解析信号的相位关于时间的导数.上述定义只对单分量信号有意义.下面分析一下将瞬时频率定义为复信号相位关于时间的导数的原因. 设一复信号c (t )=A (t )e j φ(t ),A (t )、 φ(t )分别称为信号c (t )的幅度和相位.c (t )的频谱为C (ω)=12 π∫+∞-∞c (t )e -j ωt d t c (t )的总能量E =∫+∞-∞|c (t )|2d t =∫+∞-∞ |C (ω)|2d ω 于是,归一化的函数|c (t )|2/E 和|C (ω )|2/E 可分别作为信号c (t )在时域和频域的能量密度函数,从而得到信号频谱C (ω )的平均频率: 〈ω〉=1E ∫+∞-∞ω|C (ω)|2d ω=1E ∫+∞-∞ ωC (ω)C 3(ω)d ω (3表示共轭运算)
FM调制与解调系统 一、目的 FM在通信系统中的使用非常广泛。FM广泛应用于高保真音乐广播、电视伴音信号的传输、卫星通信和蜂窝电话系统等。 本设计主要是利用MATLAB集成环境下的M文件，编写程序来实现FM调制与解调过程，并分别绘制出基带信号，载波信号，已调信号的时域波形；再进一步分别绘制出相干解调后解调基带信号的时域波形。该设计使用系统开发平台为Windows XP ,程序运行平台使用Windows XP,程序设计语言采用MATLAB，运行程序完成对FM调制和解调结果的观察。通过该本次设计，达到了实现FM信号调制和解调系统的仿真目的。 二、工作原理与计算 通信系统的作用就是将信息从信息源发送到一个或多个目的地。对于任何个通信系统，均可视为由发送端、信道和接收端三大部分组成（如图1所示）。 图1 通信系统一般模型 信息源的作用是把各种信息转换成原始信号，发送设备的作用产生适合传输的信号，信息源和发送设备统称为发送端。发送端将信息直接转换得到的较低频率的原始电信号称为基带信号。通常基带信号
不宜直接在信道中传输。因此，在通信系统的发送端需将基带信号的频谱搬移（调制）到适合信道传输的频率范围内进行传输。这就是调制的过程。信号通过信道传输后，具有将信号放大和反变换功能的接收端将已调制的信号搬移（解调）到原来的频率范围，这就是解调的过程。 调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。而解调是将位于载频的信号频谱再搬回来，并且不失真地恢复出原始基带信号。在本仿真的过程中我们选择用非相干解调方法进行解调。 2.1 FM调制原理 调制在通信系统中具有十分重要的作用。一方面，通过调制可以把基带信号的频谱搬移到所希望的位置上去，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号。另一方面，通过调制可以提高信号通过信道传输时的抗干扰能力，同时，它还和传输效率有关。具体地讲，不同的调制方式产生的已调信号的带宽不同，因此调制影响传输带宽的利用率。可见，调制方式往往决定一个通信系统的性能 2.2 FM解调原理 调制信号的解调分为相干解调和非相干解调两种。相干解调仅仅适用于窄带调频信号，且需同步信号，故应用范围受限；而非相干解调不需同步信号，且对于NBFM信号和WBFM信号均适用，因此是FM 系统的主要解调方式。。
龙源期刊网 https://www.doczj.com/doc/9712718275.html, 基于分段波形的信号瞬时频率计算方法 作者：张亢,程军圣,杨宇,邹宪军 来源：《湖南大学学报·自然科学版》2011年第11期 摘要：针对局部均值分解（Local Mean Decomposition，LMD）中乘积函数（Product Function，PF）分量的瞬时频率计算问题，引入了一种新的信号瞬时频率计算方法.该方法基于分段波形，先将信号分成若干个全波段（full wave），然后以一组递增的反正弦函数定义每个全波段的瞬时相位，进而得到信号的瞬时频率.由该方法得到的瞬时频率理论上是正的、稳定的并且能够确保信号局部特征信息的完整.应用该方法计算了仿真信号和实际齿轮故障振动信号的瞬时频率，并与其他方法求得的瞬时频率进行了对比.结果表明，本文方法非常适合求取信号的瞬时频率. 关键词：故障检测；局部均值分解；乘积函数；纯调频信号；瞬时频率；分段波形 中图分类号：TN911.7 文献标识码：A A Piece-wise Based Signal Instantaneous Frequency Computing Method ZHANG Kang, CHENG Jun-sheng, YANG Yu, ZOU Xian-jun (State key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan Univ, Changsha，Hunan 410082,China) Abstract:To address the computing instantaneous frequency of the product function (PF) in local mean decomposition (LMD), a new instantaneous frequency of a signal computing method was introduced. This method is piece-wise wave based. Firstly, a signal was separated to a number of full waves. Then, the instantaneous phase of each full wave was defined by a set of monotonic increasing arcsine functions. Therefore, the instantaneous frequency of a signal was obtained. Theoretically, the instantanoues frequency obtained in this method was positive, stable and could guarantee the characteristic information of signal integrity. This method was applied to compute the instantaneous frequency of simulated signals and actual gear fault vibration signals, and the results were compared with those obtained in other methods. It has been shown that this method is quite suitable for extracting the instantaneous frequency of a signal. Key words: fault detection;local mean decomposition; product function; pure frequency modulated signal; instantaneous frequency; piece-wise wave
P S K(D P S K)调制与解 调
实验题目——PSK（DPSK）调制与解调 一、实验目的 1、掌握绝对码、相对码的概念以及它们之间的变换关系和变换方法。 2、掌握产生PSK(DPSK)信号的方法。 3、掌握PSK(DPSK)信号的频谱特性。 二、实验内容 1、观察绝对码和相对码的波形。 2、观察PSK(DPSK)信号波形。 3、观察PSK(DPSK)信号频谱。 4、观察PSK(DPSK)相干解调器各点波形。 三、实验仪器 1、信号源模块 2、数字调制模块 3、数字解调模块 4、20M双踪示波器 5、导线若干 四、实验原理 1、2PSK(2DPSK)调制原理 2PSK信号是用载波相位的变化表征被传输信息状态的，通常规定0相位载波和π相位载波分别代表传1和传0，其时域波形示意图如图所示。
2PSK 信号是用载波的不同相位直接去表示相应的数字信号而得出的，在这种绝对移相的方式中，由于发送端是以某一个相位作为基准的，因而在接收系统也必须有这样一个固定基准相位作参考。如果这个参考相位发生变化，则恢复的数字信息就会与发送的数字信息完全相反，从而造成错误的恢复。这种现象常称为2PSK 的“倒π”现象，因此，实际中一般不采用2PSK 方式，而采用差分移相（2DPSK ）方式。 2DPSK 方式即是利用前后相邻码元的相对载波相位值去表示数字信息的一种方式。如图为对同一组二进制信号调制后的2PSK 与2DPSK 波形。 0 0 0 1 0 1 1 1 0 0 0 1 1 1 0 0 1 数字信息（绝对码）PSK 波形 DPSK 波形 相对码 从图中可以看出，2DPSK 信号波形与2PSK 的不同。2DPSK 波形的同一相位并不对应相同的数字信息符号，而前后码元相对相位的差才唯一决定信息符号。这说明，解调2DPSK 信号时并不依赖于某一固定的载波相位参考值。只要前后码元的相对相位关系不破坏，则鉴别这个关系就可以正确恢复数字信息，这就避免了2PSK 方式中的“倒π”现象发生。同时我们也可以看到，单纯从波形上看，2PSK 与2DPSK 信号是无法分辨的。这说明，一方面，只有已知移相键控方式是绝对的还是相对的，才能正确判定原信息；另一方面，相对移相信号可以看成是把数字信息序列（绝对码）变换成相对码，然后再根据相对码进行绝对移相而形成。
探究FM信号的调制与解调 方波的调制： clear all ts=0.001; %信号抽样时间间隔 t=0:ts:10-ts; %时间向量 fs=1/ts; %抽样频率 df=fs/length(t); %fft的频率分辨率 msg=square(4*t); msg2=reshape(msg.',1,length(t)); Pm=fft(msg2)/fs; %求消息信号的频谱 f=-fs/2:df:fs/2-df; subplot(2,1,1) plot(t,fftshift(abs(Pm))) title('消息信号频谱') int_msg(1)=0; %消息信号积分 for ii=1:length(t)-1 int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts; end kf=50; fc=250; %载波频率 Sfm=cos(2*pi*fc*t+2*pi*kf*int_msg); %调频信号 Pfm=fft(Sfm)/fs; % FM信号频谱 subplot(2,1,2) plot(f,fftshift(abs(Pfm))) % 画出已调信号频谱 title('FM信号频谱') Pc=sum(abs(Sfm).^2)/length(Sfm) %已调信号功率 Ps=sum(abs(msg2).^2)/length(msg2) %消息信号功率 fm=50; betaf=kf*max(msg)/fm % 调制指数 W=2*(betaf+1)*fm % 调制信号带宽
012345678910 02 4 6 消息信号频谱 -500-400-300-200-1000100200300400500 01 2 3 FM 信号频谱 已调信号的功率：Pc =0.5000 消息信号的功率：Ps =1 调制指数：betaf =1 调制信号的带宽：W =200 正弦信号的FM 调制： clear all ts=0.001; %信号抽样时间间隔 t=0:ts:10-ts; %时间向量 fs=1/ts; %抽样频率 df=fs/length(t); %fft 的频率分辨率 msg=sin(2*pi*t); msg2=reshape(msg.',1,length(t)); Pm=fft(msg2)/fs; %求消息信号的频谱 f=-fs/2:df:fs/2-df; subplot(2,1,1) plot(t,fftshift(abs(Pm))) title('消息信号频谱') int_msg(1)=0; %消息信号积分 for ii=1:length(t)-1 int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts; end kf=50;
频率和相位是周期函数的两个独立参数，想像一下两个人围着一个圆形场地跑步，离起跑点的圆弧距离是运动位置与起跑点所夹圆心角的函数，这个夹角就是相位，而一定时间所跑圈数是频率，如果两人速度相同（即频率相同），则两人之间的距离是始终不变的，也就是相位差是一定的，这个相位差大小取决于后跑者比先跑者延后起跑的时间。如果两人速度不一样，则之间距离（相位差）不断变化。所以频率不同，相位差不固定。鉴相器不管频率只比较相位，只要相位变化，就给信号给控制器对频率加以控制，使其二者频率一致。 “F(t) = sin(2πft +α)：f就是频率；2πft + α 就是相位；α是t = 0时的相位，即初相位。就是这么简单。 首先，我们通常说的“相位”这个词其实有两个含义： 一、特指周期信号的初相位 二、一般意义上的相位，即“瞬时相位” 频率和相位，一开始都是周期信号的属性，频率是单位时间内的周期数，初相位指周期信号相对所选时间原点的位置，瞬时相位则是指周期信号在任一时刻“走到了一个周期中的哪一步”。 对上面的公式，如果从数学角度理解： 频率就是相位的微分（相位的“行进速度”）或者相位是频率的积分; 这种关系，从数学上推广一步，即使f是变量也成立，再回到物理世界，就发现，不必强求“严格的”周期信号，频率和相位都可以是瞬时值。 频率不同，“初相位”之差是没有意义的，但“瞬时相位”之差仍然存在，不就是两个2πft + α 之差么？ 所谓鉴相器的“相”，指的是就是这种瞬时相位，所以自然不必局限于周期信号，当然也不必局限于“同频”信号，否则“鉴相器”就是个错误的词了。鉴相器的功能，理论上把这种瞬时相位差变换成电压值（当然实际电路总需要经过一段时间才能得出结果，不可能完全“瞬时”） 锁相环的工作原理，表面看是用鉴相器的输出控制VCO的频率，但实际是通过瞬时频率的积分达到相位控制，最终使反馈到鉴相器的瞬时相位与输入的瞬时相位之差趋于零。
密级：无 多线性调频信号瞬时频率估计迭代算法 二炮工程大学士官学院 作者 于鹏鹏 黄向阳 艾名舜 摘要：针对多线性调频信号的瞬时频率估计问题提出一种快速算法，该算法以特征子空间跟踪算法为基础，结合矩阵线性变换和多项式方程求根得到参数估计。该算法的优点是计算量小，其计算量仅与短时傅里叶变换相当；频率分辨力较高；多信号情况下不存在交叉项问题；当多信号的功率差异达到14dB 时仍能有效估计瞬时频率。由于采用了矩阵求逆的步骤，该算法在低信噪比环境下性能较差。仿真实验显示在信噪比不低于6dB 时本文算法具有明显的优越性。 关键词：线性调频 瞬时频率 时频分析 一、引言 线性调频 (Linear Frequency Modulation, LFM) 信号在雷达、声纳、通信等领域有着广泛的应用，由于瞬时频率随时间变化，LFM 信号具有非平稳特性，因此通常采用时频分析的方法对其进行分析及参数估计。短时傅里叶变换是一种简单的时频分析方法，但是时频聚集性较差；Wigner-Ville 分布 [1] （WVD ）的时频聚集性较好，但由于采用了二次型变换，在多LFM 信号情况下不可避免地存在 交叉项，为信号参数估计造成了一定的困难；在Cohen 类时频分布[2]的框架下各种核函数被设计出来用于抑制交叉项，自适应核函数[3-4]的提出进一步提高了交叉项的抑制能力，然而性能较优的时频分析方法计算量也较大，因此在一定程度上较低了此类算法的实用性。 上述方法都是描述信号功率在时频平面上的分布，即信号的功率谱，其频率分辨率受限于信号时窗长度的倒数，这个限制被称为“瑞利限”。超分辨算法利用信号特征子空间的正交性得到信号在频域上的“伪谱”，使有限长信号的频率分辨率能够突破“瑞利限”，从而获得更优的参数估计，但由于传统的超分辨频率估计算法的计算量较大，该类算法很少被用于估计非平稳信号参数。 本文提出一种基于子空间跟踪的信号瞬时频率估计算法，该算法利用数据投影实现信号特征子空间的跟踪，对特征子空间矩阵进行线性变换后得到多项式系数，进而利用多项式方程求根的方法获得信号瞬时频率的估计。本文算法得到的是信号在时频平面上的 “伪谱”，不仅具有较好的时频聚集性，而且在多LFM 信号情况下不存在交叉项的问题，更重要的是，本文算法的计算量仅与短时傅里叶变换相当，因此是一种快速算法。 二、信号模型 考虑一维时间序列S (t )由M 个调频信号线性叠加而成 1 ()()(),1,2,...,M m m m t A t t t T ==+=∑S s n (1) 这里21()exp(2())2m m m t j f t k t π=-+s ，m =1,2,…,M ， A m 、f m 和m k 分别表示第m 个信号的幅度、起 始频率和调频斜率。T 表示有限长采样点数，设采样频率为f s ，测向无模糊范围不大于1 2s f 。n (t )表 示通道噪声，这里假设为零均值高斯白噪声，设等间隔采样，将N 个连续的采样点构成的向量称为一个快拍，N > M ，忽略噪声，t 0时刻的快拍向量0()t y 可以表示为 []0 000 022(1)1111122()(),(1),...,(1)(),(),...,()[(),(),...,()]m m t t M M M j f t j f N t m m m m m m m m m t t M M t t t t t t N A t A t e A t e A t A t A t ππ=-?--?======--+??=??????=?∑∑∑y S S S s s s s s s F (2) 其中，F 是包含当前瞬时频率的矩阵，表达式为
MATLAB仿真信号的相位调制与解调 专业：通信与信息系统 姓名：赵* 学号:********* 指导老师：****教授
摘要 Psk调制是通信系统中最为重要的环节之一，Psk调制技术的改进也是通信系统性能提高的重要途径。本文首先分析了数字调制系统的基本调制解调方法，然后，运用Matlab及附带的图形仿真工具——Simulink设计了这几种数字调制方法的仿真模型。通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。另外，本文还利用Matlab的图形用户界面（GUI）功能为仿真系统设计了一个便于操作的人机交互界面，使仿真系统更加完整，操作更加方便。 关键词：数字调制；分析与仿真；Matlab；Simulink；PSK；QPSK；
1.数字调制技术 (2) 2.PSK调制系统 (3) 2.1 QPSK调制部分，原理框图如图七所示 (6) 2.2 QPSK解调部分，原理框图如图八所示： (8) 3.用Simulink实现PSK调制 (9) 3.1 2PSK仿真 (9) 3.1.1调制 (9) 3.1.2 解调仿真 (12) 3.2 QPSK仿真 (13) 3.2.1 QPSK调制框图 (13) 参考文献 (18)
1.数字调制技术 通信按照传统的理解就是信息的传输与交换。在当今信息社会，通信则与遥感，计算技术紧密结合，成为整个社会的高级“神经中枢”。没有通信，人类社会是不可想象的。一般来说，社会生产力水平要求社会通信水平与之相适应。若通信水平跟不上，社会成员之间的合作程度就受到限制。可见，通信是十分重要的。 通信传输的消息是多种多样的，可以是符号的，文字的，数据和图像的等等。各种不同的消息可以分为两类：一类称为离散消息；另一类称为连续消息。离散消息的状态是可数的或离散的，比如符号，文字或数据等。离散消息也称数字消息。而连续消息则是其状态连续变化的消息，例如，连续变化的语音，图像等。连续消息也称模拟消息。因此按照信道中传输的是模拟信号还是数字信号可以将通信系统分为模拟通信系统和数字通信系统。 数字通信有以下突出的特点：第一，数字信号传输时，信道噪声或干扰所造成的差错，原则上是可以控制的。第二，当需要保密的时候，可以有效的对基带信号进行人为的“扰乱”，即加上密码。 数字通信系统可以用下图表示： →→→→→→→→信数信信数信 信源 道 字受道源字信 息编编调 解译译信 源 码码调码码者 制 道 器 器 器 器 器 器 图一 数字通信在近20年来得到了迅速的发展，其原因是： （1） 抗干扰能力强 （2） 便于进行各种数字信号处理 （3） 易于实现集成化 （4） 经济效益正赶上或超过模拟通信 （5） 传输与交换可结合起来，传输电话与传输数据也可结合起来，成为一个 统一整体，有利于实现综合业务通信网。
郑州轻工业学院 课程设计任务书 题目利用MATLAB实现信号的AM调制与解调专业、班级电子信息工程级班学号姓名 主要内容、基本要求、主要参考资料等： 主要内容： 利用MATLAB对信号 () () ?? ? ? ?≤ = 其他 ,0 t , 100 2t t Sa t m 进行AM调制，载波信号 频率为1000Hz，调制深度为0.5。t0=0.2;首先在MATLAB中显示调制信号的波形和频谱，已调信号的波形和频谱，比较信号调制前后的变化。然后对已调信号解调，并比较解调后的信号与原信号的区别。 基本要求： 1、掌握利用MATLAB实现信号AM调制与解调的方法。 2、学习MATLAB中信号表示的基本方法及绘图函数的调用，实现对常用连续时间信号的可视化表示。 3、加深理解调制信号的变化；验证信号调制的基本概念、基本理论，掌握信号与系统的分析方法。 主要参考资料： 1、王秉钧等. 通信原理[M].北京：清华大学出版社，2006.11 2、陈怀琛.数字信号处理教程----MATLAB释义与实现[M].北京：电子工业出版社,2004. 完成期限：2014.6.9—2014.6.13 指导教师签名： 课程负责人签名： 2014年6月5日
目录 摘要 (1) 1.matlab简介 (2) 1.1matlab基本功能 (2) 1.2matlab应用 (2) 2.系统总体设计方案 (4) 2.1调制信号 (4) 2.1.1 matlab实现调制信号的波形 (4) 2.1.2 matlab实现调制信号的频谱 (4) 2.1.3 matlab实现载波的仿真 (5) 2.2信号的幅度调制 (6) 2.2.1信号的调制 (6) 2.2.2幅度调制原理 (6) 2.2.3 matlab实现双边带幅度调制 (8) 2.2.4 matlab实现已调信号的频谱图 (8) 2.2.5 幅度调制前后的比较 (9) 2.3已调信号的解调 (9) 2.3.1 AM信号的解调原理及方式 (9) 2.3.2 matlab实现已调信号的解调 (11) 2.3.3信号解调前后的比较 (12) 结论与展望 (13) 参考文献 (14) 附录 (15)
1.前言 1.1 序言 随着人类社会步入信息化社会，电子信息科学技术正以惊人的速度发展，开辟了社会发展的新纪元。从20世纪90年代开始至今，通信技术特别是移动通信技术取得了举世瞩目的成就。在通信技术日新月异的今天，学习通信专业知识不仅需要扎实的基础理论，同时需要学习和掌握更多的现代通信技术和网络技术。通信技术正向着数字化、网络化、智能化和宽带化的方向发展。全面、系统地论述了通信系统基本理沦、基本技术以及系统分析与设计中用到的基本工具和方法，并将重点放在数字通信系统上。通信系统又可分为数字通信与模拟通信。传统的模拟通信系统，包括模拟信号的调制与解调，以及加性噪声对幅度调制和角度调制模拟信号解调的影响。数字通信的基本原理，包括模数转换、基本AWGN信道中的数字调制方法、数字通信系统的信号同步方法、带限AWGN信道中的数字通信问题、数字信号的载波传输、数字信源编码以及信道编码与译码等，同时对多径信道中的数字通信、多载波调制、扩频、GSM与IS95数位蜂窝通信。随着数字技术的发展原来许多不得不采用的模拟技术部分已经可以由数字化来实现，但是模拟通信还是比较重要的 1.2 设计任务 本设计是基于MATLAB的模拟相位（PM）调制与解调仿真，主要设计思想是利用MATLAB这个强大的数学软件工具，其中的通信仿真模块通信工具箱以及M檔等，方便快捷灵活的功能实现仿真通信的调制解调设计。还借助MATLAB可视化交互式的操作，对调制解调处理，降低噪声干扰，提高仿真的准确度和可靠性。要求基于MATLAB的模拟调制与解调仿真，主要设计思想是利用MATLAB、simulink檔、M檔等，方便快捷的实现模拟通信的多种调制解调设计。基于simulink对数字通信系统的调制和解调建模。并编写相应的m檔，得出调试及仿真结果并进行分析。
%SSB信号调制解调 clear;clc; f0 = 1; %信源信号频率(Hz) E0 = 1; %信源信号振幅(V) E = 1; %载波分量振幅(V) fc = 10; %载波分量频率（Hz） t0 = 1; %信号时长 snr = 15; %解调器输入信噪比dB dt = 0.003; %系统时域采样间隔 fs = 1/dt; %系统采样频率 df = 0.001; %所需的频率分辨率 t = 0:dt:t0; Lt = length(t); %仿真过程中，信号长度 snr_lin = 10^(snr/10);%解调器输入信噪比 %-------------画出调制信号波形及频谱 %产生模拟调制信号 m = E*cos(2*pi*f0*t); L = min(abs(m));%包络最低点 R = max(abs(m));%包络最高点 %画出调制信号波形和频谱 clf; figure(1); %% %画出调制信号波形 subplot(411); plot(t,m(1:length(t))); axis([0,t0,-R-0.3,R+0.3]);%设置坐标范围 xlabel('t');title('调制信号'); set(gca,'YTick',-R:1:R); subplot(412); [M,m,df1,f] = T2F_new(m,dt,df,fs); %求出调制信号频谱 [Bw_eq] = signalband(M,df,t0); %求出信号等效带宽 f_start_low = fc - Bw_eq; %求出产生下边带信号的带通滤波器的起始频率f_cutoff_low = fc; %求出产生下边带信号的带通滤波器的截止频率f_start_high = fc; %求出产生上边带信号的带通滤波器的起始频率f_cutoff_high = fc + Bw_eq; %求出产生上边带信号的带通滤波器的截止频率 plot(f,fftshift(abs(M))); %画出调制信号频谱%M：傅里叶变换后的频谱序列