状态空间方程的求解

 

 

% 定义系统的状态矩阵、输入矩阵、输出矩阵和直流增益矩阵
A = [-1, -2; 3, -4];
B = [1, 0; 0, 1];
C = [0, 1];
D = [0, 0];

% 定义初始状态向量和输入信号向量
x0 = [1, 2]';
t = 0:0.01:10;
u = [sin(t); cos(t)];

% 求解状态方程
x = zeros(2, length(t));
x(:, 1) = x0;
for i = 2:length(t)
    dt = t(i) - t(i-1);
    x(:, i) = expm(A*dt)*x(:, i-1) + inv(A)*(expm(A*dt)-eye(2))*B*u(:, i-1);
end

% 求解输出方程
y = C*x + D*u;

% 绘制输出响应图像
plot(t, y);
xlabel('Time (s)');
ylabel('Output');
title('Output response of the system');