Unexpected NaN output from power block

Question

Quint 2024-12-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2172222-unexpected-nan-output-from-power-block

评论： Quint 2024-12-18

While building my model I came across the following integrator error:

Derivative of state '1' in block 'simpie/Motor & PI Current Controller/Integrator' at time 4.519405258285984 is not finite. The simulation will be stopped. There may be a singularity in the solution. If not, try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances)

After some investigation I found by setting the INF or NaN block output to warning that there are NaN values that are contaminating my feedback loop which in turn raises the integrator error. By setting the INF or NaN block output setting to warning I tried to figure out where the first NaN value comes from. This raises the following error:

An error occurred during simulation and the simulation was terminated

Caused by:

Block 'simpie/Controller/Power1' outputs 'NaN' for element 1 of output port 1 at minor time step 4.51941

From reading the documentation about the power block, it became clear that the block will output a NaN value from an NaN input. I assume that there is no NaN input into the power block since this is the first NaN output. I tried to use the isnan block to monitor the generated NaN value without succes. The block does not recognice any NaN output. The error occurs at a specific point in time where the value into the power block for the base come close to zero (See the included figure). The exponent is a constant.

Strangely the simulation works fine for t<1 which also has an input of 0 for the power block.

Does anyone have a clue what causes this problem and how to fix it?

Thank you in advance

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2024-12-17

I suggest that instead of testing using isnan() that you instead test ~isfinite() as that also tests for +/- inf

Quint 2024-12-17

Here is my simulink model file with corresponding data. I tried the isfinite() function without succes. Additionally, the problem arises when the base of the power block is almost zero. I kind of fixed the issue to round the value before the powerblock to zero when it almost reaches zero. There are still some NaN values in my simulation but now it does not crash anymore.

请先登录，再进行评论。

请先登录，再回答此问题。