
Differential Equations Ordinary and Partial Differential Equations Math Forum 
 LinkBack  Thread Tools  Display Modes 
March 9th, 2018, 11:23 AM  #1 
Newbie Joined: Mar 2018 From: London Posts: 1 Thanks: 0  Identification of unknown function in a DAE system
I have a system of 1st order odes given by $$ \dot{x_1}(t) = \alpha_1 f_1(x_1,t) + \beta_1 u(t) \\ \dot{x_2}(t) = \alpha_2 f_2(x_2,t) + \beta_2 u(t) $$ They are constrained by an algebraic equation $$ x_1(t) + x_2(t) = k $$ where $\left( \alpha_1,\alpha_2, \beta_1,\beta_2 , k \right) \in \mathbb{R}$ are known constants (i.e. parameters). $f_1(t)$ and $f_2(t)$ are both unknown. Starting from a rich set of inputoutput **noisefree** data available from simulating a complex proxy system, what would be the best procedure to identify (even a subset of repeatable/characteristic properties) the unknown **_possibly timevarying_** functions $f_1(x_1,t)$ and $f_2(x_2,t)?$ I am almost certain that $f_1(x_1,t)$ and $f_2(x_2,t)$ are both linear. I am looking for a greybox systemid approach that shall work well to arbitrary excitations in all future simulations. (NOT merely a curvefitting procedure to match a specific excitation inputoutput dataset) 
March 9th, 2018, 09:47 PM  #2 
Senior Member Joined: Sep 2016 From: USA Posts: 531 Thanks: 304 Math Focus: Dynamical systems, analytic function theory, numerics 
Interesting question. I have a few observations which might help get you started. Initially, lets focus on the question of determining if the equation is linear. If this is the case, then note that it is a completely uncoupled linear equation. Each of your functions has the form $f_i(x_i,t) = c_i(t)x_i(t)$ and therefore, the ODE is given by writing \[\dot x = A(t)x(t) + \mu(t) \] where $A$ is amatrix and $\mu$ is the vectorized inhomogenous term. Most importantly, the observation that it is uncoupled means that $A$ is a diagonal matrix. This leads to a fundamental matrix solution defined by the usual matrix exponential: \[x \mapsto \exp(tA(0))\cdot x \] and applying the variation of constants formula gives the fundamental solution operator (flow) \[\Phi(x,t) = \exp(tA(0)) \cdot x + \int_0^t \exp((ts)A(s)) \mu(s) \ ds \] Now, you claim you have the ability to query this flow for a bunch of initial data. In other words, you can pick a few values, $x_1,x_2$ and compute the trajectory $x_1(t) = \Phi(x_1,t), x_2(t) = \Phi(x_2,t)$ and from this data and the fact that $A$ is diagonal, hence easily invertible, you can obtain explicit values for $c_1(t),c_2(t)$. Now, compute the flow for a "bunch" of other initial data by querying the black box, and also compute those trajectories in the linear ODE defined by choosing $c_1,c_2$ to satisfy the first two trajectories, and check if these trajectories match. If they don't, then $f_1,f_2$ are not linear, and if they do, then they are linear, and moreover, you know explicitly what they are. 2. The conserved quantity $x_1 + x_2 = k$ means that a particular linear combination of solutions for the ODE will conserve this quantity. Hence, if the equations are linear, then this solution must correspond to an eigenvalue of 0 for $A$ so you only have a 1dimensional subspace of nontrivial solutions. Thus, you actually expect to be able to uniquely determine $c_1,c_2$ using only a single initial condition. 3. In the event it turns out to be nonlinear, then how much you can learn depends whether you are allow to query $f_1,f_2$ directly, rather than just the solution operator. If so, you can still uncover a great deal of information about the phase space structure. I can provide more details if this is relevant. I'll think about this question some more since it is interesting. I hope my comments so far are helpful. 
March 16th, 2018, 12:56 PM  #3 
Senior Member Joined: Dec 2015 From: Earth Posts: 316 Thanks: 42 
Try for $\displaystyle x_1 ' (t)+ x_2 ' (t)=k'=0$ Now plug to system expression above 

Tags 
control, dae, dynamicaa, function, identification, system, unknown 
Thread Tools  
Display Modes  

Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Prime Identification Game  jsmith  Elementary Math  0  November 12th, 2014 08:59 AM 
Equation system with 2 unknown, which contains quadratic eq.  miran97  Algebra  5  February 16th, 2013 06:11 AM 
Unknown Function  aaronmath  Calculus  8  October 12th, 2011 11:03 PM 
Set up linear system to solve for unknown  Len  Linear Algebra  0  November 18th, 2009 11:01 AM 
Identification protocol  trid2  Number Theory  0  April 20th, 2009 09:29 AM 