nyy 发表于 2023-8-15 18:27
@王守恩 又没有任何有效的回答，你激动个啥？
不知道是不是本论坛的人提的问题。
area :=Module[{n, m, out}, n=(b - a - Sin[(b - a)])R^2/2; m=Abs[(Cos - Cos) (Sin - Sin)]R^2/2; out = n + m]
ans = FindRoot[{area ==12, area == 20, area ==25}, {{R, 4}, {a, 10/60}, {b, 80/60}}]
S = Pi R^2 - area - area - area /. ans
Solve[{Pi*R^2/2==S+12+(a+Sin/2)R^2==20+25-(a+Sin/2)R^2==12+20+(b+Sin/2)R^2==25+S-(b+Sin/2)R^2,R>0,S>0,30/60>a>0,30/60>b>0},{a,b,R,S}]
求助:我就是不知道怎么把(4)装进(1)(2)(3)去。
或:(3)装进(2)?(1)装进(2)?
太难了，登天还难。
本帖最后由 Jack315 于 2023-8-16 06:53 编辑
31.6762 的图是这个样子的：
圆半径为 r=5.1609。
P 点坐标为 x=0.9587 和 y=-0.9587
位于 135 度线上。
（下面这个图删不掉了:L ）
https://bbs.emath.ac.cn/data/attachment/forum/202308/14/152445wslynlnztaxjtf3u.png
..
借用24#的图. 直接面积加减法.
$S_1 = S_{扇形OCD} - S_{三角形OPD} - S_{三角形OPC}$
$S_2 = S_{扇形OBC} + S_{三角形OPC} - S_{三角形OPB}$
$S_3 = S_{扇形OAB} + S_{三角形OPA} + S_{三角形OPB}$
$S_4 = S_{扇形OAD} + S_{三角形OPD} - S_{三角形OPA}$
设$OP=d,u = d \cos (\alpha ) = R \cos (x)$, $v =d \sin (\alpha ) = R \cos (y)$, 列出表达式为
$-\frac{1}{2} d R \sin (x-\alpha )+\frac{1}{2} d R \cos (\alpha +y)+\frac{1}{2} R^2 \left(x+y-\frac{\pi }{2}\right)=12$
$-\frac{1}{2} d R \sin (\alpha +x)-\frac{1}{2} d R \cos (\alpha +y)+\frac{1}{2} R^2 \left(x-y+\frac{\pi }{2}\right)=20$
$\frac{1}{2} d R \sin (\alpha +x)+\frac{1}{2} d R \cos (y-\alpha )+\frac{1}{2} R^2 \left(-x-y+\frac{3 \pi }{2}\right)=25$
$\frac{1}{2} d R \sin (x-\alpha )-\frac{1}{2} d R \cos (y-\alpha )+\frac{1}{2} R^2 \left(-x+y+\frac{\pi }{2}\right)=S$
代入 $d \cos (\alpha ) \to R \cos (x)$, $d \sin (\alpha ) \to R \cos (y)$ ,得到
$R^2[ (-\pi/4 + (x + y)/2) - (Sin + Sin)/4 + Cos Cos] = 12$,
$R^2[(\pi/4 + (x - y)/2) - (Sin - Sin)/4 - Cos Cos ]= 20$
$R^2[ (\frac{3\pi}{4} - (x + y)/2) + (Sin + Sin)/4 + Cos Cos ]= 25$
$R^2[ (\pi/4 - (x - y)/2) + (Sin - Sin)/4 - Cos Cos] = S$
sol=TrigExpand[{1/2R^2(x+y-Pi/2)-1/2R d Sin]-1/2R d Sin)]==12,
1/2R^2(x-y+Pi/2)-1/2R d Sin]+1/2R d Sin)]==20,
1/2R^2(2Pi-(x+y+Pi/2))+1/2R d Sin]+1/2R d Sin)]==25,
1/2R^2(y-x+Pi/2)+1/2R d Sin]-1/2R d Sin)]==S}]/.{dCos[\]->R Cos,dSin[\]->R Cos}
Eliminate
跟7# 和9# 的结构 都比较类似.
\[\left\{ \begin{array}{l}
64 \cos (x) \cos (y)-37 x+37 \sin (x) \cos (x)+16 \pi =0 \\
180 \cos (x) \cos (y)+74 y-74 \sin (y) \cos (y)-29 \pi =0 \\
S = \frac{17 \pi -228 \cos (x) \cos (y)}{4 \cos (x) \cos (y)+\pi } \\
\end{array} \right.\]
数值解就是
FindRoot
{x->1.479562523482503437776738698107319809676480650279014451160437712048379185605335033904215447591948991,
y->1.377272602814462853186131203706929797854803642104178086827821236519830595354315839796202069063838094,
R->4.800094607277248765851985987118642066763331826398128627436431334969912416189838911534091245956958841,
S->15.38514805508871342500004107511989278433570526627393947940780984580291479054228346327225883902552547}
画图就是
再给一个全局图
本帖最后由 Jack315 于 2023-8-16 07:44 编辑
wayne 发表于 2023-8-16 07:07
..
这个计算方法应该是比较简洁的一种。
只是在计算角度时，如果 P 点在不同象限，
公式在形式上可能会有差别。
如果用公式能得到六个解，就没有问题。
否则可能就是上面说的这个问题。
另外，能比较详细地讲解下这个帖子 [原创] 求红色线段的长度值 中 8#
关于无穷级数求和的方法吗？说明一下利用留数定理的具体思路也行。
我估计自己就卡在哪个关键点上，还望大佬点拨一下:D
再笼统的给一个一般化的方程, 除了第一个方程,每个方程的几何意义都很明显,谁能解释一下第一个方程的几何意义?
$\cos (x) \cos (y)=\frac{\pi(-S+\text{S1}-\text{S2}+\text{S3})}{4 (S+\text{S1}+\text{S2}+\text{S3})},2 x-\sin (2 x)=\frac{2 \pi(\text{S1}+\text{S2})}{S+\text{S1}+\text{S2}+\text{S3}} , 2 y-\sin (2 y)=\frac{2 \pi(S+\text{S1})}{S+\text{S1}+\text{S2}+\text{S3}},\pi R^2=S+\text{S1}+\text{S2}+\text{S3}$
Block[{S1 = 12, S2 = 20, S3 = 25},
FindRoot[{Cos Cos == (\ (-S + S1 - S2 + S3))/(4 (S + S1 + S2 + S3)), (2 \ (S1 + S2))/(S + S1 + S2 + S3) + Sin ==2 x,
(2 \ (S + S1))/(S + S1 + S2 + S3) + Sin == 2 y, \ R^2 == S + S1 + S2 + S3}, {{x, 1}, {y, 1}, {S, S1}, {R, 10}}, WorkingPrecision -> 100]]
{x->1.479562523482503437776738698107319809676480650279014451160437712048379185605335033904215447591948991,
y->1.377272602814462853186131203706929797854803642104178086827821236519830595354315839796202069063838094,
S->15.38514805508871342500004107511989278433570526627393947940780984580291479054228346327225883902552547,
R->4.800094607277248765851985987118642066763331826398128627436431334969912416189838911534091245956958841}
Block[{S1=20,S2=12,S3=25},FindRoot[{2(S1+S2)/(2a-Sin)==2(S1+S)/(2b-Sin)==(S1-S2+S3-S)/(4CosCos)==Pi*R^4/(S1+S2+S3+S)==R^2},{{R,4},{S,S1},{a,1},{b,1}},WorkingPrecision->9]]
Block[{S1=25,S2=12,S3=20},FindRoot[{2(S1+S2)/(2a-Sin)==2(S1+S)/(2b-Sin)==(S1-S2+S3-S)/(4CosCos)==Pi*R^4/(S1+S2+S3+S)==R^2},{{R,4},{S,S1},{a,1},{b,1}},WorkingPrecision->9]]
Block[{S1=12,S2=20,S3=25},FindRoot[{2(S1+S2)/(2a-Sin)==2(S1+S)/(2b-Sin)==(S1-S2+S3-S)/(4CosCos)==Pi*R^4/(S1+S2+S3+S)==R^2},{{R,4},{S,S1},{a,1},{b,1}},WorkingPrecision->9]]
Block[{S1=25,S2=20,S3=12},FindRoot[{2(S1+S2)/(2a-Sin)==2(S1+S)/(2b-Sin)==(S1-S2+S3-S)/(4CosCos)==Pi*R^4/(S1+S2+S3+S)==R^2},{{R,4},{S,S1},{a,1},{b,1}},WorkingPrecision->9]]
Block[{S1=12,S2=25,S3=20},FindRoot[{2(S1+S2)/(2a-Sin)==2(S1+S)/(2b-Sin)==(S1-S2+S3-S)/(4CosCos)==Pi*R^4/(S1+S2+S3+S)==R^2},{{R,4},{S,S1},{a,2},{b,1}},WorkingPrecision->9]]
Block[{S1=20,S2=25,S3=12},FindRoot[{2(S1+S2)/(2a-Sin)==2(S1+S)/(2b-Sin)==(S1-S2+S3-S)/(4CosCos)==Pi*R^4/(S1+S2+S3+S)==R^2},{{R,4},{S,S1},{a,2},{b,1}},WorkingPrecision->9]]
{R -> 5.52129591, 12. -> 38.7705363, a -> 1.30401157, b -> 1.75127984}
{R -> 5.52129591, 12. -> 38.7705364, a -> 1.39031282, b -> 1.83758109}
{R -> 4.80009461, 12. -> 15.3851481, a -> 1.47956252, b -> 1.37727260}
{R -> 4.80009461, 12. -> 15.3851481, a -> 1.76432005, b -> 1.66203013}
{R -> 4.59133032, 12. -> 9.22576002, a -> 1.66325722, b -> 1.28085908}
{R -> 4.59133032, 12. -> 9.22576008, a -> 1.86073357, b -> 1.47833543}
\(\D\frac{2 (S1 + S2)}{(2 a - \sin(2 a)) R^2}=\frac{ 2 (S1 + S)}{(2 b - \sin(2 b)) R^2}=\frac{ S1 - S2 + S3 - S}{4 \cos(a) \cos(b) R^2} =\frac{\pi*R^2}{S1 + S2 + S3 + S}=1\)
对不起！浪费大家的时间，我还是说几句(不是小结)。
1圆被两条相互垂直直线分为4部分, 3部分的面积是12,20,25, 求剩下部分的面积。
往大的想，4部分有4种可能：
1，4部分是相同的。可以有唯一解。
2，3部分是相同的，无解。
3，2部分相同，另2部分相同。无解。
2部分相同，另2部分不同，？？？
4，4部分都不相同，6个图形(3个R)。
4部分为S1,S2,S3,S4, 且S1≥S2≥S3≥S4,则S2*S3≥S1*S4
圆心肯定在S1。
重点：2部分相同，另2部分不同，？？？
2部分相同(肯定是对边)，另2部分不同(对边)，
4部分只能是S1≥S2≥S2≥S4, ？(怎么确定)≥S2*S2≥S1*S4
譬如：12,20,20,25没有解？？？
12,20,20,31.6762可以有解。
王守恩 发表于 2023-8-16 16:22
{R -> 5.52129591, 12. -> 38.7705363, a -> 1.30401157, b -> 1.75127984}
{R -> 5.52129591, 12. -> 3 ...
你这个求面积不需要加绝对值吗？
他是用了格林公式？
这个还是可以塞进去的。23#41#还是塞不进去。
FindRoot[{(2*S - a*b)/(2 x - Sin) == (2*12 - a*c)/(2 y - Sin) == (2*20 - c*d)/(\ - 2 x - Sin) == (2*25 - b*d)/(\ - 2 y - Sin) == R^2,
Sqrt/Sin == Sqrt/Sin == (b + c)/Sin == 2 R,a*d == b*c}, {{S, 1}, {R, 1}, {x, 1}, {y, 1}, {a, 1}, {b, 1}, {c, 1}, {d, 1}}, WorkingPrecision -> 9]
{S -> 15.3851481, R -> 4.80009461, x -> 0.734253203, y -> 0.643019400, a -> 3.85698669, b -> 5.14781319, c -> 4.27316596, d -> 5.70327613}
本帖最后由 Jack315 于 2023-8-19 14:40 编辑
【分享成果】
计算四部分面积采用的方法为：
每一部分的面积=扇形面积-扇形内三角形的面积+直角三角形面积
在 Excel 和 Mathematica 中， ArcTan 函数返回值范围为 \(-\pi<\alpha\le\pi\) 。
这导致圆心角的计算需要进行调整，相应的公式如下：
\[∠AOB=tan^{-1}\bigg(x_p,-\sqrt{r^2-x_p^2}\bigg)-tan^{-1}\bigg(-\sqrt{r^2-y_p^2},y_p\bigg)+\begin{cases} 0 & \text{if }y_p<0 \\ 2\pi & \text{if }y_p\ge0\end{cases}\]
\[∠BOC=tan^{-1}\bigg(\sqrt{r^2-y_p^2},y_p\bigg)-tan^{-1}\bigg(x_p,-\sqrt{r^2-x_p^2}\bigg)\]
\[∠COD=tan^{-1}\bigg(x_p,\sqrt{r^2-x_p^2}\bigg)-tan^{-1}\bigg(\sqrt{r^2-y_p^2},y_p\bigg)\]
\[∠DOA=tan^{-1}\bigg(-\sqrt{r^2-y_p^2},y_p\bigg)-tan^{-1}\bigg(x_p,\sqrt{r^2-x_p^2}\bigg)+\begin{cases} 2\pi & \text{if }y_p<0 \\ 0 & \text{if }y_p\ge0\end{cases}\]
由于对顶圆心角互补，即 \(∠AOB+∠COD=\pi\) 和 \(∠BOC+∠DOA=\pi\)，
所以只求 \(∠BOC\) 和 \(∠COD\)，另外两个角利用互补关系求出，相对地简化了计算公式。
下列附件中包含有 Excel 和 Mathematica 的工作文档：
（压缩文件有三个，上传受限了，最后一个明天补上）
Excel 文档可读性比较差；Mathematica 就好多了。
在 Excel 文档中，每个工作表最下面有简短的说明。
在 Mathematica 文档中定义了求面积的公式：
计算四部分面积[半径_, x值_, y值_] :=
Module[{r = 半径, x = x值, y = y值, 临时长度, 分割线长度, 临时角度, 临时圆心角, 圆心角补角, 圆心角,
扇形面积, 扇形内三角形面积, 直角三角形面积},
(*=======================================================*)
(* 计算各段分割线的长度 *)
(*=======================================================*)
(*第一个值为Sqrt，垂直分割线长度的一半*)
(*第二个值为Sqrt，水平分割线长度的一半。*)
临时长度 = Sqrt;
(*\分割线长度\\
四个值为四段分割线的长度，依次为 PA、PB、PC 和 PD 的长度。*)
分割线长度 = {x + 临时长度[], y + 临时长度[], 临时长度[] - x, 临时长度[] - y};
(*=======================================================*)
(*计算圆心角。 *)
(*=======================================================*)
(*三个值依次为tan^-1(Sqrt,y)、tan^-1(x,-Sqrt)和tan^-1(x,
Sqrt)，是与分割线与圆交点相关的极角，用于计算圆心角的中间临时变量。*)
临时角度 = ArcTan[{临时长度[], x, x}, {y, -临时长度[], 临时长度[]}];
(*两个值依次为圆心角\BOC和\COD的角度。*)
临时圆心角 = {临时角度[] - 临时角度[], 临时角度[] - 临时角度[]};
(*两个值依次为圆心角\DOA和\AOB的角度。*)
(*说明：对顶圆心角互补。*)
圆心角补角 = \ - 临时圆心角;
(*重新构造圆心角列表。*)
(*四个值依次为圆心角\AOB、\BOC、\COD和\DOA的角度。*)
圆心角 = {圆心角补角[], 临时圆心角[], 临时圆心角[], 圆心角补角[]};
(*=======================================================*)
(*计算面积。 *)
(*=======================================================*)
(*四个值依次为扇形AOB、BOC、COD和DOA的面积。*)
扇形面积 = r^2/2 圆心角;
(*四个值依次为扇形内三角形\AOB、\BOC、\
\COD和\DOA的面积。*)
扇形内三角形面积 = r^2/2 Sin[圆心角];
(*四个值依次为直角三角形\APB、\BPC、\
\CPD和\DPA的面积。*)
直角三角形面积 =
1/2 Table[分割线长度[]*分割线长度[]], {i, 1, 4}];
(*=======================================================*)
(*计算并返回四部分面积。 *)
(*=======================================================*)
(*四个值依次为 APB, BPC,CPD 和 DPA 的面积。*)
扇形面积 - 扇形内三角形面积 + 直角三角形面积
]
使用了三个不同的命令求解：Solve/NSolve、FindRoot 和 NMinimize 。
Solve/NSolve 命令求解失败；
运行效果 FindRoot 效果要好一些，数字精度高、速度也相对快一点。
最后画出了解的图形。