python - 是否有任何 Python 方法/包可用于对传递函数矩阵执行矩阵运算？

Question

对于许多工程问题，更容易将方程转换为拉普拉斯域，求解它们，然后将它们转换回时域。

在 MATLAB 中，很容易创建这些传递函数，然后使用该impulse方法进行拉普拉斯到时域的转换。这可以很容易地扩展到传递函数的线性方程组，正常求解矩阵，然后像以前一样获取脉冲。

如何在 Python 中执行这些操作？

下面是我想在 Python 中重现的 MATLAB 中的一个工作示例。

clear all; close all; clc

s = tf("s");

%% Example 1
omega = 5;
sine_s = omega/(s^2 + omega^2);
t = linspace(0, 10, 1000);
sine_t = impulse(sine_s, t);

figure
plot(t, sine_t)

%% Example 2
tf1 = 1/s;
tf2 = 3;
tf3 = -1;
tf4 = 2*s/(s+1);

A = [tf1, tf2; tf3, tf4];
b = [0; 4];
x = A\b;

x_t = impulse(x, t);

figure
plot(t, x_t)