ll.m

来自「欧拉前差公式通用程序 欧拉前差公式的算法」· M 代码 · 共 44 行

disp('分子分母同阶次时传递函数转化为状态空间的通用程序');
disp('以下a、c项数为1～n')
a=[1,16,194,640,11];
c=[4,0,160,720,567];
n=length(a)-1;
b=zeros(1,(n+1));
b(1)=c(1);
for i=0:n
    b(i+1)=c(i+1);
    
    for j=1:i
        b(i+1)=b(i+1)-a(j+1)*b(i+1-j);
    end
end
for p=1:n
    for q=1:n
        if p<n
            if q==p+1
                A(p,q)=1;
            else
                A(p,q)=0;
            end
        else
            A(p,q)=-a((n+2)-q);
        end
    end
end
A
f=1;
for e=1:n;
    B(e,f)=b(e+1);
end
B
w=1;
for l=1:n;
    if l==1
        C(w,l)=1;
    else
        C(w,l)=0;
    end 
end
C
D=b(1);
D