matlab求能观标准型

求能观标准型的步骤如下：
1. 确定系统的状态空间模型，即A、B、C、D矩阵。
2. 判断系统是否可控，若不可控则无法求解能观标准型。
3. 计算系统的能控性矩阵Mc，其中Mc的秩等于系统状态的个数n。
4. 构造变换矩阵T，使得TMcT^(-1)为对角矩阵。
5. 将状态空间模型通过变换矩阵T进行相似变换，得到新的状态空间模型，即为能观标准型。

以下是matlab代码实现：

% 确定系统的状态空间模型
A = [1 2 3; 4 5 6; 7 8 9];
B = [1; 1; 1];
C = [1 0 0];
D = 0;
sys = ss(A, B, C, D);

% 判断系统是否可控
if rank(ctrb(sys)) == size(A, 1)
    disp('The system is controllable.');
else
    disp('The system is uncontrollable.');
    return;
end

% 计算系统的能控性矩阵Mc
Mc = ctrb(A, B);

% 构造变换矩阵T
T = Mc;
for i = 1:size(A, 1)-1
    T = [T A^(i)*B];
end

% 将状态空间模型通过变换矩阵T进行相似变换，得到新的状态空间模型
sys_obs = ss2ss(sys, T);

% 输出能观标准型的状态空间模型
[A_obs, B_obs, C_obs, D_obs] = ssdata(sys_obs);
disp('The observable canonical form of the system is:');
disp('A_obs =');
disp(A_obs);
disp('B_obs =');
disp(B_obs);
disp('C_obs =');
disp(C_obs);
disp('D_obs =');
disp(D_obs);