diopha_eq.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function [X,Y]=diopha_eq(A,B,C)
AAA=[]; A1=A(:); B1=B(:);
n=length(B)-1; m=length(A)-1; k=1;
for i=1:n, AAA(i:i+length(A1)-1,k)=A1; k=k+1; end
for i=1:m, AAA(i:i+length(B1)-1,k)=B1; k=k+1; end
C1=zeros(n+m,1); C1(1:length(C))=C(:); x=AAA\C1;
X=x(1:n)'; Y=x(n+1:end)';