compan.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 10 行

function A = compan (c)
if min(size(c)) > 1
    error('Input argument must be a vector.')
end
n = length(c);
if n <= 1, A = [];
elseif n == 2, A = -c(2)/c(1); 
else, c = c(:).'; A = sym(diag(ones(1,n-2),-1)); A(1,:) = -c(2:n) ./ c(1);
end