huangkun1.m

目目目目目目目目目目目目目目目目目目目目目

clc;
clear all;
close all;
input('输入传递函数的分母的系数，从S的指数为0开始依次输入。回车开始输入。');
a=input('a=');
input('输入传递函数的分子的系数，从S的指数为0开始依次输入。回车开始输入。');
b=input('b=');

[m1,n1]=size(a);%取行数和列数
[m2,n2]=size(b);
for i=1:n1%顺序前后倒一下，为了能够输出输入的传递函数
    a1(i)=a(n1-i+1);
end
for i=1:n2
    b1(i)=b(n2-i+1);
end
G=tf(b1,a1)

b=b/a(n1);%化为首一
a=a/a(n1);

A=zeros(n1-1,n1-1);%定义A阵
i=n1-2;
for j=1:i%先弄好为1的项
    A(j,j+1)=1;
end
i=n1-1;
for j=1:i%最后一行的赋值
    A(n1-1,j)=0-a(j);
end
B=zeros(n1-1,1);%定义B阵
B(n1-1,1)=1;%B最后一行为1

C=zeros(1,n1-1);
if n1==n2%判断是不是分子分母同阶
    D=b(n2);
    for j=1:(n2-1)
        C(1,j)=b(j)-b(n2)*a(j);
    end
else   
C=zeros(1,n1-1);
D=0;
for j=1:n2
    C(1,j)=b(j);
end    
end   

input('能控标准型：');
A
B
C
D

input('能观标准型：');
A1=A'
B1=C'
C1=B'
D1=D