Represent Simulink Integrator block as Matlab Function

28 次查看(过去 30 天)
Hi.
I need to implement the following behavior :
The Integrator and my_Integrator blocks have to be equivalent I/O.
How should I write the Matlab Function ?
Thanks for any reply.

采纳的回答

Ryan G
Ryan G 2012-12-4
As this looks like a homework problem, I can't answer directly. However I will point you in the direction of persistent variables.

更多回答(4 个)

Azzi Abdelmalek
Azzi Abdelmalek 2012-12-7
编辑:Azzi Abdelmalek 2012-12-8
I don't know why do you need this, maybe if you explain exactly what you need, there is better way
  9 个评论
Azzi Abdelmalek
Azzi Abdelmalek 2012-12-8
编辑:Azzi Abdelmalek 2012-12-8
Ok, I see, If T is constant, you must then set, in model configuration parameters your fixed step time to T, and also your step block sample time to T. In this case you don't need a clock.
function y = fcn(u)
persistent uold yold
T=0.01;
if isempty(uold)
uold=0;yold=0;
end
y = u*T+yold-(u-uold)*T/2
yold=y;uold=u;
River Rock
River Rock 2012-12-8
Changing the sample time of the Step block to 0.01 removed the previous offset. Thanks

请先登录,再进行评论。


Guy Rouleau
Guy Rouleau 2012-12-5
This is not a good idea. The MATLAB function is not designed for this purpose.
  1 个评论
Edward Rodriguez
Edward Rodriguez 2020-8-10
Excuse me, so, What would be a good idea to implement numerical integration methods in blocks in Simulink?

请先登录,再进行评论。


River Rock
River Rock 2012-12-5
编辑:River Rock 2012-12-6
My main goal is to implement the differential equations of a physical system using a single Matlab Function. As the sums and gains were easy to represent, I couldn't find any alternative for the integration.
  4 个评论
Ryan G
Ryan G 2012-12-5
What you have written is close it would be more like:
y(z) = yOld+u(z)/SampleTime
You cannot use the ODE solver in the MATLAB function block.
River Rock
River Rock 2012-12-7
Any idea on how to get rid of this offset ?
Code looks like:
function y = fcn(u)
%#codegen
T=0.01;
persistent yOld;
persistent uOld;
if (isempty(yOld))
yOld = 0;
end
if (isempty(uOld))
uOld = 0;
end
y = yOld + (T/2)* (u + uOld);
%y=yOld + u*T;
yOld = y;
uOld = u;

请先登录,再进行评论。


River Rock
River Rock 2012-12-5
Can anybody suggest a better way of implementing the numerical integration ? The code has to be written inside the Matlab Function Block though.

类别

Help CenterFile Exchange 中查找有关 Event Functions 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by