Matlab 版 (精华区)

发信人: seesea (马虎学究), 信区: Matlab
标  题: [合集]有没有直接求状态转移矩阵的函数？
发信站: 哈工大紫丁香 (2003年03月06日13:53:16 星期四), 站内信件


────────────────────────────────────────
 aaaaaaaaaaai (何时做老大)            于 2000年09月20日18:26:46 星期三 说道:

请问有没有直接求状态转移矩阵的函数？

────────────────────────────────────────
 xuedy (dingyu)                       于 Mon Sep 25 01:39:24 2000) 说道:

【 在 aaaaaaaaaaai (何时做老大) 的大作中提到: 】
: 请问有没有直接求状态转移矩阵的函数？
符号运算工具箱可以求状态转移阵，但效果不佳。算二阶系统马马虎虎，
例
A=[0,1; -3 -4]
syms t
expm(-A*t)
: 
结果
ans =
[  3/2*exp(t)-1/2*exp(3*t), -1/2*exp(3*t)+1/2*exp(t)]
[  3/2*exp(3*t)-3/2*exp(t), -1/2*exp(t)+3/2*exp(3*t)]
试试三阶系统是否能解？
Mathematica解这个问题稍好，但别忘了Abel定理，一般4阶以上方程没有
解析解，故没太大用途
DY
http://xuedy.363.net

────────────────────────────────────────
 seasoul (醉过才知酒浓)               于 2000年09月27日21:30:37 星期三 说道:

不巧，我求的正是四阶的，比较难算，故求捷径！
【 在 xuedy (dingyu) 的大作中提到: 】
: 【 在 aaaaaaaaaaai (何时做老大) 的大作中提到: 】
: : 请问有没有直接求状态转移矩阵的函数？
: 符号运算工具箱可以求状态转移阵，但效果不佳。算二阶系统马马虎虎，
: 例
: A=[0,1; -3 -4]
: syms t
: expm(-A*t)
: : 
: 结果
: ans =
: [  3/2*exp(t)-1/2*exp(3*t), -1/2*exp(3*t)+1/2*exp(t)]
: [  3/2*exp(3*t)-3/2*exp(t), -1/2*exp(t)+3/2*exp(3*t)]
: 试试三阶系统是否能解？
: Mathematica解这个问题稍好，但别忘了Abel定理，一般4阶以上方程没有
: 解析解，故没太大用途
: DY
: http://xuedy.363.net
: 

────────────────────────────────────────
 xuedy (dingyu)                       于 Thu Sep 28 06:48:42 2000) 说道:

【 在 seasoul (醉过才知酒浓) 的大作中提到: 】
: 不巧，我求的正是四阶的，比较难算，故求捷径！
: 
4阶的理论上是有解析解的。但用MATLAB的符号运算工具箱好象不现实，不妨试试
Mathematica, 若A矩阵含有负号而不是全部给出的数值，则困难得出的解将有几页，
你要有心理准备啊。
定宇
: 【 在 xuedy (dingyu) 的大作中提到: 】
: : 【 在 aaaaaaaaaaai (何时做老大) 的大作中提到: 】
: : : 请问有没有直接求状态转移矩阵的函数？
: : 符号运算工具箱可以求状态转移阵，但效果不佳。算二阶系统马马虎虎，
: : 例
: : A=[0,1; -3 -4]
: : syms t
: : expm(-A*t)
: : : 
: : 结果
: : ans =
: : [  3/2*exp(t)-1/2*exp(3*t), -1/2*exp(3*t)+1/2*exp(t)]
: : [  3/2*exp(3*t)-3/2*exp(t), -1/2*exp(t)+3/2*exp(3*t)]
: : 试试三阶系统是否能解？
: : Mathematica解这个问题稍好，但别忘了Abel定理，一般4阶以上方程没有
: : 解析解，故没太大用途
: : DY
: : http://xuedy.363.net
: : 
: 

────────────────────────────────────────
 seesea (小学究)                      于 2000年09月28日20:03:14 星期四 说道:

听说Maple作符号运算的能力很突出的，
运行显示也比较友好。
【 在 xuedy (dingyu) 的大作中提到: 】
: 【 在 seasoul (醉过才知酒浓) 的大作中提到: 】
: : 不巧，我求的正是四阶的，比较难算，故求捷径！
: : 
: 4阶的理论上是有解析解的。但用MATLAB的符号运算工具箱好象不现实，不妨试试
: Mathematica, 若A矩阵含有负号而不是全部给出的数值，则困难得出的解将有几页，
: 你要有心理准备啊。
: 定宇
: : 【 在 xuedy (dingyu) 的大作中提到: 】
: : 
: 

────────────────────────────────────────