当前位置：文档之家 › 连续时间信号傅里叶级数分析MATLAB课程设计

连续时间信号傅里叶级数分析MATLAB课程设计

  • 格式：doc
  • 大小：962.00 KB
  • 文档页数：31

Matlab应用实践课程设计 1 绪论 本次课程是通过MATLAB软件来实现数字信号系统里的相关图像和相关仿真的软件。近年来，MATLAB以其强大的矩阵计算和图像视化功能逐渐为国人所知。MATLAB是mathworks公司的软件产品，MATLAB已经成为一个系列产品：MATLAB主包各种工具（toolbox）。功能丰富的工具箱大致分为两类：功能型工具箱和领域型工具箱。功能型工具箱主要用来扩充MATLAB的符号计算功能﹑图形建模仿真功能﹑文字处理功能以及与硬件实时交互功能，能用于多种学科。而领域型工具箱是专业性很强的，如控制工具（control toolbox）﹑信号处理工具箱（signal processing toolbox）等。MATLAB (MATrix LABoratory)具有用法简易、可灵活运用、程式结构强又兼具延展性[3]。 以下为其几个特色： 功能强的数值运算 - 在MATLAB环境中，有超过500种数学、统计、科学及工程方面的函数可使用，函数的标示自然，使得问 题和解答像数学式子一般简单明了，让使用者可全力发挥在解题方面，而非浪费在电脑操作上。 先进的资料视觉化功能 - MATLAB的物件导向图形架构让使用者可执行视觉数据分，并制作高品质的图形，完成科学性或工程 性图文并茂的文章。 高阶但简单的程式环境 - 做为一种直译式的程式语言，MATLAB容许使用者在短时间内写完程式，所花的时间约为用 FORTRAN 或 C 的几分之一，而且不需要编译(compile)及联结 (link) 即能执行，同时包含了更多及更容易使用的内建 功能。 开放及可延伸的架构 - MATLAB容许使用者接触它大多数的数学原使码，检视运算法，更改现存函数，甚至加入自己的函数使 MATLAB成为使用者所须要的环境。 丰富的程式工具箱 - MATLAB的程式工具箱融合了套装前软体的优点，与一个灵活的开放但容易操作之环境，这些工具箱提 供了使用者在特别应用领域所需之许多函数。现有工具箱有：符号运算（利用Maple V的计算核心执行 ）、影像处理、统计分析、讯号处理、神经网路、模拟分析、控制系统、即时控制、系统确认、强建控 制、弧线分析、最佳化、模糊逻辑、mu分析及合成、化学计量分析[4]。 Matlab应用实践课程设计

2 1 MATLAB简介

1.1 MATLAB语言功能 MATLAB功能丰富，可扩展性强。MATLAB软件包括基本部分和专业扩展两大部分的功能。基本部分包括：矩阵的运算和各种变换；代数和超越方程的求解；数据处理和傅立叶变换；数值部分等等，可以充分满足大学理工科本科的计算需要。扩展部分称为工具箱。它实际上是用MATLAB的基本语句辩称的各种子程序集，用于解决某一方面的专门问题，或实现某一类的新算法。 MATLAB 具有以下基本功能 （1）数值计算功能； （2）符号计算功能； （3）图形处理及可视化功能； （3）可视化建模及动态仿真功能[6]。

1.2 MATLAB语言特点 MATLAB 给用户带来的是最直观，最简洁的程序开发环境。它具有以下特点： （1）语言简洁紧凑，使用方便灵活，库函数极其丰富。MATLAB 程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。 （2）运算符丰富。由于MATLAB 是用C 语言编写的，MATLAB 提供了和C语言几乎一样多的运算符，灵活使用MATLAB 的运算符将使程序变得极为简短。（3）MATLAB 既具有结构化的控制语句（如for 循环，while 循环，break 语句和if 语句），又有面向对象编程的特性。 （4）程序限制不严格，程序设计自由度大。例如，在MATLAB 里，用户无需对矩阵预定义就可使用。 （5）程序的可移植性很好，基本上不做修改就可以在各种型号的计算机和操作系统上运行。 （6）MATLAB 的图形功能强大。在FORTRAN 和C 语言里，绘图都很不容易，但在MATLAB 里，数据的可视化非常简单。MATLAB 还具有较强的编辑图形界面的能力。 （7）功能强大的工具箱是MATLAB 的另一特色。MATLAB 包含两个部分： 核心部分和各种可选的工具箱。核心部分中有数百个核心内部函数。其工具箱又 分为两类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计 算功能，图示建模仿真功能，文字处理功能以及与硬件实时交互功能，而学科性工具箱是专业性比较强的，如control, toolbox, signl ，proceessing ，toolbox,commumnication toolbox 等[4]。 Matlab应用实践课程设计 3 2 傅里叶级数基本原理概要

2.1 周期信号的傅里叶分解

设有连续时间周期信号，它的周期为T，角频率，且满足狄里赫

利条件，则该周期信号可以展开成傅里叶级数，即可表示为一系列不同频率的正弦或复指数信号之和。傅里叶级数有三角形式和指数形式两种。

1. 三角形式的傅里叶级数[2]：

式中系数，称为傅里叶系数，可由下式求得： [ 2. 指数形式的傅里叶级数[2]：

式中系数称为傅里叶复系数，可由下式求得： 周期信号频谱具有三个特点[1]： （1） 离散性，即谱线是离散的； （2） 谐波性，即谱线只出现在基波频率的整数倍上； （3） 收敛性，即谐波的幅度随谐波次数的增高而减小。 周期信号的傅里叶分解用Matlab进行计算时，本质上是对信号进行数值积分运算。在Matlab中有多种进行数值积分运算的方法，我们采用quadl函数，它有两种其调用形式。 (1) y＝quadl(‘func’, a, b)。 其中func是一个字符串，表示被积函数的.m文件名（函数名）；a、b分别表示定积分的下限和上限。 (2) y＝quadl(@myfun, a, b)。其中“@”符号表示取函数的句柄，myfun表示所定义函数的文件名。[6] Matlab应用实践课程设计 4 2.3 三角形式和指数形式傅里叶级数及各系数间的关系

傅里叶级数的指数形式和三角形式是等价的，其系数可互相转换。表2-1综合了三角 形式和指数形式傅里叶级数及其系数，以及各系数间的关系。 表2-1 周期函数展开为傅里叶级数 形式 指数形式 三角函数形式

展开式

傅里叶系数

系数间的关系 Matlab应用实践课程设计

5 2.2周期信号的频谱 周期信号经过傅里叶分解可表示为一系列正弦或复指数信号之和。为了直观地表示出信号所含各分量的振幅，以频率（或角频率）为横坐标，以各谐波的振幅或虚指数函数的幅度为纵坐标，可画出幅度－频率关系图，称为幅度频谱或幅度谱。类似地，可画出各谐波初相角与频率的关系图，称为相位频谱或相位谱[2]。 Matlab应用实践课程设计

6 3 用MATLAB实现周期信号的傅立叶级数分解与综合

下面以矩形信号为例介绍用MATLAB来实现周期信号的傅立叶技术的分解与综合。 3.1 合成波形与原波形之间的关系 本文使用的连续周期矩形信号是一个周期为4，占空比为50%，幅值为1的矩形波信号，从上面的分析可以得知，这个矩形波信号可以分解成傅立叶级数也是就无数个不同频率的三角波的叠加，用MATLAB软件可以很容易的画出不同个数谐波叠加形成的合成波的形状（源代码见附录），以便将它们与原信号做比较，给理论分析以正确的实验基础。 下图为分别用1到10次谐波叠加形成的合成波与原信号的比较，可以很快地看出，当叠加的谐波数越多是，与原信号的差别就越小，这样就直接证明了傅立叶级数理论的正确性。

图3-1 不同次数的谐波的合成波与原信号的比较 Matlab应用实践课程设计

7 3.2 吉布斯现象 分析傅立叶级数的公式

满足狄里赫利条件的周期函数表示成的傅立叶级数都收敛。狄里赫利条件如下： 1. 在任何周期内，x(t)必须绝对可积； 2. 在任一有限区间中，x(t)只能取有限个最大值或最小值； 3. 在任何有限区间上，x(t)只能有有限个第一类间断点。 所谓的吉布斯现象就是：在x(t)的不可导点上，如果我们只取x(t)等式右边的无穷级数中的有限项作和X(t)，那么X(t)在这些点上会有起伏[1]。 具体现象如下图所示，以下分别为谐波次数为N=50，N=100，N=500合成波的情况。

图3.2 谐波次数为N=50合成波 图3.3 谐波次数为N=100合成波 Matlab应用实践课程设计

8 图3.4 谐波次数为N=500合成波

从上面的图像中可以看出，当N=500的时候，合成波与原来的方波拟合得非常好，但是在不可导的点上，即为x=-1.5,x=-0.5,x=0.5,x=1.5这样的点的时候，合成波会有较大的波动，这就是非常明显的吉布斯现象。 Matlab应用实践课程设计

9 4 用MATLAB实现周期信号的单边频谱及双边频谱。

用周期矩形波信号为例其中周期T=5，脉冲宽度 τ=1。在下图中，τ为一个脉冲的宽度，T为脉冲的周期，它们是方波信号的两个非常重要的参数，它们的大小关系将直接影响方波信号的形状和性质。 以下皆以方波为例来介绍方波各种参数对其频谱的影响以及单边频谱和双边频谱的关系，方波示意图如下：

图4.1 周期矩形脉冲信号 4.1 单边，双边（幅度，相位）频谱及其关系

4.1.1单边，双边（幅度，相位） 如前所述，周期信号可以分解成一系列正弦（余弦）信号或虚指数信号之和，即 NjntnnNftFe



0

11cossin2NNnnnnaantbnt



其中，1122njnnnnFAeajb

或2211||22arctannnnnnnnFAabba 幅度和相位 为了直观地表示出信号所含各分量的振幅nA或||nF，随频率的变化情况，通常以角频率为横坐标，以各次谐波的振幅nA或虚指数函数||nF的幅度为纵坐标，画出如图4.2和4.3所示的各谐波的振幅nA或||nF与角频率的关系图，称为周期信号的幅度（振幅）频谱，简称

合集下载

连续时间信号卷积运算的MATLAB实现课程设计

连续时间信号卷积运算的MATLAB实现课程设计

信号与通信综合设计项目题目：连续时间信号的卷积运算学院：电子与信息工程学院专业：通信工程班级：姓名：学号：联系方式：指导教师：报告成绩：摘要卷积积分是信号与系统的时域分析的重要方法之一。

为了对其计算进行一定的简化,本文运用MATLAB的数值计算法通过实例计算,验证了函数的有效性,并且展示了MATLAB在卷积积分方面的优越性。

关键词：MATLAB；卷积积分；信号与系统AbstractConvolution is one of the most important methods of time domain analysis of signals and systems. In order to simplify its calculation, this paper uses MATLAB's numerical calculation method to verify the validity of the function by example calculation, and shows the advantages of MATLAB in convolution integral.Keywords：MATLAB；Convolution points；Signals and systems目录摘要 (I)Abstract (II)目录 ............................................................................................................................................ I II一、设计任务及要求 (1)1.1课程设计目的 (1)1.2 设计要求 (1)1.3应解决问题 (1)二、卷积定义以及原理 (1)三、MATLAB设计卷积仿真 (1)3.1思路 (1)3.2步骤 (2)四、卷积的MATLAB设计源程序 (2)4.1方案一 (2)步骤一：函数sconv编辑如下： (2)步骤二：定义阶跃函数heaviside (3)步骤三：定义主函数 (3)步骤四：结果运行MATLAB后，图形如下： (3)4.2方案二 (4)五、总结 (6)参考文献 (7)谢辞 (8)一、设计任务及要求1.1课程设计目的（1）进一步熟悉常用MATLAB 命令；（2）提高综合运用所学的理论知识，独立分析和解决问题的能力；（3）进一步熟悉MATLAB 的正确使用方法。

连续频域分析实验报告(3篇)

连续频域分析实验报告(3篇)

第1篇一、实验目的1. 理解连续频域分析的基本概念和原理。

2. 掌握连续信号的傅里叶变换及其性质。

3. 学会使用MATLAB进行连续信号的频域分析。

4. 通过实验加深对连续信号频域特性的理解。

二、实验原理连续信号的频域分析是信号与系统分析中的重要内容，它可以将信号从时域转换到频域，便于分析和处理。

本实验主要涉及以下原理：1. 傅里叶变换：傅里叶变换是一种将信号从时域转换到频域的方法，它可以将任何连续时间信号分解为一系列不同频率的正弦波和余弦波的叠加。

傅里叶变换的数学表达式为：F(w) = ∫ f(t) e^(-jwt) dt其中，F(w)表示信号的频谱，f(t)表示信号，w表示频率，j表示虚数单位。

2. 傅里叶变换的性质：傅里叶变换具有许多性质，如时移性、频移性、尺度变换、卷积定理等，这些性质可以简化信号的频域分析。

3. MATLAB函数：MATLAB提供了丰富的函数用于连续信号的频域分析，如fourier函数用于计算信号的傅里叶变换，ifourier函数用于计算信号的傅里叶逆变换等。

三、实验内容1. 实验一：信号傅里叶变换（1）输入一段连续信号，如正弦波、方波等。

（2）使用fourier函数计算信号的傅里叶变换。

（3）绘制信号的时域波形和频谱图，观察信号的频谱特性。

2. 实验二：信号时移和频移（1）对实验一中的信号进行时移和频移操作。

（2）观察信号的时域波形和频谱图的变化，验证傅里叶变换的时移性和频移性。

3. 实验三：信号尺度变换（1）对实验一中的信号进行尺度变换操作。

（2）观察信号的时域波形和频谱图的变化，验证傅里叶变换的尺度变换性。

4. 实验四：信号卷积（1）输入两个连续信号，如矩形脉冲信号、三角脉冲信号等。

（2）使用conv函数计算两个信号的卷积。

（3）绘制信号的时域波形和频谱图，观察信号的卷积特性。

四、实验结果与分析1. 实验一：通过实验，我们得到了信号的时域波形和频谱图，可以看出信号的频谱特性与信号的时域特性密切相关。

实验四 信号的分解与合成

实验四 信号的分解与合成

实验四信号的分解与合成实验目的：1.了解信号的分解与合成原理；2.掌握连续时间信号的傅里叶级数分解公式及其应用；3.掌握离散时间信号的傅里叶变换公式及其应用。

实验原理：1.信号的分解任何信号都可以分解成若干谐波的叠加。

这是因为任何周期信号都可以表示为若干谐波的叠加。

傅里叶级数分解公式：$$x(t)=\sum_{n=-\infty}^{+\infty} C_ne^{jn\omega_0t}$$其中，$C_n$为信号的各级谐波系数，$\omega_0$为信号的基波频率。

当信号为实信号时，其傅里叶级数中只有实系数，且对称性可利用，因此实际计算中可以只计算正频率系数，即$$x(t)=\sum_{n=0}^{+\infty} A_n\cos(n\omega_0t+\phi_n)$$其中，$A_n$为信号各级谐波幅度，$\phi_n$为各级谐波相位。

若信号不是周期信号，则可以采用傅里叶变换进行分解。

2.信号的合成对于任意信号$y(t)$，都可以表示为其傅里叶系数与基波频率$\omega_0$的乘积的叠加，即$$y(t)=\sum_{n=-\infty}^{+\infty}C_ne^{jn\omega_0t}$$若$y(t)$为实信号，则其傅里叶系数中只有正频率系数，即$$y(t)=\sum_{n=0}^{+\infty}A_n\cos(n\omega_0t+\phi_n)$$实验步骤：一、连续时间信号的傅里叶级数分解1.打开Matlab软件，使用line或scatter等函数绘制出函数$f(x)=x(0<x<2\pi)$的图像。

2.使用Matlab的fft函数对f(x)进行逆傅里叶变换得到其傅里叶级数分解。

3.将得到的傅里叶级数分解与原函数的图像进行比较，分析级数中谐波幅度的变化规律。

二、离散时间信号的傅里叶变换1.使用Matlab生成一个为$sin(\pi k/4),0\le k\le 15$的离散时间信号。

通信原理实验教程(MATLAB)

通信原理实验教程(MATLAB)

实验教程目录实验一：连续时间信号与系统的时域分析-------------------------------------------------6一、实验目的及要求---------------------------------------------------------------------------6二、实验原理-----------------------------------------------------------------------------------61、信号的时域表示方法------------------------------------------------------------------62、用MATLAB仿真连续时间信号和离散时间信号----------------------------------73、LTI系统的时域描述-----------------------------------------------------------------11三、实验步骤及内容--------------------------------------------------------------------------15四、实验报告要求-----------------------------------------------------------------------------26 实验二：连续时间信号的频域分析---------------------------------------------------------27一、实验目的及要求--------------------------------------------------------------------------27二、实验原理----------------------------------------------------------------------------------271、连续时间周期信号的傅里叶级数CTFS---------------------------------------------272、连续时间信号的傅里叶变换CTFT--------------------------------------------------283、离散时间信号的傅里叶变换DTFT -------------------------------------------------284、连续时间周期信号的傅里叶级数CTFS的MATLAB实现------------------------295、用MATLAB实现CTFT及其逆变换的计算---------------------------------------33三、实验步骤及内容----------------------------------------------------------------------34四、实验报告要求-------------------------------------------------------------------------48 实验三：连续时间LTI系统的频域分析---------------------------------------------------49一、实验目的及要求--------------------------------------------------------------------------49二、实验原理----------------------------------------------------------------------------------491、连续时间LTI系统的频率响应-------------------------------------------------------492、LTI系统的群延时---------------------------------------------------------------------503、用MATLAB计算系统的频率响应--------------------------------------------------50三、实验步骤及内容----------------------------------------------------------------------51四、实验报告要求-------------------------------------------------------------------------58 实验四：调制与解调以及抽样与重建------------------------------------------------------59一、实验目的及要求--------------------------------------------------------------------------59二、实验原理----------------------------------------------------------------------------------591、信号的抽样及抽样定理---------------------------------------------------------------592、信号抽样过程中的频谱混叠----------------------------------------------------------623、信号重建--------------------- ----------------------------------------------------------624、调制与解调----------------------------------------------------------------------------------645、通信系统中的调制与解调仿真---------------------------------------------------------66三、实验步骤及内容------------------------------------------------------------------------66四、实验报告要求---------------------------------------------------------------------------75 实验五：连续时间LTI系统的复频域分析----------------------------------------------76一、实验目的及要求------------------------------------------------------------------------76二、实验原理--------------------------------------------------------------------------------761、连续时间LTI系统的复频域描述--------------------------------------------------762、系统函数的零极点分布图-----------------------------------------------------------------773、拉普拉斯变换与傅里叶变换之间的关系-----------------------------------------------784、系统函数的零极点分布与系统稳定性和因果性之间的关系------------------------795、系统函数的零极点分布与系统的滤波特性-------------------------------------------806、拉普拉斯逆变换的计算-------------------------------------------------------------81三、实验步骤及内容------------------------------------------------------------------------82四、实验报告要求---------------------------------------------------------------------------87 附录：授课方式和考核办法-----------------------------------------------------------------88实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数；2、掌握连续时间和离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程；3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质；4、掌握利用MA TLAB计算卷积的编程方法，并利用所编写的MA TLAB程序验证卷积的常用基本性质；掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线。