我工作的一个简单的模型,其中包括的DY / DX衍生物,但在Modelica的,我不能这个方程直接写,我可以使用组合的x=time和der(y),但我认为这是因为Modelica的语言限制的妥协.
我的问题是:是否有另一种更好的方法来描述 Modelica 中的导数?
这是代码:
model HowToExpressDerivative "dy/dx=5, how to describe this equation in Modelica?"
Real x,y;
equation
x = time;
der(y) = 5;
end HowToExpressDerivative;
Run Code Online (Sandbox Code Playgroud)
我也试过用der(y)/der(x)来表达dy/dx,但是当 x 等于 时出现错误time^2。
model HowToExpressDerivative "dy/dx=5, how to describe this equation in Modelica?"
Real x,y;
equation
x=time^2;
der(y)/der(x)=5;
end HowToExpressDerivative;
Run Code Online (Sandbox Code Playgroud)
Error: The following error was detected at time: 0
Model error - division by zero: (1.0) / (der(x)) = (1) / (0)
Error: Integrator failed to start model.
... "HowToExpressDerivative.mat" creating (simulation result file)
ERROR: The simulation of HowToExpressDerivative FAILED
Run Code Online (Sandbox Code Playgroud)
给定焓h和曲柄角,phi您可以替换dh/dphi=...为:
der(h)/der(phi)=...
Run Code Online (Sandbox Code Playgroud)
但是,即使正确,该公式也会在发动机静止时失效(der(phi)=0),因此并不理想。
另一种方法是重写公式。仔细观察公式似乎是:
dh/dphi=(\partial a/\partial T)*dT/dphi+...
Run Code Online (Sandbox Code Playgroud)
这表明它们可以重写为:
der(h)=(\partial a/\partial T)*der(T)+...
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
86 次 |
| 最近记录: |