Greybox nonlniear estimation model

Question

Imisi Jo 2017-7-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/348212-greybox-nonlniear-estimation-model

评论： Imisi Jo 2017-7-20

Screen Shot 2017-07-10 at 01.44.40.png

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.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Prashant Arora 2017-7-17

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/348212-greybox-nonlniear-estimation-model#answer_274516

在 MATLAB Online 中打开

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);