所以我距离成功编写此功能还有几步之遥:
这是我到目前为止所做的:
function a=myanalyzecp(f,a,b)
syms x;
v=coeffs(f(x)); % grabs function's coefficients
vertex=(-v(2))/(2*v(3)); % vertex formula
if (a<vertex && vertex<b)
if (diff(diff(f(x)))>0) % f''>0 means minima
a=1;
else
a=-1;
end
else
a=0;
end
我遇到的问题是当函数只有1或2个项时,例如x ^ 2或x ^ 2 + 4或x ^ 2 + 4 * x。因为那时我的顶点函数失败了
答案 0 :(得分:0)
由于v
可能只包含一个或两个元素,因此添加一个测试,例如
if length(v)==1,
vertex = ...
elseif length(v)==2
vertex = ...
else
vertex=(-v(2))/(2*v(3));
end
这应该取代vertex=(-v(2))/(2*v(3));
答案 1 :(得分:0)
要获取所有系数,包括零,请使用sym2poly
代替coeffs
:
v = sym2poly(f(x));
这会将结果显示为double
向量,而不是sym
。如果您需要它是sym
,您需要转换:
v = sym(sym2poly(f(x)));