求教：matlab 2015b 中没有simple函数你是怎么解决的？

% 拷贝一个

function?[r,h]?=?simple(s,varargin)

%SIMPLE?Search?for?simplest?form?of?a?symbolic?expression?or?matrix.

%?SIMPLE(S)?tries?several?different?algebraic?simplifications?of

%?S,?displays?any?which?shorten?the?length?of?S's?representation,

%?and?returns?the?shortest.?S?is?a?SYM.?If?S?is?a?matrix,?the?result

%?represents?the?shortest?representation?of?the?entire?matrix,?which?is?

%?not?necessarily?the?shortest?representation?of?each?individual?element.

%

%?[R,HOW]?=?SIMPLE(S)?does?not?display?intermediate?simplifications,

%?but?returns?the?shortest?found,?as?well?as?a?string?describing

%?the?particular?simplification.?R?is?a?SYM.?HOW?is?a?string.

%?

%?SIMPLE(S,'IgnoreAnalyticConstraints',VAL)?controls?the?level?of?

%?mathematical?rigor?to?use?on?the?analytical?constraints?while?simplifying?

%?(branch?cuts,?division?by?zero,?etc).?The?options?for?VAL?are?TRUE?or?

%?FALSE.?Specify?TRUE?to?relax?the?level?of?mathematical?rigor

%?in?the?rewriting?process.?The?default?is?FALSE.

%

%

%?Examples:

%

%SRHow

%

%cos(x)^2+sin(x)^21simplify

%2*cos(x)^2-sin(x)^23*cos(x)^2-1?simplify

%cos(x)^2-sin(x)^2cos(2*x)?simplify

%cos(x)+i*sin(x)exp(i*x)?rewrite(exp)

%(x+1)*x*(x-1)x^3-xsimplify(100)

%x^3+3*x^2+3*x+1(x+1)^3simplify

%cos(3*acos(x))?4*x^3-3*xsimplify(100)

%

%simple(asin(sin(x)))?=?asin(sin(x))

%simple(asin(sin(x)),'IgnoreAnalyticConstraints',true)=?x

%

%?See?also?SYM/SIMPLIFY,?SYM/FACTOR,?SYM/EXPAND,?SYM/COLLECT.

%?Copyright?1993-2011?The?MathWorks,?Inc.

if?builtin('numel',s)?~=?1,s?=?normalizesym(s);end

p?=?nargout?==?0;

[rsym,h]?=?mupadSimple(s,p,varargin{:});

r?=?privResolveOutput(rsym,?s);

end

function?[r,h]?=?mupadSimple(s,p,varargin)

h?=?'';

r?=?s;

x?=?symvar(s,1);

%?parse?arguments:?Look?for?options

narg?=?nargin?-?2;

args?=?varargin;

%?default:

options?=?'null()';

k?=?1;

while?k?<=?size(args,?2)

v?=?args{k};

if?ischar(v)?&&?strcmp(v,?'IgnoreAnalyticConstraints')

if?k?==?size(args,?2);

error(message('symbolic:sym:optRequiresArg',?v))

end

value?=?args{k+1};

if?value?==?true

value?=?'TRUE';

elseif?value?==?false

value?=?'FALSE';

elseif?strcmp(v,?'IgnoreAnalyticConstraints')?&&?isa(value,?'char')

if?strcmp(value,?'all')

value?=?'TRUE';

elseif?strcmp(value,?'none')

value?=?'FALSE';

else

error(message('symbolic:sym:badArgForOpt',?v))

end

else

error(message('symbolic:sym:badarg',?value))

end

options?=?[options?',?'?v?'='?value];?%#ok<AGROW>

args(k:k+1)?=?[];

else

error(message('symbolic:sym:badarg',?v))

end

end

%?Try?the?different?simplifications.

[r,h]?=?simpler('simplify',s,r,h,p,options);

[r,h]?=?simpler('radsimp',s,r,h,p,'null()');

[r,h]?=?simpler('symobj::simplify',s,r,h,p,options,'100');

[r,h]?=?simpler('combine',s,r,h,p,options,'sincos');

[r,h]?=?simpler('combine',s,r,h,p,options,'sinhcosh');

[r,h]?=?simpler('combine',s,r,h,p,options,'ln');

[r,h]?=?simpler('factor',s,r,h,p,'null()');

[r,h]?=?simpler('expand',s,r,h,p,options);

[r,h]?=?simpler('combine',s,r,h,p,options);

[r,h]?=?simpler('rewrite',s,r,h,p,'null()','exp');

[r,h]?=?simpler('rewrite',s,r,h,p,'null()','sincos');

[r,h]?=?simpler('rewrite',s,r,h,p,'null()','sinhcosh');

[r,h]?=?simpler('rewrite',s,r,h,p,'null()','tan');

[r,h]?=?simpler('symobj::mwcos2sin',s,r,h,p,'null()');

if?~isempty(x)

[r,h]?=?simpler('collect',s,r,h,p,'null()',x);

end

end

function?[r,h]?=?simpler(how,s,r,h,p,options,x)

%SIMPLER?Used?by?SIMPLE?to?shorten?expressions.

%?SIMPLER(HOW,S,R,H,P,X)?applies?method?HOW?with?optional?parameter?X

%?to?expression?S,?prints?the?result?if?P?is?nonzero,?compares?the

%?length?of?the?result?with?expression?R,?which?was?obtained?with

%?method?H,?and?returns?the?shortest?string?and?corresponding?method.

if?nargin?<?7

[t,err]?=?mupadmex('symobj::map',s.s,how,options);

elseif?ischar(x)

[t,err]?=?mupadmex('symobj::map',s.s,how,x,options);

else

[t,err]?=?mupadmex('symobj::map',s.s,how,x.s,options);

end

if?err

return;

end

if?nargin?==?7

how?=?[how?'('?char(x)?')'];

end

how?=?strrep(how,'symobj::','');

if?p?

loose?=?isequal(get(0,'FormatSpacing'),'loose');

if?loose,?disp('?'),?end

disp([how?':'])

if?loose,?disp('?'),?end

disp(t)

end

cmp?=?mupadmex('symobj::simpler',?t.s,?r.s,?0);

if?strcmp(cmp,'TRUE')?

r?=?t;

h?=?how;

end

end