习题四
习题-4-参考答案'

高分子化学与物理习题习题四 (共聚合,开环聚合部分)1. 写出两种不同形式的二元聚合微分方程,简述其应用条件。
答: 浓度表示的二元共聚微分方程式:])[][]([])[][]([][][1222211121M M r M M M r M M d M d ++= 摩尔分数表示式:222212112121112f r f f f r f f f r F +++= 应用条件:稳态,等活性,分子量大(引发和终止对共聚物的组成影响不大),低转化率(一般小于10%)2. 为什么不研究r 1>1, r 2>1的二元共聚体系?答:根据r 的定义,r>1不利于共聚;r 1>1,r 2>1的体系进行共聚合,一般情况下得到的大多数为两种单体的均聚物的混合物,也可能是含有少量嵌段聚合物(视二者的投料比和均聚合的速度而定)3. 等摩尔M 1, M 2共聚合,若r 1<1, r 2>1,分析聚合后期产物组成情况。
答:r 1<1, r 2>1时,F 1<f 1,单体M 2比单体M 1更易于进入大分子链,M 2消耗很快。
若二者等摩尔混合,聚合前期会形成M 2单体含量多的共聚物,而反应后期M 2已经消耗完,剩余单体M 1发生均聚得到均聚物。
4. 下列单体:St, MMA, V Ac (醋酸乙烯酯), VCl (氯乙烯),B (丁二烯)间,选出自由基共聚合(二元共聚)好的单体组,简要说明理由。
答:按共轭单体/共轭单体,非共轭单体/非共轭单体配对易于形成共聚合的原则来考虑,下列单体对易于共聚合:St/MMA, St/B, MMA/B, V Ac/VC 。
5. F 1-f 1图是如何做出来的?何谓恒比共聚点,此时共聚合的特点是什么? 答:根据二元共聚摩尔微分方程,已知r 1和r 2,将不同的f 1值代入方程则可以得到系列F 1值,画出F 1~ f 1图。
r 1和r 2均小于1时,F 1~ f 1图与对角线有交点,此点为恒比共聚点;此时F 1= f 1计算方法为:()212121r r r f C ---= 在恒比共聚点进行投料聚合,可以得到含量明确的共聚物。
李凡长版 组合数学课后习题答案习题4

第四章 生成函数1. 求下列数列的生成函数: (1){0,1,16,81,…,n 4,…} 解:G{k 4}=235(11111)1x x x x x +++-()(2)343,,,333n +⎧⎫⎛⎫⎛⎫⎛⎫⎨⎬ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎩⎭L 解:3n G n +⎧⎫⎛⎫⎨⎬ ⎪⎝⎭⎩⎭=41(1)x - (3){1,0,2,0,3,0,4,0,……} 解:A(x)=1+2x 2+3x 4+4x 6+…=(211x-)2. (4){1,k ,k 2,k 3,…} 解:A(x)=1+kx+k 2x 2+k 3x 3+…=11kx -. 2. 求下列和式: (1)14+24+…+n 4解:由上面第一题可知,{n 4}生成函数为A(x)=235(11111)1x x x x x +++-()=0kk k a x ∞=∑, 此处a k =k 4.令b n =14+24+…+n 4,则b n =0nk k a =∑,由性质3即得数列{b n }的生成函数为 B(x)= 0nn n b x ∞=∑=()1A x x -=34125(1111)ii i x x x x x i ∞=++++⎛⎫ ⎪⎝⎭∑. 比较等式两边x n 的系数,便得14+24+…+n 4=b n =1525354511111234n n n n n n n n -+-+-+-++++----⎛⎫⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭321(1)(691)30n n n n n =+++-(2)1·2+2·3+…+n (n +1)解:{ n (n +1)}的生成函数为A(x)=32(1)x x -=0kk k a x ∞=∑,此处a k = n (n +1).令b n =1·2+2·3+…+n (n +1),则b n =0nk k a =∑.由性质3即得数列{b n }的生成函数为B(x)=nn n b x ∞=∑=()1A x x-=42(1)xx -=032nk kk x x k =+⎛⎫⎪⎝⎭∑. 比较等式两边x n 的系数,便得1·2+2·3+…+n (n +1)= b n =2(1)(2)213n n n n n +++=-⎛⎫ ⎪⎝⎭. 3. 利用生成函数求解下列递推关系: (1)()7(1)12(2)(0)2,(1)7f n f n f n f f =---==⎧⎨⎩;解:令A(x)=0()n n f n x ∞=∑则有A(x)-f(0)-f(1)x=2()nn f n x ∞=∑=2(7(1)12(2))n nf n f n x∞=---∑=217()12()nnn n x f n x xf n x∞∞==-∑∑=7x(A(x)-f(0))-12x 2A(x).将f(0)=2,f(1)=7代入上式并整理,得22711()(34)17121314n n n x A x x x x x ∞=-==+=+-+--∑. (2)()3(1)53(0)0nf n f n f =-+⋅=⎧⎨⎩;解:令A(x)=0()nn f n x ∞=∑,则有A(x)-f(0)= 1(3(1)53)n nnf n x ∞=-+⋅∑=03()153nn n n n x f n x x x ∞∞==+∑∑=3xA(x)+15x·113x-.A(x)= 215(13)xx - (3)()2(1)(2)(0)0,(1)1f n f n f n f f =-+-==⎧⎨⎩;解:令A(x)=0()nn f n x ∞=∑,则有A(x)-f(0)-f(1)x=2(2(1)(2))n nf n f n x ∞=-+-∑=212()()nnn n x f n x xf n x∞∞==+∑∑=2x(A(x)-f(0))+x 2A(x).将f(0)=0,f(1)=1代入上式并整理,得2()12x A x x x =--.4. 设序列{n a }的生成函数为:343(1)(1)xx x x --+-,但00b a =,110b a a =-,……,1n n n b a a -=-,……,求序列{n b }的生成函数.解:由00b a =,110b a a =-,……,1n n n b a a -=-,得0nk n k b a ==∑,所以A(x)=()1B x x-.由此得B(x)=(1-x)A(x)= 3431xx x -+-,亦即序列{n b }的生成函数。
习题4 马斯京根法流量

流域图
1970年9月一次洪水流量过程
工作步骤
写出水量平衡方程和槽蓄方程 取有限差形式, 推导马斯京根法的流量演算公式
O2=C0I2+C1I1+C2O1 根据报汛要求Δt=3h, 将河段分为3段 假定ke,xe(ke,xe表示每个单元河段的参数, 假定每
个单元河段的ke,xe都是相等的);计算C0,C1,C2; 编写程序进行分段马斯京根流量演算 将计算的出流过程与实测的出流过程进行比较, 如果 拟合得较好, 则取这组ke,xe作为最后的参数;如果 拟合的不好, 则重新假定, 回到第4步重新计算。
习题四 马斯京根法流量演算
任务和要求
根据实测的水文资料,采用试错法分析确定 一次洪水的马斯京根法流量演算参数。通过 本次练习要求弄清:
马斯京根法流量演算公式的推导 分段马斯京根法流量演算的计算过程 参数ke和xe的物理意义 编写程序采用试错法确定一次洪水的ke
和xe
资料情况
沅水下游沅陵至王家河河段, 河长112 公里, 沅陵以上流域面积为76400平方 公里, 王家河以上流域面积为80500平 方公里, 区间面积为4100平方公里, 约 占总面积的5%, 流域内雨量站网的分 布见图1。
上交成果
马斯京根法流量演算公式的,xe 计算和实测的流量过程图表 谈谈对马斯京根法的认识(要求从基本
原理, 基本假定, 参数物理意义, 可能造成 的误差等方面来阐述)
习题4习题解答.docx

习题四解答4.1 设tz = (10,-8,6,-5)「,” = (2,2,0,-5),,求[3« -^,2«],||a||,||« -^||解答:[(z,a] = 102 +(-8)2 +62 +52 =225[a,)3]= 10 • 2 + (―8)-2 + 6-0 + (-5) • (-5) = 293a —” = (28,—26,18,—10)「,因此[3a — /3,2a] = 2[3a—”, a] = 2[28-10 + (-26) . (-8) + 18-6 +(-10)- (-5)] = 1292或者[3a — /32a] = 6[a, a]—2[”, a] = 6-225-2-29 = 1292M = j0,a]=喜= 15, |板『=/,切=33 因此||a —”『=0 —月& —切=00 + /,切一20,切=200 n ||a 一列=10扼4.2 设有三点A(—l,2,l),3(0,3,l),C(0,2,2),求ZBAC解答:A3 = (1,1,0),AC = (1,0,1),由向量夹角余弦公式知cos A =AB AC4.3设a,”的夹角是120°,模分别为12和6,求a,”的距离.解答:||«-=[a-/3,a-/3] = |述 +—2|a||”|cos9 = 252所以阪-”|| = 6j?4.4作为平面几何定理(平行四边形的对角线与边长的关系)的推广,证明|g|2+|H「=2(| 述+四)||or + 印=[a + /3,a + /3] = [a,a] + 2[a0 + [”,切证明:,|旧_”|「=[/-", a _"] = [a, a]-2[a,"] + [”,”] 所以诉 + "If + |a - ”『=2([«,«] + [”,”]) = 2(|述 +1”/).4.5设%=(1,2,—1产,%=(—1,3,顶,的 =(4,T,0)「,试用施密特正交化方法把这组向量正交规范化.解答:岗=% = (1,2,—1)「解答:因为角,代已经正交,所以只要让%,%正交,正交就可以了. 设a 3 = (xj,x 2,x 3)r ,由于[%,%] = O,]%,%] = 0,所以x, + x 2 + x 3 = 0<M — 2X 2 + 想=°显然,(1,0,-1)「是方程组的一个非零解,所以取«3 = (1,0,-l)r .注:由于只要求出一组数满足 叫+心+心-0就可以了,所以不必将这个齐次线性方程%] -2x, + %, = 0组的通解解出来,只要找到一组非零解即可•可以令x 2 = 0,于是就能得到一组解.4.7设A,3都是〃阶正交矩阵,证明AB 也是正交矩阵.证明:山条件知A TA = I,B TB = I,因此(AB)r (AB) = B T (A T A)B = B T B = I ,即 A3 是正交矩阵.4.8 设%= (5,3,1,1)',%=(—4,18,—2,4)',欲使向量/3 = a 2 + ka x 与向正交,求S 解答:/3 = a 2+ka { = (5k -4,3k + lS,k -2,k + 4)T ,山条件知[”,% ] = 0。
大学物理学课后习题4第四章答案

[答案:D]
4.2 填空题 (1)一质点在 X 轴上作简谐振动,振幅 A=4cm,周期 T=2s,其平衡位置
取作坐标原点。若 t=0 时质点第一次通过 x=-2cm 处且向 x 轴负方向运动,则 质点第二次通过 x=-2cm 处的时刻为__ __s。
(3) t2 5s 与 t1 1s 两个时刻的位相差;
解:(1)设谐振动的标准方程为 x Acos(t 0 ) ,相比较厚则有:
A 0.1m,
8 ,T
2
1 4
s,
0
2
/3
又
vm A 0.8 m s1 2.51 m s1
am 2 A 63.2 m s2
(2)
Fm mam 0.63N
(1) x0 A ;
(2)过平衡位置向正向运动; (3)过 x A 处向负向运动;
2
(4)过 x A 处向正向运动. 2
试求出相应的初位相,并写出振动方程.
解:因为
v
x0 A cos0 0 Asin
0
将以上初值条件代入上式,使两式同时成立之值即为该条件下的初位相.故有
1
x Acos( 2 t ) T
103
(
)2
0.17
4.2
103
N
2
方向指向坐标原点,即沿 x 轴负向.
(2)由题知, t 0 时,0 0 ,
t t时
x0
A ,且v 2
0, 故 t
3
∴
t
3
/
2
2s 3
金融学-习题4-金融市场

金融学期末考试复习资料习题四第五章金融市场一、单选题8、某债券以120元买入,9个月后以125元卖出,则投资该债券以单利计算的年收益率为。
A 6.25%B 6.67%C 4.17%D 5.56%9、某剪息票债券以120元买入,9个月后以122元卖出,中间获得利息10元,则投资该债券以单利计算的年收益率为。
A 13.33%B 12%C 10%D 8.33%10、某债券以120元买入,6个月后以125元卖出,则该债券以复利计算的年收益率为。
A 4.17%B 8.51%C 8.33%D 5.56%11、上海证券交易所的组织形式是。
A 公司制B 合作制C 会员制D 合伙制13、某股票的年收益率为0.2元/股,当前市场利率为5%,则该股票的理论价格应为。
A 1元B 4元C 8元D 10元14、一般来说,市场利率提高,股票价格。
A 上升B 下降C 波动 D不动20、金融机构之间发生的短期临时性借贷活动叫。
A贷款业务 B票据业务C同业拆借 D再贴现业务21、股票是代表股份资本所有权的证券,它是一种。
A固定资本 B可变资本 C实际资本D虚拟资本22、银行发行的金融债券属。
A直接融资B间接融资 C多边融资 D混合融资28、市场利率的变化对年期债券的价格影响最大。
A 2B 5C 10D 20二、多选题1、直接融资市场包括。
A 股票市场B企业债券市场 C 国债市场 D 存贷款市场2、金融市场的功能包括。
A 资金融通功能B 资金积累功能C 风险分散功能D 信息集散功能3、金融市场的构成要素包括。
A 交易对象B 交易场所C 交易主体D 交易工具E 交易价格4、金融市场的参与者通常包括。
A 居民B 政府C企业 D 金融机构5、按融资期限的长短,金融市场可分为。
A 现货市场B 期货市场C 货币市场D资本市场6、下列属于货币市场的有。
A 贴现市场B 同业拆借市场C 基金市场D 大额定期存单市场E 商业票据市场7、债券的发行价格可分为。
《随机过程及其在金融领域中的应用》习题四答案

第四章 习题41、对泊松过程{},0t N t ≥(1)证明:当s t <时,{}1,0,1,,kn ks t n s s P N k N n k n k t t -⎛⎫⎛⎫⎛⎫===-= ⎪⎪⎪⎝⎭⎝⎭⎝⎭(2)当2λ=时,试求:()()()112112;1,3;21P N P N N P N N ≤==≥≥(3)设顾客到达某商店是泊松事件,平均每小时以30人的速度到达。
求下列事件的概率:相继到达的两顾客的时间间隔为大于2分钟、小于2分钟、在1分钟到3分钟之间。
答:(1)证明:{}()()()()()()()()()()()()()()()()()()(),,!!!!!!!1!!s t s t s s t s s t t t t n kkt s sk n kn k nk n ktn kk n kk nP N k N n P N k N n k P N k P N n k P N k N n P N n P N n P N n t s s e ek n k s t s n k n k t t t e n n s t s n s s k t k n k t t λλλλλλλλλλ------------====-==-========-⎡⎤⎣⎦--==--⎛⎫⎛⎫⎛⎫==- ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭(2)()()()()()()()()11110121112222201211120!1!2!225P N P N N N e e e e e e e λλλλλλλ-------≤==+=+==++-=++=()()()()12121224111,31,3112224P N N P N N P N P N ee e----=====-=====()()()()()()()()()()111111121112112,122111121011311101P N N P N P N N P N P N P N P N P N e P N P N e --≥≥≥≥≥==≥≥-<-=-=-===-<-=-(3) 解法一:顾客到达事件间隔服从参数为λ的指数分布:()()()30,03030,0x x Z Z f t e x f t e x λλλ--=≥=⇒=≥①()30301111303023030106030x x P Z e dx e e e ∞∞----⎧⎫>===--=⎨⎬-⎩⎭⎰②()11303011303000230301116030x x P Z e dx e e e ----⎧⎫<===--=-⎨⎬-⎩⎭⎰ ③1131133030202022221160601330301606030x x P Z e dx e e e e e ------⎛⎫⎧⎫<<===--=-⎨⎬ ⎪-⎩⎭⎝⎭⎰解法二:()3030==0.560λ∴平均每小时有人到达人/分钟根据齐次Poisson 过程的到达时间间隔{},1,2,n X n =是独立同分布于均值为1λ的指数分布的,故可有: 相继到达的顾客的时间间隔大于2分钟的概率为:()12t n P X e e λ-->== 相继到达的顾客的时间间隔小于2分钟的概率为:()1211t n P X e e λ--<=-=-相继到达的顾客的时间间隔在1分钟到3分钟之间的概率为:()()()()1.50.50.5 1.5133111n n n P X P X P X e e e e ----<<=<-<=---=-2、{},0t N t ≥是强度为λ的泊松过程。
MATLAB 习题 4

MATLAB 习题第1章MATLAB R2010环境2、在命令窗口中输入:>>a=2.5>>b=5*6>>c=[a b]写出在命令窗口中的运行结果4、用“Format”命令设置数据输出格式,(format long 或format long g)将pi显示为3.14159265358979,(formatshort e)将pi显示为3.1416e+0006、在工作空间查看变量的变量名、数据结构、类型、大小和字节数,打开数组编辑器窗口修改第2题的变量c元素。
修改变量c元素:8、输入变量a=5.3,b=[1 2;3 4],在工作空间中使用who、whos、exist和clear命令,并用save 命令将变量存入“F:\exe0101.mat”文件。
A=[1,2;3,4];inv(A)rank(A)det(A)A^3[v,d]=eig(A)结果:>>ans =-2.0000 1.00001.5000 -0.5000ans =2ans =-2ans =37 5481 118v =-0.8246 -0.41600.5658 -0.9094d =-0.3723 00 5.3723 >>命令窗口图片:接上一张图片:工作空间窗口图片:第2章MATLAB数值计算1、选择和填空(1)下列变量名中的(A)是合法变量。
A. char-1,i,jB. x*y,a.1C. x\y,a1234D. end,1bcx(2)已知x为1个向量,计算其正弦函数的运算为(C)。
A. SIN(X)B. SIN(x)C. sin(x)D. sinx(3) 已知x为1个向量,计算ln(x)运算为(B)。
A. ln(x)B. log(x)C. Ln(x)D.log10(x)(4) 若a=2.4,使用取整函数得出3,则该取整函数名为(C)。
A. fixB. roundC. ceilD. floor(5) 已知a=0:4, b=1:5, 下面的运算表达式出错的为(D)。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第四章 控制系统的稳定性复习及习题课1 基本概念(1) 稳定性的定义 外部稳定;内部稳定;李氏稳定;渐近稳定;大范围渐近稳定;不稳定。
(2) 李氏第二法的基本判据 (3) 线性系统李氏函数的求法 2 基本要求(1)理解和掌握有关系统稳定性的基本概念和基本定理,包括外稳和内稳的定义、相关定理及两者的关系; (2)理解和掌握李氏第二法有关线性系统稳定性的判别方法和一般步骤,了解非线性系统的李氏分析方法。
(3)学会并掌握利用MATLAB 判断系统的稳定性。
注意:应从物理概念上理解稳定性的定义,不被书中的数学描述所迷惑。
可以将它的数学描述类比于高等数学中的极限定义,微分定义。
思考题1怎样判定二次函数的正定、负定、半正定、半负定? (线性代数内容复习)(赛尔维斯特判据)[]11121n 121222n 212n 111n2nn n p p p p p x x x p p x n nTij i j i j n x p V(x )p x x x Px x p ==⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦∑∑ , 根据赛尔维斯特(Sylvester),可以通过判断对称矩阵P 的定号性来确定二次型函数的定号性。
11111122212211121n 21222n 1n2nn p 0p p p p p p 0 p p n n p p det p ...p det p ∆=⎛⎫⎡⎤∆=> ⎪⎢⎥⎣⎦⎝⎭⎛⎫⎡⎤ ⎪⎢⎥ ⎪⎢⎥∆=> ⎪⎢⎥ ⎪⎢⎥ ⎪⎣⎦⎝⎭ 则(1) 矩阵P 正定的充要条件是所有的顺子主子式都是正的,即012i ,(i ,,...n )∆>=; (2) 矩阵P 负定的充要条件是00i ,i ,i >⎧∆⎨<⎩偶奇(3) 矩阵P 半正定的充要条件是01210i ,i ,,..n ,i n ≥=-⎧∆⎨==⎩;(4) 矩阵P 半负定的充要条件是00i ,i ,i ≥⎧⎪∆≤⎨⎪⎩偶奇=0,i=n2 什么是外部稳定性?什么是内部稳定性?外部稳定性和内部稳定性性有什么关系?什么是李氏稳定性?什么是平衡点?根据输入输出描述来研究系统的稳定性性属于外部稳定性分析。
对输入的不同性质可引出不同的稳定性定义。
普通应用的是有界输入有界输出(BIBO )稳定。
对于零初始状态的线性系统BIBO 稳定的充要条件是对任意有界输入,其输出是有界的。
依据状态空间描述来研究系统的稳定性属于内部稳定性分析。
如果由任意非零初始状态引起的零输入响应有界且渐近,即0o t lim x (t )→∞=成立,则称线性系统在初态时刻内部稳定。
内部稳定性可根据状态转移矩阵直接判断。
即充要条件是00t lim (t,t )→∞Φ=对于线性定常系统,若内部稳定即渐近稳定,则系统必为BIBO 稳定即外部稳定,但BIBO 稳定不能保证系统为内部稳定。
在系统完全能控和完全能观测的条件下,系统外部稳定与系统内部稳定等价。
李氏稳定性指的是系统在平衡状态下受到扰动时,经过“足够长的时间”后系统恢复到平衡状态的能力。
即李氏稳定性是相对于系统的平衡状态的。
对于线性定常系统,当A非奇异时,存在唯一的平衡点即原点。
而非线性系统的平衡点可能有多个。
平衡态可由0x(t )= 求出。
由于非零平衡点总可由坐标变换将其移到坐标原点,故大多将坐标原点作为平衡点来研究。
3 古典控制理论中的系统稳定性与李雅普诺夫意义下的稳定性有什么区别? 答:(1)古典控制理论中的稳定指的是指输入输出稳定性,与系统状态无关;而李雅普诺夫意义下的稳定性是指系统的内部稳定性,反映了系统状态在偏离平衡状态后,是否仍能保持在平衡状态附近、甚至回到平衡状态的系统能力。
(2) 对于古典控制理论中的稳定性是利用系统的传递函数定义的,因此必须要假定系统的初始条件为零。
对象是线性定常单输入单输出系统,采用的方法是判断系统的极点位置等,仅适用于线性定常系统稳定性分析。
(3) 李雅普诺夫意义稳定性理论适合于线性和非线性系统,时变和时不变系统,多变量系统。
通过分析系统能量的变化来判断系统的稳定性;李氏方法不仅适合于分析,而且更重要 的是可用于控制系统的设计,尤其是在非线性系统和自适应控制系统常常是将稳定性作为首要的控制目标。
4 请说出李雅普诺夫意义下稳定、渐近稳定、大范围渐近稳定的区别和联系,并说明它们的几何意义。
答:李雅普诺夫意义下稳定、渐近稳定、大范围渐近稳定的区别在于: (1) 李雅普诺夫意义下稳定指的是当系统初始状态在平衡状态的一个小范围内时,其之后的状态也一定在平衡状态附近。
(2) 渐近稳定性首先要求是稳定的,而且随着时间的推移,系统的状态回复到平衡状态。
(3) 大范围渐近稳定不再要求初始状态在平衡状态附近,而可以是任意的。
因此,对任意初始状态出发的轨线都逼近于平衡状态,当然,系统的平衡状态也是惟一的。
李雅普诺夫意义下稳定、渐近稳定、大范围渐近稳定的联系是:大范围渐近稳定⇒渐近稳定⇒李雅普诺夫意义下稳定。
几何意义:李雅普诺夫意义下稳定指的是,当时间无限增加时,从球域内出发的状态轨线不会超出球域;李雅普诺夫意义下的渐近稳定指的是,当时间无限增加时,从球域内出发的状态轨线不仅不会超出球域,而且最终收敛到原点;大范围渐近稳定指的是,当为整个状态空间时,系统的状态都会随着时间的推移趋向于原点。
因此,大范围渐近稳定是最高的,渐近稳定次之,李氏稳定再次之。
5 李雅普诺夫稳定性的定义是什么?它适用于哪些系统?答:考虑系统的平衡状态,如果对任意给定的,存在一个(与和初始时刻有关),使得从球域内任一初始状态出发的状态轨线始终都保持在球域内,则平衡状态称为是李雅普诺夫意义下稳定的。
进一步,如果平衡状态是李雅普诺夫意义下稳定的,并且当时,始于原点邻域中的轨线,则平衡状态称为在李雅普诺夫意义下是渐近稳定的。
它既适用于线性系统,也适用于非线性系统,既适用于时变系统,也适用于时不变系统,既适用于连续系统,也适用于离散系统。
6 李雅普诺夫稳定性定理的物理意义是什么?答: 李雅普诺夫稳定性定理的物理意义是:针对系统引入一个虚拟的能量函数(即李雅普诺夫函数),其本身要求是正定的。
该能量函数沿系统轨线关于时间的导数是负定的表明了:当系统运动时,其能量随时间的推移而持续地减少,直至消耗殆尽,则系统的状态就回到平衡状态,从而系统是渐近稳定的。
7 什么是李氏第一法? 李氏第一法又称间接法。
它通过系统状态方程的解或根来判别系统的稳定性。
李氏第二法也称直接法。
它不求出状态方程解,而是通过构造一个广义的能量函数来直接判别系统的稳定性。
李氏第一法的基本结论:(1)线性定常系统渐近稳定的充要条件是系统矩阵A的特征根均具有负实部。
(2)非线性定常系统可在平衡点用微偏线性化的方法变成线性定常系统。
微偏线性化就是进行泰勒级数展开,将展开式的一次近似式构成A(雅可比矩阵),再用线性定常系统的判据来判别。
(3)线性时变系统可用范数的概念来分析稳定性。
8 如何根据李氏定理(第二法)判别系统的稳定性?对于任何系统包括线性系统和非线性系统,首先假定能找到一个广义的能量函数,也就是李氏函数v(x,t ),如果0v(x,t )>即正定,而判定的充分条件是(1)0000v(x,t )v(x,t ),v(x,t )(x )<≤≠≠ ,渐近稳定(2)0v(x,t )≤ ,李氏稳定 (3)0000v(x,t )v(x,t ),v(x,t )(x )>≥≠≠ ,不稳定。
因此,判定系统的稳定性的关键在于能否找到或者构造一个合理的李氏函数。
这是李氏第二法的关键。
李氏定理并没有提供构造李氏函数的一般方法。
所以,尽管李氏定理原理上是简单的,但实际应用中尤其是在非线性系统中并非易事。
不过,对于线性系统和某些非线性系统,已经有了一些可行的构造李氏函数的方法。
9 构造李氏函数的方法(即李氏第二法的应用) 9.1 线性定常系统的李氏函数(掌握)对任给的一个正定实对称矩阵Q ,存在一个正定的矩阵P ,满足+=-T A P PA Q ,则标量函数v(x )=T x Px 就是要找的具有二次型形式的李氏函数,而v(x )=T x (-Q)x 。
为计算方便,常取Q =I ,+=-T A P PA I ,只要根据检验P 是否正定,就可判定系统的稳定性了。
9.2 线性时变连续系统(了解)在平衡点附近,大范围一致渐近稳定的充要条件是:对于任意给定的连续实对称、一致有界和一致正定的时变矩阵Q(t),存在一个连续的实对称、一致有界和一致正定的矩阵P(t),使得 T P(t)=-A (t)P(t)-P(t)A(t)-Q(t),则李氏函数就是[]v x(t ),t =T x (t)P(t)x(t)。
9.3线性定常离散系统(理解)在平衡点附近,大范围一致渐近稳定的充要条件是:对于任意给定的正定的实对称矩阵Q ,存在一个正定的实对称矩阵P ,满足T G P G -P =-Q ,那么李氏函数就是[]v x(k )k k =T x ()Px()。
9.4 非线性系统的李氏函数构造法(介绍)非线性不同于线性系统,没有统一的公式可用,但可针对具体问题具体分析。
针对非线性关系是解析的单值函数,阶次不太高(4阶以下),也有几种实用的工程方法。
一种是一次近似法。
即在平衡点附近线性化后,用其一次近似式,变成线性定常系统。
而变量梯度法也是一种应用性很强的方法。
10 如果一个系统的李雅普诺夫函数确实不存在,那么是否能断定此系统不稳定?为什么? 答:如果一个系统的李雅普诺夫函数确实不存在,我们也不能断定此系统不稳定,因为李雅普诺夫稳定性定理给出的稳定性条件仅仅是充分的,而非必要。
当系统是线性时,其条件就是充分必要的. 李雅普诺夫函数的形式并不是唯一的,其中最简单的形式是二次型函数,二次型函数一定适合线性系统,但对非线性系统,不一定都是这种最简单的形式。
习题四参考答案4-1试确定下列二次型是否为正定的。
(1)2221231223134262V(x )x x x x x x x x x =+++-- (2)222123122310462V(x )x x x x x x x =---++ (3)222123122313104224V(x )x x x x x x x x x =+++-- 解:(1) 二次型写成矩阵形式[]1123231 1 -1 1 4 -3-1 -3 1T x V(x )x Px x x x x x ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦,矩阵P 的三个顺序主子式分别为: 123101 1 -11 10 1 4 -301 4-1 -3 1det ,det ∆=>⎛⎫⎡⎤⎛⎫⎡⎤ ⎪⎢⎥∆=>∆=< ⎪⎢⎥ ⎪⎢⎥⎣⎦⎝⎭ ⎪⎢⎥⎣⎦⎝⎭ 根据塞尔维斯特准则,V (x )不是正定的。