|

楼主 |
发表于 2023-4-23 08:29
|
显示全部楼层
本帖最后由 天山草 于 2023-4-23 08:36 编辑
程序代码:
- Clear["Global`*"];
- (*三角形ABC各顶点坐标 \[DoubleLongRightArrow]*)
- \!\(\*OverscriptBox[\(b\), \(_\)]\) = b = 0; \!\(\*OverscriptBox[\(c\), \(_\)]\) = c; a = ((t^2 - 1) (t - c) + 2 I t (t - c))/(1 - c t + t^2);
- \!\(\*OverscriptBox[\(a\), \(_\)]\) = ((t^2 - 1) (t - c) - 2 I t (t - c))/(1 - c t + t^2);
- (*\[EmptyUpTriangle]ABC各顶角平分线与对边的交点坐标 D、E、F \[DoubleLongRightArrow]*)
- d = ((t^2 + 1) (c - t))/(c t - t^2 + 1); \!\(\*OverscriptBox[\(d\), \(_\)]\) = d; e = ( 2 c (t + I) (t - c))/(-c^2 + t^2 + 1);
- \!\(\*OverscriptBox[\(e\), \(_\)]\) = ( 2 c (t - I) (t - c))/(-c^2 + t^2 + 1); f = (c (t + I)^2)/( 2 c t - t^2 - 1);
- \!\(\*OverscriptBox[\(f\), \(_\)]\) = (c (t - I)^2)/(2 c t - t^2 - 1);
- (*\[EmptyUpTriangle]ABC的内心坐标 \[DoubleLongRightArrow]*) i = t + I; \!\(\*OverscriptBox[\(i\), \(_\)]\) = t - I;
- (*过ABC三点的圆的圆心坐标:*)
- WX[a_, b_, c_] := (a \!\(\*OverscriptBox[\(a\), \(_\)]\) (c - b) + b \!\(\*OverscriptBox[\(b\), \(_\)]\)(a - c) + c \!\(\*OverscriptBox[\(c\), \(_\)]\)(b - a) )/(\!\(\*OverscriptBox[\(a\), \(_\)]\)(c - b) + \!\(\*OverscriptBox[\(b\), \(_\)]\)(a - c) + \!\(\*OverscriptBox[\(c\), \(_\)]\)(b - a));
- \!\(\*OverscriptBox[\(WX\), \(_\)]\)[a_, b_, c_] := -((a \!\(\*OverscriptBox[\(a\), \(_\)]\)(\!\(\*OverscriptBox[\(c\), \(_\)]\) - \!\(\*OverscriptBox[\(b\), \(_\)]\)) + b \!\(\*OverscriptBox[\(b\), \(_\)]\)(\!\(\*OverscriptBox[\(a\), \(_\)]\) - \!\(\*OverscriptBox[\(c\), \(_\)]\)) + c \!\(\*OverscriptBox[\(c\), \(_\)]\)(\!\(\*OverscriptBox[\(b\), \(_\)]\) - \!\(\*OverscriptBox[\(a\), \(_\)]\)))/(\!\(\*OverscriptBox[\(a\), \(_\)]\)(c - b) + \!\(\*OverscriptBox[\(b\), \(_\)]\)(a - c) + \!\(\*OverscriptBox[\(c\), \(_\)]\)(b - a)));
- n = Simplify@ WX[d, e, f]; \!\(\*OverscriptBox[\(n\), \(_\)]\) = Simplify@ \!\(\*OverscriptBox[\(WX\), \(_\)]\)[d, e, f];
- (*复斜率定义:*) k[a_, b_] := (a - b)/(\!\(\*OverscriptBox[\(a\), \(_\)]\) - \!\(\*OverscriptBox[\(b\), \(_\)]\));
- Simplify@Solve[{(n - d) (\!\(\*OverscriptBox[\(n\), \(_\)]\) - \!\(\*OverscriptBox[\(d\), \(_\)]\)) == (n - x) (\!\(\*OverscriptBox[\(n\), \(_\)]\) - \!\(\*OverscriptBox[\(x\), \(_\)]\)), k[b, c] == k[b, x], x != d}, {x, \!\(\*OverscriptBox[\(x\), \(_\)]\)}];
- x = (c (c^2 (5 t^2 - 3) + c (2 t - 6 t^3) + (t^2 + 1)^2))/( 2 (2 c^3 t - c^2 (t^2 + 1) - 2 c (t^3 + t) + (t^2 + 1)^2));
- \!\(\*OverscriptBox[\(x\), \(_\)]\) = ( c (c^2 (5 t^2 - 3) + c (2 t - 6 t^3) + (t^2 + 1)^2))/( 2 (2 c^3 t - c^2 (t^2 + 1) - 2 c (t^3 + t) + (t^2 + 1)^2));
- Simplify@Solve[{(n - d) (\!\(\*OverscriptBox[\(n\), \(_\)]\) - \!\(\*OverscriptBox[\(d\), \(_\)]\)) == (n - y) (\!\(\*OverscriptBox[\(n\), \(_\)]\) - \!\(\*OverscriptBox[\(y\), \(_\)]\)), k[a, c] == k[c, y], y != e}, {y, \!\(\*OverscriptBox[\(y\), \(_\)]\)}];
- y = (c^3 t (t^2 + 8 I t + 1) - c^2 (t^4 + 13 I t^3 - 5 t^2 - 3 I t + 2) - c (t^5 - 4 I t^4 + 8 t^3 + 6 I t^2 - t + 2 I) + t (t - I)^2 (t + I)^3)/(2 (2 c^3 t^2 + c^2 (-5 t^3 + 2 I t^2 + t) + c (4 t^4 - 3 I t^3 - t^2 + I t - 1) - t^5 + I t^4 + t - I));
- \!\(\*OverscriptBox[\(y\), \(_\)]\) = ( c^3 t (t^2 - 8 I t + 1) - c^2 (t^4 - 13 I t^3 - 5 t^2 + 3 I t + 2) + c (-t^5 - 4 I t^4 - 8 t^3 + 6 I t^2 + t + 2 I) + t (t - I)^3 (t + I)^2)/(2 (2 c^3 t^2 + c^2 (-5 t^3 - 2 I t^2 + t) + c (4 t^4 + 3 I t^3 - t^2 - I t - 1) - t^5 - I t^4 + t + I));
- Simplify@Solve[{(n - d) (\!\(\*OverscriptBox[\(n\), \(_\)]\) - \!\(\*OverscriptBox[\(d\), \(_\)]\)) == (n - z) (\!\(\*OverscriptBox[\(n\), \(_\)]\) - \!\(\*OverscriptBox[\(z\), \(_\)]\)), k[b, a] == k[b, z], z != f}, {z, \!\(\*OverscriptBox[\(z\), \(_\)]\)}];
- z = -(((t + I) (c^3 (3 - 5 t^2) + c^2 t (11 t^2 - 5) - c (7 t^4 + 1) + t (t^2 + 1)^2))/( 2 (t - I) (c^3 t - c^2 (t^2 - 1) - c (t^3 + t) + t^4 - 1)));
- \!\(\*OverscriptBox[\(z\), \(_\)]\) = -(((t - I) (c^3 (3 - 5 t^2) + c^2 t (11 t^2 - 5) - c (7 t^4 + 1) + t (t^2 + 1)^2))/(2 (t + I) (c^3 t - c^2 (t^2 - 1) - c (t^3 + t) + t^4 - 1)));
- DX = Simplify[x - d] // Factor;
- EY = Simplify[Sqrt[(e - y) (\!\(\*OverscriptBox[\(e\), \(_\)]\) - \!\(\*OverscriptBox[\(y\), \(_\)]\))]] // Factor;
- EY = -(1/2) (7 c^4 t^3 - c^4 t - 20 c^3 t^4 - 10 c^3 t^2 + 2 c^3 + 18 c^2 t^5 + 12 c^2 t^3 - 6 c^2 t - 4 c t^6 + 2 c t^4 +
- 8 c t^2 + 2 c - t^7 - 3 t^5 - 3 t^3 - t)/(2 c^4 t^2 - 3 c^3 t^3 + c^3 t - c^2 t^4 - 2 c^2 t^2 - c^2 + 3 c t^5 + 2 c t^3 - c t - t^6 - t^4 + t^2 + 1);
- FZ = Simplify[Sqrt[(f - z) (\!\(\*OverscriptBox[\(f\), \(_\)]\) - \!\(\*OverscriptBox[\(z\), \(_\)]\))]] // Factor;
- FZ = 1/2 (8 c^4 t^3 - 8 c^4 t - 25 c^3 t^4 + 8 c^3 t^2 + c^3 + 27 c^2 t^5 + 10 c^2 t^3 - c^2 t - 11 c t^6 - 13 c t^4 - c t^2 + c + t^7 + 3 t^5 + 3 t^3 + t)/(2 c^4 t^2 - 3 c^3 t^3 + c^3 t - c^2 t^4 - 2 c^2 t^2 - c^2 + 3 c t^5 + 2 c t^3 - c t - t^6 - t^4 + t^2 + 1);
- FullSimplify[DX] // Factor
- FullSimplify[EY + FZ] // Factor
- Simplify[DX == EY + FZ]
复制代码 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
|