simulink is giving unexpected simulation answer

Question

abdo salm 2024-9-6，19:55

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2150879-simulink-is-giving-unexpected-simulation-answer

评论： Sam Chak 2024-9-7，6:43

采纳的回答： Sam Chak

I am trying to simulate the following block diagram

which simulates the equation :

with the second integrator having the initial condition of 1.

where this equation can have the solution :

but upon simulation I got the following result:

I can notice the mathematical solution doesn't follow actual simulation where the simulation will always produce sinusoidal trying to reach stability state which is the oscillation between (inital_condition / 2) and (-inital_condition / 2). so why doesn't the mathematical solution match the simulation.

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

Sam Chak 2024-9-6，21:40

Can you you show us the scope when wo is set to 1?

abdo salm 2024-9-6，21:44

@Sam Chak, I set wo to 1 and I got a nearly stable cosine wave

请先登录，再进行评论。

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

Answer 1

Sam Chak 2024-9-6，22:22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2150879-simulink-is-giving-unexpected-simulation-answer#answer_1512504

Nope, not stable. The reason you got the high oscillations exponential decay response is because there are not enough samples for the scope to plot it out correctly.

Solution: Use a smaller time step, preferably (2π/wo)/100.

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

Sam Chak 2024-9-7，6:43

在 MATLAB Online 中打开

Hi @abdo salm

You set the time step at the Model Configuration Parameters > Solver Pane. The simulation time step size is specified using the Max step size parameter, if the variable-step solver is selected.

Or, you can set the value from the Command Window. In the following example, the Simulink model 'vdp.slx' is opened. Use the 'get_param()' command check the current value, and then the 'set_param()' command to set the desired value.

openExample('simulink_general/VanDerPolOscillatorExample');
% Check the Maximum time step size for variable-step solver
get_param('vdp', 'MaxStep')
ans = 'auto'
% Set the Maximum time step size to 0.01, preferably using this formula (2π/wo)/100
set_param('vdp', 'MaxStep', '0.01')
% Check the Maximum time step size again
get_param('vdp', 'MaxStep')
ans = '0.01'