LTI连续系统零状态响应求解

方法一

利用经典解方法,通解+特解
求通解,列写特征方程,根据特征根写出方程齐次解,存在未知系数
求特解,根据输入信号f(t)写出方程特解,特解中含未知数,代入原方程求出未知数
列写零状态响应=通解+特解,代入0+时刻值零状态响应及其一阶导数值,求出未知系数

方法二

先求出仅由f(t)引起的零状态响应(求解方法参考经典解方法),再利用LTI连续系统齐次性和可加性,求出整个系统的零状态响应


0+时刻值求解方法

对于任意LTI连续系统,0-时刻,零状态响应及其一阶导数等于0,由于输入信号的存在,在0时刻的值可能会发生跳变。

根据方程两边最高阶(一般为二阶)导数的系数相等,设写出零状态响应。即可由输入信号最高阶设写出最高阶零状态响应(两者最高阶可不相等),然后对最高阶零状态进行积分运算,写出次高阶、……、一阶、原函数。

求0+时刻值:已知0-时刻值,0+减去0-等于对一阶零状态响应在(0-,0+)求积分。

求0+一阶导数值:同上,对二阶零状态响应在(0-,0+)求积分。

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
在Matlab中,可以使用控制系统工具箱中的函数`ss`将系统转化为状态空间形式,然后使用矩阵运算求解。具体步骤如下: 1. 定义系统的系数矩阵和输入输出向量: ```matlab A = [0 1; -1 -1]; % 系数矩阵 B = [0; 1]; % 输入矩阵 C = [1 0]; % 输出矩阵 D = 0; % 直通矩阵 ``` 2. 将系数矩阵和输入输出向量传入`ss`函数中,得到系统状态空间模型: ```matlab sys = ss(A, B, C, D); ``` 3. 求解输入响应: ```matlab x0 = [1; 0]; % 初始状态向量 t = 0:0.1:10; % 时间向量 [yzi, t] = initial(sys, x0, t); plot(t, yzi); ``` 4. 求解状态响应: ```matlab u = sin(t); % 输入向量 [yzs, t] = lsim(sys, u, t); plot(t, yzs); ``` 5. 求解响应: ```matlab y = yzi + yzs; plot(t, y); ``` 完整代码如下: ```matlab % 定义系统的系数矩阵和输入输出向量 A = [0 1; -1 -1]; % 系数矩阵 B = [0; 1]; % 输入矩阵 C = [1 0]; % 输出矩阵 D = 0; % 直通矩阵 % 将系数矩阵和输入输出向量传入ss函数中,得到系统状态空间模型 sys = ss(A, B, C, D); % 求解输入响应 x0 = [1; 0]; % 初始状态向量 t = 0:0.1:10; % 时间向量 [yzi, t] = initial(sys, x0, t); plot(t, yzi); xlabel('Time'); ylabel('Zero Input Response'); % 求解状态响应 u = sin(t); % 输入向量 [yzs, t] = lsim(sys, u, t); plot(t, yzs); xlabel('Time'); ylabel('Zero State Response'); % 求解响应 y = yzi + yzs; plot(t, y); xlabel('Time'); ylabel('Total Response'); ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值