Nested functions procedure error
2 次查看(过去 30 天)
显示 更早的评论
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 个评论
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?
采纳的回答
Walter Roberson
2013-9-9
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 个评论
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) ?
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!