clc;
clear all;
close all;
input('輸入微分方程左側y前的系數(傳遞函數分母s前的系數），從y的0階導（從s的0次方）開始依次輸入');
input('(最高階的導數系數（s最高次冪系數）為1,不要輸入），構成行數為1的矩陣?；剀囬_始輸入。');
y1=input('y1=');
input('輸入微分方程右側u前的系數(傳遞函數分子s前的系數），從u的0階導（從s的0次方）開始依次輸入');
input(',構成行數為1的矩陣。回車開始輸入。');
u1=input('u1=');

[m1,n1]=size(y1);
[m1 n2]=size(u1);
for i=1:n1
    a1(i)=y1(n1-i+1);
end
for i=1:n2
    b1(i)=u1(n2-i+1);
end
G=tf(b1,a1)
u=zeros(1,n1-1);
for i=1:n2
    u(i)=u1(i);
end
u=u/y1(n1);
y1=y1/y1(n1);
y=zeros(1,n1-1);
for i=1:(n1-1)
    y(i)=y1(i);
end
n1=n1-1;
A=zeros(n1,n1);
C=zeros(n1,n1-1);

input('輸入y的初始條件，從y的0階導數值開始輸入，構成一個列數為1的矩陣。回車開始輸入。');
Y=input('Y=');
input('輸入u的初始條件，從u的0階導數值開始輸入，構成一個列數為1的矩陣?；剀囬_始輸入。');
U=input('U=');

for i=1:n1
    A(i,n1+1-i)=1;
end

for i=1:(n1-1)
    k=i;
    for j=1:(n1-i)
        k=k+1;
        A(i,j)=y(k);
    end
end

for i=1:(n1-1)
    k=i;
    for j=1:(n1-i)
        k=k+1;
        C(i,j)=0-u(k);
    end
end

input('狀態初值為：');
x=A*Y+C*U