[求助] 线性系统问题求助

» 猜你喜欢

谈谈两天一夜的“延安行” 已经有6人回复
博士申请都是内定的吗？已经有12人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助线性系统中的一些概念已经有4人回复
【求助】非线性方程组的求解问题已经有6人回复
求教：线性系统与线性微分方程的关系。非线性系统能用传递函数来描述吗？已经有7人回复
求助时间序列的最大李雅普诺夫指数问题已经有14人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
求助投稿期刊杂志已经有6人回复
急切求助！一个mathematica求解一个简单一元非线性方程的问题已经有4人回复
求助液相方法线性考查问题！急！！！！已经有27人回复
求助有关Matlab有约束非线性最优化问题已经有9人回复
求助SAS 多元非线性回归拟合编程已经有3人回复
【求助】质量控制分析方法验证指导原则中线性与范围问题的讨论已经有4人回复
【求助】求助《信号与线性系统》已经有11人回复
【求助】线性极化曲线的若干问题。已经有14人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】有关气相色谱检测线性范围的问题已经有3人回复
【求助】请教：液相精密度、线性问题？已经有7人回复
【讨论】线性时不变系统的脉冲响应已经有7人回复
【求助】HPLC方法学考察具体内容，如何做？已经有19人回复

在绝望中寻找希望，人生终将辉煌。

1楼2013-08-11 13:41:12

红星爱锅头

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 121.9
帖子: 54
在线: 8.8小时
虫号: 2096284
注册: 2012-10-30
专业: 控制理论与方法

【答案】应助回帖

感谢参与，应助指数 +1

我觉得应该可以从从稳定矩阵的性质和矩阵范数的定义两个角度综合考虑，希望能够提供一定的思路。

2楼2013-08-11 16:12:36

dfh3919

金虫 (小有名气)

应助: 6 (幼儿园)
金币: 1235.4
帖子: 117
在线: 106.5小时
虫号: 701453
注册: 2009-02-15
专业: 控制理论与方法

【答案】应助回帖

★
感谢参与，应助指数 +1
lipeng0327: 金币+1 2013-08-11 22:13:39

将F矩阵对角化，然后利用F与对角矩阵（其中主对角线元素为F矩阵的特征值）的等价性就可以证明了。

3楼2013-08-11 21:49:38

liuyang5264

铁虫 (初入文坛)

应助: 1 (幼儿园)
金币: 455.4
帖子: 20
在线: 6.9小时
虫号: 2589207
注册: 2013-08-09
专业: 控制理论与方法

【答案】应助回帖

感谢参与，应助指数 +1

将F化为Jordan标准型，按照子块的几种不同类型分别讨论就可以了。

4楼2013-08-12 09:45:37

gaofeng79

银虫 (正式写手)

应助: 19 (小学生)
金币: 5016.6
红花: 4
帖子: 943
在线: 32.8小时
虫号: 1766233
注册: 2012-04-19
性别: GG
专业: 交通工程

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
lipeng0327: 金币+2 2013-08-12 13:11:16

将F分解为F=U*S*V'，U和V为单位矩阵，S为F的奇异值构成的对角阵。
这样就有|e^Ft|=|U*(e^St)*V|<=|U|*|e^St|*|V|。由于U和V为单位阵，所以有|e^Ft|<=|e^St|。
S是对角阵，剩下的应该明白了吧。

发动机控制;快速原型;硬件在环测试;汽车主动安全;

5楼2013-08-12 12:16:10

lipeng0327

木虫 (正式写手)

水手

应助: 15 (小学生)
金币: 3965.6
散金: 662
红花: 3
帖子: 777
在线: 584.6小时
虫号: 837280
注册: 2009-08-31
性别: GG
专业: 控制理论与方法

引用回帖:

5楼: Originally posted by gaofeng79 at 2013-08-12 12:16:10
将F分解为F=U*S*V'，U和V为单位矩阵，S为F的奇异值构成的对角阵。
这样就有|e^Ft|=|U*(e^St)*V|<=|U|*|e^St|*|V|。由于U和V为单位阵，所以有|e^Ft|<=|e^St|。
S是对角阵，剩下的应该明白了吧。

|e^Ft|=|U*(e^St)*V|此式不成立的

在绝望中寻找希望，人生终将辉煌。

6楼2013-08-12 18:41:26

lipeng0327

木虫 (正式写手)

水手

应助: 15 (小学生)
金币: 3965.6
散金: 662
红花: 3
帖子: 777
在线: 584.6小时
虫号: 837280
注册: 2009-08-31
性别: GG
专业: 控制理论与方法

顶顶

在绝望中寻找希望，人生终将辉煌。

7楼2013-08-12 23:27:54

gaofeng79

银虫 (正式写手)

应助: 19 (小学生)
金币: 5016.6
红花: 4
帖子: 943
在线: 32.8小时
虫号: 1766233
注册: 2012-04-19
性别: GG
专业: 交通工程

【答案】应助回帖

★
lipeng0327: 金币+1, ★有帮助, 利用级数展开时，UV~=I，只是化为约旦型可以， F=Q*J*inv(Q), J为约旦块 2013-08-13 16:45:33

引用回帖:

6楼: Originally posted by lipeng0327 at 2013-08-12 18:41:26
|e^Ft|=|U*(e^St)*V|此式不成立的...

请看一下自动控制原理里面关于状态方程求解和变换部分的内容。

赞一下(1人)

发动机控制;快速原型;硬件在环测试;汽车主动安全;

8楼2013-08-13 16:03:44

lipeng0327

木虫 (正式写手)

水手

应助: 15 (小学生)
金币: 3965.6
散金: 662
红花: 3
帖子: 777
在线: 584.6小时
虫号: 837280
注册: 2009-08-31
性别: GG
专业: 控制理论与方法

我的证明思路是这样的：
将F化为 F=P*J*inv(P), J为标准的约当块，P为F的广义特征向量组成的变换阵，e^Ft = P * M * inv(P)，M是可以直接由J列写的，M为上三角或下三角矩阵。以2范数为例来证明：
|| e^(Ft)|| = || P * M * inv(P)|| <=||P||*||M||*|| inv(P)||
因为||P||和|| inv(P)|| 为常数，
则  || e^(Ft)||<=k*||M||，此处 k = ||P||*|| inv(P)||
||M||=(max_ lambda(M*M'))^(1/2)  （此处lambda(.) 表示.的特征值）
而M*M'为正定对称阵，其所有的特征值都大于零，
所以有 max_ lambda(M*M')<trace(M*M'),  trace(.)表示方阵.的迹。
trace(M*M')< n*(M*M'的最大主对角元)，此处n为M*M'的维数。
M*M'的最大主对角元应该是（e^(bt))^2，此处b为F的所有特征值中具有的最大的负实部。
综合得到
|| e^(Ft)||<=k*||M||

=k*(max_ lambda(M*M'))^(1/2)
=k*(trace(M*M'))^(1/2)
<=k*(n*(M*M'的最大主对角元))^(1/2)
=k*n^(1/2)*e^(bt))
=m*e^(-alpha*t)    此处（m = k*n^(1/2), alpha = -b）

在绝望中寻找希望，人生终将辉煌。

9楼2013-08-13 17:13:01

lipeng0327

木虫 (正式写手)

水手

应助: 15 (小学生)
金币: 3965.6
散金: 662
红花: 3
帖子: 777
在线: 584.6小时
虫号: 837280
注册: 2009-08-31
性别: GG
专业: 控制理论与方法

我的结果可能很粗糙，应为下面的这一步缩放太大了
max_ lambda(M*M')<trace(M*M'), trace(.)表示方阵.的迹

在绝望中寻找希望，人生终将辉煌。

10楼2013-08-13 17:16:07