Nested functions procedure error

2013 9 9
1 个回答
回答已采纳
3 次查看（30 天）

0 个投票

Hi,
I created some nested function but it reports the error. In the first function, I calculated values for temperatures x1 and x2. Then, I defined function result that represents thermostat which should determine should thermostat be on or off in each iteration. In case that thermostat does not work in each iteration, temperature raise. Here are the functions:
function xprime = heat_sys(t,x);
params...
xprime = (system of diff equat)
function result = thermostat(x);
on = true;
off = false;
for i=1:n(end)
if(x(i,1)>=18 && x(i,1)<=20)
result(i)=on;
else if(x(i,1)<18)
result(i)=on;
else if(x(i,1)>20)
result(i,1)=off;
end;
end;
end;
end;
By the way, is this right approach to do this or I could use something different also?
Thanks in advance!

8 个评论
显示 6更早的评论 隐藏 6更早的评论

Arthur
Arthur 2013-9-9
What is the error you get?
Nikola
Nikola 2013-9-9
Maybe I expressed wrong. It is not really error but I got values that are so high. That indicate that thermostat function didn't act in each iteration. My goal is to keep x values within bounds given in if else statement.
Arthur
Arthur 2013-9-9
How do you call the two functions?
Nikola
Nikola 2013-9-9
I called only heat_sys function since thermostat function is nested inside. I want to check these temperature conditions in each iteration in order to prevent temperature to raise.
Walter Roberson
Walter Roberson 2013-9-9
Nesting a function does not cause the function to be called. MATLAB would not know what arguments to call the function on, or what variable to write the output to. You need to have an explicit call to thermostat in your code.
Nikola
Nikola 2013-9-9
编辑：Nikola 2013-9-9
Aha ok. Thus, how I suppose to write a procedure which execute that? What should be change in existed procedure that I wrote here?
Walter Roberson
Walter Roberson 2013-9-9
The xprime you are returning does not look to involve values calculated by thermostat(). If it actually does, then you would have the calls in the way you defined your system of equations... e.g.,
xprime = [ .... t^5+thermostat(x), ....]
If, though, the differential equation matrix is not calculated in terms of thermostat calls, then some other layer needs to receive the result of the thermostat calls. What layer is that?
Nikola
Nikola 2013-9-9
I think that I also thought about the same thing today but I couldn't figure a way how to solve the system of equations separately and then relate to the thermostat function.
Thermostat input should be xprime output, correct? Is that answer to your question?

请先登录，再进行评论。

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

 采纳的回答

Walter Roberson
Walter Roberson 2013-9-9

0 个投票

I am going to guess that you have missed out on telling us that your heat_sys function is being invoked by ode45 or similar. I am going to further guess that heat_sys is your main function, and that your thermostat function is needed to calculate non-linear constraints. If I am correct, then a handle to your thermostat function would need to be passed as an additional parameter to the ode call, and it would also have to be modified to return information for both equality constraints and inequality constraints (I cannot tell which you intend in your existing code.)
If I am correct so far, please have another look at the documentation for the ode function you are using.

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

Nikola
Nikola 2013-9-9
Right, I used ode45 to solve heat_sys function and it is main function because I want to obtain x1 and x2 values. Thermostat is necessary to check if x values are within given bounds and should be switch on or off.
To handle thermostat function into ode, is this correct: [t,x] = ode45(@heat_sys,@thermostat,x0,tspan)
Since I am not an experienced Matlab user, I need a help to make corrections in procedure that I wrote above.
Walter Roberson
Walter Roberson 2013-9-9
Okay so how do the thermostat values actually get used? What takes them as inputs ? Or are they purely constraints (which it doesn't quite sound to me that they are) ?
Nikola
Nikola 2013-9-10
Output thermostat values are on and off. It means if temperature is below 18, then is on. Also, I want to keep temperature between 18 and 20 but if goes above 20, then thermostat should be switched off. Therefore, inputs of thermostat function are output values of xprime because in every another iteration, we have new value of xprime, or x1 and x2.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by