8. MATLAB在信号与系统中的应用.pptVIP

例6.21 连续系统状态方程的解 线性连续系统在输入信号为零时的状态方程表示式为 ，其中x为n×1的列向量，A为n×n阶的方阵，求其在初始条件作用下的解。 解：建模 ：按线性方程理论，此齐次方程的解为 其中 此处要把方阵At作为一个整体求指数函数，就需要调用MATLAB矩阵指数函数expm。注意它与exp函数不同。 例6.21 连续系统状态方程的解 如果A是n×n阶，则expm(At)也是n×n阶，这时如果要算一系列的t值所对应的expm(At)，就不可能像标量指数那样用元素群运算方法了，必须用for循环。不仅如此，如果t的长度为nt，则状态转移矩阵F(t)=expm(At)将是一个n×n×nt的三维矩阵，需要用MATLAB中高维矩阵的概念来解决问题。举以下数字例来说明。 设A = [-2，1;-17，-4]，x0 = [3 ;4]; 求上述线性方程的状态转移矩阵F(t)及x的解。 例6.22 离散系统状态方程的响应 离散系统的状态方程为 初始条件 输入信号为阶跃函数，求前10步的解x。 解：已知x(0)和u(0)就可求出k=1时的x(1)，再依次递推求x(2)，…，x(10)。在程序中可用循环语句来完成。离散系统的计算比连续系统简单得多，只要作四则矩阵运算即可。 第六章  MATLAB在信号与系统中的应用 例6.1 连续信号的MATLAB