|
本帖最后由 TSC999 于 2023-12-18 13:47 编辑 这个题目,用复斜率解析几何很容易做。
两个弧长相等等价于两个弦长相等。把 F 点的复数坐标(是一个实数)作为未知数,可求出含有未知数的 G 和 G' 的坐标。最后由 CA=GG' 的条件,可求出 F 点的坐标。于是可求得三角形 DEF 的面积精确到 50 位的数字值为: 8.0820177829901484072083229106793268984084477052172....
面积的理论表达式见后面的图片。
- Clear["Global`*"];
- \!\(\*OverscriptBox[\(b\), \(_\)]\) = b = 0; c = -2 Sqrt[2] I; \!\(\*OverscriptBox[\(c\), \(_\)]\) = 2 Sqrt[2] I; d = -2 + 2 I;
- \!\(\*OverscriptBox[\(d\), \(_\)]\) = -2 - 2 I; e = -Sqrt[2] I; \!\(\*OverscriptBox[\(e\), \(_\)]\) = Sqrt[2] I; R = 2 Sqrt[2];
- \!\(\*OverscriptBox[\(f\), \(_\)]\) = f;
- k[a_, b_] := (a - b)/(\!\(\*OverscriptBox[\(a\), \(_\)]\) - \!\(\*OverscriptBox[\(b\), \(_\)]\)); (*复斜率定义*)
- W1 = {g, \!\(\*OverscriptBox[\(g\), \(_\)]\)} /. Simplify@Solve[{(b - g) (\!\(\*OverscriptBox[\(b\), \(_\)]\) - \!\(\*OverscriptBox[\(g\), \(_\)]\)) == R^2, k[e, f] == k[e, g]}, {g, \!\(\*OverscriptBox[\(g\), \(_\)]\)}] // Flatten;
- g = Part[W1, 3]; \!\(\*OverscriptBox[\(g\), \(_\)]\) = Part[W1, 4];
- W2 = {g1, \!\(\*OverscriptBox[\(g1\), \(_\)]\)} /. Simplify@Solve[{(b - g1) (\!\(\*OverscriptBox[\(b\), \(_\)]\) - \!\(\*OverscriptBox[\(g1\), \(_\)]\)) == R^2, k[d, f] == k[d, g1]}, {g1, \!\(\*OverscriptBox[\(g1\), \(_\)]\)}] // Flatten;
- g1 = Part[W2, 1]; \!\(\*OverscriptBox[\(g1\), \(_\)]\) = Part[W2, 2];
- W3 = {a, \!\(\*OverscriptBox[\(a\), \(_\)]\)} /. Simplify@Solve[{(b - a) (\!\(\*OverscriptBox[\(b\), \(_\)]\) - \!\(\*OverscriptBox[\(a\), \(_\)]\)) == R^2, k[d, e] == k[d, a]}, {a, \!\(\*OverscriptBox[\(a\), \(_\)]\)}] // Flatten;
- a = Part[W3, 1]; \!\(\*OverscriptBox[\(a\), \(_\)]\) = Part[W3, 2];
- W = {f} /. Simplify@Solve[{(c - a) (\!\(\*OverscriptBox[\(c\), \(_\)]\) - \!\(\*OverscriptBox[\(a\), \(_\)]\)) == (g - g1) (\!\(\*OverscriptBox[\(g\), \(_\)]\) - \!\(\*OverscriptBox[\(g1\), \(_\)]\)), k[b, f] == 1}, {f}] // Flatten;
- f = Part[W, 2]; \!\(\*OverscriptBox[\(f\), \(_\)]\) = f;
- Print["F = ", f, " \[TildeTilde] ", N@Re[f]];
- S[a_, b_, c_] := (\!\(\*OverscriptBox[\(a\), \(_\)]\) (b - c) + \!\(\*OverscriptBox[\(b\), \(_\)]\) (c - a) + \!\(\*OverscriptBox[\(c\), \(_\)]\) (a - b))/(4 I);(*已知三角形各顶点坐标,求其有向面积,ABC逆时针向环绕*)
- Print["\[EmptyUpTriangle]DEF 的面积 = ", Simplify@S[d, e, f], " \[TildeTilde] ", N[Re@S[d, e, f], 50]];
复制代码
运行结果为:
- F = 1/69 (-84+72 Sqrt[2]-(2 I 3^(2/3) (-887+113 Sqrt[2]) (-I+Sqrt[3]))/Power[145566-143667 Sqrt[2]+23 I Sqrt[-83346270+66834876 Sqrt[2]], (3)^-1]+(1-I Sqrt[3]) Power[436698-431001 Sqrt[2]+69 I Sqrt[-83346270+66834876 Sqrt[2]], (3)^-1]) \[TildeTilde] 3.90591
- \[EmptyUpTriangle]DEF 的面积 = (2 Power[3, (6)^-1] (1548+661 Sqrt[2]) (3 I+Sqrt[3])+6 (-4+33 Sqrt[2]) Power[145566-143667 Sqrt[2]+23 I Sqrt[-83346270+66834876 Sqrt[2]], (3)^-1]+Power[3, (3)^-1] (2+Sqrt[2]) (1-I Sqrt[3]) (145566-143667 Sqrt[2]+23 I Sqrt[-83346270+66834876 Sqrt[2]])^(2/3))/(138 Power[145566-143667 Sqrt[2]+23 I Sqrt[-83346270+66834876 Sqrt[2]], (3)^-1]) \[TildeTilde] 8.0820177829901484072083229106793268984084477052172
复制代码
运行结果的图片:
上面一行是 F 点坐标的理论表达式及大约的数字值。理论表达式应该能够进一步简化的!请高人简化。
下面一行是三角形 DEF 面积的理论表达式及精确到50位的数字值。面积的理论表达式也应该能够简化的!请高人简化。 |
|