线性化非线性感应电机模型的正向欧拉计算gydF4y2Ba

你好,gydF4y2Ba

我拼命地试图通过线性化(通过计算雅可比矩阵)来模拟感应电机模型(非线性)，然后通过欧拉正向方法将其离散化。您可以在asm_discrete temodelmatlab .m文件中找到该代码gydF4y2Ba

不知怎么的，我得不到任何正确的结果。我可以用ode45模拟非线性模型，得到正确的结果。gydF4y2Ba

我用TOOLS_jacobian文件计算线性化的状态系统矩阵。在ASM_discreteModelMatlab.m中进行模拟时，在每个时间步骤更新它gydF4y2Ba

我不知道我在线性化模型时是否做错了什么，因为线性化被定义为(从gydF4y2Ba https://de.mathworks.com/help/slcontrol/ug/linearizing-nonlinear-models.html):gydF4y2Ba

”“gydF4y2Ba

将线性化的概念扩展到动态系统，你可以写出这样的连续时间非线性微分方程:gydF4y2Ba

x”(gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba fgydF4y2Ba （gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba ，gydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba ，gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba ggydF4y2Ba （gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba ，gydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba ，gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba 。gydF4y2Ba

在这些方程中，gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba )表示系统状态，gydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba )表示系统的输入，和gydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba )表示系统的输出。gydF4y2Ba

该系统的线性化模型在工作点附近的小范围内是有效的gydF4y2Ba tgydF4y2Ba =gydF4y2Ba tgydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba 0gydF4y2Ba ）=gydF4y2Ba xgydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba 0gydF4y2Ba ）=gydF4y2Ba ugydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba 0gydF4y2Ba ）=gydF4y2Ba ggydF4y2Ba （gydF4y2Ba xgydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba ugydF4y2Ba 0gydF4y2Ba ，gydF4y2Ba tgydF4y2Ba 0gydF4y2Ba ）=gydF4y2Ba ygydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba

为了表示线性化模型，定义以工作点为中心的新变量:gydF4y2Ba

δgydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba −gydF4y2Ba xgydF4y2Ba 0gydF4y2Ba

δgydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba −gydF4y2Ba ugydF4y2Ba 0gydF4y2Ba

δgydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba −gydF4y2Ba ygydF4y2Ba 0gydF4y2Ba

用δ表示的线性化模型gydF4y2Ba xgydF4y2Ba ,δgydF4y2Ba ugydF4y2Ba ，和δgydF4y2Ba ygydF4y2Ba 当这些变量的值很小时有效:gydF4y2Ba

δgydF4y2Ba x”(gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba 一个gydF4y2Ba δgydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba +gydF4y2Ba BgydF4y2Ba δgydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba δgydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba =gydF4y2Ba CgydF4y2Ba δgydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba +gydF4y2Ba DgydF4y2Ba δgydF4y2Ba ugydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba

”“gydF4y2Ba

对我来说，不清楚该如何处理gydF4y2Ba δgydF4y2Ba x(t)因为我在x0处线性化了x0 = x(k)当我想设置的时候gydF4y2Ba δgydF4y2Ba X (k) = X (k) - x0结果是0。这就是我刚刚设置的原因gydF4y2Ba δgydF4y2Ba X (k) = X (k)正如你所看到的，我确实在理解线性化方法方面有一些困难。gydF4y2Ba