Matlab里如何求解含参数二元函数的极值？

函数形式为
M/(t*n*d)+K/(t*d)+V*h/(t*T*(r+1)*(r+2))*((t+T)^(r+2)-T^(r+2)-t^(r+2)
其中 T和t 是两个变量，其它的都是参数

现在要让函数对T和t求偏导等于零，解出函数的最优解

1 如果含参的话，使用哪个命令？
2 如果求数值解的话应该怎么办？（参数有赋值的情况下）

按道理说，方程不复杂的话，这样既可：

syms M t n d K V h r T
f=M/(t*n*d)+K/(t*d)+V*h/(t*T*(r+1)*(r+2))*((t+T)^(r+2)-T^(r+2)-t^(r+2))
dfdT=diff(f,T)
dfdt=diff(f,t)
solve(dfdT,dfdt,'T,t')

但是由于方程太复杂，没有解析解。
Warning: Warning, solutions may have been lost
Warning: Explicit solution could not be found.

ans =

[ empty sym ]

syms a b s; s=30*sqrt(1-((20-30*sin(b))/30)^2)-30*cos(b); diff(s,b); value_a=solve('30*sin(b) - (30*cos(b)*(sin(b) - 2/3))/(1 - (sin(b) - 2/3)^2)^(1/2)=0','b'); value_b=solve('30*sin(a)+30*sin(value_a)=20','a'); value_a%解得的a的值 value_b%解得的b的值