在 Matlab 中使用循环从表达式中查找最大值

Question

当x介于0和1之间时，我想使用表达式的二阶导数找到最大值。 换句话说，我两次取cox(x^2)的导数得到二阶导数，结果是- 2*sin(x^2) - 4*x^2*cos(x^2) ，然后我想在x = 0到x = 1处计算这个二阶导数，并显示填充值的maximum 。

我有：

syms x
f = cos(x^2);
secondD = diff(diff(f));

for i = 0:1
y = max(secondD(i))
end

有人可以帮忙吗？

Answer 1

您可以通过subs和double轻松完成：

syms x
f = cos(x^2);
secondD = diff(diff(f));

% instead of the for loop
epsilon = 0.01;
specified_range = 0:epsilon:1;
[max_val, max_ind] = max(double(subs(secondD, specified_range)));

请注意，这是一种寻找最大值的数值方法，并且返回的答案并非始终完全正确。 但是，通过增加epsilon ，您通常可以期待更好的结果（同样在某些情况下它并不完全正确）。