matlab-自控原理 已知A矩阵求线性定常系统的矩阵指数函数

原创
2019/02/23 18:35
阅读数 688

     matlab : R2018a 64bit
      OS : Windows 10 x64
typesetting : Markdown
       blog : my.oschina.net/zhichengjiu

code

clear
clc

% 制作时间:2016-09-26
% 制作原因:自控原理留了作业,写完了后,想验证一下写的对不对,于是写了这个小代码
% 这个代码,没有捕捉异常的工序,所以请使用的时候,小心点哦。

syms s;
A=input('请输入题目已知的A矩阵,如果不会输入,请参考模板:一个三行三列[1 2 3;4 5 6;7 8 9]:\n');
fprintf('您输入的矩阵是:');
A
sizeOfA=size(A);
fprintf('s*I-A的结果:');
sI_A=s*eye(sizeOfA(1))-A
fprintf('对s*I-A求逆矩阵:\n');
fprintf('行列式为:')
hOfsI_A=det(sI_A)
fprintf('伴随矩阵为:')
bOfsI_A=det(sI_A)*inv(sI_A)
fprintf('逆矩阵为:')
nOfsI_A=inv(sI_A)
fprintf('对sI_A的逆矩阵进行拉普拉斯逆变换的结果是:');
lnOfsI_A=ilaplace(nOfsI_A)
fprintf('感谢您的使用程序结束\n');

result

请输入题目已知的A矩阵,如果不会输入,请参考模板:一个三行三列[1 2 3;4 5 6;7 8 9]:
[0 1;-2 -3]
您输入的矩阵是:
A =

     0     1
    -2    -3

s*I-A的结果: 
sI_A =
 
[ s,    -1]
[ 2, s + 3]
 
对s*I-A求逆矩阵:
行列式为: 
hOfsI_A =
 
s^2 + 3*s + 2
 
伴随矩阵为: 
bOfsI_A =
 
[ s + 3, 1]
[    -2, s]
 
逆矩阵为: 
nOfsI_A =
 
[ (s + 3)/(s^2 + 3*s + 2), 1/(s^2 + 3*s + 2)]
[      -2/(s^2 + 3*s + 2), s/(s^2 + 3*s + 2)]
 
对sI_A的逆矩阵进行拉普拉斯逆变换的结果是: 
lnOfsI_A =
 
[   2*exp(-t) - exp(-2*t),   exp(-t) - exp(-2*t)]
[ 2*exp(-2*t) - 2*exp(-t), 2*exp(-2*t) - exp(-t)]
 
感谢您的使用程序结束
>> 

resource

  • [文档] ww2.mathworks.cn/help/matlab
  • [文档] ww2.mathworks.cn/help/simulink
  • [平台] www.oschina.net
  • [平台] gitee.com


感谢帮助 志成就 的人们。
matlab优秀,值得学习。基础知识 + 专业知识 + matlab = ?
Simulink,用于仿真和基于模型的设计,值得学习。
该博文仅可用于测试与参考。

展开阅读全文
加载中
点击引领话题📣 发布并加入讨论🔥
打赏
0 评论
0 收藏
0
分享
返回顶部
顶部