gusucode.com > symbolic工具箱matlab源码程序 > symbolic/@sym/minpoly.m
function P = minpoly(A,varargin) %MINPOLY Minimal polynomial of a matrix. % P = MINPOLY(A) returns the coefficients of the minimal polynomial % of the matrix A. If A is a SYM object, the vector returned is a % SYM object, too. Otherwise a vector with doubles is returned. % % P = MINPOLY(A,x) returns the minimal polynomial of the matrix A in % terms of the variable x. Here x must be a free symbolic variable. % % Examples: % syms x; % A = sym([1 1 0; 0 1 0; 0 0 1]); % % P = minpoly(A) % returns [ 1, -2, 1] % % P = sym2poly(minpoly(A,x)) % returns 1 -2 1 % % P = minpoly(A,x) % returns x^2 - 2*x + 1 % % P = poly2sym(minpoly(A),x) % returns x^2 - 2*x + 1 % % See also CHARPOLY, SYM/POLY2SYM, SYM/SYM2POLY, SYM/JORDAN, SYM/EIG, % SOLVE. % Copyright 2012-2014 The Mathworks, Inc. oldDigits = digits(16); cleanupObj = onCleanup(@() digits(oldDigits)); p = inputParser; p.addRequired('A', @checkA); p.addOptional('x', sym([]), @checkx); p.parse(A,varargin{:}); A = p.Results.A; x = p.Results.x; if ~isa(A,'sym') Asym = sym(A); else Asym = A; end if ~isa(x,'sym'), x = sym(x); end if builtin('numel',Asym) ~= 1, Asym = normalizesym(Asym); end if builtin('numel',x) ~= 1, x = normalizesym(x); end if any(reshape(~isfinite(Asym),[],1)) error(message('symbolic:sym:InputMustNotContainNaNOrInf')); end if isempty(x) P = privUnaryOp(Asym,'symobj::minpoly'); if ~isa(A,'sym') P = cast(P,'like',A); end else P = privBinaryOp(Asym,x,'symobj::minpoly'); end end function checkA(x) if isempty(x) error(message('symbolic:sym:NonEmptySymExpected')); end end function checkx(x) expr = sym(x); if ~sym.isVariable(expr) error(message('symbolic:sym:SymVariableExpected')); end end