持续时间TB码元频域波形

多路子载波频谱密度

北京理工大学信号与系统实验报告5-连续时间系统的复频域分析

北京理工大学信号与系统实验报告5-连续时间系统的复频域分析

实验5连续时间系统的复频域分析 （综合型实验） 一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换定义为 (s)(t)e st X x dt +∞ --∞ = ? （1） 拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ -∞ =? （2） MATLAB 中相应函数如下： (F) L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 () F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 (,) F ilaplace L x =用x 替换结果中的变量t 。

的连续时间系统，其系统函数为s 的有理函数 110 110 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++= +++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下： r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。 求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。 [p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。 还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。调用格

(实验三)连续时间LTI系统的频域分析汇总

实验三 连续时间LTI 系统的频域分析 一、实验目的 1、掌握系统频率响应特性的概念及其物理意义； 2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用； 3、学习和掌握幅度特性、相位特性以及群延时的物理意义； 4、掌握用MA TLAB 语言进行系统频响特性分析的方法。 基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理解LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。 二、实验原理及方法 1 连续时间LTI 系统的频率响应 所谓频率特性，也称为频率响应特性，简称频率响应（Frequency response ），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。 上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：)(*)()(t h t x t y =，由傅里叶变换的时域卷积定理可得到： )()()(ωωωj H j X j Y = 3.1 或者： ) () ()(ωωωj X j Y j H = 3.2 )(ωj H 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即 ? ∞ ∞ --= dt e t h j H t j ωω)()( 3.3 由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说 是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数，

时域抽样与频域抽样

实验三时域抽样与频域抽样 一、实验目的 1.加深理解连续时间信号的离散化过程中的数学概念和物理概念，掌握时域抽样定理（奈奎斯特采样定理）的基本内容。 2.加深对时域取样后信号频谱变化的认识。掌握由抽样序列重建原连续信号的基本原理与实现方法，理解其工程概念。 3.加深理解频谱离散化过程中的数学概念和物理概念，掌握频域抽样定理的基本内容。 二、实验原理 1.时域抽样。 时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件：信号抽样频率f s 大于等于2倍的信号最高频率f m，即f s≥ 2f m。时域抽样先把连续信号x(t)变成适合数字系统处理的离散信号x[k]；然后根据抽样后的离散信号x[k]恢复原始连续时间信号x(t)完成信号重建。信号时域抽样（离散化）导致信号频谱的周期化，因此需要足够的抽样频率保证各周期之间不发生混叠；否则频谱的混叠将会造成信号失真，使原始时域信号无法准确恢复。 2.频域抽样。 非周期离散信号的频谱是连续的周期谱，计算机在分析离散信号的频谱时，必须将其连续频谱离散化。频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件：频域采样点数N 大于等于序列长度M，即N≥M。频域抽样把非周期离散信号x(n)的连续谱X(e jω)变成适合数字系统处理的离散谱X(k)；要求可由频域采样序列X(k)变换到时域后能够不失真地恢复原信号x(n)。

三、实验内容 1.已知模拟信号，分别以T s =0.01s 、0.05s 、0.1s 的采样间隔采样得到x (n )。 （1）当T=0.01s 时，采样得到x(n),所用程序为： %产生连续信号x （t ） t=0:0.001:1; x=sin(20*pi*t); subplot(4,1,1) plot(t,x,'r') hold on title('原信号及抽样信号') %信号最高频率fm 为10 Hz %按100 Hz 抽样得到序列 fs=100; n=0:1/fs:1; y=sin(20*pi*n); subplot(4,1,2) stem(n,y) 对应的图形为： ()sin(20),01a x t t t =π≤≤

连续时间LTI系统的频率特性及频域分析

实验报告 实验项目名称：运用Matlab进行连续时间信号卷积运算 （所属课程：信号与系统） 学院：电子信息与电气工程学院 专业： 10电气工程及其自动化 姓名： xx 学号： 201002040077 指导老师： xxx

一、实验目的 1、学会运用MATLAB 分析连续系统的频率特性。 2、掌握相关函数的调用。 二、实验原理 1、一个连续LTI 系统的数学模型通常用常系数线性微分方程描述，即 )()()()()()(01 )(01)(t e b t e b t e b t r a t r a t r a m m n n +'++=+'++ (1) 对上式两边取傅里叶变换，并根据FT 的时域微分性质可得： )(])([)(])([0101ωωωωωωE b j b j b R a j a j a m m n n +++=+++ 101)()()()()(a j a j a b j b j b j E j R j H n n m m ++++++==ωωωωωωω H ( j ω )称为系统的频率响应特性，简称系统频率响应或频率特性。一般H ( j ω )是复函数，可表示为： )()()(ω?ωωj e j H j H = 其中， )(ωj H 称为系统的幅频响应特性，简称为幅频响应或幅频特性；)(ω?称为系统的相频响应特性，简称相频响应或相频特性。H ( j ω )描述了系统响应的傅里叶变换与激励的傅里叶变换间的关系。H ( j ω )只与系统本身的特性有关，与激励无关，因此它是表征系统特性的一个重要参数。 MATLAB 信号处理工具箱提供的freqs 函数可直接计算系统的频率响应的数值解，其语句格式为：H=freqs(b,a,w)其中，b 和a 表示H ( j ω )的分子和分母多项式的系数向量；w 为系统频率响应的频率范围，其一般形式为w1:p:w2，w1 为频率起始值，w2 为频率终止值，p 为频率取值间隔。 H 返回w 所定义的频率点上系统频率响应的样值。注意，H 返回的样值可能为包含实部和虚部的复数。因此，如果想得到系统的幅频特性和相频特性，还需要利用abs 和angle 函数来分别求得。

北京理工大学信号与系统实验实验5连续时间系统地复频域分析报告报告材料

实验5 连续时间系统的复频域分析 一、实验目的 1.掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握连续时间系统系统函数的定义及复频域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号)(t x 的拉普拉斯变换定义为 )1.....(..........)()(dt e t x s X st ? +∞ ∞ --= 拉普拉斯反变换定义为 )2....(..........)(21)(ds e s X j t x j j st ?∞ +∞ -=σσπ 在MATLAB 中，可以采用符号数学工具箱的laplace 函数和ilaplace 函数进行拉氏变换和反拉氏变换。 L=laplace(F)符号表达式F 的拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 L=laplace(F,t)用t 替换结果中的变量s 。 F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 F=ilaplace(L,x)用x 替换结果中的变量t 。 除了上述ilaplace 函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下： 当 X (s )为有理分式时，它可以表示为两个多项式之比： )3.(..........)()()(0 110 11a s a s a b s b s b s D s N s X N N N N M M M M +?+++?++==---- 式（3）可以用部分分式法展成一下形式 )4.....(.............)(2211N N p s r p s r p s r s X -++-+-= 通过查常用拉普拉斯变换对，可以由式（1-2）求得拉普拉斯逆变换。 利用 MATLAB 的residue 函数可以将 X (s )展成式（1-2）所示的部分分式展开式，该 函数的调用格式为：[r,p,k] = residue(b,a) 其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数

频域采样定理

频域采样定理 （2）频域采样理论的验证。 给定信号如下： ?? ???≤≤-≤≤+=其它02614271301)(n n n n n x 编写程序分别对频谱函数()FT[()]j X e x n ω=在区间]2,0[π上等间隔采样32 和16点，得到)()(1632k X k X 和： 32232()() , 0,1,2,31j k X k X e k ωπω=== 16216()() , 0,1,2,15 j k X k X e k ωπω=== 再分别对)()(1632k X k X 和进行32点和16点IFFT ，得到)()(1632n x n x 和： 323232()I F F T [()] , 0,1,2,,31 x n X k n == 161616()I F F T [()] , 0,1,2,,15x n X k n == 分别画出()j X e ω、)()(1632k X k X 和的幅度谱，并绘图显示x (n)、)()(1632n x n x 和的波形， 进行对比和分析，验证总结频域采样理论。 提示：频域采样用以下方法容易变程序实现。 ① 直接调用MATLAB 函数fft 计算3232()FFT[()]X k x n =就得到()j X e ω 在]2,0[π的32点频率域采样 ② 抽取32()X k 的偶数点即可得到()j X e ω在]2,0[π的16点频率域采样16()X k ，即1632()(2) , 0,1,2,,15X k X k k == 。 ○ 3 当然也可以按照频域采样理论，先将信号x (n)以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT(FFT),得到的就是()j X e ω 在]2,0[π的16点频率域采样16()X k 。 2 频域采样理论的验证程序清单 %频域采样理论验证程序exp2b.m M=27;N=32;n=0:M; %产生M 长三角波序列x(n)

实验八 系统的复频域分析

实验八系统的复频域 分析

一、实验目的 1、掌握系统的复频域分析方法。 2、掌握测试系统的频率响应的方法。 二、预习内容 1、系统频响的方法。（见第四章波特图的介绍） 三、实验原理 1. N 阶系统系统的传递函数 用微分方程描述的N 阶系统为： 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 则系统传递函数可表达为： 用差分方程描述的N 阶系统为： 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 则系统传递函数可表达为： 2.根据系统传递函数的零极点图分析系统 零点：传递函数分子多项式的根。 极点：传递函数分母多项式的根。 根据零极点图的不同分布分析系统。 3.涉及到的Matlab 函数 (1)freqz 函数：实验六中出现过，可用来求单位圆上的有理z 变换的值。调用格式：同实验六 (2)zplane 函数：得到有理z 变换的零极点图。 调用格式：zplane(num,den)

其中，num和 den是按z ?1 的升幂排列的、z 变换分子分母多项式系数的行向量。 (3)roots 函数：求多项式的根。 调用格式：r=roots(c), c 为多项式系数向量；r 为根向量。 四、实验内容 1.系统零极点的求解 (1)求解系统和的零极点，验 证下面程序的运行结果，根据系统零极点图分析系统性质。 b=[1,0,-1]; a=[1,2,3,2]; zr=roots(b); pr=roots(a); plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2); grid; legend('零点','极点'); figure; zplane(b,a); (2)参考上述程序，绘制系统和 的零极点图，并分析系统性质。与用zplane 函数直接绘制系统零极点图（注:圆心的圆圈并非系统的零点）做比较。

时域采样理论与频域采样定理验证

实验4时域采样理论与频域采样定理验证 一 一、实验目的 1时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。 二、实验原理及方法 时域采样定理的要点是： (a)对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(?Ωj X 是原模拟信号频谱()a X j Ω以采样角频率s Ω（T s /2π=Ω）为周期进行周期延拓。公 式为： )](?[)(?t x FT j X a a =Ω )(1∑∞ -∞ =Ω-Ω=n s a jn j X T （b ）采样频率s Ω必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的 频谱不产生频谱混叠。 利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。 理想采样信号)(?t x a 和模拟信号)(t x a 之间的关系为： ∑∞ -∞ =-=n a a nT t t x t x )()()(?δ 对上式进行傅立叶变换，得到： dt e nT t t x j X t j n a a Ω-∞∞ -∞ -∞ =?∑ -=Ω])()([)(?δ dt e nT t t x t j n a Ω-∞ -∞ =∞ ∞ -∑? -)()( δ＝ 在上式的积分号内只有当nT t =时，才有非零值，因此： 课程名称 实验成绩 指导教师 实 验 报 告 院系 班级 学号 姓名 日期

∑∞ -∞ =Ω-=Ωn nT j a a e nT x j X )()(? 上式中，在数值上)(nT x a ＝)(n x ，再将T Ω=ω代入，得到： ∑ ∞ -∞ =-=Ωn n j a e n x j X ω)()(? 上式的右边就是序列的傅立叶变换)(ωj e X ，即 T j a e X j X Ω==Ωωω)()(? 上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变 量ω用T Ω代替即可。 频域采样定理的要点是： a) 对信号x(n)的频谱函数X(e j ω )在[0，2π]上等间隔采样N 点，得到 2()() , 0,1,2,,1j N k N X k X e k N ωπω===- 则N 点IDFT[()N X k ]得到的序列就是原序列x(n)以N 为周期进行周期延拓后的主值区序列，公式为： ()IDFT[()][ ()]()N N N N i x n X k x n iN R n ∞ =-∞ ==+∑ (b)由上式可知，频域采样点数N 必须大于等于时域离散信号的长度M(即N ≥M)，才能使时域不产生混叠，则N 点IDFT[()N X k ]得到的序列()N x n 就是原序列x(n),即()N x n =x(n)。如果N>M ，()N x n 比原序列尾部多N-M 个零点；如果N

理工大学信号与系统实验报告连续时间系统的复频域分析

理工大学信号与系统实验报告连续时间系统的 复频域分析 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】

实验5连续时间系统的复频域分析 （综合型实验） 一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换定义为(s)(t)e st X x dt +∞ --∞ =? （1） 拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ - ∞ = ? （2） MATLAB 中相应函数如下： (F)L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 ()F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量 为t 的结果表达式。 (,)F ilaplace L x =用x 替换结果中的变量t 。 拉氏变换还可采用部分分式法，当(s)X 为有理分式时，它可以表示为两个多项式之比： 110 1 10 ...(s)(s)(s)...M M M M N N N N b s b s b N X D a s a s a ----+++==+++ （3）

上式可以采用部分分式法展成以下形式 1212(s)...N N r r r X s p s p s p = +++--- （4） 再通过查找常用拉氏变换对易得反变换。 利用residue 函数可将X(s)展成（4）式形式，调用格式为： [r,p,k]residue(b,a)=其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分 别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数 连续时间系统的系统函数是指系统单位冲激响应的拉氏变换 (s)(t)e st H h dt +∞ --∞ = ? （5） 连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。 (s)(s)/X(s)H Y = （6） 单位冲激响应(t)h 反映了系统的固有性质，而(s)H 从复频域反映了系统的固有性质。由（6）描述的连续时间系统，其系统函数为s 的有理函数 110 1 10 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++=+++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下：

连续时间信号的频域分析.

课程设计任务书 题目 专业、班级电信1班学号姓名 主要内容、基本要求、主要参考资料等： 基于钟表设计的常识，给出时、分、秒的设计思路，并利用硬件编程语言VHDL或者Verilog-HDL来实 现。要求具有基本功能如调整时间对表、闹铃、计时器等，给出完成控制电路所需要的设计模块；给出硬 件编程语言的实现，并进行仿真；给出下载电路的设计，设计为2种下载方法，其中一种必须为JTAG；同 时设计者报告不允许雷同。 参考资料： 1、潘松、黄继业《EDA技术及其应用》（第四版）科学出版社 2009 2、樊昌信《通信原理》电子出版社 完成期限： 指导教师签名： 课程负责人签名： 年月日

目录 摘要…………………………………………………………………………………II

ABSTRACT……………………………………………………………………………III 绪论…………………………………………………………………………………III 1傅里叶变换原理概述 (1) 1.1 傅里叶变换及逆变换的MATLAB实现 (2) 2 用MATLAB实现典型非周期信号的频域分析 (3) 2.1 单边指数信号时域波形图、频域图 (3) 2.2 偶双边指数信号时域波形图、频域图 (4) 2.3 奇双边指数信号时域波形图、频域图 (4) 2.4 直流信号时域波形图、频域图 (5) 2.5 符号函数信号时域波形图、频域图 (5) 2.6 单位阶跃信号时域波形图、频域图 (6) 2.7 单位冲激信号时域波形图、频域图 (6) 2.8 门函数信号时域波形图、频域图 (7) 3 用MATLAB实现信号的幅度调制 (8) 3.1 实例1 (8) 3.2 实例2 (10) 4 实现傅里叶变换性质的波形仿真 (11) 4.1 尺度变换特性 (11) 4.2 时移特性 (14) 4.3 频移特性 (16) 4.4 时域卷积定理 (18) 4.5 对称性质 (20) 4.6 微分特性 (22) 心得体会 (25) 参考文献 (26) 附录 (27)

连续系统的复频域分析

实验四：连续系统的复频域分析 一、实验目的： 1、掌握连续与离散时间系统的正反复频域与Z域变换 2、掌握利用MATLAB进行零极点分析，进一步了解零极点对整个系统的影响 3、掌握simulink环境下系统建模与仿真以及系统求解。 二、实验内容： 1、已知某连续系统的系统函数为： （1）利用[r, p, k]=residue(num, den),求H（s）的极零点以及多项式系数; （2）画出系统的零极点分布图，判断系统得稳定性。 （3）求h(t)，判断系统得稳定性。 2、已知某离散系统的系统函数为：， （1）利用[r, p, k]=residuez(num, den)求H（z）的极零点以及多项式系数; （2）画出零极点分布图，判断系统得稳定性。 （3）求单位函数响应用impz(b, a)，判断系统是否稳定； 3、已知线性时不变微分方程 在Simulink环境下搭建起系统的仿真模型，并查看仿真结果曲线。（1）写出传递函数H（s），绘出系统模拟框图； （2）当f(t)分别为，，的零状态响应；且当与课本P81的结果进行比较（3）方程的初值为, ,求全响应； 4、已知某信号，n(t)为正态噪声干扰且服从N(0,0.22)分布，对此信号进行采样，采样间隔为0.001s，之后对此信号进行Botterworth低通滤波，从信号中过滤10HZ的输出信号，试对系统进行建模与仿真。 三、实验数据处理与结果分析： 第一题：题1_1：

>> num=[2,5]; den=[1,1,3,2]; [r,p,k]=residue(num,den) r = -0.5750 - 0.7979i -0.5750 + 0.7979i 1.1499 p =-0.1424 + 1.6661i -0.1424 - 1.6661i -0.7152 k =[]

带通采样定理

3.1.3 带通抽样定理 实际中遇到的许多信号是带通型信号，这种信号的带宽往往远小于信号中心频率。若带通信号的上截止频率为，下截止频率为，这时并不需要抽样频率高于两倍上截止频率，可按照带通抽样定理确定抽样频率。 [定理3-2] 带通抽样定理：一个频带限制在内的时间连续信号，信号带宽，令，这里为不大于的最大正整数。如果抽样频率满足条件 ， （3.1-9） 则可以由抽样序列无失真的重建原始信号。 对信号以频率抽样后，得到的采样信号的频谱是的频谱经过周期延拓而成，延拓周期为，如图3-3所示。为了能够由抽样序列无失真的重建原始信号，必须选择合适的延拓周期（也就是选择采样频率），使得位于和的频带分量不会和延拓分量出现混叠，这样使用带通滤波器就可以由采样序列重建原始信号。 由于正负频率分量的对称性，我们仅考虑的频带分量不会出现混叠的条件。 在抽样信号的频谱中，在频带的两边，有着两个延拓频谱分量：和。为了避免混叠，延拓后的频带分量应满足 （3.1-10） （3.1-11） 综合式（3.1-10）和式（3.1-11）并整理得到 （3.1-12） 这里是大于等于零的一个正数。如果取零，则上述条件化为 （3.1-13） 这时实际上是把带通信号看作低通信号进行采样。 取得越大，则符合式（3.1-12）的采样频率会越低。但是有一个上限，因为，而为了避免混叠，延拓周期要大于两倍的信号带宽，即。 因此 （3.1-14） 由于为不大于的最大正整数，因此不大于的最大正整数为，故有 综上所述，要无失真的恢复原始信号，采样频率应满足 ， （3.1-15）H f L f H f ),(H L f f )(t x L H f f B -=N B f M H -=/N B f H /s f m f f m f L s H 212≤≤+10-≤≤N m )(t x )(t x s f )(s nT x )(t x s f )(t x ),(H L f f ),(L H f f --),(H L f f ),(H L f f ),(s L s H mf f mf f +-+-))1(,)1((s L s H f m f f m f ++-++-L s L f mf f ≤+-H s H f f m f ≥++-)1(m f f m f L s H 212≤≤+m m H s f f 2≥m m m f f L s 2≤B f s 2≥B f B f f f m L L s L =≤≤222N B f H /B f L /1-N 10-≤≤N m )(t x s f m f f m f L s H 212≤≤+10-≤≤N m

实验4：连续系统的频域分析

实验4：连续系统的频域分析 一、实验目的 （1）掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。 （2）掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。 二、实验原理 1.周期信号的分解 根据傅里叶级数的原理，任何周期信号都可以分解为三角级数的组合——称为 ()f t 的傅里叶级数。在误差确定的前提下，可以由一组三角函数的有限项叠加而得到。 例如一个方波信号可以分解为： 11114111 ()sin sin 3sin 5sin 7357E f t t t t t ωωωωπ?? = ++++ ??? 合成波形所包含的谐波分量越多，除间断点附近外，它越接近于原波形，在间断点附近，即使合成的波形所含谐波次数足够多，也任存在约9%的偏差，这就是吉布 斯现象（Gibbs ）。 2.连续时间信号傅里叶变换的数值计算 由傅里叶变换的公式： ()()lim ()j t j n n F j f t e dt f n e ωωττωττ∞ ∞ ---∞ →=-∞ ==∑ ? 当 ()f t 为时限信号时，上式中的n 取值可以认为是有限项N ，则有： ()(),0k N j n n F k f n e k N ωτττ-==≤≤∑，其中2k k N π ωτ = 3.系统的频率特性 连续LTI 系统的频率特性称为频率响应特性，是指在正弦信号激励作用下稳态响应随激励信号频率的变化而变化的情况，表示为 () ()() Y H X ωωω= 三、实验内容与方法 1.周期信号的分解 【例1】用正弦信号的叠加近似合成一个频率为50Hz 的方波。 MATLAB 程序如下： clear all; fs=10000; t=[0:1/fs:0.1]; f0=50;sum=0; subplot(211) for n=1:2:9 plot(t,4/pi*1/n*sin(2*pi*n*f0*t),’k ’); hold on; end title(‘信号叠加前’); subplot(212) for n=1:2:9;

连续系统的时域、频域分析

学生实验报告实验课程:信号与 系统E D A 实验地点:东1教 414 学院: 专业: 学号 : 姓名 :

2．信号卷积,根据PPT 中的实验2、2与2、3内容完成课堂练习,写出程序及运行结果。 用Matlab 实现卷积运算)(*)(t h t f ,其中 )()()],2()([2)(t e t h t t t f t εεε-=--=,)2 ()(2t h t h =;对比说明信号)( t f 分别输入系统)(和)(2t h t h 时的输出有什么区别并分析原因。 >> p=0、01; nf=0:p:4; f=2*(heaviside(nf)-heaviside(nf-2)); nh=0:p:6; h=exp(-nh)、*(nh>0); y=conv(f,h);

t=0:length(y)-1; subplot(3,1,1),stairs(nf,f);title('f(t)');axis([0 6 0 2、1]); subplot(3,1,2),plot(nh,h);title('h(t)');axis([0 6 0 1、1]); subplot(3,1,3),plot(0、01*t,y); title('y(t)=f(t)*h(t)'); >> p=0、01; nf=0:p:4; f=2*(heaviside(nf)-heaviside(nf-2)); nh=0:p:6; h=exp(-2*nh)、*(2*nh>0); y=conv(f,h); t=0:length(y)-1; subplot(3,1,1),stairs(nf,f);title('f(t)');axis([0 6 0 2、1]);

连续信号的频域分析

第四章 连续信号的频域分析 将信号分解为若干不同频率的正弦信号或虚指数信号，实质上是将信号在频率域上进行分解，因此根据这种基本思想对信号和系统的分析称为频域分析。这种分解过程是通过傅里叶级数和傅里叶变换这一数学工具来实现的。 本章首先介绍连续信号的傅里叶级数和傅里叶变换，熟悉信号频谱的概念。 4.1 基本要求 1．基本要求 ? 了解傅里叶级数和傅里叶变换的定义及其物理含义； ? 掌握信号频谱和频谱密度的概念； ? 了解连续谱和离散谱的特点和区别； ? 掌握傅里叶变换的常用性质； ? 掌握周期信号傅里叶变换的求解方法。 2．重点和难点 ? 傅里叶变换的性质及其应用 4.2 知识要点 1．周期信号的傅里叶级数 （1）傅里叶级数展开式 三角形式：∑∑∞ =∞=+Ω+=Ω+Ω+=1010)cos(2)]sin()cos([2)(n n n n n n t n A A t n b t n a a t f ?（4-1） 指数形式： ∑∑∞ -∞ =+Ω∞ -∞ =Ω= =n t n n n t n n n F F t f )j(j e e )(? （4-2） 其中 ? +Ω= T t t n t t n t f T a 00 d cos )(2 ，n =0，1，2，? （4-3） ? +Ω= T t t n t t n t f T b 00 d sin )(2，n =1，2，? （4-4） 且

n n n n n n a b b a A a A arctg , ,2 200-=+==? （4-5） ?+Ω-= T t t t n n t t f T F 00 d e )(1j （4-6） （2）两种形式之间的转换关系 0)( e 2 1 j ≥=n A F n n n ? （4-7） 并且|F n |为偶函数，?n 为奇函数，即 ||||n n F F -=，||||n n -=?? （4-8） （3）傅里叶级数的物理含义 通过傅里叶级数可以将任意周期信号f (t )分解为若干个正弦信号（三角形式）或复简谐信号（指数形式）的叠加。每个正弦信号分量的频率为周期信号基波频率的n 倍（n ?0），即n ?，而幅度为A n 或者2|F n |，相位为?n ，将其称作第n 次谐波分量。特别地，将频率为0（即n =0）的分量称为直流分量，幅度为A 0/2或者F 0；频率等于基波频率?（即n =1）的分量称为基波分量。 2．周期信号的频谱 通过傅里叶级数可以将时域中的周期信号分解为直流分量、基波分量和各次谐波分量之和，傅里叶级数展开式中的A n 、?n 或傅里叶系数F n 分别代表了各分量的幅度和相位随谐波次数n （从而频率n ?）的变化关系，称为周期信号的频谱，其中A n 或|F n |称为幅度谱，?n 称为相位谱。 A n 或|F n |、?n 都是关于整型变量n 的实函数，分别以其为纵轴，以n （或者n ?）为横轴，得到的图形称为周期信号的幅度谱图和相位谱图，合称为周期信号的频谱图。 但是，在三角形式的傅里叶级数中，A n 和?n 的自变量n 只能取非负的整数，因此称为单边频谱，而在F n 中，n 可以为任意的整数，相应地将F n 称为双边频谱。对同一个周期信号，其单边和双边频谱可以通过式（4-7）进行相互转换。 所有周期信号的频谱都具有离散性，因此称为离散谱。 3．非周期信号的傅里叶变换及其频谱密度 非周期信号的傅里叶变换及傅里叶反变换的定义为 ?∞ ∞--=t t f F t d e )()j (j ωω （4-9） ?∞ ∞ -= ωωωd )e (j 2π1)(j t F t f （4-10） 其中正变换用于根据信号的时域表达式求其频谱表达式，反变换用于根据其频谱表达式求时域表达式。 通过傅里叶变换可以将信号分解为不同频率的复简谐信号的叠加，而信号的傅里叶变换F (j ?)反映了信号中各分量的幅度和相位随其频率? 的变化关系，称为信号的频谱密度，又称为频谱密度函数或频谱函数。 教材表4-1中列出了一些基本信号的傅里叶变换，在求解复杂信号的傅里叶变换和频谱密度时经常用到。 4．傅里叶变换的性质

信号与系统报告 实验5 连续系统的复频域分析实验

信号与系统 实验报告 实验五连续系统的复频域分析 实验五连续系统的复频域分析 一、实验目的 1. 深刻理解拉普拉斯变换、逆变换的定义，掌握用MATLAB实现拉普拉斯变换、逆变换的方法。 2会求几种基本信号的拉氏变换。 3 掌握用MATLAB绘制连续系统零、极点的方法。 4 求解系统函数H(s)。 二

1已知连续时间信号f（t）=sin（t）u（t）、求出该信号的拉普拉斯变换，并用MATLAB 绘制拉普拉斯变换的曲面图。 syms t; ft=sin(t)*heaviside(t); Fs=Laplace(ft); a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); c=a+i*b; d=ones(size(a)); c=c.*c; c=c+d; c=1./c; c=abs(c); mesh(a,b,c); surf(a,b,c) axis([-0.5,0.5,-2,2,0,10]) colormap(hsv

) 2求[(1-e^(-at))]/t的拉氏变换。 syms t s a f1=(1-exp(-a*t))/t; F=laplace(f1,t,s) F = log(s+a)-log(s) 3求F（s）=-log（s）+ log（s+a）的拉氏逆变换syms t s a F =log(s+a)-log(s); f1=ilaplace(F,s,t) f1 = (1-exp(-a*t))/t

4已知某连续系统的系统函数为： H（s）=(s^2+3s+2)/(8s^4+2s^3+3s^2+5)试用MATLAB求出该系统的零极点，画出零极点分布图。 b=[1 3 2]; a=[8 2 3 0 5]; zs=roots(b); ps=roots(a); hold on plot(real(zs),imag(zs),'o'); plot(real(ps),imag(ps),'x'); grid axis([-2.5,1,-1,1]) 5已知H(s)=（s+1）/（s^2+s+1）,绘制阶跃响应图形，冲激响应图形，频率激响应图形。 syms t s H=(s+1)/(s^2+s+1); f1=ilaplace(H,s,t); f2=heaviside(t);

连续系统的频域分析

第三章傅立叶变换 时域分析：f(t) y f(t)=h(t)*f(t) ↓分解↑ 基本信号δ(t)→LTI →h(t) 频域分析： f(t) ye jωt =h(t)* H(jω)Fe jωt ↓分解↑ 基本信号 sinωt →LTI →H(jω)e jωt e jωt H(jω)：系统的频域响应函数，是信号角频率ω的函数，与t无关. 主要内容： 一、信号的分解为正交函数。 二、周期信号的频域分析?付里叶级数（求和），频谱的特点。信号 三、非周期信号的频域分析?付里叶变换（积分），性质。分析 四、LTI系统的频域分析：频域响应H(jω)；y(jω)= H(jω)?F(jω). (系统分析) 五、抽样定理：连续信号→离散信号.

§3.1 信号分解为正交函数 一、正交： 两个函数满足φ1(t)φ2(t)dt=0,称φi(t),φj(t)在区间(t1 ,t2)正交。 二、正交函数集：几个函数φi(t)φi(t)dt= 0 当i≠j; K i 当i=j. 三、完备正交函数集：在{φ1(t)…φn(t)}之外， 不存在ψ(t)满足ψ (t)φi(t)dt= 0 (i=1,2,…n). 例、三角函数集:{1,cosΩt,cos2Ωt,… ,cosmΩt,…,sinΩt, sin2Ωt,…sin(nΩt),…}区间:（t0,t0+T）,t=2π/Ω为周期. 满足： cosmΩtcosnΩtdt= 0 m≠n T/2 m=n≠0 T m=n=0 sin（mΩt）sin(nΩt)dt= 0 m≠n T/2 m=n≠0 sin(mΩt)cos(nΩt)dt= 0. 所有的m和n. 结论:三角函数集是完备正交集。 推导： cosmΩtcosnΩtdt =(1/2) [cos(m+n) Ωt+cos(m-n) Ωt]dt =(1/2)sin(m+n)Ωt +(1/2)sin(m-n)Ωt =(1/2)[sin(m+n) Ω(t0+T)-sin(m+n)Ωt0] +(1/2)[sin(m-n) Ω(t0+T)-sin(m-n)Ωt0] =0 当m≠n时.

实验三-时域及频域采样定理

广州大学学生实验报告 开课学院及实验室:电子楼３１7 ?２０13年月日 ｎ）以Ｎ为周期进行周期延拓后的主值区序列，

-- 三、实验用ＭＡTL ＡＢ函数介绍 １． fft 功能:一维快速傅立叶变换(FFT)。 Xk=fft(x ｎ,N)：采用FFT 算法计算序列向量x ｎ的Ｎ点DFT ，缺省N 时，ｆft 函数自动按xn 的长度计算xn 的ＤFT 。当N 为2的整数次幂时，ff ｔ按基2算法计算,否则用混合基算法。 2. i ｆf ｔ 功能:一维快速逆傅立叶变换（IFFT)。 调用格式：与f ｆt 相同。 四、实验内容和步骤 （一） 时域采样定理实验 １. 给定模拟信号如下: 0()sin()() at a x t Ae t u t -=Ω 假设式中A=4４4.１28,250π=a , 2500π=Ωｒad ／s ，将这些参数代入上式中,对 () a x t 进行傅立叶变换,得到()a X j Ω,画出它的幅频特性()~a X jf f ,如图３．1所示。根据该曲线可以选 择采样频率。 图3.１ () a x t 的幅频特性曲线 实验过程及原始数据： ｃlf ； A=4４4.128；a=50*ｐi ＊sqrt(２);ｗ０＝5０＊pi ＊sqrt （2）; ｆs=１000; %采样频率10０0HZ T ＝1/fs; n=０：０.05＊fs －1; ％产生的长度区间n x ｔ=A*ｅxp(－ａ*n*T).*sin(w0＊ｎ*Ｔ）; %产生采样序列xt(n) f ＝f ｆt(xt,l ｅng ｔh(n)); %采样序列xt(n)的FFT 变换 k1=0:ｌe ｎg ｔh （f)-1； ｆk1=k1/0.0５; %设置xt(ｎ）的频谱的横坐标的取值 ｓubplot （1，２,1) st ｅｍ(n,xt,＇.＇) %xt 离散图 tit ｌe('(a)fs=1000Hz'); xlab ｅl （'n ＇);ylab ｅl('x(n ）'）; subplot （１,２，2) plot(fk1,ab ｓ（ｆ)) title(＇(a ） FT[x(ｎＴ）]，Fs ＝100０Hz ＇)； xlabel('ｆ(H ｚ)');ylabel('幅度'） 2. 按照选定的采样频率对模拟信号进行采样,得到时域离散信号()x n : 0()()sin()() anT a x n x nT Ae nT u nT ==Ω 这里给定采样频率如下：1s f kHz =，300Ｈz ，20０Hz 。分别用这些采样频率形成时域离散信号, 按顺序分别用 1() x n 、 2() x n 、3() x n 表示。选择观测时间 50p T ms =。 实验过程及原始数据： c ｌf; A=444.12８;a=50＊ｓqrt(2)*pi;ｗ０=５０*sqrt(2)*p ｉ; Ｔｐ=50/１０00;F1=1000;F2=3０0;Ｆ3=200; Ｔ１=１/F1;T2=１／F ２;T3=1/F ３； n1＝0:0．5： T ｐ*F1-1;n2＝0:0.5:T ｐ*F2-1;n3=0:0．５:T ｐ＊Ｆ3－1; x1=A*exp （-ａ＊n1*T1).＊s ｉn （w0*n １*T １)； ｘ2=A*ｅxp （－ａ*n2＊T ２).*ｓin(w0*n ２*T2）； x3=A*exp （-a*n3*T3).*ｓｉn(w ０*ｎ３*T3）； f1=fft （x １，length （n1)); f2＝fft(x2，len ｇth(n2））； ｆ3＝fft(x3，l ｅng ｔｈ(n3)); k １=0:ｌｅｎgth （f1)-1; f ｋ1=ｋ1/Ｔp; ｋ２=0：length(f2）－1； fk2=ｋ2/Tp; k3＝0：le ｎg ｔh(ｆ3)-１;