【线性系统理论作业：传递函数转为状态空间算法】

一、题目要求

编程作业1：

根据I/O描述到状态空间描述的两种不同算法，编程实现并给出能控、能观标准型

输入为传递函数系数（a，b），输出为能控、能观标准型（A,B,C,D)

要求：

程序的通用性，任意阶次系统均可检验；语言不限；

二、编程解答

本次解答采用matlab语言编写，原创作品，仿冒必究。

%**************************************************************************%
% 编程作业1：根据I/O描述到状态空间描述的两种不同算法，编程实现并给出能控、能观标准型 
%            输入为传递函数系数（a，b），输出为能控、能观标准型（A,B,C,D)
% 要求：程序的通用性，任意阶次系统均可检验；
%**************************************************************************%

% 测试样例
% a=[1,14,56,160];
% b=[0,0,1,0];

% 输入参数
clc
clear
% 输入传递函数时务必同样长度，长度不足用0补齐
a=input('输入矩阵a:');
b=input('输入矩阵b:');
clc
a
b

num=b;
den=a;
bn=num(1,1);
m=size(num,2)-1;
n=size(den,2)-1;

num0=num(:,2:n+1);
den0=den(:,2:n+1);

numf=fliplr(num0);
denf=fliplr(den0);
denx=denf.*bn;

n0=zeros(1,n-1);
A=eye(n-1);

% 以下为能控标准型计算
Ac=[n0' A;-denf];
Bc=[n0,1]';
Cc=numf-denx;
Dc=bn;

% 求解能控标准型的未知数
Bs=[bn];
for i=1:n
   Bs0=num(1,i+1)-den(1,2:i+1)*Bs;
   Bs=[Bs;Bs0];
end

% 以下为能观标准型计算
As=[n0' A;-denf];
Bs=Bs(2:n+1,:);
Cs=[1,n0];
Ds=bn;

% 程序运行结果
if m>n
    disp("分子阶次大于分母阶次，系统不稳定")
else
    disp("输出：")
    disp("能控标准型状态空间表示：")
    Ac
    Bc
    Cc
    Dc
    disp("能观标准型状态空间表示：")
    As
    Bs
    Cs
    Ds
end