matlab : R2018a 64bit
OS : Windows 10 x64
typesetting : Markdown
blog : my.oschina.net/zhichengjiu
gitee : gitee.com/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,用于仿真和基于模型的设计,值得学习。
该博文仅可用于测试与参考。