第九章线性系统的状态空间分析与综合习题

第九章线性系统的状态空间分析与综合

9-1 已知电枢控制的直流司服电机的微分方程组及传递函数为

b a a a a a E t d di L i R u ++=，t d d K E m b b θ=，a m m i C M =，t d d f t d d J M m

m m m m θθ+=2

2； )]

()([)()(2

m b m a a m m a m a m

a m C K f R s R J f L s J L s C s U s ++++=Θ。 ⑴ 设状态变量m x θ=1，m x θ&=2，m x θ&&=3，输出量m

y θ=，试建立其动态方程； ⑵ 设状态变量a i x =1，m x θ=2，m x θ&=3，输出量m

y θ=，试建立其动态方程； ⑶ 设x T x =，确定两组状态变量间的变换矩阵T 。

解：⑴ 由传递函数得 a m a m m a m b m a m a u C x R J f L x C K f R x J L ++-+-=323)()(&，动态方程为

[]x

y u x x x x x x 001100010001032121321=??????????+????????????????????--=??????????αα&&&，其中)/()()/()()/(21m a a m m a m a m b m a m a a m J L R J f L J L C K f R J L u C u +=+==αα； ⑵ 由微分方程得

133

2311x f x C x J x x u x K x R x L m m m a b a a -==---=&&&，即 []x

y u x x x a a a a x x x a 0200010100032133311311321=????

?????+??????????????????

??=??????????&&&，其中 m

m m m a b a a J f a J C a L K a L R a ////33311311-==-=-=；

⑶ 由两组状态变量的定义，直接得到????

???????????????=??????????3213331

321010001

0x x x a a x x x 。

9-2 设系统的微分方程为

u x x x =++23&&&

其中u 为输入量，x 为输出量。

⑴ 设状态变量x x =1，x x &=2，试列写动态方程；

⑵ 设状态变换211x x x +=，2122x x x --=，试确定变换矩阵T 及变换后的动态方程。

解：⑴ u x x x x ??????+????????????--=??????1032102121&&，[]??

????=2101x x y ； ⑵ ??????=??????2121x x T x x ，??????--=2111T ；??

????--=-11121

T ；AT T A 1-=，B T B 1-=，CT C =；得，?

?????--=2111T ；u x x x x ??

????-+????????????-=??????1110012121&&，[]??????=2111x x y 。 9-3 设系统的微分方程为

u y y y y 66116=+++&&&&&&

其中u 、y 分别系统为输入、输出量。试列写可控标准型(即A 为友矩阵)及可观标准型(即A 为友矩

阵转置)状态空间表达式，并画出状态变量图。

解：可控标准型和可观标准型状态空间表达式依次为，

[]x y u x x 00610061161

00010=??????????+??????????---=&；[]x

y u x x 100

006610

1101600=????

?????+???????

???---=&；可控标准型和可观标准型的状态变量图依次为，

9-4 已知系统结构图如图所示，其状态变量为1x 、2x 、3x 。试求动态方程，并画出状态变量图。

解：由图中信号关系得，31x x =&，u x x x 232212+--=&，32332x x x -=&，1x y =。动态方程为

u x x ???

?????+??????????---=020*********&，[]x y 001；状态变量图为

9-5 已知双输入双-输出系统状态方程和输出方程

32132

1321

216

1162u x x x x u u x x u x x +---=-+=+=&&&，3

2122

112x x x y x x y -+=-=，写出其向量-矩阵形式并画出状态变量图。

解：状态方程 u x x ????

?????-+??????????---=1012016116100010&，

x y ??????--=112011；状态变量图为

9-6 已知系统传递函数为

6)(22++++=s s s s s G ，

试求出可控标准型(A 为友矩阵)、可观标准型(A 为友矩阵转置)、对角型(A 为对角阵)动态方程。

解：135

.015.113

452)(2++++=++++=s s s s s s G ；可控标准型、可观标准型和对角型依次为

[]u x y u x x +=??????+??????--=25104310&；[]u x y u x x +=??????+??????--=10254130&；[]u

x y u x x +=??

???+??????--=115.05.13001&。

9-7 已知系统传递函数为

)

2()1(5

)(2

++=

s s s G ，试求约当型(A 为约当阵)动态方程。

解：2

)1(5)1(525)(+++-++=s s s s G ；

u x x ????

?????-+??????????---=555100110002&，[]x y 011=。 9-8 已知矩阵

???

???=0001

100001000010

A , 试求A 的特征方程、特征值、特征向量，并求出变换矩阵将A 约当化。

解：特征方程0)I det(=-A s ，即014

=-s ；特征值11=λ、12-=λ、j =3λ、j -=4

λ；

特征向量依次对应矩阵1

-T 的列，所求变换矩阵为T ；

????????????------=

-j j j j T 111111*********；????????????------=j j j j T 111

11111111121；??

???

???????--==-j j TAT A 00

0000001000

11。 9-9 已知矩阵

???-=1001A ，试用幂级数法及拉普拉斯变换法求出矩阵指数(即状态转移矩阵)。解：幂级数法求解，

??????-=k k k

A 100)1(；??

???===Φ-∞=∑t t k k k t

A e e t A k e t 00!1)(0；拉普拉斯变换法求解，

??????-+=??????-+=---)1/(100)1/(11001)I (1

s s s s A s ；??

???=-=Φ---t t e e A s L s 0

0])I [()(11。 9-10 求下列状态方程的解：

x x ??????????---=30002000

1&。解：????????

?=Φ---t t t

e e e t 320

000

00)(，得到 ?????

??????????????=---)0()0()0(00000023132x x x e e e x t t t 。 9-11 已知系统的状态方程为

u x x ??

???+??????=111101&，初始条件为1)0(1=x ，0)0(2=x 。试求系统在单位阶跃输入作用下的响应。解法1：

?????

?=???? ????????---=Φ--t t t e te e s s L t 01101)(1

1； ??

????-=??????-+??????=???

?????????-+????????????=?---t t t t t t t t t t t t t te e te e te e d e e t e e te e x 212111)(00100τττττ。解法2：

???--=??????--+??

????--=

+-=-s s s s s s s s s s x s Bu A s s x 21)1(1

11)1(1

1)1(1

)}0()({)I ()(22

； ?

????-==-t t te e s x L x 212)]([1

。 9-12 已知系统的状态转移矩阵

???+-+---=Φ--------t t t

t t t t

t e e e e e e e e t 222232332223)(，

试求该系统的状态阵A 。

解：?

????--=Φ==4321)

t t A &。(注：原题给出的)(t Φ不满足A =Φ)0(&及A t t A t )()()(Φ=Φ=Φ&。) 9-13 已知系统动态方程

u x x ????

?????+??????????---=210311032010&，[]x y 100=，试求传递函数)(s G 。解：B A s C s G 1

)I ()(--=，

[][]??????

???????????

???++----+-----=?????

???????

??????--+-=-210231503620396710021031103201

100)(22231

s s s s s s s s s s s s s s s G ； 6

7372)(32

--++=s s s s s G 。

9-14 试求所示系统的传递函数矩阵。

u x x ??

?????-+??????????---=1012016116100010&，x y ??????--=112011。解：????

??????---+-++++++=?????

?????+--=---22

2231

1611666161166116161161001)I (s s s s s s s s s s s s s s s A s ； ?????

????-??????????---+-+++??????-+++=1012016116661611611201161161)(22

23s s s s s s s s s s s s s G ； ??

????-+-++--+++==442845429461161

)(2

223s s s s s s s s s s G 。 9-15 已知差分方程

)(3)1(2)(2)1(3)2(k u k u k y k y k y ++=++++，

试列写可控标准型(A 为友矩阵)离散动态方程，并求出1)(=k u 时的系统响应。给定0)0(=y ，1)1(=y 。

解：系统的脉冲传递函数为

2332)(2

+++=

z z z z G ，1)(-=z z

z U ；)(10)(3210)1(k u k x k x ??

????+??????--=+，[])(23)(k x k y =。 )2(32)1(2)1(65)2)(1)(1(32)}0(3)1()0()(){()(232

++++-=++-+=+++=z z z z z z z z z z z y z y z y z U z G z Y ；

322)1(65)(1

++-+=k k k y 。

9-16 已知连续系统动态方程为

u x x ??

???+??????=102010&，[]x y 01=，设采样周期s T 1=，试求离散化动态方程。解：设)()(k u t u =，T k t kT )1(+<≤；

??????--=??????--=---)2/(10)]2(/[1/1201)(1

s s s s s s A sI ，??

????-=Φt

t e e t 220)1(5.01)(；

??????-=Φ=Φ220)1(5.01)(e e T ，??

????--=???

???-Φ=Γ?)1(5.0)3(25.010)(2

20e e t d t T T ； )()1(5.0)3(25.0)(0)1(5.01)1(2

222k u e e k x e e k x ??

???--+??????-=+，[])(01)(k x k y =。 9-17 判断下列系统的状态可控性：

⑴ u x x ??

?????+??????????----=100041020122

&；

⑵ u x x ??

???

?????+??????????=010********

1&； ⑶ u x x ??

???

?????+??????????=011000010010011&； ⑷ u x x ??

?????+??????????--=121100040004&；

⑸ u x x ????????????+????????????=11100000000000012111λλλλ&； ⑹ u x x ????

??????+????????????=11000000000100012111λλλλ&。解：

⑴ ????

??????--=101000210U ，n U <=2rank ；状态不完全可控；

⑵ ??????????=210111210U ，n U <=2rank ；

状态不完全可控；

⑶ ????

??????=1101101010001U ，3rank 1=U ；状态完全可控；

⑷ ????

??????--=11132821641U ，n U <=2rank ；状态不完全可控；

⑸ ??????

?????

???=3222231

21131

21121

1111

3210λλλλλλλλλλλU ，n U <=3rank ；

状态不完全可控；

⑹ ??

???

???=322223

121121

11113210

3100λλλλλλλλλU ，4rank =U ；

状态完全可控；

9-18 已知bc ad =，试计算=??

???100

d c b a ？

解：矩阵A 的特征方程为 0)()(2

=+-=s d a s s α，据凯莱哈密尔定理得知：

0)(2=+-A d a A ，k k A d a A )(1+=+；A d a A 99100)(+=；

???+=??

????d c b a d a d c b a 99100

)(。

9-19 设系统状态方程为

u b x a x ??

???+??????-=1110&，且状态完全可控。试求a 、b 。解：??

??-=11ab b b U ，01det 2

≠--=b ab U ，只需b b a 1+≠。 9-20 设系统传递函数为

147)(23++++=

s s s a

s s G ，

且状态完全可控。试求a 。

解：可控标准型实现的系统，无论a 取何值，系统状态完全可控。在可观标准型实现中

u a x x ??????????+??????????---=017101401800&，[]x y 100=；????

??????---=71014180a a a U ， 08147det 2

3≠-+-=a a a U ；只需1≠a 、2≠a 且4≠a 。

注：由)(s G 分子和分母的多项式互质条件，同样得到081472

3≠-+-a a a 。

9-21 判断下列系统的输出可控性：

⑴ u x d c b a x

?????

??????+???????

?????----=1100000

00000000

0&，[]x y 0001=。 ⑵ u x x ??

?????+??????????---=100611

6100010

&，[]x y 001=；

解：输出可控性判别矩阵CU B A AB B C B CA CAB CB S n n o

===--][][11Λ

Λ。

⑴ ??

???

???----=32

3211000

00000d d d

c c c U ，]0000[=o S ，q S o <=0rank ，系统的输出不可控。

⑵ ????

??????=001010100U ，]100[=o S ，q S o ==1rank ，系统的输出可控； 9-22 判断下列系统的可观测性：

⑴u x x ??????????+??????????-----=102101110221&，[]x y 011=； ⑵x x ?????

???=130

020002

&，[]x y 111=； ⑶x x ????????????----=200

012000010001

1&，x y ??????-=01000001；⑷x x ????

??????-=30002001

2&，[]x y 110=。

解：应用可观测性判别矩阵。

⑴ ??

????????---=252131011

V ，3rank =V ；

系统完全可观测；

⑵ ????

??????=1134152111V ，3rank =V ；系统完全可观测；

⑶ ????

?????

???---=12000011010000012V ，4rank =V ；系统完全可观测；

⑷ ????

??????-=940320110V ，n V <=2rank ；

系统不完全可观测；

9-23 试确定使下列系统可观测的a 、b ：

x b a x ??

???=01&，[]x y 11-=。

解：??

??--=b a V 111，01det ≠+-=a b V ，只需1-≠b a 。

9-24 已知系统各矩阵为

??????????=100240231A ，??

?????=010010B ，??????=100001C ，

试用传递函数矩阵判断系统的可控性、可观测性。

解：????

???

???-----=??????????------=---400210234

)4)(1(1100240231)I (1

s s s s s s s s A s ，传递函数矩阵为 ??????----=04

42)4)(1(1

)(s s s s s G ； ????

??????=010*********U ，n U ==3rank ；??

???

????

???=100

231

100001

V ，n V ==3rank ；该实现是完全可控且完全可观测的。 9-25 将下列状态方程化为可控标准型

u x x ??

????+??????-=114321&。解： x T x =；1

-=TAT A ，TB B =；

65)I det(2+-=-s s A s ，??????-=Λ0115，??????-=7111U ；??

????-=Λ=-12161U T ，???

???-=621181T ； u x x ??

????+??????-=1051010&。注：若不要求计算变换矩阵，可根据特征多项式直接列写可控标准型。 9-26 已知系统传递函数为

)(2+++=

s s s s G ,

试写出系统可控不可观测、可观测不可控、不可控不可观测的动态方程。

解：系统传递函数的分子和分母多项式中有公因式)1(+s ，任何2维动态方程不可能是既完全可控又完全可观测的。

可控不可观测动态方程 u x x ???

???+??????--=103210&， []x y 11=；可观测不可控动态方程 u x x ??????+??????--=113120&， []x y 10=；不可控不可观测动态方程 u x x ??

????+??????--=011002&， []x y 01=。

9-27 已知系统各矩阵为

?????

???---=4023032600200001

A ，?????

??????=2301b ，[]1134--=c ，试求可控子系统、不可控子系统的态方程。解：1-=TAT A ，TB B =，1

-=CT

C ；

?????????

???=19147112333300001111U ，2rank =U ；选取????????????=-001121033010000111T ，?

???

??????---=09027

0090100210011

91T ； u x x ?????

??????+????????????---=00013200020009/25109/240&，[]x y 13101-=；可控子系统动态方程： u x x c c ??

????+??????-=015140&，[]c c x y 101=；不可控子系统动态方程： c c x x ?

??-=3202&，[]c c x y 13-=。 9-28 系统各矩阵同习题9-27，试求可观测子系统、不可观测子系统的态方程。

解：1

-=TAT A ，TB B =，1-=CT C ；

?????????

???--------=6427118191693413431071134V 初等变换成?????

???????---0000000120129113

4，n V <=3rank ，选取变换矩阵，?????????

???------=10001693413431071134T ，????

????????------=-120002419971081252324017121211T ； u x x

?????

??????+????????????--=246101212/712/25108191201000010&，[]x y 0001=；不可观测子系统动态方程： u x x o o 22+=&。

可观测子系统动态方程：

u x x o o ?????

?????+??????????-=4610181912100010&，[]o o x y 001=； 9-29 设被控系统状态方程为

u x x ????

?????+??????????--=10001010110010&，可否用状态反馈任意配置闭环极点？求状态反馈矩阵，使闭环极点位于}31,10{j ±--，并画出状

态变量图。

解：????

??????=99010010901001000

U ，3rank =U ，系统完全可控，可用状态反馈任意配置闭环极点。

期望的特征多项式为 402412)42)(10()(2

32+++=+++=s s s s s s s K α；

待定参数特征多项式为 1322

3310)91010()910()(k s k k s k s s +-++-+=α；解得， []1.22.14=K 。状态变量图如下：

9-30 设被控系统动态方程为

u x x ??

????+??????=100010&，[]x y 01=，试设计全维状态观测器，使其闭环极点位于}2,{r r --，)0(>r ，并画出状态变量图。

解：期望的观测器特征多项式为 2

223)2)(()(r rs s r s r s s L ++=++=α；

待定系数的特征多项式为 212

)I det()(l s l s LC A s s ++=+-=α；

????=223r r L ；

u y r r z r r z ??

???+?

?????+??????--=1023021322&，z x =?。状态变量图如右图所示。

9-31 设被控系统动态方程为

u x x ??

??????--+??????????-=102102

310101500&，[]x y 100=，试检查被控系统的可控性、可观测性；求输出至输入的反馈矩阵，使闭环极点位于3.122.0j ±-，

57.0-，并画出状态变量图。

解：??

??????----=511012215002

1U ，3rank 1=U ，可控性判别矩阵满秩；动态方程是可观测标准型；

被控系统是完全可控且完全可观测的；期望的特征多项式为

9909.09892.101.1)7384.144.0)(57.0()(232+++=+++=s s s s s s s K α；

选取状态反馈矩阵 ???

???=32

1000k k k K ；则待定参数特征多项式为

)5105()1255()32()(123212233

--+---++-+=k k s k k k s k k s s α

解得 ?

????---=6068.23084.07533.0000

K ；构造全维状态观测器，其极点选为}5,3,2{---；则，

303110)5)(3)(2()(23+++=+++=s s s s s s s L α， )5()1()3()(12233l s l s l s s +-++-+++=α；

即??

??????=73235L ；y u z z ?????

?????+??????????--+??????????---=73235102102101031013000&，z x =?；

9-32 已知系统的传递函数为

)

3)(2)(1()

2)(1()(+-++-=

s s s s s s G ，

能否利用状态反馈将传递函数变成

)

3)(2(1

)(++-=

s s s s G k 。

若有可能，求出一个满足要求的状态反馈阵K ，并画出状态变量图。(提示：状态反馈不改变原传递函数零点。)

解：系统能控规范形是完全可控的，完全可控系统可利用状态反馈任意配置闭环极点。将闭环极点配

置在}3,2,1{---上，即可满足要求。取 ][321k k k K =；原系统和状态反馈系统的状态实现分别为：

[]x y u x x 112100256100010-=??????????+??????????-=&； []x

y r x x 112,10071612100010-=??????????+??????????---=&；

即，1261-=-k ；1652-=-k ；721-=--k ；

解得：[]52118=K 。

9-33 已知系统动态方程各矩阵为

??????????-=020113021A ，??

?????=100b ，[]111-=c ，

试检查可观测性，设计)(q n -维观测器，并使所有极点配置在4-。解：3=n ，1=q ；

???

?????---=171112111V ，n V ==3rank ，该系统完全能观；选取变换矩阵

??????????-=100010111P ， ??????????-==-1000101111Q P ；则 ??????????--=020413312A ；??

?????=101b ；??????=21l l L ；

u b L b y A L A L A L A z A L A z )()]()[()(12112112221222-+-+-+-=&，z Q y L Q Q x

221)(?++=； 1682++=s s L α，)86()13()det()(212121222-++-++=+-=l l s l l s A L A sI s α；解得，

???=??????=7647.17059.317/3017/63L ；)(q n -维观测器方程如下： y u z z ??????+??????--+??????---=5294.34118.47647.07059.32941.52353.01176.77059.2&，z y x ????

?????+??????????=1001117647.17059.34706.4?。

9-34 试用李亚普诺夫第二法判断下列系统平衡态的稳定性：

211x x x +-=&，21232x x x -=&。解：李亚普诺夫方程Q PA P A -=+T

，

其中系统矩阵为 ??????--=3211

A ；取??????=22121211p p p p P ，??????=221100q q Q ，060422

11>=>=q q ；解得 089920>??

???=P ，系统的平衡态是渐近稳定的；

（或采用李亚普诺夫方程I T -=+PA P A ，解得 050.175.175.100.4>?

???=P 。） 9-35 已知系统的状态方程为

u x x ????

?????+??????????---=01200115.0001035.02&，当I =Q 时，=P ？若选Q 为半正定矩阵，=Q ？对应=P ？判断系统稳定性。解：系统稳定性与所选取的矩阵Q 无关。

当I =Q 时，由李亚普诺夫方程Q PA P A -=+T

得到

=+-+--=--=-+--=+-=++-=0

5.025.031260231

205.05.01433231312331323131123221213121111p p p p p p p p p p p p p p p p ，解得 ????

??????-----=535135141013107281P ，由011

可取]001[diag =Q ，（必须011>q ），经检验知},{Q A 完全可观测；解得?????

???

??-----=459

9151115117

281P ，P 非正定，系统不稳定。 9-36 设线性定常离散系统状态方程为

)(002/100010)1(k x a k x ??

?????=+，0>a ，

试求使系统稳定的a 值范围。

解1：离散系统渐近稳定的充要条件是所有特征值均在单位圆内。

由05.0)I det(3

=-=-a z A z ，得2||

解2：选取I P =；由Q P P -=-ΦΦT

，计算得4/1211a q -=，其余0=ij q ，因},{Q Φ完全可观测，

只需2||

线性系统的时域分析法(第七讲)

第三章线性系统的时域分析法 3.1 引言分析控制系统的第一步是建立模型，数学模型一旦建立，第二步分析控制性能，分析有多种方法，主要有时域分析法，频域分析法，根轨迹法等。每种方法，各有千秋。均有他们的适用范围和对象。本章先讨论时域法。实际上，控制系统的输入信号常常是不知的，而是随机的。很难用解析的方法表示。只有在一些特殊的情况下是预先知道的，可以用解析的方法或者曲线表示。例如，切削机床的自动控制的例子。在分析和设计控制系统时，对各种控制系统性能得有评判、比较的依据。这个依据也许可以通过对这些系统加上各种输入信号比较它们对特定的输入信号的响应来建立。许多设计准则就建立在这些信号的基础上，或者建立在系统对初始条件变化（无任何试验信号）的基础上，因为系统对典型试验信号的响应特性，与系统对实际输入信号的响应特性之间，存在着一定的关系；所以采用试验信号来评价系统性能是合理的。 3.1.1 典型试验信号经常采用的试验输入信号： ① 实际系统的输入信号不可知性； ② 典型试验信号的响应与系统的实际响应，存在某种关系； ③ 电压试验信号是时间的简单函数，便于分析。突然受到恒定输入作用或突然的扰动。如果控制系统的输入量是随时间逐步变化的函数，则斜坡时间函数是比较合适的。（单位）阶跃函数（Step function ） 0,)(1≥t t 室温调节系统和水位调节系统（单位）斜坡函数（Ramp function ）速度 0,≥t t ∝ （单位）加速度函数（Acceleration function ）抛物线 0,2 12 ≥t t （单位）脉冲函数（Impulse function ） 0,)(=t t δ 正弦函数（Simusoidal function ）Asinut ，当输入作用具有周期性变化时。通常运用阶跃函数作为典型输入作用信号，这样可在一个统一的基础上对各种控制系统的特性进行比较和研究。本章讨论系统非周期信号（Step 、Ramp 、对正弦试验信号相应，将在第五章频域分析法，第六章校正方法中讨论）作用下系统的响应。 3.1.2 动态过程和稳态过程

线性系统的时域分析方法

第三章线性系统的时域分析方法教学目的：通过本章学习，熟悉控制系统动态性能指标定义，掌握线性系统稳定的充要条件和劳斯判椐的应用，以及稳态误差计算方法，掌握一阶、二阶系统的时域分析方法。教学重点：掌握系统的动态性能指标，能熟练地应用劳斯判椐判断系统稳定性，二阶系统的动态响应特性分析。教学难点：高阶系统的的动态响应特性分析。本章知识结构图：系统结构图闭环传递函数一阶标准式二阶标准式特征方程稳定性、稳定域代数判据误差传递函数误差象函数终值定理稳态误差开环传递函数系统型别、开环增益公式静态误差系数第九讲

3.1 系统时间响应的性能指标一、基本概念 1、时域分析方法：根据系统的数学模型求出系统的时间响应来直接分析和评价系统的方法。（1）响应函数分析方法：建立数学模型→确定输入信号→求出输出响应→ 根据输出响应→系统分析。（2）系统测试分析方法：系统加入扰动信号→测试输出变化曲线→系统分析。系统举例分析：举例：原料气加热炉闭环控制系统 2、分析系统的三大要点（1）动态性能(快、稳) （2）稳态性能（准）（3）稳定性（稳）二、动态性能及稳态性能 1、动态过程（过渡过程）：在典型信号作用下，系统输出从初始状态到最终状态的响应过程。（衰减、发散、等幅振荡） 2、稳态过程：在典型信号作用下，当t → ∞ 系统输出量表现的方式。表征输出量最终复现输入量的程度。（稳态误差描述） 3、动态稳态性能指标图3-1温度控制系统原理图（1）上升时间tr ：从稳态值的10%上升到稳态值的90%所需要的时间。（2）峰值时间tp ：从零时刻到达第一个峰值h(tp)所用的时间。（3）超调量δ%：最大峰值与稳态值的差与稳态值之比的百分数。（稳）（3-1） %100)(()(%?∞∞-= h h t h p ） δ

线性系统状态空间分析报告与运动解

【实验地点】课外(宿舍) 【实验目的】 1、学会利用MATLAB 实现离散系统传递函数模型的生成 2、学会利用MATLAB 将连续系统离散化【实验设备与软件】 1、MATLAB/Simulink 数值分析软件 2、计算机一台【实验原理】 1、求矩阵特征值和特征向量命令格式[V J]=eig （A ） Cv=eig(A) 说明：V 特征向量，J 是Jordan 型，cv 是特征值列向量 2、求运动的方法（1）利用Laplace 逆变换----适合于连续/离散线性系统采用ilaplace/iztrans 对传递函数求逆，这种方法一般是零输入情况下求响应。（2）用连续（离散）状态转移矩阵表示系统解析解----适合于线性定常系统对连续定常系统有：假设初始时刻为零，LTI 系统的解析解为dt Bu e e x e t x t At At At ??+=0 )()0()(τ。若u （t ）是单位阶跃输入，则上述解可写成dtBu e e x e t x t At At At ? ?+=0 )()0()(τ。进一步简化为： Bu A Bu A x e t x At 11))0(()(---+= 对离散线性定常系统有： ∑---+ =1 1 )()0()(k i k k i Hu G x G k x

（3）状态方程的数值分析方法----适合于连续线性系统和非线性系统采用直接数值积分很容易的处理各种定常/时变和线性/非线性系统。有很多数值积分方法，其中有一类预测-修正数值积分方法+自适应步长调整的算法比较有效。在MATLAB/Simulink 中包含的多种有效的、适用于不同类型的ODE 求解算法，典型的是Runge-Ktuta 算法，其通常使用如下的函数格式： [t,x]=ode45(odefun,[ti,tf],x0,options)----采用四阶、五阶Runge-Ktuta 算法 [t,x]=ode23(odefun,[ti,tf],x0,options)----采用二阶、三阶Runge-Ktuta 算法说明：a.这两个函数是求解非刚性常微分方程的函数。 b.参数options 为积分的误差设置，取值为相对误差‘reltol ’和绝对误差‘abstol ’;[ti,tf]求解的时间围；x0是初值是初值向量；[t,x]是解。（4）利用CotrolToolBox 的离散化求解函数----适合于TLI 系统用step （）/impulse()函数求取阶跃输入/冲激输入时系统的状态响应：当系统G 是连续的情况下：调用[y,t,x]=step/impulse(G )会自动对连续系统G 选取采样时间围和周期；调用[y,t,x]=step/impulse(G ，ti:Ts:tf)由用户自己定义对连续系统G 的样时间围和周期；当系统G 是离散的情况下: 调用[y,t,x]=step/impulse(G )会按离散系统G 给出的采样周期计算；调用[y,t,x]=step/impulse(G ，ti:Ts:tf)是Ts 必须与离散系统G 的采样时间围和周期一致。另外lsim()函数调用格式：[y,x,t]=lsim(G,u,ti,TS,tf,x0) 零输入响应调用函数initial （）,格式：[y,x,t]=(G,x0) (5)利用simulink 环境求取响应----适用于所有系统求取响应使用simulink 求取线性或非线性系统的响应，调用格式如下： [t,x,y]=sim(‘XX.mdl ’,ti:Ts:tf,options,u) 【实验容】已知线性系统：]) (201)() (2 10)(404040202119201921)(t x t y t u t x t x +-----? 已知线性系统 1、利用Matlab 求零状态下的阶跃响应（包括状态和输出），生成两幅图：第一幅绘制各状态响应曲线并标注；第二幅绘制输出响应曲线。

信号与系统期末考试试题(有答案的)

信号与系统期末考试试题一、选择题（共10题，每题3分，共30分，每题给出四个答案，其中只有一个正确的） 1、卷积f 1(k+5)*f 2(k-3) 等于。（A ）f 1(k)*f 2(k) （B ）f 1(k)*f 2(k-8)（C ）f 1(k)*f 2(k+8)（D ）f 1(k+3)*f 2(k-3) 2、积分 dt t t ? ∞ ∞ --+)21()2(δ等于。（A ）1.25（B ）2.5（C ）3（D ）5 3、序列f(k)=-u(-k)的z 变换等于。（A ） 1-z z （B ）-1-z z （C ）11-z （D ）1 1--z 4、若y(t)=f(t)*h(t),则f(2t)*h(2t)等于。（A ） )2(41t y （B ）)2(21t y （C ）)4(41t y （D ）)4(2 1 t y 5、已知一个线性时不变系统的阶跃相应g(t)=2e -2t u(t)+)(t δ,当输入f(t)=3e —t u(t)时，系统的零状态响应y f (t)等于（A ）(-9e -t +12e -2t )u(t) （B ）(3-9e -t +12e -2t )u(t) （C ）)(t δ+(-6e -t +8e -2t )u(t) （D ）3)(t δ +(-9e -t +12e -2t )u(t) 6、连续周期信号的频谱具有（A ）连续性、周期性（B ）连续性、收敛性（C ）离散性、周期性（D ）离散性、收敛性 7、周期序列2)455.1(0 +k COS π的周期N 等于（A ） 1（B ）2（C ）3（D ）4 8、序列和 ()∑∞ -∞ =-k k 1δ等于（A ）1 (B) ∞ (C) ()1-k u (D) ()1-k ku 9、单边拉普拉斯变换()s e s s s F 22 12-+= 的愿函数等于 ()()t tu A ()()2-t tu B ()()()t u t C 2- ()()()22--t u t D 10、信号()()23-=-t u te t f t 的单边拉氏变换()s F 等于 ()A ()()()232372+++-s e s s ()() 2 23+-s e B s

自考信号与线性系统分析内部题库含答案

单项选择题。 1. 已知序列3()cos( )5 f k k π=为周期序列，其周期为（） A ． 2 B. 5 C. 10 D. 12 2. 题2图所示 () f t 的数学表示式为（）图题2 A ．()10sin()[()(1)]f t t t t πεε=+- B. ()10sin()[()(1)]f t t t t πεε=-- C. ()10sin()[()(2)] f t t t t πεε=-- D. ()10sin()[()(2)] f t t t t πεε=+- 3.已知sin() ()()t f t t dt t πδ∞ -∞=? ，其值是（） A ．π B. 2π C. 3π D. 4π 4.冲激函数()t δ的拉普拉斯变换为（） A ． 1 B. 2 C. 3 D. 4 5.为了使信号无失真传输，系统的频率响应函数应为（） A ． ()d jwt H jw e = B. ()d jwt H jw e -= C. ()d jwt H jw Ke = D. ()d jwt H jw Ke -= 1 f( t 0 10 正弦函数

6.已知序列1()()()3 k f k k ε=,其z 变换为（） A ． 1 3 z z + B. 1 3 z z - C. 1 4 z z + D. 1 4 z z - 7.离散因果系统的充分必要条件是（ A ） A ．0,0)(<=k k h B. 0,0)(>=k k h C. ,0)(<>k k h 8.已知()f t 的傅里叶变换为()F jw ，则(3)f t +的傅里叶变换为（） A ．()jw F jw e B. 2()j w F jw e C. 3()j w F jw e D. 4()j w F jw e 9.已知)()(k k f k εα=，)2()(-=k k h δ，则()()f k h k *的值为（） A ．) 1(1 --k k εα B. ) 2(2--k k εα C. ) 3(3--k k εα D. ) 4(4--k k εα 10.连续时间系统的零输入响应的“零”是指（ A ） A. 激励为零 B. 系统的初始状态为零 C. 系统的冲激响应为零 D. 系统的阶跃响应为零 11. 已知序列k j e k f 3 )(π=为周期序列，其周期为（） A ． 2 B. 4 C. 6 D. 8 12. 题2 图所示 () f t 的数学表示式为（）

实验25线性系统状态空间分析和运动解

广西大学实验报告纸【实验时间】2014年06月15日【实验地点】(课外) 【实验目的】 1、掌握线性系统状态空间的标准型、解及其模型转换。【实验设备与软件】 1、MATLAB数值分析软件【实验原理】 Matlab提供了非常丰富的线性定常连续系统的状态空间模型求解(即系统运动轨迹的计算)的功能,主要的函数有 ①、阶跃响应函数step()可用于计算在单位阶跃输入和零初始状态(条件)下传递函数模型的输出响应,或状态空间模型的状态和输出响应，其主要调用格式为 step(sys,t) [y,t] = step(sys,t) [y,t,x] = step(sys,t) ②、脉冲激励下的仿真函数impulse()可用于计算在脉冲刺激输入下传递函数模型的输出响应,或状态空间模型的状态和输出响应,其主要调用格式为 impulse(sys,t) [y,t] = impulse(sys,t) [y,t,x] = impulse(sys,t) ③、任意输入激励下的仿真函数lsim()可用于计算在给定的输入信号序列(输入信号函数的采样值)下传递函数模型的输出响应，其主要调用格式为 lsim(sys,u,t,x0) [y,t,x] = lsim(sys,u,t,x0) 【实验内容、方法、过程与分析】已知线性系统 1、利用Matlab求零状态下的阶跃响应（包括状态和输出），生成两幅图：第一幅绘制各状态响应曲线并标注；第二幅绘制输出响应曲线。状态响应曲线： A=[-21 19 -20;19 -21 20;40 -40 -40]; B=[0;1;2]; C=[1 0 2]; D=[0]; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵 X0=[0;0;0]; % 输入初始状态 sys=ss(A,B,C,D); %构造传递函数 [y,x,t]=step(sys); % 绘以时间为横坐标的状态响应曲线图 plot(t,x); grid;

《信号与线性系统》试题与答案5

综合测试(三) 一、选择题（本题共6小题，每小题3分，共18分） 1、若想使连续时间信号在通过线性非时变系统传输时，波形不会产生失真，而仅仅是延时一段时间输出，则要求系统的单位冲激响应必须满足（） A. B. C. D. 2、序列和等于（） A. 1 B. C. D. 3、连续时间信号的单边拉普拉斯变换为（） A. B. C. D. 4、下列各式中正确的是（） A． B. C.D． 5、单边Z变换对应的原时间序列为（） A．B． C．D． 6．请指出是下面哪一种运算的结果？（）

A ．左移6 B. 右移6 C ．左移2 D. 右移2 三、描述某系统的微分方程为 y ”(t) + 4y ’(t) + 3y(t) = f(t) 求当f(t) = 2e -2t ，t ≥0；y(0)=2，y ’(0)= -1时的解；（ 15分）解: (1) 特征方程为λ2 + 4λ+ 3 = 0 其特征根λ1= –1，λ2= –2。齐次解为 y h (t) = C 1e -t + C 2e -3t 当f(t) = 2e –2 t 时，其特解可设为 y p (t) = Pe -2t 将其代入微分方程得 P*4*e -2t + 4(–2 Pe -2t ) + 3Pe -t = 2e -2t 解得 P=2 于是特解为 y p (t) =2e -t 全解为： y(t) = y h (t) + y p (t) = C 1e -t + C 2e -3t + 2e -2t 其中待定常数C 1,C 2由初始条件确定。 y(0) = C 1+C 2+ 2 = 2， y ’(0) = –2C 1 –3C 2 –1= –1 解得 C 1 = 1.5 ，C 2 = –1.5 最后得全解 y(t) = 1.5e – t – 1.5e – 3t +2 e –2 t , t ≥0 三、描述某系统的微分方程为 y ”(t) + 5y ’(t) + 6y(t) = f(t) 求当f(t) = 2e -t ，t ≥0；y(0)=2，y ’(0)= -1时的解；（ 15分）解: (1) 特征方程为λ2 + 5λ+ 6 = 0 其特征根λ1= –2，λ2= –3。齐次解为 y h (t) = C 1e -2t + C 2e -3t 当f(t) = 2e – t 时，其特解可设为 y p (t) = Pe -t 将其代入微分方程得 Pe -t + 5(– Pe -t ) + 6Pe -t = 2e -t 解得 P=1 于是特解为 y p (t) = e -t 全解为： y(t) = y h (t) + y p (t) = C 1e -2t + C 2e -3t + e -t 其中待定常数C 1,C 2由初始条件确定。 y(0) = C 1+C 2+ 1 = 2， y ’(0) = –2C 1 –3C 2 –1= –1 解得 C 1 = 3 ，C 2 = – 2 最后得全解 y(t) = 3e – 2t – 2e – 3t + e – t , t ≥0 四、如图信号f(t)的拉氏变换F(s) = ，试观 )e e 1(e 2s s s s s -----)e e 1(e 2 s s s s s -----

信号与系统期末考试试题

重庆大学信号与线性系统期末考试试题一、填空题：（30分，每小题3分） 1. =-? ∞ ∞ -dt t t )()5cos 2(δ 。 2. ()dt t e t 12-?+∞ ∞ --δ= 。 3. 已知 f (t )的傅里叶变换为F (j ω), 则f (2t -3)的傅里叶变换为。 4. 已知 6 51 )(2 +++= s s s s F ，则=+)0(f ; =∞)(f 。 5. 已知 ω ωπδεj t FT 1 )()]([+=，则=)]([t t FT ε 。 6. 已知周期信号 )4sin()2cos()(t t t f +=，其基波频率为 rad/s ；周期为 s 。 7. 已知 )5(2)2(3)(-+-=n n k f δδ，其Z 变换 =)(Z F ；收敛域为。 8. 已知连续系统函数1342 3)(2 3+--+= s s s s s H ，试判断系统的稳定性：。 9．已知离散系统函数1.07.02 )(2+-+=z z z z H ，试判断系统的稳定性：。 10．如图所示是离散系统的Z 域框图，该系统的系统函数H(z)= 。二．(15分)如下方程和非零起始条件表示的连续时间因果LTI 系统，

?????==+=++-- 5 )0(',2)0() (52)(4522y y t f dt df t y dt dy dt y d 已知输入 )()(2t e t f t ε-=时，试用拉普拉斯变换的方法求系统的零状态响应 )(t y zs 和零输入响应)(t y zi ，0≥t 以及系统的全响应),(t y 0≥t 。三．（14分） ① 已知2 36 62)(22++++=s s s s s F ,2]Re[->s ,试求其拉氏逆变换f (t )； ② 已知) 2(2 35)(2>+-=z z z z z X ，试求其逆Z 变换)(n x 。四（10分）计算下列卷积： 1. }1,0,6,4,3{}4,1,2,1{)()(21--*=*k f k f ； 2． )(3)(23t e t e t t εε--* 。

线性系统的状态空间分析与综合

第九章线性系统的状态空间分析与综合一、教学目的与要求：通过本章内容的学习，使学生建立起状态变量和状态空间的概念，掌握线性定常系统状态空间模型的建立方法，状态空间表达式的线性变换，状态完全能控或状态完全能观测的定义，及其多种判据方法，状态转移矩阵的求法，传递函数矩阵与状态空间表达式的关系。二、授课主要内容： 1．线性系统的状态空间描述 2．线性系统的可控性与可观测性 3．线性定常系统的状态反馈与状态观测器（详细内容见讲稿）三、重点、难点及对学生的要求（掌握、熟悉、了解、自学） 1.重点掌握线性定常系统状态空间模型的建立方法与其他数学描述（微分方程、传递函数矩阵）之间的关系。 2.掌握采用状态空间表述的系统运动分析方法，状态转移矩阵的概念和求解。 3.掌握系统基本性质——能控性和能观测性的定义、有关判据及两种性质之间的对偶性。 4.理解状态空间表达式在线性变换下的性质，对于完全能控或能观测系统，构造能控、能观测标准形的线性变换方法，对于不完全能控或不完全能观测系统，基于能控性或能观测性的结构分解方法。 5.掌握单变量系统的状态反馈极点配置和全维状态观测器设计方法，理解分离定理，带状态观测器的状态反馈控制系统的设计。重点掌握线性系统的状态空间描述和求解，线性系统的可控性与可观测性及状态反馈与状态观测器。四、主要外语词汇线性系统 linear system 状态空间 state space 状态方程 state equation

状态向量 state vector 传递函数矩阵 translation function matrix 状态转换矩阵 state-transition matrix 可观测标准形 observational standard model 可控标准形 manipulative standard model 李亚普诺夫方程Lyaponov equation 状态观测器 state observation machine 对偶原理 principle of duality 五、辅助教学情况（见课件）六、复习思考题 1．什么是系统的状态空间模型？状态空间模型中的状态变量、输入变量、输出变量各指什么？ 2．通过机理分析法建立系统状态空间模型的主要步骤有哪些？ 3．何为多变量系统？如何用传递矩阵来描述多变量系统的动态特性？在多变量系统中，环节串联、并联、反馈连接时，如何求取总的传递矩阵？4．试简述数学模型各种表达式之间的对应关系。 5．用非奇异矩阵P对状态方程式进行线性状态变换后，与原状态方程式之间存在什么关系？ 6．试简述系统能控性与能观性两个概念的含义及意义。 7．试述能控性和能观性定义。 8．试述系统能控性和能观性常用判据。 9．何谓对偶系统和对偶原理？ 10．什么是状态方程的线性变换？ 11．试述系统状态方程规范型变换的条件、特点及变换的基本方法。 12．试述状态能控性与能观性和系统传递函数（阵）的关系。七、参考教材（资料） 1．《自动控制原理与系统》上、下册清华大学吴麒等国防工业出版社

《信号与线性系统》期末试卷

2006-2007学年第二学期《信号与线性系统》（课内）试卷A 卷一、计算题（共45分） 1.（5分）计算积分dt t t t )6 ()sin (π δ- +?+∞ ∞-的值。 2.（5分）绘出函数)1()]1()([-+--t u t u t u t 的波形图。 3.（6分）已知)2()()(),1()()(21--=--=t u t u t f t u t u t f ，求卷积)()(21t f t f *。 4.（6分）若)(t f 的傅里叶变换已知，记为)(ωF ，求)1()1(t f t --对应的傅里叶变换。

5.（6分）如下图所示信号，已知其傅里叶变换，记为)(ωF ，求：（1）)0(F ；（2）?+∞ ∞ -ωωd F )(。 6.（5分）已知)(t f 对应的拉氏变换为)(s F ，求)/(/a t f e a t -（0>a ）对应的拉氏变换。 7.（6分）已知)(t f 对应的拉氏变换2 3)(2 +-=-s s e s F s ，求)(t f

8．（6分）线性时不变系统的单位样值响应为)(n h ，输入为)(n x ，且有 )4()()()(--==n u n u n x n h ，求输出)(n y ，并绘图示出)(n y 。二、综合题（共计55分） 1、（10分）系统如图所示，已知t t x 2000cos )(=，t t t f 2000cos 100cos )(=，理想低通滤波器)300()300()(--+=ωωωu u H ，求滤波器的响应信号)(t y 。 x(t) y(t) f(t)

信号与线性系统分析习题答案

1 / 257 信号与线性系统课后答案第一章信号与系统（一） 1-1画出下列各信号的波形【式中)() (t t t r ε=】为斜升函数。（2）∞<<-∞=- t e t f t ,)( （3）)()sin()(t t t f επ= （4）)(sin )(t t f ε= （5）)(sin )(t r t f = （7）)(2)(k t f k ε= （10）)(])1(1[)(k k f k ε-+= 解：各信号波形为（2）∞<<-∞=-t e t f t ,)( （3）) ()sin()(t t t f επ=

2 / 257 （4）)(sin )(t t f ε= （5）) (sin )(t r t f =

3 / 257 （7）)(2)(k t f k ε= （10）) (])1(1[)(k k f k ε-+=

4 / 257 1-2 画出下列各信号的波形[式中)()(t t t r ε=为斜升函数]。（1）)2()1(3)1(2)(-+--+=t t t t f εεε （2）)2()1(2)()(-+--=t r t r t r t f （5） )2()2()(t t r t f -=ε （8）)]5()([)(--=k k k k f εε （11）)]7()()[6 sin()(--=k k k k f εεπ （12） )]()3([2)(k k k f k ---=εε 解：各信号波形为（1） ) 2()1(3)1(2)(-+--+=t t t t f εεε

5 / 257 （2） )2()1(2)()(-+--=t r t r t r t f （5） ) 2()2()(t t r t f -=ε

自动控制原理》实验2(线性系统时域响应分析

实验二线性系统时域响应分析一、实验目的 1．熟练掌握step( )函数和impulse( )函数的使用方法，研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2．通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。二、基础知识及MATLAB 函数（一）基础知识时域分析法直接在时间域中对系统进行分析，可以提供系统时间响应的全部信息，具有直观、准确的特点。为了研究控制系统的时域特性，经常采用瞬态响应（如阶跃响应、脉冲响应和斜坡响应）。本次实验从分析系统的性能指标出发，给出了在MATLAB 环境下获取系统时域响应和分析系统的动态性能和稳态性能的方法。用MATLAB 求系统的瞬态响应时，将传递函数的分子、分母多项式的系数分别以s 的降幂排列写为两个数组num 、den 。由于控制系统分子的阶次m 一般小于其分母的阶次n ，所以num 中的数组元素与分子多项式系数之间自右向左逐次对齐，不足部分用零补齐，缺项系数也用零补上。 1．用MATLAB 求控制系统的瞬态响应 1）阶跃响应求系统阶跃响应的指令有： step(num,den) 时间向量t 的范围由软件自动设定，阶跃响应曲线随即绘出 step(num,den,t) 时间向量t 的范围可以由人工给定（例如 t=0:0.1:10） [y ，x]=step(num,den) 返回变量y 为输出向量，x 为状态向量在MATLAB 程序中，先定义num,den 数组，并调用上述指令，即可生成单位阶跃输入信号下的阶跃响应曲线图。考虑下列系统： 25 425 )()(2++=s s s R s C 该系统可以表示为两个数组，每一个数组由相应的多项式系数组成，并且以s

第九章线性系统的状态空间分析与综合习题

第九章线性系统的状态空间分析与综合 9-1 已知电枢控制的直流司服电机的微分方程组及传递函数为 b a a a a a E t d di L i R u ++=，t d d K E m b b θ=，a m m i C M =，t d d f t d d J M m m m m m θθ+=2 2； )] ()([)()(2 m b m a a m m a m a m a m C K f R s R J f L s J L s C s U s ++++=Θ。 ⑴ 设状态变量m x θ=1，m x θ&=2，m x θ&&=3，输出量m y θ=，试建立其动态方程； ⑵ 设状态变量a i x =1，m x θ=2，m x θ&=3，输出量m y θ=，试建立其动态方程； ⑶ 设x T x =，确定两组状态变量间的变换矩阵T 。解：⑴ 由传递函数得 a m a m m a m b m a m a u C x R J f L x C K f R x J L ++-+-=323)()(&，动态方程为 []x y u x x x x x x 001100010001032121321=??????????+????????????????????--=??????????αα&&&，其中)/()()/()()/(21m a a m m a m a m b m a m a a m J L R J f L J L C K f R J L u C u +=+==αα； ⑵ 由微分方程得 3 133 2311x f x C x J x x u x K x R x L m m m a b a a -==---=&&&，即 []x y u x x x a a a a x x x a 0200010100032133311311321=???? ? ?????+?????????????????? ??=??????????&&&，其中 m m m m a b a a J f a J C a L K a L R a ////33311311-==-=-=； ⑶ 由两组状态变量的定义，直接得到???? ? ???????????????=??????????3213331 321010001 0x x x a a x x x 。 9-2 设系统的微分方程为 u x x x =++23&&& 其中u 为输入量，x 为输出量。 ⑴ 设状态变量x x =1，x x &=2，试列写动态方程； ⑵ 设状态变换211x x x +=，2122x x x --=，试确定变换矩阵T 及变换后的动态方程。解：⑴ u x x x x ??????+????????????--=??????1032102121&&，[]?? ????=2101x x y ； ⑵ ??????=??????2121x x T x x ，??????--=2111T ；?? ????--=-11121 T ；AT T A 1-=，B T B 1-=，CT C =；得，? ?????--=2111T ；u x x x x ?? ????-+????????????-=??????1110012121&&，[]??????=2111x x y 。 9-3 设系统的微分方程为 u y y y y 66116=+++&&&&&& 其中u 、y 分别系统为输入、输出量。试列写可控标准型(即A 为友矩阵)及可观标准型(即A 为友矩阵转置)状态空间表达式，并画出状态变量图。

信号与系统期末试题与答案

课程名称信号与线性系统A 考试学期 08-07 得分适用专业微电、物理、考试形式闭卷考试时间 120分钟姓名班级学号一、选择题（每小题可能有一个或几个正确答案，将正确的题号填入[ ]内） 1．f （5-2t ）是如下运算的结果————————（ C ）（A ）f （-2t ）右移5 （B ）f （-2t ）左移5 （C ）f （-2t ）右移 2 5 （D ）f （-2t ）左移25 2．已知)()(),()(21t u e t f t u t f at -==，可以求得=)(*)(21t f t f —————（ C ）（A ）1-at e - （B ）at e - （C ）)1(1at e a -- （D ）at e a -1 3．线性系统响应满足以下规律————————————（AD ）（A ）若起始状态为零，则零输入响应为零。（B ）若起始状态为零，则零状态响应为零。（C ）若系统的零状态响应为零，则强迫响应也为零。（D ）若激励信号为零，零输入响应就是自由响应。 4．若对f （t ）进行理想取样，其奈奎斯特取样频率为f s ，则对)23 1 (-t f 进行取样，其奈奎斯特取样频率为————————（B ）（A ）3f s （B ） s f 31 （C ）3（f s -2）（D ）)2(3 1 -s f 5．理想不失真传输系统的传输函数H （jω）是 ————————（B ）（A ）0j t Ke ω- （B ）0 t j Ke ω- （C ）0 t j Ke ω-[]()()c c u u ωωωω+-- （D ）00 j t Ke ω- （00,,,c t k ωω为常数） 6．已知Z 变换Z 1 311 )]([--= z n x ，收敛域3z >，则逆变换x （n ）为——（ A ）（A ）)(3n u n （C ）3(1)n u n - （B ）)(3n u n -- （D ）)1(3----n u n

信号与线性系统分析习题答案-(吴大正-第四版--高等教育出版社)

第一章信号与系统（二） 1-1画出下列各信号的波形【式中)()(t t t r ε=】为斜升函数。（2）∞<<-∞=-t e t f t ,)( （3）)()sin()(t t t f επ= （4）)(sin )(t t f ε= （5）)(sin )(t r t f = （7）)(2)(k t f k ε= （10）)(])1(1[)(k k f k ε-+= 解：各信号波形为（2）∞<<-∞=-t e t f t ,)(

（3）) ()sin()(t t t f επ= （ 4）)(sin )(t t f ε=

（5）) t f= r ) (sin (t （7）) f kε = t ) ( 2 (k

（10）)(])1(1[)(k k f k ε-+= 1-2 画出下列各信号的波形[式中)()(t t t r ε=为斜升函数]。（1）)2()1(3)1(2)(-+--+=t t t t f εεε （2） )2()1(2)()(-+--=t r t r t r t f （5）)2()2()(t t r t f -=ε （8）)]5()([)(--=k k k k f εε

（11）)]7()()[6 sin()(--=k k k k f εεπ （12）)]()3([2)(k k k f k ---=εε 解：各信号波形为（1）)2()1(3)1(2)(-+--+=t t t t f εεε （2）)2()1(2)()(-+--=t r t r t r t f

（5）)2()2()(t t r t f -=ε （8）)]5()([)(--=k k k k f εε

线性系统的时域分析习题答案

第3章线性系统的时域分析学习要点 1控制系统时域响应的基本概念，典型输入信号及意义； 2控制系统稳定性的概念、代数稳定判据及应用； 3控制系统的时域指标，一阶二阶系统的阶跃响应特性与时域指标计算； 4高阶系统时域分析中主导极点和主导极点法； 5 控制系统稳态误差概念、计算方法与误差系数，减小稳态误差的方法。思考与习题祥解题思考与总结下述问题。（1）画出二阶系统特征根在复平面上分布的几种情况，归纳ξ值对二阶系统特征根的影响规律。（2）总结ξ和n ω对二阶系统阶跃响应特性的影响规律。（3）总结增加一个零点对二阶系统阶跃响应特性的影响规律。（4）分析增加一个极点可能对二阶系统阶跃响应特性有何影响（5）系统误差与哪些因素有关试归纳减小或消除系统稳态误差的措施与方法。（6）为减小或消除系统扰动误差，可采取在系统开环传递函数中增加积分环节的措施。请问，该积分环节应在系统结构图中如何配置，抗扰效果是否与扰动点相关答：（1）二阶系统特征根在复平面上分布情况如图所示。图二阶系统特征根在复平面上的分布当0ξ=，二阶系统特征根是一对共轭纯虚根，如图中情况①。当01ξ<<，二阶系统特征根是一对具有负实部的共轭复数根，变化轨迹是以n ω为半径的圆弧，如图中情况②。当1ξ=，二阶系统特征根是一对相同的负实根，如图中情况③。当1ξ>，二阶系统特征根是一对不等的负实根，如图中情况④。

（2）ξ和n ω是二阶系统的两个特征参量。 ξ是系统阻尼比，描述了系统的平稳性。当0ξ=，二阶系统特征根是一对共轭纯虚根，二阶系统阶跃响应为等幅振荡特性，系统临界稳定。当01ξ<<，二阶系统特征根是一对具有负实部的共轭复数根，二阶系统阶跃响应为衰减振荡特性，系统稳定。ξ越小，二阶系统振荡性越强，平稳性越差； ξ越大，二阶系统振荡性越弱，平稳性越好。因此，二阶系统的时域性能指标超调量由ξ值唯一确定，即001_ 100%2 ?=-π ξξ σe 。在工程设计中，对于恒值控制系统，一般取 ξ＝～；对于随动控制系统ξ＝～。 n ω是系统无阻尼自然振荡频率，反映系统的快速性。当ξ一定，二阶系统的时域性能指标调节时间与n ω值成反比，即34 s n t ξω≈:。（3）二阶系统增加一个零点后，增加了系统的振荡性，将使系统阶跃响应的超调量增大，上升时间和峰值时间减小。所增加的零点越靠近虚轴，则上述影响就越大；反之，若零点距离虚轴越远，则其影响越小。（4）二阶系统增加一个极点后，减弱了系统的振荡性，将使系统阶跃响应的超调量减小，上升时间和峰值时间减小；所增加的极点越靠近虚轴，则上述影响就越大；反之，若极点距离虚轴越远，则其影响越小。（5）系统误差与系统的误差度（开环传递函数所含纯积分环节的个数或系统型别）、开环放大系数，以及作用于系统的外部输入信号有关。如果是扰动误差还与扰动作用点有关。因此，减小或消除系统稳态误差的措施与方法有：增大开环放大系数，增加系统开环传递函数中的积分环节，引入按给定或按扰动补偿的复合控制结构。无论采用何种措施与方法减小或消除系统稳态误差，都要注意系统须满足稳定的条件。（6）采取在系统开环传递函数中增加积分环节的措施来减小或消除系统扰动误差时，所增加的积分环节须加在扰动作用点之前。若所增加的积分环节加在扰动作用点之后，则该积分环节无改善抗扰效果作用。这一点可以通过误差表达式分析得到。题系统特征方程如下，试判断其稳定性。（a ）0203.002.023=+++s s s ；（b ）014844122345=+++++s s s s s ；（c ）025266.225.11.0234=++++s s s s 解：（a ）稳定；（b ）稳定；（c ）不稳定。题系统结构如题图所示。控制器)1 1()(s T K s G i p c + =，为使该系统稳定，控制器参数p K 、i T 应满足什么关系

《信号与线性系统》期末试卷要点

2012-2013学年第二学期《信号与线性系统》（课内）试卷A 卷一、计算题（共45分） 1.（5分）计算积分dt t t t )6 ()sin (π δ- +? +∞ ∞ -的值。 2.（5分）绘出函数)1()]1()([-+--t u t u t u t 的波形图。 3.（6分）已知)2()()(),1()()(21--=--=t u t u t f t u t u t f ，求卷积)()(21t f t f *。 4.（6分）若)(t f 的傅里叶变换已知，记为)(ωF ，求)1()1(t f t --对应的傅里叶变换。

5.（6分）如下图所示信号，已知其傅里叶变换，记为)(ωF ，求：（1）)0(F ；（2）? +∞ ∞ -ωωd F )(。 6.（5分）已知)(t f 对应的拉氏变换为)(s F ，求)/(/a t f e a t -（0>a ）对应的拉氏变换。 7.（6分）已知)(t f 对应的拉氏变换2 3)(2+-=-s s e s F s ，求)(t f

8．（6分）线性时不变系统的单位样值响应为)(n h ，输入为)(n x ，且有 )4()()()(--==n u n u n x n h ，求输出)(n y ，并绘图示出)(n y 。二、综合题（共计55分） 1、（10分）系统如图所示，已知t t x 2000 cos )(=，t t t f 2000cos 100cos )(=，理想低通滤波器)300()300()(--+=ωωωu u H ，求滤波器的响应信号)(t y 。 y(t) f(t)

线性系统的状态空间描述

第一章线性系统的状态空间描述 1.内容系统的状态空间描述化输入—输出描述为状态空间描述由状态空间描述导出传递函数矩阵线性系统的坐标转换组合系统的状态空间方程与传递函数矩阵 2.基本概念系统的状态和状态变量状态：完全描述系统时域行为的一个最小变量组状态变量：构成系统状态的变量状态向量设系统状态变量为X i(t),X2(t)厂，X n(t)写成向量形式称为状态向量，记为 _X i (t) x(t)= _X n(t) 状态空间状态空间：以状态变量为坐标轴构成的n维空间状态轨迹：状态变量随时间推移而变化，在状态空间中形成的一条

轨迹。

3. 状态空间表达式设系统r 个输入变量：U i (t ),u 2(t )^ ,u r (t ) m 个输出：yQM), ,y m (t) n 个状态变量：X i (t),X 2(t), ,X n (t) 例：图示RLC 电路，建立状态空间描述 i L C 电容C 和电感L 两个独立储能元件，有两个状态变量, 方程为如图中所注, L di L (t) dt Ri L (t) U c (t) =u(t) C 沁 “L (t) dt X i (t)二 L(t), X 2(t)二 U c (t) 二 LX i (t) RX i (t) X 2(t)二 u(t) Cx (t)二 X (t) N(t) - R/L 殳⑴门1/C 0 匚X 2(— O u(t) U c

输出方程一般定义状态方程：状态变量与输入变量之间的关系 dX i (t) dt = X i (t)二 f i 〔X i (t),X 2(t), ,X n (t)；U i (t),U 2(t), ,U r (t);tl dX 2(t) dt = X 2(t)二 f 2'X i (t),X 2(t)^ ,X n (t)；U i (t),U 2(t), ,U r (t)；t 】 dX n (t) dt 二 X n (t)二 f n 〔X i (t),X 2(t), ^⑴小⑴心⑴,,U 「(t)；t 】用向量表示，得到一阶的向量微分方程 x(t)二 f 'X(t),u(t), t 1 其中 X i (t) U ](t) fQ) “、 X 2(t) - U 2(t) . f 2(?)?Qn X(t) - c R ,u(t)戶；c R , f (?) ^^ ： c R N(t) 一 JU r (t) 一 -f n (叽输出方程：系统输出变量与状态变量、输入变量之间的关系，即 %(t)二 g i X i (t),X 2(t), ,X n (t)；U i (t),U 2(t), ,U r (t)；t ] y 2(t)二 g 2 X i (t), X 2(t), ,X n (t)；U i (t),U 2(t), ,U r (t)；t 〔 y(t)二％(t)二 1 01 X i (t) 殳(t).

第九章 线性系统的状态空间分析与综合习题

第九章线性系统的状态空间分析与综合习题