[{"@context":"http:\/\/schema.org\/","@type":"BlogPosting","@id":"https:\/\/wiki.edu.vn\/en\/wiki24\/tp-model-transformation-in-control-theory\/#BlogPosting","mainEntityOfPage":"https:\/\/wiki.edu.vn\/en\/wiki24\/tp-model-transformation-in-control-theory\/","headline":"TP model transformation in control theory","name":"TP model transformation in control theory","description":"Baranyi and Yam proposed the TP model transformation[1][2][3][4][5][6][7] as a new concept in quasi-LPV (qLPV) based control, which plays a","datePublished":"2019-12-28","dateModified":"2019-12-28","author":{"@type":"Person","@id":"https:\/\/wiki.edu.vn\/en\/wiki24\/author\/lordneo\/#Person","name":"lordneo","url":"https:\/\/wiki.edu.vn\/en\/wiki24\/author\/lordneo\/","image":{"@type":"ImageObject","@id":"https:\/\/secure.gravatar.com\/avatar\/c9645c498c9701c88b89b8537773dd7c?s=96&d=mm&r=g","url":"https:\/\/secure.gravatar.com\/avatar\/c9645c498c9701c88b89b8537773dd7c?s=96&d=mm&r=g","height":96,"width":96}},"publisher":{"@type":"Organization","name":"Enzyklop\u00e4die","logo":{"@type":"ImageObject","@id":"https:\/\/wiki.edu.vn\/wiki4\/wp-content\/uploads\/2023\/08\/download.jpg","url":"https:\/\/wiki.edu.vn\/wiki4\/wp-content\/uploads\/2023\/08\/download.jpg","width":600,"height":60}},"image":{"@type":"ImageObject","@id":"https:\/\/wikimedia.org\/api\/rest_v1\/media\/math\/render\/svg\/51d289eb3d5ec0cdc26b7e09406d558f266a2b20","url":"https:\/\/wikimedia.org\/api\/rest_v1\/media\/math\/render\/svg\/51d289eb3d5ec0cdc26b7e09406d558f266a2b20","height":"","width":""},"url":"https:\/\/wiki.edu.vn\/en\/wiki24\/tp-model-transformation-in-control-theory\/","about":["Wiki"],"wordCount":15246,"articleBody":"Baranyi and Yam proposed the TP model transformation[1][2][3][4][5][6][7] as a new concept in quasi-LPV (qLPV) based control, which plays a central role in the highly desirable bridging between identification and polytopic systems theories. It is also used as a TS (Takagi-Sugeno) fuzzy model transformation. It is uniquely effective in manipulating the convex hull of polytopic forms (or TS fuzzy models), and, hence, has revealed and proved the fact that convex hull manipulation is a necessary and crucial step in achieving optimal solutions and decreasing conservativeness[8][9][2] in modern linear matrix inequality based control theory. Thus, although it is a transformation in a mathematical sense, it has established a conceptually new direction in control theory and has laid the ground for further new approaches towards optimality.For details please visit: TP model transformation.TP-tool MATLAB toolboxA free MATLAB implementation of the TP model transformation can be downloaded at [1] or an old version of the toolbox is available at  MATLAB Central [2]. Be careful, in the MATLAB toolbox the assignments of the dimensions of the core tensor is in the opposite way in contrast to the notation used in the related literature. In some variants of the ToolBox, the first two dimension of the core tensor is assigned to the vertex systems. In the TP model literature the last two. A simple example is given below.clearM1=20; \t\t% Grid densityM2=20;omega1=[-1,1]; \t\t%Interval omega2=[-1,1];domain=[omega1; omega2];for m1=1:M1    for m2=1:M2        p1=omega1(1)+(omega1(2)-omega1(1))\/M1*(m1-1); \t%sampling grid        p2=omega2(1)+(omega2(2)-omega2(1))\/M2*(m2-1);        SD(m1,m2,1,:)=[1 0];\t\t\t\t% SD is the discretized system matrix        SD(m1,m2,2,:)=[(-1-0.67*p1*p1) (1.726*p2*p2)];    end    end[S,U, sv]=hosvd(SD,[1,1,0,0],1e-12);\t\t% Finding the TP structureUA{1}=U{1};\t\t\t\t\t% This is the HOSVD based canonical formUA{2}=U{2};ns1 = input('Results of SNNN TS fuzzy model');UC=genhull(UA,'snnn'); \t\t\t\t% snnn weightinf functionsUCP{1}=pinv(UC{1});UCP{2}=pinv(UC{2});SC=tprods(SD,UCP);  \t\t\t\t%This is to find  the core tensor H(:,:)=SC(1,1,:,:)                              %This is to show the vertices of the TP modelH(:,:)=SC(1,2,:,:)H(:,:)=SC(2,1,:,:)H(:,:)=SC(2,2,:,:)figure(1) hold all plothull(U{1}, omega1)                          %Draw the waiting functions of p1title('Weighting functions for p_{1}');xlabel('p_{1}')ylabel('Weighting functions')grid onbox onfigure(2) hold all plothull(UC{2}, omega2)                         %Show the waiting functions of p2title('Weighting functions for p_{2}');xlabel('p_{2}')ylabel('Weighting functions')grid onbox onns2 = input('Results of CNO TS fuzzy model');UC=genhull(UA,'cno');                          %Create CNO type waiting functions UCP{1}=pinv(UC{1});UCP{2}=pinv(UC{2});SC=tprods(SD,UCP);                             %Find the cortensorH(:,:)=SC(1,1,:,:)                             %Show the vertices of the TP modelH(:,:)=SC(1,2,:,:)H(:,:)=SC(2,1,:,:)H(:,:)=SC(2,2,:,:)figure(1) hold all plothull(U{1}, omega1)                         %Show the waiting functions of p1title('Weighting functions for p_{1}');xlabel('p_{1}')ylabel('Weighting functions')grid onbox onfigure(2) hold all plothull(UC{2}, omega2)                        %Show the waiting functions of p2title('Weighting functions for p_{2}');xlabel('p_{2}')ylabel('Weighting functions')\tOnce you have the feedback vertexes derived to each vertexes of the TP model then you may want to calculate the controller over the same polytope (see PDC design by Tanaka)W = queryw1(UC,domain,p); \u00a0% computing the weighting values over the parameter vectorF = tprods(K,W);\t\t% calculating the parameter dependent feedback F(p)F = shiftdim(F)U=-F*x\t\t\t\t% calculate the control value.Key features for control analysis and design[edit]The TP model transformation transforms a given qLPV model into a (tensor product type) polytopic form, irrespective of whether the model is given in the form of analytical equations resulting from physical considerations, or as an outcome of soft computing based identification techniques (such as neural networks or fuzzy logic based methods, or as a result of a black-box identification).Further the TP model transformation is capable of manipulating the convex hull defined by the polytopic form that is a necessary step in polytopic qLPV model-based control analysis and design theories.Related definitions[edit]Linear Parameter-Varying (LPV) state-space model(x\u02d9(t)y(t))=S(p(t))(x(t)u(t)),{displaystyle {begin{pmatrix}{mathbf {dot {x}} }(t)\\{mathbf {y} }(t)end{pmatrix}}={mathbf {S} }({mathbf {p} }(t)){begin{pmatrix}{mathbf {x} }(t)\\{mathbf {u} }(t)end{pmatrix}},}with input u(t){displaystyle {mathbf {u} }(t)}, output y(t){displaystyle {mathbf {y} }(t)} and statevector x(t){displaystyle {mathbf {x} }(t)}. The system matrix S(p(t))\u2208RL1\u00d7L2{displaystyle {mathbf {S} }({mathbf {p} }(t))in mathbb {R} ^{L_{1}times L_{2}}} is a parameter-varying object, where p(t)\u2208\u03a9{displaystyle {mathbf {p} }(t)in Omega } is a time varying N{displaystyle N}-dimensional parameter vector which is an element ofclosed hypercube \u03a9=[a1,b1]\u00d7[a2,b2]\u00d7\u22ef\u00d7[aN,bN]\u2282RN{displaystyle Omega =[a_{1},b_{1}]times [a_{2},b_{2}]times cdots times [a_{N},b_{N}]subset mathbb {R} ^{N}}. As a matter of fact, further parameter dependent channels can be inserted to  S(p(t)){displaystyle {mathbf {S} }({mathbf {p} }(t))} that represent various control performance requirements.quasi Linear Parameter-Varying (qLPV) state-space modelp(t){displaystyle {mathbf {p} }(t)} in the above LPV model can also include some elements of the state vectorx(t){displaystyle {mathbf {x} }(t)}, and, hence this model belongs to the class of non-linear systems, and is also referred to as a quasi LPV (qLPV) model.TP type polytopic Linear Parameter-Varying (LPV) state-space model(x\u02d9(t)y(t))=S\u22a0n=1Nwn(pn(t))(x(t)u(t)),{displaystyle {begin{pmatrix}{mathbf {dot {x}} }(t)\\{mathbf {y} }(t)end{pmatrix}}={mathcal {S}}boxtimes _{n=1}^{N}mathbf {w} _{n}(p_{n}(t)){begin{pmatrix}{mathbf {x} }(t)\\{mathbf {u} }(t)end{pmatrix}},}with input u(t){displaystyle {mathbf {u} }(t)}, output y(t){displaystyle {mathbf {y} }(t)} and statevector x(t){displaystyle {mathbf {x} }(t)}. The system matrix S(p(t))=S\u22a0n=1Nw(pn(t))\u2208RL1\u00d7L2{displaystyle {mathbf {S} }({mathbf {p} }(t))={mathcal {S}}boxtimes _{n=1}^{N}mathbf {w} (p_{n}(t))in mathbb {R} ^{L_{1}times L_{2}}} is a parameter-varying object, where p(t)\u2208\u03a9{displaystyle {mathbf {p} }(t)in Omega } is a time varying N{displaystyle N}-dimensional parameter vector which is an element ofclosed hypercube \u03a9=[a1,b1]\u00d7[a2,b2]\u00d7\u22ef\u00d7[aN,bN]\u2282RN{displaystyle Omega =[a_{1},b_{1}]times [a_{2},b_{2}]times cdots times [a_{N},b_{N}]subset mathbb {R} ^{N}}, and the weighting functions wn,in(pn(t))\u2208[0,1]{displaystyle w_{n,i_{n}}(p_{n}(t))in [0,1]} are the elements of vector wn(pn(t)){displaystyle mathbf {w} _{n}(p_{n}(t))}. Core tensor contains elements Si1,i2,\u2026,iN{displaystyle mathbf {S} _{i_{1},i_{2},ldots ,i_{N}}} which are the vertexes of the system.As a matter of fact, further parameter dependent channels can be inserted to  S(p(t)){displaystyle {mathbf {S} }({mathbf {p} }(t))} that represent various control performance requirements.Here\u2200n:\u2211in=1Inwn,in(pn(t))=1.{displaystyle forall n:sum _{i_{n}=1}^{I_{n}}w_{n,i_{n}}(p_{n}(t))=1.} and wn,in(pn(t))\u2208[0,1].{displaystyle w_{n,i_{n}}(p_{n}(t))in [0,1].}This means that S(p(t)){displaystyle mathbf {S} ({mathbf {p} }(t))} is within the vertexes Si1,i2,\u2026,iN{displaystyle mathbf {S} _{i_{1},i_{2},ldots ,i_{N}}} of the system (within the convex hull defined by the vertexes) for all p(t)\u2208\u03a9{displaystyle mathbf {p} (t)in Omega }.Note that the TP type polytopic model can always be given in the formS(p(t))=\u2211r=1RSrwr(p(t)),{displaystyle mathbf {S} (mathbf {p} (t))=sum _{r=1}^{R}mathbf {S} _{r}w_{r}(mathbf {p} (t)),}where the vertexes are the same as in the TP type polytopic form and the multi variable weighting functions are the product of the one variable weighting functions according to the TP type polytopic form, and r is the linear index equivalent of the multi-linear indexing i1,i2,\u2026iN{displaystyle i_{1},i_{2},ldots i_{N}}.TP model transformation for qLPV modelsAssume a given qLPV model S(p(t)){displaystyle mathbf {S} (mathbf {p} (t))}, where p(t)\u2208\u03a9\u2282RN{displaystyle mathbf {p} (t)in Omega subset R^{N}}, whose TP polytopic structure may be unknown (e.g. it is given by neural networks). The TP model transformation determines its TP polytopic structure asS(p(t))=S\u22a0n=1Nwn(pn(t)){displaystyle mathbf {S} (mathbf {p} (t))={mathcal {S}}boxtimes _{n=1}^{N}mathbf {w} _{n}(p_{n}(t))},namely it generates core tensor S{displaystyle {mathcal {S}}} and weighting functions of wn(pn(t)){displaystyle mathbf {w} _{n}(p_{n}(t))} for all n=1\u2026N{displaystyle n=1ldots N}. Its free MATLAB implementation is downloadable at [3] or at  MATLAB Central [4].If the given model does not have (finite element) TP polytopic structure, then the TP model transformation determines its approximation:S(p(t))\u2248S\u22a0n=1Nwn(pn(t)),{displaystyle mathbf {S} (mathbf {p} (t))approx {mathcal {S}}boxtimes _{n=1}^{N}mathbf {w} _{n}(p_{n}(t)),}where trade-off is offered by the TP model transformation between complexity (number of vertexes stored in the core tensor or the number of weighting functions) and the approximation accuracy.[10] The TP model can be generated according to various constrains. Typical TP models generated by the TP model transformation are:HOSVD canonical form of qLPV models,Various kinds of TP type polytopic form (this feature is very important in control performance optimization).TP model based control design[edit]Key methodologySince the TP type polytopic model is a subset of the polytopic model representations, the analysis and design methodologies developed for polytopic representations are applicable for the TP type polytopic models as well.One typical way is to search the nonlinear controller in the form:u=\u2212F(p(t))x(t)=\u2212\u2211r=1RFrwr(p(t))x(t),{displaystyle u=-mathbf {F} (mathbf {p} (t))mathbf {x} (t)=-sum _{r=1}^{R}F_{r}w_{r}(mathbf {p} (t))mathbf {x} (t),}where the vertexes Fr{displaystyle mathbf {F} _{r}} of the controller is calculated from Sr{displaystyle mathbf {S} _{r}}. Typically, the vertexes Sr{displaystyle mathbf {S} _{r}} are substituted into Linear Matrix Inequalities in order to determine Fr{displaystyle mathbf {F} _{r}}.In TP type polytopic form the controller is:u=\u2212F(p(t))x(t)=\u2212F\u22a0n=1Nwn(pn(t))x(t),{displaystyle u=-mathbf {F} (mathbf {p} (t))mathbf {x} (t)=-{mathcal {F}}boxtimes _{n=1}^{N}mathbf {w} _{n}(p_{n}(t))mathbf {x} (t),}where the vertexes Fi1,i2,\u2026,iN{displaystyle mathbf {F} _{i_{1},i_{2},ldots ,i_{N}}} stored in the core tensor F{displaystyle {mathcal {F}}} are determined from the vertexes Si1,i2,\u2026,iN{displaystyle mathbf {S} _{i_{1},i_{2},ldots ,i_{N}}} stored in S{displaystyle {mathcal {S}}}. Note that the polytopic observer or other components can be generated in similar way, such as these vertexes are also generated from S{displaystyle {mathcal {S}}}.Convex hull manipulation based optimizationThe polytopic representation of a given qLPV model is not invariant. I.e. a given S(p(t)){displaystyle mathbf {S} (mathbf {p} (t))} has Z=\u221e{displaystyle Z=infty } number of different representation as:S(p(t))=Sz\u22a0n=1Nwz,n(pn(t)),{displaystyle mathbf {S} (mathbf {p} (t))={mathcal {S}}_{z}boxtimes _{n=1}^{N}mathbf {w} _{z,n}(p_{n}(t)),}where z=1\u2026Z{displaystyle z=1ldots Z}. In order to generate an optimal control of the given model S(p(t)){displaystyle mathbf {S} (mathbf {p} (t))} we apply, for instance LMIs. Thus, if we apply the selected LMIs to the above polytopic model we arrive at:uz=\u2212Fz(p(t))x(t)=\u2212Fz\u22a0n=1Nwz,n(pn(t))x(t).{displaystyle u_{z}=-mathbf {F} _{z}(mathbf {p} (t))mathbf {x} (t)=-{mathcal {F}}_{z}boxtimes _{n=1}^{N}mathbf {w} _{z,n}(p_{n}(t))mathbf {x} (t).}Since the LMIs realize a non-linear mapping between the vertexes in S{displaystyle {mathcal {S}}} and F{displaystyle {mathcal {F}}} we may find very different controllers for each z=1\u2026Z{displaystyle z=1ldots Z}. This means that we have Z{displaystyle Z} different number of “optimal” controllers to the same system S(p(t)){displaystyle mathbf {S} (mathbf {p} (t))}. Thus, the question is: which one of the “optimal” controllers is really the optimal one. The TP model transformation let us to manipulate the weighting functions systematically that is equivalent to the manipulation of the vertexes. The geometrical meaning of this manipulation is the manipulation of the convex hull defined by the vertexes. We can easily demonstrate the following facts:Tightening the convex hull typically decreases the conservativeness of the solution, so as may lead to better control performance. For instance, if we have a polytopic representationS(p(t))=S\u22a0n=1Nwn(pn(t)),{displaystyle mathbf {S} (mathbf {p} (t))={mathcal {S}}boxtimes _{n=1}^{N}mathbf {w} _{n}(p_{n}(t)),}of a given model S(p(t)){displaystyle mathbf {S} (mathbf {p} (t))}, then we can generate a controller asu=\u2212F\u22a0n=1Nwn(pn(t))x(t),{displaystyle u=-{mathcal {F}}boxtimes _{n=1}^{N}mathbf {w} _{n}(p_{n}(t))mathbf {x} (t),}then we solved the control problem of all systems k=1\u2026K=\u221e{displaystyle k=1ldots K=infty } that can be given by the same vertexes, but with different weighting functions as:Sk(p(t))=S\u22a0n=1Nwk,n(pn(t)),{displaystyle mathbf {S} _{k}(mathbf {p} (t))={mathcal {S}}boxtimes _{n=1}^{N}mathbf {w} _{k,n}(p_{n}(t)),}whereuk=\u2212F\u22a0n=1Nwk,n(pn(t))x(t).{displaystyle u_{k}=-{mathcal {F}}boxtimes _{n=1}^{N}mathbf {w} _{k,n}(p_{n}(t))mathbf {x} (t).}If one of these systems are very hardly controllable (or even uncontrollable) then we arrive at a very conservative solution (or unfeasible LMIs). Therefore, we expect that during tightening the convex hull we exclude such problematic systems.It can also be easily demonstrated that the observer design is typically needs large convex hull. So, as when we design controller and observer we need to find the optimal convex hull between the tight one and the large one. Same papers also demonstrate that using different convex hulls (if the separation principal is applicable) for observer and controller may lead to even better solution.Properties of the TP model transformation in qLPV theories[edit]It can be executed uniformly (irrespective of whether the model is given in the form of analytical equations) resulting from physical considerations, or as an outcome of soft computing based identification techniques (such as neural networks or fuzzy logic based methods, or as a result of a black-box identification), without analytical interaction, within a reasonable amount of time. Thus, the transformation replaces the analytical and in many cases complex and not obvious conversions to numerical, tractable, straightforward operations that can be carried out in a routine fashion.It generates the HOSVD-based canonical form of qLPV models, which is a unique representation. This form extracts the unique structure of a given qLPV model in the same sense as the HOSVD does for tensors and matrices, in a way such that:the number of LTI components are minimized;the weighting functions are one variable functions of the parameter vector in an orthonormed system for each parameter (singular functions);the LTI components (vertex components) are also in orthogonal positions;the LTI systems and the weighting functions are ordered according to the higher-order singular values of the parameter vector;it has a unique form (except for some special cases);introduces and defines the rank of the qLPV model by the dimensions of the parameter vector;The core step of the TP model transformation was extended to generate different types of convex polytopic models, in order to focus on the systematic (numerical and automatic) modification of the convex hull instead of developing new LMI equations for feasible controller design (this is the widely adopted approach). It is worth noting that both the TP model transformation and the LMI-based control design methods are numerically executable one after the other, and this makes the resolution of a wide class of problems possible in a straightforward and tractable, numerical way.Based on the higher-order singular values (which express the rank properties of the given qLPV model, see above, for each element of the parameter vector in L2{displaystyle L_{2}} norm), the TP model transformation offers a trade-off between the complexity of the TP model (polytopic form),[10] hence, the LMI design and the accuracy of the resulting TP model.The TP model transformation is executed before utilizing the LMI design. This means that when we start the LMI design we already have the global weighting functions and during control we do not need to determine a local weighting of the LTI systems for feedback gains to compute the control value at every point of the hyperspace the system should go through. Having predefined continuous weighting functions also ensures that there is no friction in the weighting during control.References[edit]^ Baranyi, P. (2004). “TP Model Transformation as a Way to LMI-Based Controller Design”. IEEE Transactions on Industrial Electronics. 51 (2): 387\u2013400. doi:10.1109\/TIE.2003.822037. S2CID\u00a07957799.^ a b Baranyi, P\u00e9ter (2016). TP-Model Transformation-Based-Control Design Frameworks. doi:10.1007\/978-3-319-19605-3. ISBN\u00a0978-3-319-19604-6.^ Baranyi, P\u00e9ter; Tikk, Domonkos; Yam, Yeung; Patton, Ron J. (2003). “From differential equations to PDC controller design via numerical transformation”. Computers in Industry. 51 (3): 281\u2013297. doi:10.1016\/S0166-3615(03)00058-7.^ Baranyi, Peter (2014). “The Generalized TP Model Transformation for T\u2013S Fuzzy Model Manipulation and Generalized Stability Verification”. IEEE Transactions on Fuzzy Systems. 22 (4): 934\u2013948. doi:10.1109\/TFUZZ.2013.2278982.^ P. Baranyi; Y. Yam; P. V\u00e1rlaki (2013). Tensor Product model transformation in polytopic model-based control. Boca Raton FL: Taylor & Francis. p.\u00a0240. ISBN\u00a0978-1-43-981816-9.^ Baranyi, Peter (February 2022). “How to Vary the Input Space of a T\u2013S Fuzzy Model: A TP Model Transformation-Based Approach”. IEEE Transactions on Fuzzy Systems. 30 (2): 345\u2013356. doi:10.1109\/TFUZZ.2020.3038488. ISSN\u00a01063-6706.^ Baranyi, Peter (March 2020). “Extracting LPV and qLPV Structures From State-Space Functions: A TP Model Transformation Based Framework”. IEEE Transactions on Fuzzy Systems. 28 (3): 499\u2013509. doi:10.1109\/TFUZZ.2019.2908770. ISSN\u00a01063-6706.^ Szollosi, Alexandra; Baranyi, Peter (2016). “Influence of the Tensor Product Model Representation of QLPV Models on the Feasibility of Linear Matrix Inequality”. Asian Journal of Control. 18 (4): 1328\u20131342. doi:10.1002\/asjc.1238.^ Szollosi, Alexandra; Baranyi, Peter (2017). “Improved control performance of the 3-DoF aeroelastic wing section: A TP model based 2D parametric control performance optimization”. Asian Journal of Control. 19 (2): 450\u2013466. doi:10.1002\/asjc.1418. S2CID\u00a0126195046.^ a b D. Tikk, P.Baranyi, R. J. Patton (2007). “Approximation Properties of TP Model Forms and its Consequences to TPDC Design Framework”. Asian Journal of Control. 9 (3): 221\u2013331. doi:10.1111\/j.1934-6093.2007.tb00410.x. S2CID\u00a0121716136.{{cite journal}}:  CS1 maint: multiple names: authors list (link)"},{"@context":"http:\/\/schema.org\/","@type":"BreadcrumbList","itemListElement":[{"@type":"ListItem","position":1,"item":{"@id":"https:\/\/wiki.edu.vn\/en\/wiki24\/#breadcrumbitem","name":"Enzyklop\u00e4die"}},{"@type":"ListItem","position":2,"item":{"@id":"https:\/\/wiki.edu.vn\/en\/wiki24\/tp-model-transformation-in-control-theory\/#breadcrumbitem","name":"TP model transformation in control theory"}}]}]