Question

matlab中对函数求导并求极值

z=exp(4.1399-0.0235*x)
y=0.000157*z*(1-z)*(1+z)^3/((1+z)^4+0.00004476*z^2)^(3/2)
求y的一阶导数的极值，x的取值范围为[1,225]
求一个代码，谢谢！

Answer 1

用matlab对函数求导及其极值 ，可以先用diff（）函数求出y'，然后令y'=0，用vpasolve（）解出x值。实现代码如下：

syms x

z=exp(4.1399-0.0235*x);

y=0.000157*z*(1-z)*(1+z)^3/((1+z)^4+0.00004476*z^2)^(3/2);

dy=diff(y,1)

x=vpasolve(dy==0)

运行结果

x=120.125

追问

这个函数的极值应该不止一个，请问如何都计算出来呢？

追答

是软件的问题。如用R14可以求得全解。

x=solve（dy）

实数解为232.206，120.125