Greybox nonlniear estimation model

1 次查看(过去 30 天)
I am trying to use the nonlinear greybox model to estimate the parameter R0 in my state space model (attached). However, I have an extra variable Voc which is a column vector and varies depending on the time. How can I write this in my function under modelfile name (ie nlm = idnlgrey('batteryeq',order,parameters,initial_states,Ts)? I get an error saying "Subscript indices must either be real positive integers or logicals." with my current code:
function [dx,y] = batteryeq(t,x,u,T1,T2,C1,C2,R0,OCV,varargin)
% Output equation.
y = [-x(1);-x(2)]-R0*u+OCV(t); % output equation
% State equations.
dx = [-1/T1+(1/C1)*u;-1/T2+(1/C2)*u];
end
Thank you.

采纳的回答

Prashant Arora
Prashant Arora 2017-7-17
Hi Imisi,
If you already know the dependence of OCV on time, you can definitely pass it as a parameter. However, as time is a floating value, you'll need to "index" into OCV with the corresponding index of time, rather than the actual decimal value of time. (OCV(t) in your code)
For example, let's say OCV = 100X1 where OCV defines the Open-circuit voltage for time = 1 to 100 seconds. Given that you know the value of current time in the function as variable, t, you can compute the OCV at the time using linear interpolation:
interp1(1:100,OCV,t);
  3 个评论
Prashant Arora
Prashant Arora 2017-7-19
This is probably happening because the variable "t" spans outside the domain of "time_OCV", which results in interp1 to produce NaNs.
You can specify interp1 to use extrapolation outside the domain using
interp1(time_OCV,OCV,t,'linear','extrap')
Imisi Jo
Imisi Jo 2017-7-20
Ok I will try that, thank you! I assume there is no way for me to limit the time duration of 't'?

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Battery Pack Modeling 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by