How can I implement the following blocks and code in function block ?

1 次查看（过去 30 天）
Willim
Willim 2020-4-8
评论： Walter Roberson 2020-4-10
I would use the following integral equation as function block in my simulation. I know I can do it as what I attached but I need a way to implement it by code inside function block.
  2 个评论
Willim
Willim 2020-4-8
I have read it but you know the simlink will run the blocks at each instant but in the link the result will give you over one period. Sorry I=the link cannot help me.
Just consider M=4 ; N=2 and I have a clock start from t=0:0.1:3, I need the y(t) to have all value at each t=0:0.1:3;

请先登录，再进行评论。

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

回答（1 个）

Mehmed Saad
Mehmed Saad 2020-4-8
编辑：Mehmed Saad 2020-4-8
This is an example of matlab given in integral help
% Without For loop
fun = @(x,c) 1./(x.^3-2*x-c);
q_f = integral(@(x)fun(x,5),0,2);
q_f = -0.4605
Now i just change the time of integration and sum them
% With For loop
t = 0:0.1:2;
q =0;
for i =2:length(t)
q =q+ integral(@(x)fun(x,5),t(i-1),t(i));
end
q = -0.4605
  6 个评论
显示 4更早的评论
Willim
Willim 2020-4-10
I am using solver type ( Variable-step ) solver ode45. I am using this to have current time results
Walter Roberson
Walter Roberson 2020-4-10
With variable-step solvers, the block will be called at irregular intervals that might not match the interval that you want solutions at. The variable-step solvers can also go backwards in time (I encountered a case in which this really happened within the last few days.) ode45() invokes the function 6 times per step, using slightly different times and slightly different boundary conditions in order to make projections and check the error measurement.
In order to get results for a fixed time interval, because the function could be called for times spanning multiple intervals, you would have to emit a variable number of results -- and you would have to be prepared to have the values "revoked" if ode45 decided that the step was a failure. ode45 does not always go backwards in time when it decides that a step is a failure.
What is the reason that you are trying to replace those blocks with a MATLAB Function Block ? Are you trying to move to a discrete time system, using MATLAB Function Block to replace the ode45() work over a known interval?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 General Applications 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by