第五章 线性系统的频域分析法习题

501

第五章 线性系统的频域分析法

5-1 设闭环系统稳定，闭环传递函数为)(s Φ，试根据频率特性的定义证明：系统输入信号为余弦函数)cos()(φω+=t A t r 时，系统的稳态输出为

)](cos[|)(|)(ωφωωj t j A t c ss Φ∠++Φ=。

证明：根据三角定理，输入信号可表示为 )90sin()( ++=φωt A t r ，

根据频率特性的定义，有 ]90)(sin[|)(|)( +Φ∠++Φ=ωφωωj t j A t c ss ， 根据三角定理，得证： )](cos[|)(|)(ωφωωj t j A t c ss Φ∠++Φ=。

5-2 若系统的单位阶跃响应

t t e e t c 948.08.11)(--+-=，

试确定系统的频率特性。 解：s s s s C 1

361336)(2++=

，36

1336)(2++=s s s G ，)9)(4(36)(ωωωj j j G ++=；

2

/122/12)

81()16(36

|)(|ωωω++=j G ，9arctan 4arctan )(ωωω--=∠j G 。 或：)(2.7)()(94t t e e t c

t g ---== ；36

1336

)]([)(2

++==s s t g L s G ； 5-3 设系统如下图所示，试确定输入信号

)452cos()30sin()(

--+=t t t r

作用下，系统的稳态误差)(t e ss 。

解：2

1)(++=Φs s s e ； )452sin()30sin()(

+-+=t t t r

6325.0|)(|=Φj e ， 4.186.2645)(=-=Φ∠j ；

7906.0|)2(|=Φj e ， 4.18454.63)2(=-=Φ∠j ； 答案：)4.632sin(7906.0)4.48sin(6325.0)( +-+=t t t e ss 。

5-4 典型二阶系统的开环传递函数

)

2()(2

n n

s s s G ωζω+=

， 当取t t r sin 2)(=时，系统的稳态输出为

)45sin(2)( -=t t c ss ，

试确定系统参数n ω和ζ。

解：2

222)(n

n n

s s s ωζωω++=Φ； 1]

4)1[(2

2222=+-n n n

ωζωω，

451

2arctan

2

-=--n n ωζω； 122

-=n n ωζω，

答案：414.12==n ω，3536.04/2==ζ。

502

5-5 已知系统开环传递函数

)

1()

1()(2

++=

Ts s s K s G τ，0,,>T K τ， 试分析并绘制T >τ和τ>T 情况下的概略幅相曲线。 解：

其中

2/12231)(2-+=ττT K A ；T K A ττ=2；2/1223)(τ+=T KT A ；

)/arctan(451τφT -= ；]))((5.0arctan[2/1--=T T m ττφ；

参考：ωτω

ωωωτωωωτωωτωωω)

()1()()1()1()1()1()(2

022222222-+-+-++--=+-+=?→T K j K T T K j T T K jT j K j G 。 5-6 已知系统开环传递函数

)

2)(1(1

)(++=

s s s s G v

， 试分别绘制4,3,2,1=v 时的概略开环幅相曲线。

解：∞=|)0(|j G ， 90)0(?-=∠v j G ；0|)(|=∞j G ， 90)2()(?+-=∞∠v j G ；

2/122/12)4()1(|)(|---++=ωωωωv j G 和ωωω5.0arctan arctan 90)(--?-=∠ v j G 都是递减函数。所有幅相曲线的终止相角均小于起始相角180o ，以 90)2(?+-

v 趋于原点。

当1=v 时，有22

=x ω，204.0|)(|=x j G ω，与负实轴有交点)0,204.0(j -。

5-7

已知系统开环传递函数

)

1()

1(

)(12

++-=

s T s s T K

s G ，0,,21>T T K ，

当取1=ω时， 180)(-=∠ωj G ，5.0|)(|=ωj G 。当输入为单位速度信号时，系统的稳态误差为 0.1。试写出)(ωj G 的表达式。

解：据题义有下列结果，

503

10=K ； 180arctan 90arctan 12-=---T T ；2/1212/122)1(5.0)1(10T T +=+；

90)]1/()arctan[(2121=-+T T T T ，121=T T ；201=T ，05.02=T 。

所求的表达式为 )

201()

05.01(10)(ωωωωj j j j G +-=

。 5-8 已知系统开环传递函数

)

15.0)(12(10

)(2

+++=

s s s s s G ， 试分别计算5.0=ω和2=ω时，开环频率特性的幅值|)(|ωj G 和相位)(ωj G ∠。 解：5.0=ω，

89.17791

.0414.15.010

|)(|=??=ωj G ， 4.1534.184590)(-=---=∠ωj G ；

2=ω，

383.0162

.3123.4210

|)(|=??=

ωj G ， 6.3274.181800.7690)(-=+---=∠ωj G 。

5-9 已知系统开环传递函数

)

125.0)(1(10

)(2

++=

s s s s G ， 试绘制系统的概略开环幅相曲线。 解：{参考：

5

.22)

(2

+?→ωωj j G }

5-10 已知系统开环传递函数

?

?

? ??++??? ??++=

131911211

)(2s s s s s s G ，

选择频率点，列表计算对应的幅值与相位，绘制对数幅频特性曲线和相频特性曲线。 解：(过于烦琐，绘制渐近幅频特性)

504

5-11 绘制下列开环传递函数的对数渐近幅频特性曲线：

(1))

18)(12(2

)(++=

s s s G ；

(2))

110)(1(200

)(2

++=

s s s s G ； (3))

12/)(1()11.0/(8)(2

++++=s s s s s s G ； (4))11.0/)(1()

110/400/(10)(2++++=s s s s s s G 。 解：(1) 125.01=ω，5.02=ω； (2) 1.01=ω，12=ω；

(3) 1.01=ω，12=ω，23=ω； (4) 1.01=ω，12=ω，203=ω；

5-12 已知最小相位系统的对数渐近幅频特性如下试确定系统的开环传递函数。

解：(a) )1)(1()

1()(312+++=s T s T s T K s G ；01.01003==T K ，100

1.012==T T ；

)101.0)(1100()

11.0(100)(+++=s s s s G 。

(b) )

1()1()(22

1++=

s T s s T K s G ；

10

100

1==ωK ,101002=ω，

00316

.0316

.021==T T ；

)

100316.0()

1316.0(100)(2++=s s s s G 。

(c) )

1)(12()(212212

+++=s T s T s T Ks s G ζ；05.010==ζK ，1.0121==T T ；

)

11.0)(11.0(10)(22

+++=s s s s s G 。

505

5-13 试用Nyquist 稳定判据判断题5-5、5-6系统的稳定性。 解：题5-5中，0=P ；T >τ时，

Nyquist 曲线Γ

G 不包围临界点，系统稳定； τ>T 时， Nyquist 曲线Γ

G 包围临界点，系统不稳定。 题5-6中，0=P ；1=v 时， Nyquist 曲线Γ

G 不包围临界点，系统稳定； 4,3,2=v 时， Nyquist 曲线Γ

G 包围临界点，系统不稳定。

5-14 已知下列系统的开环传递函数(所有参数均大于0)

(1) )1)(1)(1()(321+++=s T s T s T K s G ； (2) )

1)(1()(21++=s T s T s K

s G ；

(3) )1()(2+=Ts s K

s G ； (4) )

1()1()(221++=s T s s T K s G ；

(5) 3)(s

K

s G =； (7) )1)(1)(1)(1()1)(1()(432165++++++=s T s T s T s T s s T s T K s G ；

(6) 3

21)1)(1()(s

s T s T K s G ++=； (8) 1)(-=Ts K

s G ； (9) 1

)(+--=

Ts K

s G ；

(10) )

1()(-=

Ts s K

s G 。

及其对应的幅相曲线分别如下图所示，应用Nyquist 稳定判据判断各系统的稳定性，若闭环系统不稳定指出系统在S 平面右半部的闭环极点数。 解：(1)0=P ，2-=R ，2=Z ； 不稳定；

(2)0=P ，0=R ； 稳定； (3)0=P ，2=R ，2=Z ； 不稳定； (4)0=P ，0=R ； 稳定； (5)0=P ，2-=R ，2=Z ； 不稳定； (6)0=P ，0=R ； 稳定； (7)0=P ，0=R ；

稳定； (8)1=P ，1=R ；

稳定；

(9)1=P ，0=R ，1=Z ；

不稳定；

(10)1=P ，1-=R ，2=Z ； 不稳定。

注：第(6)小题的幅相曲线未包围临界点。应用劳斯稳定判据能够说明闭环系统是稳定的：图中)(ωj G 曲线与负实轴交点处2/1211)(-=T T ω，且1|)(|1>ωj G ，得到1)(2121>+T T T KT 。 5-15 试用Nyquist 稳定判据判断题5-9系统的稳定性。 解：0=P ，2-=R ，闭环系统不稳定。 5-16 已知系统开环传递函数

)

1)(1()(++=

s Ts s K

s G ，0,>T K ，

506

试用Nyquist 稳定判据判断系统闭环稳定条件：

(1)2=T 时，K 值的范围； (2)10=K 时，T 值的范围； (3)K 、T 值的范围。

解：(计算)(ωj G 与虚轴的交点是解该题的要点，即计算临界稳定条件)

180arctan arctan 90)(-=---=∠x x x T j G ωωω，1|)(|=x j G ω；

2/1-=T x ω，1)1/(=+T KT ；

(1)2=T 时，5.10<

5-17 试用对数稳定判据判定题5-10系统的闭环稳定性。

解：采用对数频率稳定判据判，0=P ，且在0|)(|log 20>ωj G 区，相频曲线未穿越

180-线，闭

环系统稳定。

采用稳定裕度判断，

19

/)9/1(25.011|)(|22

222

=+-++=

c

c

c

c c c j G ωωωωωω；

7.1=c ω， 3.69)(180=∠+=c j G ωγ；

180)]81/(27arctan[5.0arctan 90arctan |)(2

-=----=∠g g g g g j G ωωωωω；

90106162)134)(29162(arctan 4

222=-+--g g g

g g ωωωωω，解得，36855.1050688.1072

==g

g ωω， 0118.048.1140.737.1042.10|)(|=??=

g j G ω；56.380118

.01

log 20log 20==g K ；

最小相位系统， 0>γ且0log 20>g K ，闭环系统稳定。

5-18 已知两个最小相位系统开环对数相频特性曲线如图所示。试分别确定系统的稳定性。鉴于改变

系统开环增益可使系统剪切频率变化，试确定闭环系统稳定时，剪切频率

ω的范围。

解：右图：闭环系统稳定；10ωω<

1

)(8.0+=-s Ke

s G s

，

试确定系统稳定时的K 值范围。

507

解：计算临界点，πωωω-=--=∠c c c j G arctan 8.0)(，1)1/(|)(|2

/12=+=c c c K j G ωω；

44822.2=c ω，64458.2=c K ；

使闭环系统稳定的K 值范围：64.20<

4

2)

1(5)(+=-s e s s G s

τ， 试确定闭环系统稳定时的τ值范围。

解：计算临界点，1)1/(5|)(|2

22=+=c c c j G ωωω，πωωτω-=--=∠c c c c j G arctan

4)(； 618.0=c ω，)618.1(=c ω， 5.1=c τ。

使闭环系统稳定时的τ值范围：5.10<<τ。

5-21 设单位反馈系统的传递函数为

21

)(s

s a s G +=

， 试确定相角裕度为45o 时参数a 的值。

解： 45180arctan 180=-+=c a ωγ，1=c a ω；

1/)1(22/122=+c c a ωω，2/12

2=c ω；

8409.0=a 。

5-22 对于典型二阶系统，已知参数3=n ω，7.0=ζ，试确定剪切频率c ω和相角裕度γ。

解：典型二阶系统的开环传递函数)

2()(2

n n

s s s G ζωω+=；

据1|)(|=c j G ω，01)(4)(2

24=-+n

c n c ωωζωω；

2/124)214(ζζωω-+=n c ；})214(2arctan{2/124ζζζγ++=；

答案：94455.1=c ω， 16.65=γ。

5-23 对于典型二阶系统，已知%15=σ，s t s 3=，试计算剪切频率c ω和相角裕度γ。 解：15.0])1/(exp[2/12=--=ζζπσ，517.0=ζ；3)/(3==n s t ζω，934.1=n ω；

应用5-22题的结果：2/124)214(ζζωω-+=n c ；})214(2arctan{2

/124ζζζγ++=；

答案：497.1=c ω， 18.53=γ。

5-24 根据题5-11所绘渐近幅频特性曲线，近似确定剪切频率c ω，并由此确定相角裕度γ的近似值。 解：(1)0)125.0/log(202log 20=-c ω，25.0=c ω， 90=γ；

精确计算，196.0=c ω， 11.101=γ；

(2)0)1/log(60)001.0/200log(20=-c ω，115.2=c ω， 98.151-=γ；

精确计算，06.2=c ω， 33.151-=γ；

(3)0)2/log(60)1/2log(40)1.0/8log(20=--c ω，429.5=c ω， 60-=γ；

精确计算，34.5=c ω， 54.60-=γ；

508

(4) 0)1.0/log(40)1.0/10log(20=-c ω，1=c ω， 57.33-=γ；

精确计算，87.0=c ω， 49.29-=γ。

5-25 某一控制系统，若要求%17=σ，s t s 5=，试由近似公式确定频域指标要求c ω和γ。 解：(1) 最简处理方案，按二阶系统处理。

17.0])1/(exp[2/12=--=ζζπσ，491.0=ζ；5)/(3==n s t ζω，222.1=n ω；

2/124)214(ζζωω-+=n c ；})214(2arctan{2/124ζζζγ++=；

968.0=c ω， 1.52=γ。

(2) 按高阶系统估算。无相应的超调量估算公式。

)1(4.016.0-+≈r M σ；γsin /1≈r M ；c r ωω≈；

})1(5.2)1(5.12{2-+-+=r r c

s M M t ωπ；

17.0)1(4.016.0=-+r M ，25.1=r M ；59.1=c ω， 1.53=γ。

第5章频域分析法习题解答

第5章频域分析法 学习要点 1 频率特性的概念，常用数学描述与图形表示方法； 2 典型环节的幅相频率特性与对数频率特性表示及特点； 3 系统开环幅相频率特性与对数频率特性的图示要点； 4 应用乃奎斯特判据判断控制系统的稳定性方法； 5 对数频率特性三频段与系统性能的关系； 6 计算频域参数与性能指标； 思考与习题祥解 题判断下列概念的正确性 ω的正弦信号加入线性系统，这个系统的稳态输出也将是同 (1) 将频率为 一频率的。 M仅与阻尼比ξ有关。 (2) 对于典型二阶系统，谐振峰值 p (3) 在开环传递函数中增加零点总是增加闭环系统的带宽。 (4) 在开环传递函数中增加极点通常将减少闭环系统的带宽并同时降低稳定性。 (5) 对于最小相位系统，如果相位裕量是负值，闭环系统总是不稳定的。 (6) 对于最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (7) 对于最小相位系统，如果幅值裕量是负分贝值，闭环系统总是不稳定的。 (8) 对于非最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (9) 对于非最小相位系统，须幅值裕量大于1且相位裕量大于0，闭环系统才是稳定的。 (10) 相位穿越频率是在这一频率处的相位为0。 (11) 幅值穿越频率是在这一频率处的幅值为0dB。 (12) 幅值裕量在相位穿越频率处测量。 (13) 相位裕量在幅值穿越频率处测量。 (14) 某系统稳定的开环放大系数25 K<，这是一个条件稳定系统。 (15) 对于（-2/ -1/ -2)特性的对称最佳系统，具有最大相位裕量。 (16) 对于（-2/ -1/ -3)特性的系统，存在一个对应最大相位裕量的开环放大系数值。 (17) 开环中具有纯时滞的闭环系统通常比没有时滞的系统稳定性低些。 (18) 开环对数幅频特性过0分贝线的渐近线斜率通常表明了闭环系统的相对稳定性。 M和频带宽BW (19) Nichols图可以用于找到一个闭环系统的谐振峰值 p 的信息。

大作业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 即

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

实验名称： 控制系统的频域分析 实验类型：________________同组学生姓名：__________ 一、实验目的和要求 用计算机辅助分析的方法，掌握频率分析法的三种方法，即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 图，并判断稳定性。 （三）实验要求

第五章 频域分析法

第五章 频域分析法 时域分析法具有直观、准确的优点。如果描述系统的微分方程是一阶或二阶的，求解后可利用时域指标直接评估系统的性能。然而实际系统往往都是高阶的，要建立和求解高阶系统的微分方程比较困难。而且，按照给定的时域指标设计高阶系统也不是一件容易的事。 本章介绍的频域分析法，可以弥补时域分析法的不足。因为频域法是基于频率特性或频率响应对系统进行分析和设计的一种图解方法，故其与时域分析法相比有较多的优点。首先，只要求出系统的开环频率特性，就可以判断闭环系统是否稳定。其次，由系统的频率特性所确定的频域指标与系统的时域指标之间存在着一定的对应关系，而系统的频率特性又很容易和它的结构、参数联系起来。因而可以根据频率特性曲线的形状去选择系统的结构和参数，使之满足时域指标的要求。此外，频率特性不但可由微分方程或传递函数求得，而且还可以用实验方法求得。对于某些难以用机理分析方法建立微分方程或传递函数的元件(或系统)来说，具有重要的意义。因此，频率法得到了广泛的应用，它也是经典控制理论中的重点内容。 5.1 频率特性 对于线性定常系统，若输入端作用一个正弦信号 t U t u ωsin )(= (5—1) 则系统的稳态输出y(t)也为正弦信号，且频率与输人信号的频率相同，即 ) t Y t y ?ω+=sin()( (5—2) u(t)和y(t)虽然频率相同，但幅值和相位不同，并且随着输入信号的角频率ω的改变，两者之间的振幅与相位关系也随之改变。这种基于频率ω的系统输入和输出之间的关系称之为系统的频率特性。 不失一般性，设线性定常系统的传递函数G(s)可以写成如下形式 ) () () () () ())(() ()()()(1 21s A s B p s s B p s p s p s s B s U s Y s G n j j n = +=+++== ∏=Λ (5—3) 式中B(s)——传递函数G(s)的m 阶分子多项式，s 为复变量； A(s)——传递函数G(s)的n 阶分母多项式 (n ≥m)； n p p p ---,,,21Λ—传递函数G(s)的极点，这些极点可能是实数，也可能是复数，对稳定的系统采说，它们都应该有负的实部。 由式(5—1)，正弦输入信号u(t)的拉氏变换为(查拉氏变换表) ) )(()(22ωωω ωωj s j s U s U s U -+=+= （5—4）

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

第五章 线性系统的频域分析法习题

501 第五章 线性系统的频域分析法 5-1 设闭环系统稳定，闭环传递函数为)(s Φ，试根据频率特性的定义证明：系统输入信号为余弦函数)cos()(φω+=t A t r 时，系统的稳态输出为 )](cos[|)(|)(ωφωωj t j A t c ss Φ∠++Φ=。 证明：根据三角定理，输入信号可表示为 )90sin()( ++=φωt A t r ， 根据频率特性的定义，有 ]90)(sin[|)(|)( +Φ∠++Φ=ωφωωj t j A t c ss ， 根据三角定理，得证： )](cos[|)(|)(ωφωωj t j A t c ss Φ∠++Φ=。 5-2 若系统的单位阶跃响应 t t e e t c 948.08.11)(--+-=， 试确定系统的频率特性。 解：s s s s C 1 361336)(2++= ，36 1336)(2++=s s s G ，)9)(4(36)(ωωωj j j G ++=； 2 /122/12) 81()16(36 |)(|ωωω++=j G ，9arctan 4arctan )(ωωω--=∠j G 。 或：)(2.7)()(94t t e e t c t g ---== ；36 1336 )]([)(2 ++==s s t g L s G ； 5-3 设系统如下图所示，试确定输入信号 )452cos()30sin()( --+=t t t r 作用下，系统的稳态误差)(t e ss 。 解：2 1)(++=Φs s s e ； )452sin()30sin()( +-+=t t t r 6325.0|)(|=Φj e ， 4.186.2645)(=-=Φ∠j ； 7906.0|)2(|=Φj e ， 4.18454.63)2(=-=Φ∠j ； 答案：)4.632sin(7906.0)4.48sin(6325.0)( +-+=t t t e ss 。 5-4 典型二阶系统的开环传递函数 ) 2()(2 n n s s s G ωζω+= ， 当取t t r sin 2)(=时，系统的稳态输出为 )45sin(2)( -=t t c ss ， 试确定系统参数n ω和ζ。 解：2 222)(n n n s s s ωζωω++=Φ； 1] 4)1[(2 2222=+-n n n ωζωω， 451 2arctan 2 -=--n n ωζω； 122 -=n n ωζω， 答案：414.12==n ω，3536.04/2==ζ。

第5章频域分析法习题解答

第5章频域分析法 5.1 学习要点 1 频率特性的概念，常用数学描述与图形表示方法； 2 典型环节的幅相频率特性与对数频率特性表示及特点； 3 系统开环幅相频率特性与对数频率特性的图示要点； 4 应用乃奎斯特判据判断控制系统的稳定性方法； 5 对数频率特性三频段与系统性能的关系； 6 计算频域参数与性能指标； 5.2 思考与习题祥解 题5.1 判断下列概念的正确性 ω的正弦信号加入线性系统，这个系统的稳态输出也将是同 (1) 将频率为 一频率的。 M仅与阻尼比ξ有关。 (2) 对于典型二阶系统，谐振峰值 p (3) 在开环传递函数中增加零点总是增加闭环系统的带宽。 (4) 在开环传递函数中增加极点通常将减少闭环系统的带宽并同时降低稳定性。 (5) 对于最小相位系统，如果相位裕量是负值，闭环系统总是不稳定的。 (6) 对于最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (7) 对于最小相位系统，如果幅值裕量是负分贝值，闭环系统总是不稳定的。 (8) 对于非最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (9) 对于非最小相位系统，须幅值裕量大于1且相位裕量大于0，闭环系统才是稳定的。 (10) 相位穿越频率是在这一频率处的相位为0。 (11) 幅值穿越频率是在这一频率处的幅值为0dB。 (12) 幅值裕量在相位穿越频率处测量。 (13) 相位裕量在幅值穿越频率处测量。 (14) 某系统稳定的开环放大系数25 K<，这是一个条件稳定系统。 (15) 对于（-2/ -1/ -2)特性的对称最佳系统，具有最大相位裕量。 (16) 对于（-2/ -1/ -3)特性的系统，存在一个对应最大相位裕量的开环放大系数值。 (17) 开环中具有纯时滞的闭环系统通常比没有时滞的系统稳定性低些。 (18) 开环对数幅频特性过0分贝线的渐近线斜率通常表明了闭环系统的相对稳定性。 M和频带宽BW (19) Nichols图可以用于找到一个闭环系统的谐振峰值 p 的信息。

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

理工大学信号与系统实验报告连续时间系统的 复频域分析 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）分子分母根的求解，调用格式如下：

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

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

系统频域分析课程设计报告

系统频域分析课程设计 报告 Company number【1089WT-1898YT-1W8CB-9UUT-92108】

《综合仿真》课程设计报告 姓名 学号 同组成员 指导教师 时间 11周至14周

系统的频域分析 【目的】 (1) 加深对系统频域分析基本原理和方法的理解。 (2) 加深对信号幅度调制与解调基本原理和方法的理解。 (3) 锻炼学生综合利用所学理论和技术，分析与解决工程实际 问题的能力。 【研讨内容】 题目1．幅度调制和连续信号的Fourier 变换 本题研究莫尔斯码的幅度调制与解调。本题中信号的形式为 )π2sin()()π2sin()()π2cos()()(132211t f t m t f t m t f t m t x ++= 其中信号x (t )由文件定义，可用命令Load ctftmod 将文件定义的变量装入系统内存。运行命令Load ctftmod 后，装入系统的变量有 af bf dash dot f1 f2 t x 其中 bf af ： 定义了一个连续系统H (s )的分子多项式和分母多项式。可利用freqs(bf,af,w)求出该系统的频率响应，也可用sys=tf(bf,af)得到系统的模型，从而用lsim 求出信号通过该系统的响应。 dash dot ： 给出了莫尔斯码中的基本信号dash 和dot 的波形 f1 f2： 载波频率 t: 信号x (t )的抽样点 x: 信号x (t )的在抽样点上的值 信号x (t )含有一段简单的消息。Agend 007的最后一句话是

The future of technology lies in ··· 还未说出最后一个字，Agend 007就昏倒了。你(Agend 008)目前的任务就是要破解Agend 007的最后一个字。该字的信息包含在信号x (t )中。信号x (t )具有式(1)的形式。式中的调制频率分别由变量f1和f2给出，信号m 1(t )，m 2(t )和m 3(t )对应于字母表中的单个字母，这个字母表已用国际莫尔斯码进行编码，如下表所示： (1)字母B 可用莫尔斯码表示为b=[dash dot dot dot]，画出字母B 莫尔 斯码波形； (2) 用freqs(bf,af,w)画出系统的幅度响应； (3) 利用lsim 求出信号dash 通过由sys=tf(bf,af)定义的系统响应，解释你所获得的结果； (4)用解析法推导出下列信号的Fourier 变换 )π2cos()π2cos()(21t f t f t m )π2sin()π2cos()(21t f t f t m

第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) 或

连续系统的频域分析

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

实验三线性系统的频域分析

自动控制理论 上 机 实 验 报 告 学院：机电工程学院 班级：13级电信一班

： 学号： 实验三 线性系统的频域分析 一、实验目的 1．掌握用MATLAB 语句绘制各种频域曲线。 2．掌握控制系统的频域分析方法。 二、基础知识及MATLAB 函数 频域分析法是应用频域特性研究控制系统的一种经典方法。它是通过研究系统对正弦信号下的稳态和动态响应特性来分析系统的。采用这种方法可直观的表达出系统的频率特性，分析方法比较简单，物理概念明确。 1．频率曲线主要包括三种:Nyquist 图、Bode 图和Nichols 图。 1）Nyquist 图的绘制与分析 MATLAB 中绘制系统Nyquist 图的函数调用格式为: nyquist(num,den) 频率响应w 的围由软件自动设定 nyquist(num,den,w) 频率响应w 的围由人工设定 [Re,Im]= nyquist(num,den) 返回奈氏曲线的实部和虚部向量， 不作图 例4-1:已知系统的开环传递函数为2 526 2)(2 3++++=s s s s s G ，试绘制Nyquist 图，并判断系统的稳定性。

num=[2 6]; den=[1 2 5 2]; [z,p,k]=tf2zp(num,den); p nyquist(num,den) 极点的显示结果及绘制的Nyquist 图如图4-1所示。由于系统的开环右根数P=0，系统的Nyquist 曲线没有逆时针包围（-1，j0）点，所以闭环系统稳定。 p = -0.7666 + 1.9227i -0.7666 - 1.9227i -0.4668 若上例要求绘制)10,10(32-∈ω间的Nyquist 图，则对应的MATLAB 语句为: num=[2 6]; den=[1 2 5 2]; w=logspace(-1,1,100); 即在10-1和101之间，产生100个等距 离的点 nyquist(num,den,w) 2）Bode 图的绘制与分析 系统的Bode 图又称为系统频率特性的对数坐标图。Bode 图有两图，分别绘制开环频率特性的幅值和相位与角频率ω的关系曲线，称为对数幅频特性曲线和对数相频特性曲线。 MATLAB 中绘制系统Bode 图的函数调用格式为: bode(num,den) 频率响应w 的围由软件自动设定 bode(num,den,w) 频率响应w 的围由人工设定 图4-1 开环极点的显示结果及Nyquist 图

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

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

控制系统时域与频域性能指标的联系 经典控制理论中，系统分析与校正方法一般有时域法、复域法、频域法。时域响应法是一种直接法，它以传递函数为系统的数学模型，以拉氏变换为数学工具，直接可以求出变量的解析解。这种方法虽然直观，分析时域性能十分有用，但是方法的应用需要两个前提，一是必须已知控制系统的闭环传递函数，另外系统的阶次不能很高。 如果系统的开环传递函数未知，或者系统的阶次较高，就需采用频域分析法。频域分析法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法，而且当系统的数学模型未知时，还可以通过实验的方法建立。此外，大量丰富的图形方法使得频域分析法分析高阶系统时，分析的复杂性并不随阶次的增加而显著增加。 在进行控制系统分析时，可以根据实际情况，针对不同数学模型选用最简洁、最合适的方法，从而使用相应的分析方法，达到预期的实验目的。 系统的时域性能指标与频域性能指标有着很大的关系，研究其内在联系在工程中有着很大的意义。 一、系统的时域性能指标 延迟时间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% 二、系统频率特性的性能指标 采用频域方法进行线性控制系统设计时，时域内采用的诸如超调量，调整时间等描述系统性能的指标不能直接使用，需要在频域内定义频域性能指标。

信号与系统报告 实验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);

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

实验四 专业 自动化 班号 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)

第三章连续系统的频域分析

习题三 31证明题图囂所示矩形函数/址）与冷“尬帆为整数｝在区间（0.2^ ）上正応 J 1 /W P -7T r 卷也J 3.2设了①的正交展开式为 /0 =養恥 是iiE 明f ⑴和护o ”6呵￡ ｝是11对应关系E [1 (「1)＜于＜2 0其他 II 试问函数组 苗⑦務②焉②爲②｝在（山4）区间上是否为正交函数值，是否为归一 牝正交函數组，是否为完备正交画数爼「并用它们的线t 删合精确的表示题團玄2所示函数 “） 9 /(i) 題要1 3 2 M4证明下列函数集在匕心*— 匡间上是正交函数集右肯任意一个正实数? \ 叫丿 （1）｛ cos^ivof, sinMw e f | M - 0,±1,±2^.,,）； ⑵｛*叫1沪蚣…｝h

3.5试求题因3.3所示信号的三角形傅立叶级数展开式，并画出频谙因。 1/w A A n 1,n[ 1 , ：J72?T t KS 3.3 3.6试求题图34所示周期信号的指数形傅立叶级数系鹽，并画出它的幅度谙。 3?己知剛函数前四分之一的周期的波形女廳图？.5所示.根据下列各恬况的要求，画出/(/)在一个周期(0*T)的波形? (1)/(f)是偶函数'只含有偶次谐波： (2)/(f)是偶函数，只含有奇次谐波； (3)/(f)是偶函数，含有偶次和奇次谐波； (4)/(f)是奇函数，只含有偶次谐波； (5)/(f)是奇函数，只含有奇次谐波； (6)/(f)是奇函数，含有偶次和奇次谐波.

3.8设是满足以下两个条件的周期函数：条件1 ： /(0 = -/(~0 ; 条件2：/a± j)= -/(o ? 试证明/(◎中只含有奇次谐波的正弦分星。 3.9设周期信号/(f)的指城傅立叶级数系数为尺，试证明缪的指数形傅立叶级数系 at 数为感(式中叫=亨). 3.10设有一周期信号/O) ＞其奇波频率为w。= X ,且/(f)的指数形傅立叶级数为 这里，丘“；阿|"/4 ；|^|=1/2 ; |^|= 1/3 o 试写出的三角形傅立叶级数表达式? 3.L1求题图3.6所示信号的傅立叶变换? 题图3.6

第四章 连续系统的频域分析例题详解

第四章 连续系统的频域分析例题详解 1．一带限信号的频谱图如下图1所示，若次信号通过图2所示系统，请画出A 、B 、C 三点处的信号频谱。理想低通滤波器的频率函数为 )15()15()(--+=ωεωεωj H ，如图3所示。 解：设A 处的信号为：A f ，B 处的信号为：B f ，C 处的信号为：C f )30cos()(t t f f A = )30cos(t f f A B = )]] 30([)]30([[2 1 )()]]30([)]30([[21 )(++-=++-= w j F w j F jw F w j F w j F jw F A A B A

1. 如图2（a ）所示的系统，带通滤波器的频率响应如图2（b ）所示，其相频特 性()0?ω=，若输入 sin(2) (),()cos(1000)2t f t s t t t π==，求输出信号()y t 。 f () H j ω()0 ?ω=1/(.) rad s ω--1001 -999 0 999 10011 -1000 1000 图（b ） 图2

解 4sin(2)1 ()[ ]()22 t F j F g t ωωπ== [cos(1000)][(1000)(1000)]F t πδωδω=++- 441 [()cos(1000)][()][cos(1000)]21 [(1000)(1000)] 4 F f t t F f t F t g g πωω= ?*=++- 则系统输出信号的傅里叶变换为 ()[()cos(1000)]()Y j F f t t H j ωω= 由()H j ω的波形图及相频特性可得 22()(1000)(1000)H j g g ωωω=++- 所以可得 2221 ()[(1000)(1000)] 4 1 ()[(1000)(1000)]4 Y j g g g ωωωωδωδω=++-=*++- 由此可得输出信号为 1 ()()cos(1000)2y t Sa t t π = 3．一理想低通滤波器的频率响应如图3示，其相频特性φ(ω)=0。若输入信号 t t t f ππ) sin()(= ，求输出信号的频谱函数，并画出其频谱图。 图 3 解：信号t t t f ππ) sin()(= 的频域表达式为 )(2)(2ωπωg j F =