问题标签 [transfer-function]

我正在尝试使用 python 绘制奈奎斯特图。但是，我无法正确处理。

这是我正在使用的代码：

代码输出此图：

但情节应该是这样的（来自 wolfram）：

为什么控制包中的奈奎斯特图是错误的？我怎样才能正确？

非常感谢。

python-control我试图演示如何使用系统传递函数和模块的定义来“解决”（模拟解决方案）微分方程初始值问题（IVP） 。事实上，我真的是一个关于控制的新手。

我以这个简单的微分为例：y'' - 4y' + 13y = 0，具有这些初始条件：y(0) = 1和y'(0) = 0。

我手动实现了这个传递函数： Y(s) = (s - 4)/(s^2 - 4*s + 13).

因此，在 Python 中，我正在编写这段代码（请注意，这是这里看到y_ans的差分 IVP 的答案）：

这段代码让我得到这个图表：

但是当我在前面插入一个负号时yout，我得到了一个我想要的匹配：

我究竟做错了什么？Python 控制文档对我来说不是很清楚。另外，我不知道我对X0参数 for的解释control.forced_response是否正确。是否有可能按照我的意图做到这一点？

欢迎任何能对这个主题有所了解的人做出贡献。

编辑

设置X0 = [0,0]给了我这个图表：

你好，希望你们一切都好。

我有两个分母不同的符号分数，我想知道它们相减的结果。我尝试了简化和简化分数函数；但是他们没有产生我想要的结果。我希望结果是 A/B 的形式。

这是我得到的分数和结果：

你能告诉我我需要使用哪个具体功能来获得A/B形式的T2吗？

我不明白为什么scipy.signal.ss2tf()andscipy.signal.StateSpace().to_tf()会给出（相同的）意外结果。

例子：

scipy 的结果是

在 Matlab 中，结果是

似乎分母总是对的，我试过其他例子，但是传递函数的分子不对应。

我是否以不正确的方式使用 scipy？

（另一个例子）

我正在编写动态系统模拟（固定步骤，非实时；它在我的桌面上运行），我想通过 scipy.signal 提供的工具（例如dlsim）。对于这些组件，我知道它们以经典形式的传递函数表示。

使用 scipy.signal，“静态”模拟传递函数的输出非常简单直接，即在已知时间和输入向量的情况下；另一方面，我找不到在每个模拟步骤中计算它的方法。我的模拟器还包括一些闭环控制器，因此随着模拟的进行，输出会发生动态变化。

有任何想法吗？

PS我发现this thread似乎很相似，但我必须承认我不明白作者给出的解决方案......：如何在python中模拟一个传递函数的一步

我测量了系统输入（粉红噪声）和输出的时间信号。我需要计算一个系统的 TF 和它的脉冲响应。

为了计算 TF，我可以直接在 Matlab 中执行TF = tfestimate(in, out)或计算fft两者然后执行TF = fft(out)./fft(in)。在这一点上，一切对我来说都很清楚和预期，但是当我通过这样做计算系统的脉冲响应时，ifft(fft(out)./fft(in))我得到了最后有尾巴的脉冲响应。我在想我的测量可能有问题，并决定在 matlab 中重复实验，合成指数啁啾，并通过滤波器通过以下脉冲响应计算。然后我在结尾得到了“相同”的故事：结尾 带有“尾巴”的冲动反应

请看下面的代码

1. 所以我想知道这个故事的本质是什么？也许我的代码有错误？我是谷歌并寻找这个问题/解释，但没有找到任何东西

2. 它可以与我如何定义 nfft 以及点数的正确选择有关吗？我们可以说如果 nfft 定义不正确，脉冲响应将不正确吗？

3. 当我试图恢复从 获得的频谱的对称部分时tfestimate，我是否正确地做：？

谢谢！

我正在尝试绘制一些注释以配合我的阶跃函数图。我目前有这些图表，我一直在试图弄清楚如何绘制指向垂直线的水平箭头。我还需要指向水平线的标记垂直线。

我附上了一张图片，显示（红色）我的意思。我已经尝试过 annotation() 函数，但是将箭头放在我想要的位置确实很痛苦。如果有人不介意解释如何使用该功能或​​我想要实现的替代方法，那将是惊人的！

编辑：有没有办法编辑箭袋箭头大小？

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

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

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

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

我得到了一个阶跃响应图，并从中获得了传递函数。该循环包括一个工厂和一个 PID 控制器，我知道产生输出图的 PID 值。我通过应用阶跃输入验证了我找到的传递函数是正确的，并且输出图与我看到的相匹配。

我的代码如下所示：

由此，我需要找到工厂传递函数，以便我可以使用它来找到最佳 PID 值，而不是现在使用的那些。

我正在做一个附带项目，一旦光电传感器检测到运动，直流电机 (2x) 将打开 5 秒钟然后关闭。如果运动停止，则直流电机将关闭。如何确定 LDR 传感器“光/暗传感器”的传递函数或计算其传递函数？我需要光电池控制窗帘项目的帮助，它是一个闭环反馈系统，我创建了一个检测光并驱动直流电机朝一个方向转动的电路和另一个检测黑暗驱动直流电机朝一个方向转动的电路