(电气)自动控制原理21(频域分析)

大作业1(机电控制系统时域频域分析)

《机电系统控制基础》大作业一 基于MATLAB的机电控制系统响应分析 哈尔滨工业大学 2013年11月4日

1 作业题目 1. 用MATLAB 绘制系统2 ()25()() 425 C s s R s s s Φ== ++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2 ()25 ()()425 C s s R s s s Φ==++的单位阶跃响应性能指标：上升时间、峰值时间、调节时间和超调量。 3. 数控直线运动工作平台位置控制示意图如下： X i 伺服电机原理图如下： L R (1)假定电动机转子轴上的转动惯量为J 1，减速器输出轴上的转动惯量为J 2，减速器减速比为i ，滚珠丝杠的螺距为P ，试计算折算到电机主轴上的总的转动惯量J ； (2)假定工作台质量m ，给定环节的传递函数为K a ，放大环节的传递函数为K b ，包括检测装置在内的反馈环节传递函数为K c ，电动机的反电势常数为K d ，电动机的电磁力矩常数为K m ，试建立该数控直线工作平台的数学模型，画出其控制系统框图； (3)忽略电感L 时，令参数K a =K c =K d =R=J=1，K m =10，P/i =4π，利用MATLAB 分析kb 的取值对于系统的性能的影响。

2 题目1 单位脉冲响应曲线 单位阶跃响应曲线

源代码 t=[0:0.01:1.6]; %仿真时间区段和输入 nC=[25]; dR=[1,4,25]; fi=tf(nC,dR); %求系统模型 [y1,T]=impulse(fi,t); [y2,T]=step(fi,t); %系统响应 plot(T,y1); xlabel('t(sec)'),ylabel('x(t)'); grid on; plot(T,y2); xlabel('t(sec)'),ylabel('x(t)'); grid on; %生成图形 3 题目2 借助Matlab，可得： ans = 0.4330 0.6860 25.3826 1.0000 即

(实验三)连续时间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 ω)通常是复数，

连续时间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 函数来分别求得。

控制系统的频域分析实验报告

实验名称： 控制系统的频域分析 实验类型：________________同组学生姓名：__________ 一、实验目的和要求 用计算机辅助分析的方法，掌握频率分析法的三种方法，即Bode 图、Nyquist 曲线、Nichols 图。 二、实验内容和原理 （一）实验原理 1．Bode(波特)图 设已知系统的传递函数模型： 1 1211121)(+-+-+???+++???++=n n n m m m a s a s a b s b s b s H 则系统的频率响应可直接求出： 1 1211121)()()()()(+-+-+???+++???++=n n n m m m a j a j a b j b j b j H ωωωωω MATLAB 中，可利用bode 和dbode 绘制连续和离散系统的Bode 图。 2．Nyquist(奈奎斯特)曲线 Nyquist 曲线是根据开环频率特性在复平面上绘制幅相轨迹，根据开环的Nyquist 线，可判断闭环系统的稳定性。 反馈控制系统稳定的充要条件是，Nyquist 曲线按逆时针包围临界点(-1，j0)p 圈，为开环传递函数位于右半s 一平面的极点数。在MATLAB 中，可利用函数nyquist 和dnyquist 绘出连续和离散系统的乃氏曲线。 3．Nicho1s(尼柯尔斯)图 根据闭环频率特性的幅值和相位可作出Nichols 图，从而可直接得到闭环系统的频率特性。在 MATLAB 中，可利用函数nichols 和dnichols 绘出连续和离散系统的Nichols 图。 （二）实验内容 1．一系统开环传递函数为 ) 2)(5)(1(50)(-++=s s s s H 绘制系统的bode 图，判断闭环系统的稳定性，并画出闭环系统的单位冲击响应。 2．一多环系统 ) 10625.0)(125.0)(185.0(7.16)(+++=s s s s s G 其结构如图所示 试绘制Nyquist 频率曲线和Nichols 图，并判断稳定性。 （三）实验要求

2018年自动控制原理期末考试题[附答案解析]

. 2017 年自动控制原理期末考试

卷与答案 一、填空题（每空1分，共20分） 1、对自动控制系统的基本要求可以概括为三个方面，即：稳定性、快速性和准确性。 2、控制系统的输出拉氏变换与输入拉氏变换在零初始条件下的比值称为传递函数。 3、在经典控制理论中，可采用劳斯判据(或：时域分析法)、根轨迹法或奈奎斯特判据( 或：频域分析 法) 等方法判断线性控制系统稳定性。 4、控制系统的数学模型，取决于系统参数, 与外作用及初始条件无关。和结构 A( )L( )lg) 或：，横坐标为( 、线性系统的对数幅频特性，纵坐标取值为5。20lg

，其中P 是指开环传函中具有正实部的极点的个数，6、奈奎斯特稳定判据中，Z = P - RZ是指闭环 传函中具有正实部的极点的个数，R 指奈氏曲线逆时针方向包围(-1, j0 )整圈数。 定义为调整时间。%是超调量。、在二阶系统的单位阶跃响应图中，7t s K)A(22 (T)1K)，则其开环幅频特性为(T8、设系统的开环传递函数为1，相12 s(Ts 1)(T s 1)21110) (T ) tg (Ttg 。频特性为()9021 9、反馈控制又称偏差控制，其控制作用是通过给定值与反馈量的差值进行的。 0.5t0.2 t，则该系统的传递函数G(s) 为、若某系统的单位脉冲响应为10。5e g (t) 10e 510 0.2 sss 0.5s

11、自动控制系统有两种基本控制方式，当控制装置与受控对象之间只有顺向作用而无反向联系时，称 为开环控制系统；当控制装置与受控对象之间不但有顺向作用而且还有反向联系时，称为闭环控制系 统；含有测速发电机的电动机速度控制系统，属于闭环控制系统。 、根轨迹起始于开环极点，终止于开环零点。12 、稳定是对控制系统最基本的要求，若一个控制系统的响应曲线为衰减振荡，则该系统13稳定。判断 页脚

自动控制原理-线性系统的频域分析实验报告

自动调节系统频域分析 班级11081801 学号1108180135 姓名王佳炜 日期2014.1.5

线性系统的频域分析 一、实验目的 1．掌握用MATLAB 语句绘制各种频域曲线。 2．掌握控制系统的频域分析方法。 二、实验内容 1．典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω，1.0=ζ，0.3，0.5，0.8，2的bode 图，记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w) bode(num,den3,w) bode(num,den4,w) bode(num,den5,w)

-100-80-60-40-200 20M a g n i t u d e (d B )10 -2 10 -1 10 10 1 10 2 10 3 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析：随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2．系统的开环传递函数为 ) 5)(15(10 )(2 +-= s s s s G ) 106)(15() 1(8)(22++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100); nyquist(num1,den1,w)

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

理工大学信号与系统实验报告连续时间系统的 复频域分析 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)

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

学生实验报告实验课程:信号与 系统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]);

自动控制原理线性系统的频域分析实验四

武汉工程大学实验报告专业电气自动化班号指导教师 姓名同组者无

M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Frequency (rad/sec) 当3.0=ζ时，程序如下： num=[0 0 36];den=[1 3.6 36];w=logspace(-2,3,100);bode(num,den,w) grid M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Bode Diagram Frequency (rad/sec) 当5.0=ζ时，程序如下： num=[0 0 36];den=[1 6 36];w=logspace(-2,3,100);bode(num,den,w) grid

M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Frequency (rad/sec) 当8.0=ζ时，程序如下： num=[0 0 36];den=[1 9.6 36];w=logspace(-2,3,100);bode(num,den,w) grid M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 当2=ζ时，程序如下： num=[0 0 36];den=[1 24 36];w=logspace(-2,3,100);bode(num,den,w) grid

连续信号的频域分析

第四章 连续信号的频域分析 将信号分解为若干不同频率的正弦信号或虚指数信号，实质上是将信号在频率域上进行分解，因此根据这种基本思想对信号和系统的分析称为频域分析。这种分解过程是通过傅里叶级数和傅里叶变换这一数学工具来实现的。 本章首先介绍连续信号的傅里叶级数和傅里叶变换，熟悉信号频谱的概念。 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章_用MATLAB进行控制系统频域分析

第5章 用MATLAB 进行控制系统频域分析 一、基于MATLAB 的线性系统的频域分析基本知识 （１）频率特性函数)(ωj G 。 设线性系统传递函数为： n n n n m m m m a s a s a s a b s b s b s b s G ++???++++???++=---1101110)( 则频率特性函数为： n n n n m m m m a j a j a j a b j b j b j b jw G ++???++++???++=---)()()()()()()(1101110ωωωωωω 由下面的MATLAB 语句可直接求出G(jw)。 i=sqrt(-1) % 求取-1的平方根 GW=polyval(num ，i*w)./polyval(den ，i*w) 其中（num ，den ）为系统的传递函数模型。而w 为频率点构成的向量，点右除（./）运算符表示操作元素点对点的运算。从数值运算的角度来看，上述算法在系统的极点附近精度不会很理想，甚至出现无穷大值，运算结果是一系列复数返回到变量GW 中。 （２）用MATLAB 作奈魁斯特图。 控制系统工具箱中提供了一个MATLAB 函数nyquist( )，该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时，nyquist （）函数仅在屏幕上产生奈氏图，命令调用格式为： nyquist(num,den) nyquist(num,den,w) 或者 nyquist(G) nyquist(G,w) 该命令将画出下列开环系统传递函数的奈氏曲线： ) () ()(s den s num s G = 如果用户给出频率向量w,则w 包含了要分析的以弧度/秒表示的诸频率点。在这些频率点上，将对系统的频率响应进行计算，若没有指定的w 向量，则该函数自动选择频率向量进行计算。 w 包含了用户要分析的以弧度/秒表示的诸频率点,MATLAB 会自动计算这些点的频率响应。 当命令中包含了左端的返回变量时，即: [re,im,w]=nyquist(G) 或

自动控制原理实验六线性系统的频域分析

实验六 线性系统的频域分析 一. 实验目的 （1）熟练掌握使用MA TLAB 命令绘制控制系统Nyquist 图的方法； （2）能够分析控制系统Nyquist 图的基本规律； （3）加深理解控制系统乃奎斯特稳定性判据的实际应用； （4）学会利用奈氏图设计控制系统； （5）熟练掌握运用MA TLAB 命令绘制控制系统伯德图的方法； （6）了解系统伯德图的一般规律及其频域指标的获取方法； （7）熟练掌握运用伯德图分析控制系统稳定性的方法； （8）设计超前校正环节并绘制Bode 图； （9）设计滞后校正环节并绘制Bode 图。 二. 实验原理及内容 1、频率特性函数)(ωj G 。 频率特性函数为： n n n n m m m m a j a j a j a b j b j b j b jw G ++???++++???++= ---)()()()()()()(1101110ωωωωωω 由下面的MATLAB 语句可直接求出G(jw)。 i=sqrt(-1) % 求取-1的平方根 GW=polyval(num ，i*w)./polyval(den ，i*w) 2、用MATLAB 作奈魁斯特图。 控制系统工具箱中提供了一个MA TLAB 函数nyquist( )，该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时，nyquist （）函数仅在屏幕上产生奈氏图，命令调用格式为： nyquist(num,den) ; 作Nyquist 图， nyquist(num,den,w); 作开环系统的奈氏曲线， 3、奈奎斯特稳定性判据（又称奈氏判据） 反馈控制系统稳定的充分必要条件是当ω从-∞变到∞时，开环系统的奈氏曲线不穿过点（-1，j0）且逆时针包围临界点（-1，j0）点的圈数R 等于开环传递函数的正实部极点数。 4、用MATLAB 作伯德图 控制系统工具箱里提供的bode()函数可以直接求取、绘制给定线性系统的伯德图。 命令的调用格式为： [mag,phase,w]=bode(num,den) [mag,phase,w]=bode(num,den,w) 由于伯德图是半对数坐标图且幅频图和相频图要同时在一个绘图窗口中绘制，因此，要用到半对数坐标绘图函数和子图命令。 （1） 对数坐标绘图函数 利用工作空间中的向量x ，y 绘图，要调用plot 函数，若要绘制对数或半对数坐标图，只需要用相应函数名取代plot 即可，其余参数应用与plot 完全一致。 （2） 子图命令

连续系统的频域分析

第三章傅立叶变换 时域分析：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时.

线性系统的频域分析-自动控制

实验三·线性系统的频域分析 一、实验目的 1．掌握用MATLAB 语句绘制各种频域曲线。 2．掌握控制系统的频域分析方法。 二、实验内容 1.典型二阶系统 2 22 ()2n n n G s s s ωζωω=++ 绘制出6n ω=，0.1ζ =，0.3，0.5，0.8，2的bode 图，记录并分析ζ对系统bode 图的影响。 2.系统的开环传递函数为 210 ()(51)(5)G s s s s =-+ 228(1) ()(15)(610) s G s s s s s += +++ 4(/31) ()(0.021)(0.051)(0.11) s G s s s s s += +++ 绘制系统的Nyquist 曲线、Bode 图和Nichols 图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。 3.已知系统的开环传递函数为21()(0.11) s G s s s += +。求系统的开环截止频率 穿越频率、幅值裕度和相位裕度。应用频率稳定判据判定系统的稳定性。 三、实验内容及分析 1. 系统1：2 22 ()2n n n G s s s ωζωω=++中6n ω=，（1）0.1ζ=时 Matlab 文本如下： num=[36 0 0]; den=[1 1.2 36]; w=logspace(-2,3,100); bode(num,den,w) Grid 得到图像：

同理，得到其他值情况下的波特图：ξ=0.3时 ξ=0.5时 ξ=0.8时

ξ=2时 从上面的图像中可以看出：随着ξ的不断增大，波特图中震荡的部分变得越来越平滑。而且，对幅频特性曲线来说，其上升的斜率越来越慢；对相频特性曲线来说，下降的幅度也在变缓。 2. 开环传递函数1：210 ()(51)(5) G s s s s = -+ 奈奎斯特图函数及图像如下： num=[0 10]; den=[conv([5,-1],[1,5]),0,0]; [z,p,k]=tf2zp(num,den); p

实验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;

实验二_连续时间信号的频域分析

实验二 连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义； 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质； 5、学习掌握利用Matlab 语言编写计算CTFS 、CTFT 和DTFT 的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT 、DTFT 的若干重要性质。 基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用Matlab 编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS 分析 任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。 三角傅里叶级数为： ∑∞ =++=1000)]sin()cos([)(k k k t k b t k a a t x ωω 2.1 或： ∑∞ =++ =1 00)cos()(k k k t k c a t x ?ω 2.2 其中1 02T π ω= ，称为信号的基本频率（Fundamental frequency ），k k b a a ，和，0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度，k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率0ωk 的函数，绘制出它们与0ωk 之间的图像，称为信号的频谱图（简称“频谱”），k c －0ωk 图像为幅度谱，k ?－0ωk 图像为相位谱。 三角形式傅里叶级数表明，如果一个周期信号x(t)，满足狄里克利条件，就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的正弦信号所组成，其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component)，其幅度（amplitude ）为k c 。

控制系统时域与频域性能指标的联系

控制系统时域与频域性能指标的联系 经典控制理论中，系统分析与校正方法一般有时域法、复域法、频域法。时域响应法是一种直接法，它以传递函数为系统的数学模型，以拉氏变换为数学工具，直接可以求出变量的解析解。这种方法虽然直观，分析时域性能十分有用，但是方法的应用需要两个前提，一是必须已知控制系统的闭环传递函数，另外系统的阶次不能很高。 如果系统的开环传递函数未知，或者系统的阶次较高，就需采用频域分析法。频域分析法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法，而且当系统的数学模型未知时，还可以通过实验的方法建立。此外，大量丰富的图形方法使得频域分析法分析高阶系统时，分析的复杂性并不随阶次的增加而显著增加。 在进行控制系统分析时，可以根据实际情况，针对不同数学模型选用最简洁、最合适的方法，从而使用相应的分析方法，达到预期的实验目的。 系统的时域性能指标与频域性能指标有着很大的关系，研究其内在联系在工程中有着很大的意义。 一、系统的时域性能指标 延迟时间t d 阶跃响应第一次达到终值h (∞)的50%所需的时间 上升时间 t r 阶跃响应从终值的10%上升到终值的90%所需的时间；对有振荡的系 统，也可定义为从0到第一次达到终值所需的时间 峰值时间t p 阶跃响应越过终值h (∞)达到第一个峰值所需的时间 调节时间 t s 阶跃响应到达并保持在终值h (∞)的±5%误差带内所需的最短时间 超调量%σ 峰值h( t p )超出终值h (∞)的百分比，即 %σ= () ()() ∞∞-h h h t p ?100% 二、系统频率特性的性能指标 采用频域方法进行线性控制系统设计时，时域内采用的诸如超调量，调整时间等描述系统性能的指标不能直接使用，需要在频域内定义频域性能指标。

自动控制原理线性系统的频域分析实验报告

实验四 专业 自动化 班号 03班 指导教师 陈艳飞 姓名 胡波 实验名称 线性系统的频域分析 实验日期 第 次实验 一、实验目的 1．掌握用MATLAB 语句绘制各种频域曲线。 2．掌握控制系统的频域分析方法。 二、实验内容 1．典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω，1.0=ζ，0.3，0.5，0.8，2的bode 图，记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w) bode(num,den3,w) bode(num,den4,w) bode(num,den5,w)

-100-80-60-40-200 20M a g n i t u d e (d B )10 -2 10 -1 10 10 1 10 2 10 3 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析：随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2．系统的开环传递函数为 ) 5)(15(10 )(2+-= s s s s G ) 106)(15() 1(8)(22++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图和Nichols 图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100); nyquist(num1,den1,w)

连续时间信号的频域分析(信号与系统课设).

福建农林大学计算机与信息学院 信息工程类 课程设计报告 课程名称：信号与系统 课程设计题目：连续时间信号的频域分析 姓名： 系：电子信息工程 专业：电子信息工程 年级：2008 学号： 指导教师： 职称： 2011 年 1 月10 日

福建农林大学计算机与信息学院信息工程类 课程设计结果评定

目录 1课程设计的目的 (1) 2课程设计的要求 (1) 3课程设计报告内容.....................................................................1-13 3.1连续信号的设计..................................................................1-11 3.2验证傅里叶变换的调制定理 (11) 3.3周期信号及其频谱 (12) 4总结 (13) 参考文献 (14)

连续时间信号的频域分析 1.课程设计的目的 （1）熟悉MATLAB语言的编程方法及MATLAB指令； （2）掌握连续时间信号的基本概念； （3）掌握门函数、指数信号和抽样信号的表达式和波形； （4）掌握连续时间信号的傅里叶变换及其性质； （5）掌握连续时间信号频谱的概念以及幅度谱、相位谱的表示； （6）掌握利用MATLAB进行信号的傅里叶变换以及时域波形和频谱的表示；（7）通过连续时间信号的频域分析，更深刻地理解了连续时间信号的时域和频域间的关系，加深了对连续时间信号的理解。 2.课程设计的要求 （1）自行设计以下连续信号：门函数、指数信号和抽样信号。要求：（a）画出以上信号的时域波形图； （b）实现以上信号的傅里叶变换，画出以上信号的幅度谱及相位谱，并对相关结果予以理论分析； （c）对其中一个信号进行时移和尺度变换，分别求变换后信号的傅里叶变换，验证傅里叶变换的时移和尺度变换性质。 （2）自行设计信号，验证傅里叶变换的调制定理。 （3）自行设计一个周期信号，绘出该信号的频谱，并观察周期信号频谱的特点。 3.课程设计报告内容 3.1（a）①门函数(矩形脉冲)： MATLAB中矩形脉冲信号用rectpuls函数表示： y=rectpuls (t,width) %width缺省值为1 >> t=-2:0.001:2; T=2; yt=rectpuls (t,T); plot(t,yt); axis([-2,2,0,1.5]); grid on; %显示格线