我需要帮助找到系统h的脉冲响应在时间t (t)

6视图(30天)
任务:为例 E = , 脉冲响应 的系统是defifined h ( t )= C e 一个 t B + D δ ( t ) t 0,
,代表了系统输出 y ( t 当输入 u ( t )= δ ( t )。
使用命令 expm () 计算矩阵指数,然后画出脉冲响应
对于一个等距的网格 t = linspace (0500500) 。( 提示: 你只需要计算
矩阵指数一次。)解释为什么脉冲响应 h ( t ) 0 t →∞
我已经考虑矩阵的值如下:
= (3 5 1;1 0 0 0 0 1);
B = [1;0;0),
C = (0 0 1)
D = 1;
E =眼睛(3);

答案(1)

山姆翟
山姆翟 2022年11月2日
我使用了 冲动() 方法。但是你可以试着使用 expm () 方法。查找下面的例子,验证与脉冲响应。
一个= [3 5 1;1 0 0;0 1 0];
B = [1;0;0);
C = (0 0 1);
D = 1;
sys = ss (A, B, C, D)
sys = = (x1, x2) x3 x1 3 5 1 x2 0 0 x3 0 1 0 B = u1 x1 1 x2 0 x3 0 C = (x1, x2) x3 y1 0 0 1 D = u1 y1 1连续时间状态空间模型。
冲动(系统)
tFinal = 500;
冲动(sys tFinal)
帮助expm
EXPM矩阵指数。EXPM (A)的矩阵指数计算和使用比例和平方Pade逼近算法。虽然它不是这样计算,如果有一组完整的与相应的特征值特征向量V D然后[V D] = EIG (A)和EXPM (A) = V *诊断接头(exp(诊断接头(D))) / V。EXP()计算的指数中的元素。参见LOGM、SQRTM FUNM。文档expm doc expm其他用途的expm gpuArray / expm sym / expm象征性/ expm

s manbetx 845


释放

R2021a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!