The controllable canonical form for multiple input system is not unique. That is, there may be multiple possible controllable canonical form possible for a multiple input system. This statement is explained in a very nice manner, using the concept of controllability matrix and controllability index in the IEEE paper:
D. Luenberger, "Canonical forms for linear multivariable systems," in IEEE Transactions on Automatic Control, vol. 12, no. 3, pp. 290-293, June 1967, doi: 10.1109/TAC.1967.1098584. keywords: {MIMO;Matrices;Control systems;Controllability;Control design;Design engineering;Differential equations;Contracts},
The above sited paper discribes the possible ways to convert a multiple input system into the controllable canonical form.
The file exchange link attached by you in the comments section, works on the similiar principle as discussed in the above paper.
Comming to your question:
" I have found a code file but I can't find the logic behind of this and how It can gives controllable canonical form for multi-inputs system and I don't know If this works correctly or not ? "
After reviewing the code, I believe it effectively transforms the multiple input system into a controllable canonical form, in accordance with one of the methodologies outlined in the aforementioned IEEE paper.
I am summarising the things I collected from the paper, which will help you to understand the algorithm used in the file-exchange code.
consider the state-space representation of the system:
Assume a transformation matrix "S" that will take the system in its controllable canonical form as demonstrated below:
For a single input system, this transformation matrix is nothing but the controllability matrix "C" as shown below:
Suppose that the number of states are n and it is only a single input system. Under that case, the controllability matrix is of order (n x n). Therefore, if the rank of "C" or "S" matrix is equal to "n" ,it is invertible and we can convert the system into the controllable canonical form. If the rank of "C" matrix is less than "n", it is not invertible and we say that the system cannot be converted to its controllable canonical form.
Now comming to the multiple input case. Suppose that we have "r" number of inputs. Therefore, the controllability matrix "C" will be (n x nr). As it is a rectangular matrix, it is not invertible. Therefore, in this case we cannot take :
We have to make a choice of matrix "S" that will take the system to the controllable canonical form. This choice is not unique. There can be many possible ways to choose "S" matrix. Therefore, the above cited papaer describes some approaches to find the matrix "S" that transforms the sytem into its controllable canonical form.
The approach used in the file exchange code takes:
here, row vectors 
are noting but the linearly independent row vectors appearing in the inverse of "P" matrix. (Refer to the above cited research paper for the defination of "P" matrix and to get clearity regarding the terms appearing in the above equation.)
are noting but the linearly independent row vectors appearing in the inverse of "P" matrix. (Refer to the above cited research paper for the defination of "P" matrix and to get clearity regarding the terms appearing in the above equation.)According to it , the transformed system has "A" and "B" matrixes of the form:
I tried to run the file-exchange code on your system of equations. I am running it here, in order to demonstrate the results of the transformation.
% Original state-space Plant
A = [-0.978707531551846 2.37478961356718e-14 7.42345117560053e-14 0.318698761698607;
-2.63927929244335e-14 -2.00000000000000 -3 -4.93782914225613e-16;
8.09939070469741e-19 -3.33066907387547e-16 -1.00000000000000 5.83814072314816e-19;
-0.135043844787908 -3.13829005324231e-15 -2.83836564755302e-14 -3.02129246844816];
B = [-0.993655438699628 -0.136196562445474;
-2.31027060974740e-14 -1.78885438199983;
1.75613973071755e-14 -0.894427190999916;
-0.360965272521026 0.863803437554497];
[T,Ad,Bd,ni]=MI_ctrb(A,B) %Using the file-exchange method
You can observe that the structure of the transformed matrices "Ad" and "Bd" matches exactly as per the above attached images.
In conclusion, MATLAB does not natively support the transformation of multiple input systems into a controllable canonical form. To achieve this, one should consult research papers that discuss various approaches. If you come across third-party MATLAB code or File Exchange contributions, it is advisable to first verify them against the original resources and methodologies to ensure accuracy and reliability.
I hope you find the information useful !!!
