有没有什么化简的方法，计算值只要得到大概的一个值就可以啦
整个mathmatic的代码如下：
\[Alpha] =.(*Subscript[r, 0]'*)
\[Beta] =.(*Subscript[r, 0]*)
\[Zeta] =.(*Subscript[\[Theta], h]*)
\[Epsilon] =.(*Subscript[\[Theta], 0]*)
T = \[Alpha]/\[Beta]
h = 20
Subscript[h, sr] = 1
\[CurlyPhi] = 10*Pi/180
c = 10
\[Gamma] = 19
\[Psi] = h/Subscript[h, sr]
r = \[Beta]*E^((\[Theta] - \[Epsilon])*Tan[\[CurlyPhi]])
\[CurlyRho] = \[Alpha]*
E^(-(\[Theta] - \[Epsilon])*Tan[\[CurlyPhi]])(*r'*)
\[CurlyTheta] = (r + \[CurlyRho])/2(*Subscript[r, m]*)
R = (r - \[CurlyRho])/2
a = Sin[\[Epsilon]]/Sin[\[Theta]]*\[Beta] - \[CurlyTheta]
d = Cos[\[Zeta]]/Cos[\[Theta]]*\[Beta]*
E^((\[Zeta] - \[Epsilon])*Tan[\[CurlyPhi]]) - \[CurlyTheta]
\[Omega] = Sqrt[R^2 - a^2](*Subscript[x, 1]*)
\[Tau] = Sqrt[R^2 - d^2](*Subscript[x, 2]*)
\[Sigma] = Sqrt[R^2 - x^2](*Subscript[y, 1]*)
Subscript[f, 01] = \[Psi]*1/(
Sin[\[Zeta]]*E^((\[Zeta] - \[Epsilon])*Tan[\[CurlyPhi]]) -
Sin[\[Epsilon]])
Subscript[f, 02] =
Sin[\[Zeta] - \[Epsilon]]/Sin[\[Zeta]] -
Cos[\[Zeta]]/
Sin[\[Zeta]]*(Sin[\[Zeta]]*
E^((\[Zeta] - \[Epsilon])*Tan[\[CurlyPhi]]) - Sin[\[Epsilon]])
\[Eta] = ArcTan[Sin[\[Epsilon]]/(
Cos[\[Zeta]]*
E^((\[Zeta] - \[Epsilon])*
Tan[\[CurlyPhi]]))](*Subscript[\[Theta], B]=...*)
\[Lambda] =
ArcCos[(Subscript[f, 01]*Cot[\[Eta]]*
Sin[\[Epsilon]])/(\[Sqrt](1 +
Subscript[f, 01]*
Sin[\[Epsilon]]*(2 +
Subscript[f, 01]*(Csc[\[Eta]])^2*
Sin[\[Epsilon]])))](*Subscript[\[Theta], sr]=...*)
FindMinimum[{1/(2*\[Gamma]*(Integrate[(\[CurlyTheta] + y)^2*
Cos[\[Theta]], {\[Theta], \[Epsilon], \[Eta]}, {x,
0, \[Omega]}, {y, a, \[Sigma]}] +
Integrate[(\[CurlyTheta] + y)^2*
Cos[\[Theta]], {\[Theta], \[Eta], \[Zeta]}, {x,
0, \[Tau]}, {y, d, \[Sigma]}])) (c*
Cot[\[CurlyPhi]]*(\[Beta])^2*(E^(
2*(\[Zeta] - \[Epsilon])*Tan[\[CurlyPhi]])*(Cos[\[Zeta]])^2*
Integrate[
Sin[\[Theta]]/(Cos[\[Theta]])^3*\[Tau], {\[Theta], \[Eta], \
\[Zeta]}] - (Sin[\[Epsilon]])^2*
Integrate[
Cos[\[Theta]]/(Sin[\[Theta]])^3*\[Omega], {\[Theta], \
\[Epsilon], \[Eta]}]) +
Subscript[\[Gamma],
sr]*((\[Beta])^3*E^(
3*(\[Zeta] - \[Epsilon])*Tan[\[CurlyPhi]])*(Cos[\[Zeta]])^3*
Integrate[(2*\[Tau]*
Tan[\[Theta]]*(Tan[\[Theta]] -
Tan[\[Lambda]]))/(Cos[\[Theta]])^3, {\[Theta], \
\[Lambda], \[Zeta]}])), Pi/2 > \[Zeta] > \[Epsilon] > 0,
1 > T > 0}, {\[Epsilon], \[Zeta], T}]

13q1RUIzdMlpnTKnwhWk7tA

补充一下

笔记本文件我就不看了……目测把Integrate用NIntegrate加上NumericQ改写可能还有点希望，但是也够呛。