Synthesis of a smoothing bicubic spline for differentiating experimental data of nonparametric identification

Over the past two decades, so-called Voltaire series have been used to describe the dynamics of nonlinear systems in terms of input-output. Nonparametric identification of models using Voltaire series consists in constructing estimates for impulse transition functions (IPFs) depending on two or more arguments, which naturally makes identification algorithms much more complicated than in the one-dimensional case. So, in order to identify the two-dimensional IPF (corresponding to the quadratic term of the Voltaire series), it is necessary to calculate the second-order mixed derivatives of the output two-dimensional signal of the system, when a series of rectangular pulses of different amplitudes at different times are fed to its input. Everyone knows, the problem of differentiation is an ill-posed problem and one of the manifestations of incorrectness is poor resistance to errors in the initial data. It is proposed to use two-dimensional smoothing cubic (bicubic) spline (abbreviated SBS) to overcome this problem. The two tasks that constitute SBS synthesis: assignment and implementation of different types of boundary conditions at the border of the rectangular region where SBS is determined; optimal values estimation of two smoothing parameters due to the different “smoothness” of IPF for different two arguments. An acceptable solution to this synthesis problem is proposed in the paper. Our performed computational experiment showed the efficiency of the proposed algorithm for calculating second-order mixed derivative from noisy initial data.