coprime.m

用MATLAB的脚本语言M语言编写的

function [N,M,X,Y]=coprime(G,K1,K2)
G=ss(G); a=G.a; b=G.b; c=G.c; d=G.d;
if length(K1)==1, K1=K1*ones(size(c)); end
if length(K2)==1, K2=K2*ones(size(c)); end
F=-acker(a,b,K1); H=-acker(a',c',K2)';
M=ss(a+b*F,b,F,1);N=ss(a+b*F,b,c+d*F,d);
X=ss(a+H*c,H,F,0);Y=ss(a+H*c,-b-H*d,F,1);