weng1.m

来自「运用MATLAB对控制系统进行仿真操作」· M 代码 · 共 53 行

function [A,B,C,D]=weng(N,M,q,value_n,value_m)

if N(1,1)==1
   A(value_n,:)=-fliplr(N(2:end));
   A(1:value_n-1,1)=zeros(value_n-1,1);
   A(1:value_n-1,2:l)=eye(value_n-1,value_n-1);
   A
   B(1:value_n-1,1)=zeros(value_n-1,1);
   B(value_n,1)=eye(1,1);
   B 
if value_m<value_n   
   C(1,1:(value_m+1))=fliplr(M(1:end));
   C(1,(value_m+2):(value_n-value_m))=zeros(1,value_n-value_m-1);
   C
   
   D=0;
   D
end;
if value_m==value_n
    D=fliplr(M(2:end))-M(1,1)*fliplr(N(2:end));
        D
        E(1,1)=M(1,1);
        E
end
end


if N(1,1)~=1
    z=N(1,1);
    N1=N/z;
    M1=M/z;
   A(value_n,:)=-fliplr(N1(2:end));
   A(1:value_n-1,1)=zeros(value_n-1,1);
   A(1:value_n-1,2:value_n)=eye(value_n-1,value_n-1);
   A
   B(1:value_n-1,1)=zeros(value_n-1,1);
   B(value_n,1)=eye(1,1);
   B 
 if value_m<value_n 
   C(1,1:(value_m+1))=fliplr(M1(1:end));
   C(1,(value_m+2):(value_n-value_m))=zeros(1,value_n-value_m-1);
   C
   
   D=0;
   D
end;
 if value_m==value_n
        D=fliplr(M1(2:end))-M1(1,1)*fliplr(N1(2:end));
        D
        E(1,1)=M1(1,1);
        E
    end
end