hujunhua 发表于 2014-8-27 10:01:12

看到17#方程组,我猛然觉醒,发现自己又犯錯了。

1、不是任意的四面体都存在等角中心,即使四面体的四个立体角都小于 `\pi`。从12#的方程组中消去{x,y,z,w}, 应可得到存在等角中心的四面体的六棱长约束条件。
      不幸的是,我用Mathematica10试了一下,一直Runing, 不出结果。
2、对于存在等角中心的四面体,等角中心应该既是1#所要的费马点,即到四顶点距离之和最小的点,也是6#所说的6个三角形面积之和最小的点。
3、对于不存在等角中心的四面体,6楼所说的6个三角形面积之和最小的点的几何特征需要另行确定,但从该点所张的6个三角形之结构应该不是肥皂膜实验的结果(最小曲面)。

数学星空 发表于 2014-8-27 19:52:36

关于\(w\),重新消元计算得到:


\(27(3a^2-3b_1^2+2b_1c_1-3c_1^2)^2(3a^2-3b_1^2-2b_1c_1-3c_1^2)^2+(12960a^6-40608a^4b_1^2-40608a^4c_1^2+42336a^2b_1^4+69696a^2b_1^2c_1^2+42336a^2c_1^4-14688b_1^6-31392b_1^4c_1^2-31392b_1^2c_1^4-14688c_1^6)w^2+(33408a^4-63744a^2b_1^2-63744a^2c_1^2+33408b_1^4+46848b_1^2c_1^2+33408c_1^4)w^4+(34816a^2-30720b_1^2-30720c_1^2)w^6+12288w^8=0\)

\(27(3a_1^2+2a_1c_1-3b^2+3c_1^2)^2(3a_1^2-2a_1c_1-3b^2+3c_1^2)^2+(-14688a_1^6+42336a_1^4b^2-31392a_1^4c_1^2-40608a_1^2b^4+69696a_1^2b^2c_1^2-31392a_1^2c_1^4+12960b^6-40608b^4c_1^2+42336b^2c_1^4-14688c_1^6)w^2+(33408a_1^4-63744a_1^2b^2+46848a_1^2c_1^2+33408b^4-63744b^2c_1^2+33408c_1^4)w^4+(-30720a_1^2+34816b^2-30720c_1^2)w^6+12288w^8=0\)

\(27(3a_1^2-2a_1b_1+3b_1^2-3c^2)^2(3a_1^2+2a_1b_1+3b_1^2-3c^2)^2+(-14688a_1^6-31392a_1^4b_1^2+42336a_1^4c^2-31392a_1^2b_1^4+69696a_1^2b_1^2c^2-40608a_1^2c^4-14688b_1^6+42336b_1^4c^2-40608b_1^2c^4+12960c^6)w^2+(33408a_1^4+46848a_1^2b_1^2-63744a_1^2c^2+33408b_1^4-63744b_1^2c^2+33408c^4)w^4+(-30720a_1^2-30720b_1^2+34816c^2)w^6+12288w^8=0\)



然后将上面两两消元,得到独立的约束条件(存在费马点的条件)

2187*a^16+21384*a^14*a1^2-3888*a^14*b1^2-20088*a^14*c^2-17496*a^14*c1^2+83700*a^12*a1^4-28944*a^12*a1^2*b1^2-163656*a^12*a1^2*c^2-148392*a^12*a1^2*c1^2+15552*a^12*b1^4+12528*a^12*b1^2*c^2+19440*a^12*b1^2*c1^2+83700*a^12*c^4+141912*a^12*c^2*c1^2+61236*a^12*c1^4+175032*a^10*a1^6-86448*a^10*a1^4*b1^2-527160*a^10*a1^4*c^2-498456*a^10*a1^4*c1^2+79488*a^10*a1^2*b1^4+105312*a^10*a1^2*b1^2*c^2+129888*a^10*a1^2*b1^2*c1^2+540360*a^10*a1^2*c^4+975216*a^10*a1^2*c^2*c1^2+441288*a^10*a1^2*c1^4-86400*a^10*b1^4*c^2-86400*a^10*b1^4*c1^2-14256*a^10*b1^2*c^4-36000*a^10*b1^2*c^2*c1^2-34992*a^10*b1^2*c1^4-189000*a^10*c^6-498456*a^10*c^4*c1^2-429624*a^10*c^2*c1^4-122472*a^10*c1^6+221826*a^8*a1^8-142864*a^8*a1^6*b1^2-893256*a^8*a1^6*c^2-877224*a^8*a1^6*c1^2+181056*a^8*a1^4*b1^4+297232*a^8*a1^4*b1^2*c^2+325648*a^8*a1^4*b1^2*c1^2+1368204*a^8*a1^4*c^4+2612520*a^8*a1^4*c^2*c1^2+1240524*a^8*a1^4*c1^4-379776*a^8*a1^2*b1^4*c^2-379776*a^8*a1^2*b1^4*c1^2-164592*a^8*a1^2*b1^2*c^4-312352*a^8*a1^2*b1^2*c^2*c1^2-217584*a^8*a1^2*b1^2*c1^4-942408*a^8*a1^2*c^6-2651352*a^8*a1^2*c^4*c1^2-2427960*a^8*a1^2*c^2*c1^4-729000*a^8*a1^2*c1^6+205632*a^8*b1^4*c^4+400512*a^8*b1^4*c^2*c1^2+205632*a^8*b1^4*c1^4+5616*a^8*b1^2*c^6-27120*a^8*b1^2*c^4*c1^2-10224*a^8*b1^2*c^2*c1^4+19440*a^8*b1^2*c1^6+246402*a^8*c^8+918360*a^8*c^6*c1^2+1240524*a^8*c^4*c1^4+722520*a^8*c^2*c1^6+153090*a^8*c1^8+175032*a^6*a1^10-142864*a^6*a1^8*b1^2-877224*a^6*a1^8*c^2-893256*a^6*a1^8*c1^2+234240*a^6*a1^6*b1^4+420672*a^6*a1^6*b1^2*c^2+420672*a^6*a1^6*b1^2*c1^2+1769712*a^6*a1^6*c^4+3550752*a^6*a1^6*c^2*c1^2+1769712*a^6*a1^6*c1^4-713472*a^6*a1^4*b1^4*c^2-713472*a^6*a1^4*b1^4*c1^2-415328*a^6*a1^4*b1^2*c^4-716608*a^6*a1^4*b1^2*c^2*c1^2-447584*a^6*a1^4*b1^2*c1^4-1796880*a^6*a1^4*c^6-5367216*a^6*a1^4*c^4*c1^2-5214192*a^6*a1^4*c^2*c1^4-1651536*a^6*a1^4*c1^6+748800*a^6*a1^2*b1^4*c^4+1448448*a^6*a1^2*b1^4*c^2*c1^2+748800*a^6*a1^2*b1^4*c1^4+131904*a^6*a1^2*b1^2*c^6+182464*a^6*a1^2*b1^2*c^4*c1^2+207040*a^6*a1^2*b1^2*c^2*c1^4+150336*a^6*a1^2*b1^2*c1^6+918360*a^6*a1^2*c^8+3649056*a^6*a1^2*c^6*c1^2+5239824*a^6*a1^2*c^4*c1^4+3232800*a^6*a1^2*c^2*c1^6+722520*a^6*a1^2*c1^8-269568*a^6*b1^4*c^6-762624*a^6*b1^4*c^4*c1^2-762624*a^6*b1^4*c^2*c1^4-269568*a^6*b1^4*c1^6+5616*a^6*b1^2*c^8+131904*a^6*b1^2*c^6*c1^2+249760*a^6*b1^2*c^4*c1^4+150336*a^6*b1^2*c^2*c1^6+19440*a^6*b1^2*c1^8-189000*a^6*c^10-942408*a^6*c^8*c1^2-1796880*a^6*c^6*c1^4-1651536*a^6*c^4*c1^6-729000*a^6*c^2*c1^8-122472*a^6*c1^10+83700*a^4*a1^12-86448*a^4*a1^10*b1^2-498456*a^4*a1^10*c^2-527160*a^4*a1^10*c1^2+181056*a^4*a1^8*b1^4+325648*a^4*a1^8*b1^2*c^2+297232*a^4*a1^8*b1^2*c1^2+1240524*a^4*a1^8*c^4+2612520*a^4*a1^8*c^2*c1^2+1368204*a^4*a1^8*c1^4-713472*a^4*a1^6*b1^4*c^2-713472*a^4*a1^6*b1^4*c1^2-447584*a^4*a1^6*b1^2*c^4-716608*a^4*a1^6*b1^2*c^2*c1^2-415328*a^4*a1^6*b1^2*c1^4-1651536*a^4*a1^6*c^6-5214192*a^4*a1^6*c^4*c1^2-5367216*a^4*a1^6*c^2*c1^4-1796880*a^4*a1^6*c1^6+1089408*a^4*a1^4*b1^4*c^4+2102016*a^4*a1^4*b1^4*c^2*c1^2+1089408*a^4*a1^4*b1^4*c1^4+249760*a^4*a1^4*b1^2*c^6+401504*a^4*a1^4*b1^2*c^4*c1^2+401504*a^4*a1^4*b1^2*c^2*c1^4+249760*a^4*a1^4*b1^2*c1^6+1240524*a^4*a1^4*c^8+5239824*a^4*a1^4*c^6*c1^2+8002632*a^4*a1^4*c^4*c1^4+5239824*a^4*a1^4*c^2*c1^6+1240524*a^4*a1^4*c1^8-762624*a^4*a1^2*b1^4*c^6-2137344*a^4*a1^2*b1^4*c^4*c1^2-2137344*a^4*a1^2*b1^4*c^2*c1^4-762624*a^4*a1^2*b1^4*c1^6-27120*a^4*a1^2*b1^2*c^8+182464*a^4*a1^2*b1^2*c^6*c1^2+401504*a^4*a1^2*b1^2*c^4*c1^4+207040*a^4*a1^2*b1^2*c^2*c1^6-10224*a^4*a1^2*b1^2*c1^8-498456*a^4*a1^2*c^10-2651352*a^4*a1^2*c^8*c1^2-5367216*a^4*a1^2*c^6*c1^4-5214192*a^4*a1^2*c^4*c1^6-2427960*a^4*a1^2*c^2*c1^8-429624*a^4*a1^2*c1^10+205632*a^4*b1^4*c^8+748800*a^4*b1^4*c^6*c1^2+1089408*a^4*b1^4*c^4*c1^4+748800*a^4*b1^4*c^2*c1^6+205632*a^4*b1^4*c1^8-14256*a^4*b1^2*c^10-164592*a^4*b1^2*c^8*c1^2-415328*a^4*b1^2*c^6*c1^4-447584*a^4*b1^2*c^4*c1^6-217584*a^4*b1^2*c^2*c1^8-34992*a^4*b1^2*c1^10+83700*a^4*c^12+540360*a^4*c^10*c1^2+1368204*a^4*c^8*c1^4+1769712*a^4*c^6*c1^6+1240524*a^4*c^4*c1^8+441288*a^4*c^2*c1^10+61236*a^4*c1^12+21384*a^2*a1^14-28944*a^2*a1^12*b1^2-148392*a^2*a1^12*c^2-163656*a^2*a1^12*c1^2+79488*a^2*a1^10*b1^4+129888*a^2*a1^10*b1^2*c^2+105312*a^2*a1^10*b1^2*c1^2+441288*a^2*a1^10*c^4+975216*a^2*a1^10*c^2*c1^2+540360*a^2*a1^10*c1^4-379776*a^2*a1^8*b1^4*c^2-379776*a^2*a1^8*b1^4*c1^2-217584*a^2*a1^8*b1^2*c^4-312352*a^2*a1^8*b1^2*c^2*c1^2-164592*a^2*a1^8*b1^2*c1^4-729000*a^2*a1^8*c^6-2427960*a^2*a1^8*c^4*c1^2-2651352*a^2*a1^8*c^2*c1^4-942408*a^2*a1^8*c1^6+748800*a^2*a1^6*b1^4*c^4+1448448*a^2*a1^6*b1^4*c^2*c1^2+748800*a^2*a1^6*b1^4*c1^4+150336*a^2*a1^6*b1^2*c^6+207040*a^2*a1^6*b1^2*c^4*c1^2+182464*a^2*a1^6*b1^2*c^2*c1^4+131904*a^2*a1^6*b1^2*c1^6+722520*a^2*a1^6*c^8+3232800*a^2*a1^6*c^6*c1^2+5239824*a^2*a1^6*c^4*c1^4+3649056*a^2*a1^6*c^2*c1^6+918360*a^2*a1^6*c1^8-762624*a^2*a1^4*b1^4*c^6-2137344*a^2*a1^4*b1^4*c^4*c1^2-2137344*a^2*a1^4*b1^4*c^2*c1^4-762624*a^2*a1^4*b1^4*c1^6-10224*a^2*a1^4*b1^2*c^8+207040*a^2*a1^4*b1^2*c^6*c1^2+401504*a^2*a1^4*b1^2*c^4*c1^4+182464*a^2*a1^4*b1^2*c^2*c1^6-27120*a^2*a1^4*b1^2*c1^8-429624*a^2*a1^4*c^10-2427960*a^2*a1^4*c^8*c1^2-5214192*a^2*a1^4*c^6*c1^4-5367216*a^2*a1^4*c^4*c1^6-2651352*a^2*a1^4*c^2*c1^8-498456*a^2*a1^4*c1^10+400512*a^2*a1^2*b1^4*c^8+1448448*a^2*a1^2*b1^4*c^6*c1^2+2102016*a^2*a1^2*b1^4*c^4*c1^4+1448448*a^2*a1^2*b1^4*c^2*c1^6+400512*a^2*a1^2*b1^4*c1^8-36000*a^2*a1^2*b1^2*c^10-312352*a^2*a1^2*b1^2*c^8*c1^2-716608*a^2*a1^2*b1^2*c^6*c1^4-716608*a^2*a1^2*b1^2*c^4*c1^6-312352*a^2*a1^2*b1^2*c^2*c1^8-36000*a^2*a1^2*b1^2*c1^10+141912*a^2*a1^2*c^12+975216*a^2*a1^2*c^10*c1^2+2612520*a^2*a1^2*c^8*c1^4+3550752*a^2*a1^2*c^6*c1^6+2612520*a^2*a1^2*c^4*c1^8+975216*a^2*a1^2*c^2*c1^10+141912*a^2*a1^2*c1^12-86400*a^2*b1^4*c^10-379776*a^2*b1^4*c^8*c1^2-713472*a^2*b1^4*c^6*c1^4-713472*a^2*b1^4*c^4*c1^6-379776*a^2*b1^4*c^2*c1^8-86400*a^2*b1^4*c1^10+12528*a^2*b1^2*c^12+105312*a^2*b1^2*c^10*c1^2+297232*a^2*b1^2*c^8*c1^4+420672*a^2*b1^2*c^6*c1^6+325648*a^2*b1^2*c^4*c1^8+129888*a^2*b1^2*c^2*c1^10+19440*a^2*b1^2*c1^12-20088*a^2*c^14-163656*a^2*c^12*c1^2-527160*a^2*c^10*c1^4-893256*a^2*c^8*c1^6-877224*a^2*c^6*c1^8-498456*a^2*c^4*c1^10-148392*a^2*c^2*c1^12-17496*a^2*c1^14+2187*a1^16-3888*a1^14*b1^2-17496*a1^14*c^2-20088*a1^14*c1^2+15552*a1^12*b1^4+19440*a1^12*b1^2*c^2+12528*a1^12*b1^2*c1^2+61236*a1^12*c^4+141912*a1^12*c^2*c1^2+83700*a1^12*c1^4-86400*a1^10*b1^4*c^2-86400*a1^10*b1^4*c1^2-34992*a1^10*b1^2*c^4-36000*a1^10*b1^2*c^2*c1^2-14256*a1^10*b1^2*c1^4-122472*a1^10*c^6-429624*a1^10*c^4*c1^2-498456*a1^10*c^2*c1^4-189000*a1^10*c1^6+205632*a1^8*b1^4*c^4+400512*a1^8*b1^4*c^2*c1^2+205632*a1^8*b1^4*c1^4+19440*a1^8*b1^2*c^6-10224*a1^8*b1^2*c^4*c1^2-27120*a1^8*b1^2*c^2*c1^4+5616*a1^8*b1^2*c1^6+153090*a1^8*c^8+722520*a1^8*c^6*c1^2+1240524*a1^8*c^4*c1^4+918360*a1^8*c^2*c1^6+246402*a1^8*c1^8-269568*a1^6*b1^4*c^6-762624*a1^6*b1^4*c^4*c1^2-762624*a1^6*b1^4*c^2*c1^4-269568*a1^6*b1^4*c1^6+19440*a1^6*b1^2*c^8+150336*a1^6*b1^2*c^6*c1^2+249760*a1^6*b1^2*c^4*c1^4+131904*a1^6*b1^2*c^2*c1^6+5616*a1^6*b1^2*c1^8-122472*a1^6*c^10-729000*a1^6*c^8*c1^2-1651536*a1^6*c^6*c1^4-1796880*a1^6*c^4*c1^6-942408*a1^6*c^2*c1^8-189000*a1^6*c1^10+205632*a1^4*b1^4*c^8+748800*a1^4*b1^4*c^6*c1^2+1089408*a1^4*b1^4*c^4*c1^4+748800*a1^4*b1^4*c^2*c1^6+205632*a1^4*b1^4*c1^8-34992*a1^4*b1^2*c^10-217584*a1^4*b1^2*c^8*c1^2-447584*a1^4*b1^2*c^6*c1^4-415328*a1^4*b1^2*c^4*c1^6-164592*a1^4*b1^2*c^2*c1^8-14256*a1^4*b1^2*c1^10+61236*a1^4*c^12+441288*a1^4*c^10*c1^2+1240524*a1^4*c^8*c1^4+1769712*a1^4*c^6*c1^6+1368204*a1^4*c^4*c1^8+540360*a1^4*c^2*c1^10+83700*a1^4*c1^12-86400*a1^2*b1^4*c^10-379776*a1^2*b1^4*c^8*c1^2-713472*a1^2*b1^4*c^6*c1^4-713472*a1^2*b1^4*c^4*c1^6-379776*a1^2*b1^4*c^2*c1^8-86400*a1^2*b1^4*c1^10+19440*a1^2*b1^2*c^12+129888*a1^2*b1^2*c^10*c1^2+325648*a1^2*b1^2*c^8*c1^4+420672*a1^2*b1^2*c^6*c1^6+297232*a1^2*b1^2*c^4*c1^8+105312*a1^2*b1^2*c^2*c1^10+12528*a1^2*b1^2*c1^12-17496*a1^2*c^14-148392*a1^2*c^12*c1^2-498456*a1^2*c^10*c1^4-877224*a1^2*c^8*c1^6-893256*a1^2*c^6*c1^8-527160*a1^2*c^4*c1^10-163656*a1^2*c^2*c1^12-20088*a1^2*c1^14+15552*b1^4*c^12+79488*b1^4*c^10*c1^2+181056*b1^4*c^8*c1^4+234240*b1^4*c^6*c1^6+181056*b1^4*c^4*c1^8+79488*b1^4*c^2*c1^10+15552*b1^4*c1^12-3888*b1^2*c^14-28944*b1^2*c^12*c1^2-86448*b1^2*c^10*c1^4-142864*b1^2*c^8*c1^6-142864*b1^2*c^6*c1^8-86448*b1^2*c^4*c1^10-28944*b1^2*c^2*c1^12-3888*b1^2*c1^14+2187*c^16+21384*c^14*c1^2+83700*c^12*c1^4+175032*c^10*c1^6+221826*c^8*c1^8+175032*c^6*c1^10+83700*c^4*c1^12+21384*c^2*c1^14+2187*c1^16=0



2187*a^16+21384*a^14*a1^2-20088*a^14*b^2-17496*a^14*b1^2-3888*a^14*c1^2+83700*a^12*a1^4-163656*a^12*a1^2*b^2-148392*a^12*a1^2*b1^2-28944*a^12*a1^2*c1^2+83700*a^12*b^4+141912*a^12*b^2*b1^2+12528*a^12*b^2*c1^2+61236*a^12*b1^4+19440*a^12*b1^2*c1^2+15552*a^12*c1^4+175032*a^10*a1^6-527160*a^10*a1^4*b^2-498456*a^10*a1^4*b1^2-86448*a^10*a1^4*c1^2+540360*a^10*a1^2*b^4+975216*a^10*a1^2*b^2*b1^2+105312*a^10*a1^2*b^2*c1^2+441288*a^10*a1^2*b1^4+129888*a^10*a1^2*b1^2*c1^2+79488*a^10*a1^2*c1^4-189000*a^10*b^6-498456*a^10*b^4*b1^2-14256*a^10*b^4*c1^2-429624*a^10*b^2*b1^4-36000*a^10*b^2*b1^2*c1^2-86400*a^10*b^2*c1^4-122472*a^10*b1^6-34992*a^10*b1^4*c1^2-86400*a^10*b1^2*c1^4+221826*a^8*a1^8-893256*a^8*a1^6*b^2-877224*a^8*a1^6*b1^2-142864*a^8*a1^6*c1^2+1368204*a^8*a1^4*b^4+2612520*a^8*a1^4*b^2*b1^2+297232*a^8*a1^4*b^2*c1^2+1240524*a^8*a1^4*b1^4+325648*a^8*a1^4*b1^2*c1^2+181056*a^8*a1^4*c1^4-942408*a^8*a1^2*b^6-2651352*a^8*a1^2*b^4*b1^2-164592*a^8*a1^2*b^4*c1^2-2427960*a^8*a1^2*b^2*b1^4-312352*a^8*a1^2*b^2*b1^2*c1^2-379776*a^8*a1^2*b^2*c1^4-729000*a^8*a1^2*b1^6-217584*a^8*a1^2*b1^4*c1^2-379776*a^8*a1^2*b1^2*c1^4+246402*a^8*b^8+918360*a^8*b^6*b1^2+5616*a^8*b^6*c1^2+1240524*a^8*b^4*b1^4-27120*a^8*b^4*b1^2*c1^2+205632*a^8*b^4*c1^4+722520*a^8*b^2*b1^6-10224*a^8*b^2*b1^4*c1^2+400512*a^8*b^2*b1^2*c1^4+153090*a^8*b1^8+19440*a^8*b1^6*c1^2+205632*a^8*b1^4*c1^4+175032*a^6*a1^10-877224*a^6*a1^8*b^2-893256*a^6*a1^8*b1^2-142864*a^6*a1^8*c1^2+1769712*a^6*a1^6*b^4+3550752*a^6*a1^6*b^2*b1^2+420672*a^6*a1^6*b^2*c1^2+1769712*a^6*a1^6*b1^4+420672*a^6*a1^6*b1^2*c1^2+234240*a^6*a1^6*c1^4-1796880*a^6*a1^4*b^6-5367216*a^6*a1^4*b^4*b1^2-415328*a^6*a1^4*b^4*c1^2-5214192*a^6*a1^4*b^2*b1^4-716608*a^6*a1^4*b^2*b1^2*c1^2-713472*a^6*a1^4*b^2*c1^4-1651536*a^6*a1^4*b1^6-447584*a^6*a1^4*b1^4*c1^2-713472*a^6*a1^4*b1^2*c1^4+918360*a^6*a1^2*b^8+3649056*a^6*a1^2*b^6*b1^2+131904*a^6*a1^2*b^6*c1^2+5239824*a^6*a1^2*b^4*b1^4+182464*a^6*a1^2*b^4*b1^2*c1^2+748800*a^6*a1^2*b^4*c1^4+3232800*a^6*a1^2*b^2*b1^6+207040*a^6*a1^2*b^2*b1^4*c1^2+1448448*a^6*a1^2*b^2*b1^2*c1^4+722520*a^6*a1^2*b1^8+150336*a^6*a1^2*b1^6*c1^2+748800*a^6*a1^2*b1^4*c1^4-189000*a^6*b^10-942408*a^6*b^8*b1^2+5616*a^6*b^8*c1^2-1796880*a^6*b^6*b1^4+131904*a^6*b^6*b1^2*c1^2-269568*a^6*b^6*c1^4-1651536*a^6*b^4*b1^6+249760*a^6*b^4*b1^4*c1^2-762624*a^6*b^4*b1^2*c1^4-729000*a^6*b^2*b1^8+150336*a^6*b^2*b1^6*c1^2-762624*a^6*b^2*b1^4*c1^4-122472*a^6*b1^10+19440*a^6*b1^8*c1^2-269568*a^6*b1^6*c1^4+83700*a^4*a1^12-498456*a^4*a1^10*b^2-527160*a^4*a1^10*b1^2-86448*a^4*a1^10*c1^2+1240524*a^4*a1^8*b^4+2612520*a^4*a1^8*b^2*b1^2+325648*a^4*a1^8*b^2*c1^2+1368204*a^4*a1^8*b1^4+297232*a^4*a1^8*b1^2*c1^2+181056*a^4*a1^8*c1^4-1651536*a^4*a1^6*b^6-5214192*a^4*a1^6*b^4*b1^2-447584*a^4*a1^6*b^4*c1^2-5367216*a^4*a1^6*b^2*b1^4-716608*a^4*a1^6*b^2*b1^2*c1^2-713472*a^4*a1^6*b^2*c1^4-1796880*a^4*a1^6*b1^6-415328*a^4*a1^6*b1^4*c1^2-713472*a^4*a1^6*b1^2*c1^4+1240524*a^4*a1^4*b^8+5239824*a^4*a1^4*b^6*b1^2+249760*a^4*a1^4*b^6*c1^2+8002632*a^4*a1^4*b^4*b1^4+401504*a^4*a1^4*b^4*b1^2*c1^2+1089408*a^4*a1^4*b^4*c1^4+5239824*a^4*a1^4*b^2*b1^6+401504*a^4*a1^4*b^2*b1^4*c1^2+2102016*a^4*a1^4*b^2*b1^2*c1^4+1240524*a^4*a1^4*b1^8+249760*a^4*a1^4*b1^6*c1^2+1089408*a^4*a1^4*b1^4*c1^4-498456*a^4*a1^2*b^10-2651352*a^4*a1^2*b^8*b1^2-27120*a^4*a1^2*b^8*c1^2-5367216*a^4*a1^2*b^6*b1^4+182464*a^4*a1^2*b^6*b1^2*c1^2-762624*a^4*a1^2*b^6*c1^4-5214192*a^4*a1^2*b^4*b1^6+401504*a^4*a1^2*b^4*b1^4*c1^2-2137344*a^4*a1^2*b^4*b1^2*c1^4-2427960*a^4*a1^2*b^2*b1^8+207040*a^4*a1^2*b^2*b1^6*c1^2-2137344*a^4*a1^2*b^2*b1^4*c1^4-429624*a^4*a1^2*b1^10-10224*a^4*a1^2*b1^8*c1^2-762624*a^4*a1^2*b1^6*c1^4+83700*a^4*b^12+540360*a^4*b^10*b1^2-14256*a^4*b^10*c1^2+1368204*a^4*b^8*b1^4-164592*a^4*b^8*b1^2*c1^2+205632*a^4*b^8*c1^4+1769712*a^4*b^6*b1^6-415328*a^4*b^6*b1^4*c1^2+748800*a^4*b^6*b1^2*c1^4+1240524*a^4*b^4*b1^8-447584*a^4*b^4*b1^6*c1^2+1089408*a^4*b^4*b1^4*c1^4+441288*a^4*b^2*b1^10-217584*a^4*b^2*b1^8*c1^2+748800*a^4*b^2*b1^6*c1^4+61236*a^4*b1^12-34992*a^4*b1^10*c1^2+205632*a^4*b1^8*c1^4+21384*a^2*a1^14-148392*a^2*a1^12*b^2-163656*a^2*a1^12*b1^2-28944*a^2*a1^12*c1^2+441288*a^2*a1^10*b^4+975216*a^2*a1^10*b^2*b1^2+129888*a^2*a1^10*b^2*c1^2+540360*a^2*a1^10*b1^4+105312*a^2*a1^10*b1^2*c1^2+79488*a^2*a1^10*c1^4-729000*a^2*a1^8*b^6-2427960*a^2*a1^8*b^4*b1^2-217584*a^2*a1^8*b^4*c1^2-2651352*a^2*a1^8*b^2*b1^4-312352*a^2*a1^8*b^2*b1^2*c1^2-379776*a^2*a1^8*b^2*c1^4-942408*a^2*a1^8*b1^6-164592*a^2*a1^8*b1^4*c1^2-379776*a^2*a1^8*b1^2*c1^4+722520*a^2*a1^6*b^8+3232800*a^2*a1^6*b^6*b1^2+150336*a^2*a1^6*b^6*c1^2+5239824*a^2*a1^6*b^4*b1^4+207040*a^2*a1^6*b^4*b1^2*c1^2+748800*a^2*a1^6*b^4*c1^4+3649056*a^2*a1^6*b^2*b1^6+182464*a^2*a1^6*b^2*b1^4*c1^2+1448448*a^2*a1^6*b^2*b1^2*c1^4+918360*a^2*a1^6*b1^8+131904*a^2*a1^6*b1^6*c1^2+748800*a^2*a1^6*b1^4*c1^4-429624*a^2*a1^4*b^10-2427960*a^2*a1^4*b^8*b1^2-10224*a^2*a1^4*b^8*c1^2-5214192*a^2*a1^4*b^6*b1^4+207040*a^2*a1^4*b^6*b1^2*c1^2-762624*a^2*a1^4*b^6*c1^4-5367216*a^2*a1^4*b^4*b1^6+401504*a^2*a1^4*b^4*b1^4*c1^2-2137344*a^2*a1^4*b^4*b1^2*c1^4-2651352*a^2*a1^4*b^2*b1^8+182464*a^2*a1^4*b^2*b1^6*c1^2-2137344*a^2*a1^4*b^2*b1^4*c1^4-498456*a^2*a1^4*b1^10-27120*a^2*a1^4*b1^8*c1^2-762624*a^2*a1^4*b1^6*c1^4+141912*a^2*a1^2*b^12+975216*a^2*a1^2*b^10*b1^2-36000*a^2*a1^2*b^10*c1^2+2612520*a^2*a1^2*b^8*b1^4-312352*a^2*a1^2*b^8*b1^2*c1^2+400512*a^2*a1^2*b^8*c1^4+3550752*a^2*a1^2*b^6*b1^6-716608*a^2*a1^2*b^6*b1^4*c1^2+1448448*a^2*a1^2*b^6*b1^2*c1^4+2612520*a^2*a1^2*b^4*b1^8-716608*a^2*a1^2*b^4*b1^6*c1^2+2102016*a^2*a1^2*b^4*b1^4*c1^4+975216*a^2*a1^2*b^2*b1^10-312352*a^2*a1^2*b^2*b1^8*c1^2+1448448*a^2*a1^2*b^2*b1^6*c1^4+141912*a^2*a1^2*b1^12-36000*a^2*a1^2*b1^10*c1^2+400512*a^2*a1^2*b1^8*c1^4-20088*a^2*b^14-163656*a^2*b^12*b1^2+12528*a^2*b^12*c1^2-527160*a^2*b^10*b1^4+105312*a^2*b^10*b1^2*c1^2-86400*a^2*b^10*c1^4-893256*a^2*b^8*b1^6+297232*a^2*b^8*b1^4*c1^2-379776*a^2*b^8*b1^2*c1^4-877224*a^2*b^6*b1^8+420672*a^2*b^6*b1^6*c1^2-713472*a^2*b^6*b1^4*c1^4-498456*a^2*b^4*b1^10+325648*a^2*b^4*b1^8*c1^2-713472*a^2*b^4*b1^6*c1^4-148392*a^2*b^2*b1^12+129888*a^2*b^2*b1^10*c1^2-379776*a^2*b^2*b1^8*c1^4-17496*a^2*b1^14+19440*a^2*b1^12*c1^2-86400*a^2*b1^10*c1^4+2187*a1^16-17496*a1^14*b^2-20088*a1^14*b1^2-3888*a1^14*c1^2+61236*a1^12*b^4+141912*a1^12*b^2*b1^2+19440*a1^12*b^2*c1^2+83700*a1^12*b1^4+12528*a1^12*b1^2*c1^2+15552*a1^12*c1^4-122472*a1^10*b^6-429624*a1^10*b^4*b1^2-34992*a1^10*b^4*c1^2-498456*a1^10*b^2*b1^4-36000*a1^10*b^2*b1^2*c1^2-86400*a1^10*b^2*c1^4-189000*a1^10*b1^6-14256*a1^10*b1^4*c1^2-86400*a1^10*b1^2*c1^4+153090*a1^8*b^8+722520*a1^8*b^6*b1^2+19440*a1^8*b^6*c1^2+1240524*a1^8*b^4*b1^4-10224*a1^8*b^4*b1^2*c1^2+205632*a1^8*b^4*c1^4+918360*a1^8*b^2*b1^6-27120*a1^8*b^2*b1^4*c1^2+400512*a1^8*b^2*b1^2*c1^4+246402*a1^8*b1^8+5616*a1^8*b1^6*c1^2+205632*a1^8*b1^4*c1^4-122472*a1^6*b^10-729000*a1^6*b^8*b1^2+19440*a1^6*b^8*c1^2-1651536*a1^6*b^6*b1^4+150336*a1^6*b^6*b1^2*c1^2-269568*a1^6*b^6*c1^4-1796880*a1^6*b^4*b1^6+249760*a1^6*b^4*b1^4*c1^2-762624*a1^6*b^4*b1^2*c1^4-942408*a1^6*b^2*b1^8+131904*a1^6*b^2*b1^6*c1^2-762624*a1^6*b^2*b1^4*c1^4-189000*a1^6*b1^10+5616*a1^6*b1^8*c1^2-269568*a1^6*b1^6*c1^4+61236*a1^4*b^12+441288*a1^4*b^10*b1^2-34992*a1^4*b^10*c1^2+1240524*a1^4*b^8*b1^4-217584*a1^4*b^8*b1^2*c1^2+205632*a1^4*b^8*c1^4+1769712*a1^4*b^6*b1^6-447584*a1^4*b^6*b1^4*c1^2+748800*a1^4*b^6*b1^2*c1^4+1368204*a1^4*b^4*b1^8-415328*a1^4*b^4*b1^6*c1^2+1089408*a1^4*b^4*b1^4*c1^4+540360*a1^4*b^2*b1^10-164592*a1^4*b^2*b1^8*c1^2+748800*a1^4*b^2*b1^6*c1^4+83700*a1^4*b1^12-14256*a1^4*b1^10*c1^2+205632*a1^4*b1^8*c1^4-17496*a1^2*b^14-148392*a1^2*b^12*b1^2+19440*a1^2*b^12*c1^2-498456*a1^2*b^10*b1^4+129888*a1^2*b^10*b1^2*c1^2-86400*a1^2*b^10*c1^4-877224*a1^2*b^8*b1^6+325648*a1^2*b^8*b1^4*c1^2-379776*a1^2*b^8*b1^2*c1^4-893256*a1^2*b^6*b1^8+420672*a1^2*b^6*b1^6*c1^2-713472*a1^2*b^6*b1^4*c1^4-527160*a1^2*b^4*b1^10+297232*a1^2*b^4*b1^8*c1^2-713472*a1^2*b^4*b1^6*c1^4-163656*a1^2*b^2*b1^12+105312*a1^2*b^2*b1^10*c1^2-379776*a1^2*b^2*b1^8*c1^4-20088*a1^2*b1^14+12528*a1^2*b1^12*c1^2-86400*a1^2*b1^10*c1^4+2187*b^16+21384*b^14*b1^2-3888*b^14*c1^2+83700*b^12*b1^4-28944*b^12*b1^2*c1^2+15552*b^12*c1^4+175032*b^10*b1^6-86448*b^10*b1^4*c1^2+79488*b^10*b1^2*c1^4+221826*b^8*b1^8-142864*b^8*b1^6*c1^2+181056*b^8*b1^4*c1^4+175032*b^6*b1^10-142864*b^6*b1^8*c1^2+234240*b^6*b1^6*c1^4+83700*b^4*b1^12-86448*b^4*b1^10*c1^2+181056*b^4*b1^8*c1^4+21384*b^2*b1^14-28944*b^2*b1^12*c1^2+79488*b^2*b1^10*c1^4+2187*b1^16-3888*b1^14*c1^2+15552*b1^12*c1^4=0

数学星空 发表于 2014-8-27 20:01:47

更难以接受的是:

将\(\{a = 1, b = x, c = y, a1 = 1, b1 = x, c1 = y\}\) 代入上面两个约束式:

\(262144(x-1)^6(x+1)^6(3x^4-2x^2y^2+3y^4-6x^2-2y^2+3)=0\)

\(262144(y-1)^6(y+1)^6(3x^4-2x^2y^2+3y^4-2x^2-6y^2+3)=0\)

其中\(3x^4-2x^2y^2+3y^4-2x^2-6y^2+3=0\)的函数曲线为:



我们需要确认的是上面两个条件约束式是否正确?能否进一步简化?

数学星空 发表于 2014-8-27 20:17:44

当然我们也可以利用下面的体积公式:

\

对四个四面体\(F-ABC,F-ABD,F-BCD,F-ACD\),使用上面的体积公式可以得到:

\

\[12V=\sqrt
{-a^4a_1^2-a^2a_1^4+a^2a_1^2b^2+a^2a_1^2b_1^2+a^2a_1^2c^2+a^2a_1^2c_1^2+a^2b^2b_1^2-a^2b^2c^2-a^2b_1^2c_1^2+a^2c^2c_1^2+a_1^2b^2b_1^2-a_1^2b^2c_1^2-a_1^2b_1^2c^2+a_1^2c^2c_1^2-b^4b_1^2-b^2b_1^4+b^2b_1^2c^2+b^2b_1^2c_1^2+b^2c^2c_1^2+b_1^2c^2c_1^2-c^4c_1^2-c^2c_1^4}\]


参与上面的消元计算,或许得到的结论要简洁一些?

数学星空 发表于 2014-8-28 09:57:22

对于楼上得到的两个约束条件:

1.对于正三棱锥,若令\(\{a = x, b = x, c = x, a_1 = y, b_1 = y, c_1 = y\}\),则恒为\(0\)


2.对于等四面体,若令\(\{a = x, b = y, c = z, a_1 = x, b_1 = y, c_1 = z\}\),则简化为:

    \(262144(x-y)^6(x+y)^6(3x^4-6x^2y^2-2x^2z^2+3y^4-2y^2z^2+3z^4)=0\)

    \(262144(x-z)^6(x+z)^6(3x^4-2x^2y^2-6x^2z^2+3y^4-2y^2z^2+3z^4)=0\)


3.对于垂心四面体,若令\(a = \sqrt{m^2+n^2}, a_1 =\sqrt{s^2+t^2}, b =\sqrt{m^2+s^2}, b_1 =\sqrt{n^2+t^2}, c = \sqrt{m^2+t^2}, c_1 =\sqrt{n^2+s^2}\),则简化为:

    \(768(n-s)^4(n+s)^4(m-t)^4(m+t)^4=0\)

    \(768(n-t)^4(n+t)^4(m-s)^4(m+s)^4=0\)

hujunhua 发表于 2014-8-28 15:58:52

从17#的方程组中任选5个消去{x,y,z,w}, 都可以得到一个5棱长约束,共计6个五棱长约束。
这6个五棱长约束不是或的关系,而是与的关系。

hujunhua 发表于 2014-8-29 02:18:56

计算了一下,不存在非正四面体的等面四面体具有等角中心,垂心四面体具有等角中心的只有正三棱锥。

倪举鹏 发表于 2014-8-29 20:02:50

找两点,这两点所在的直线会通过四面体两个对棱,这线段与,这点与两个顶点连线段都成120度

数学星空 发表于 2015-8-10 19:43:50

其实2#~3#已给出了求解方程:

对于四面体\(D-ABC\)的费马点\(F\),有\[\angle BFC=\angle AFD=\alpha,\angle AFC=\angle BFD=\beta,\angle AFB=\angle CFD=\gamma\]并且每对相等的角被两者所在平面的交线平分。即:
\[\left\{\begin{split}\frac{y^2+z^2-a^2}{2yz}=\frac{x^2+w^2-a_1^2}{2xw}=\cos(\alpha)\\ \frac{x^2+z^2-b^2}{2xz}=\frac{y^2+w^2-b_1^2}{2yw}=\cos(\beta)\\ \frac{x^2+y^2-c^2}{2xy}=\frac{z^2+w^2-c_1^2}{2zw}=\cos(\gamma)\end{split}\right.\]

\[\begin{vmatrix}
0 & c^2 & b^2 & a_1^2 & x^2 & 1\\
c^2 &0 &a^2 & b_1^2 & y^2 & 1\\
b^2 &a^2 & 0 &c_1^2 &z^2 & 1\\
a_1^2 &b_1^2 &c_1^2 &0 & w^2 & 1\\
x^2 &y^2 &z^2 &w^2 &0 & 1\\
1 & 1 &1 & 1 & 1&0
\end{vmatrix}=0\]

通过四个方程及\(a,b,c,a_1,b_1,c_1\) 可能求解四个未知数\(x,y,z,w\)

数学星空 发表于 2015-8-10 19:58:17

对于特殊四面体:\(a=a_1,b=b_1,c=c_1\) 费马点即为四面体的外心.

对于一般的四面体,我可以求解楼上的方程算出费马点

例如: 取\(a=3,b=4,c=5,a_1=6,b_1=6,c_1=6\)

为了方便输入:记\(\cos(\alpha)=k_1,\cos(\beta)=k_2,\cos(\gamma)=k_3\)

\(-a^4a_1^4+2a^4a_1^2w^2+2a^4a_1^2x^2-a^4w^4+2a^4w^2x^2-a^4x^4+2a^2a_1^4y^2+2a^2a_1^4z^2+2a^2a_1^2b^2b_1^2-2a^2a_1^2b^2w^2-2a^2a_1^2b^2y^2-2a^2a_1^2b_1^2x^2-2a^2a_1^2b_1^2z^2+2a^2a_1^2c^2c_1^2-2a^2a_1^2c^2w^2-2a^2a_1^2c^2z^2-2a^2a_1^2c_1^2x^2-2a^2a_1^2c_1^2y^2+4a^2a_1^2w^2x^2-2a^2a_1^2w^2y^2-2a^2a_1^2w^2z^2-2a^2a_1^2x^2y^2-2a^2a_1^2x^2z^2+4a^2a_1^2y^2z^2-2a^2b^2b_1^2w^2-2a^2b^2b_1^2x^2+4a^2b^2c^2w^2+2a^2b^2w^4-2a^2b^2w^2x^2-2a^2b^2w^2y^2+2a^2b^2x^2y^2+4a^2b_1^2c_1^2x^2-2a^2b_1^2w^2x^2+2a^2b_1^2w^2z^2+2a^2b_1^2x^4-2a^2b_1^2x^2z^2-2a^2c^2c_1^2w^2-2a^2c^2c_1^2x^2+2a^2c^2w^4-2a^2c^2w^2x^2-2a^2c^2w^2z^2+2a^2c^2x^2z^2-2a^2c_1^2w^2x^2+2a^2c_1^2w^2y^2+2a^2c_1^2x^4-2a^2c_1^2x^2y^2-a_1^4y^4+2a_1^4y^2z^2-a_1^4z^4-2a_1^2b^2b_1^2y^2-2a_1^2b^2b_1^2z^2+4a_1^2b^2c_1^2y^2-2a_1^2b^2w^2y^2+2a_1^2b^2w^2z^2+2a_1^2b^2y^4-2a_1^2b^2y^2z^2+4a_1^2b_1^2c^2z^2+2a_1^2b_1^2x^2y^2-2a_1^2b_1^2x^2z^2-2a_1^2b_1^2y^2z^2+2a_1^2b_1^2z^4-2a_1^2c^2c_1^2y^2-2a_1^2c^2c_1^2z^2+2a_1^2c^2w^2y^2-2a_1^2c^2w^2z^2-2a_1^2c^2y^2z^2+2a_1^2c^2z^4-2a_1^2c_1^2x^2y^2+2a_1^2c_1^2x^2z^2+2a_1^2c_1^2y^4-2a_1^2c_1^2y^2z^2-b^4b_1^4+2b^4b_1^2w^2+2b^4b_1^2y^2-b^4w^4+2b^4w^2y^2-b^4y^4+2b^2b_1^4x^2+2b^2b_1^4z^2+2b^2b_1^2c^2c_1^2-2b^2b_1^2c^2w^2-2b^2b_1^2c^2z^2-2b^2b_1^2c_1^2x^2-2b^2b_1^2c_1^2y^2-2b^2b_1^2w^2x^2+4b^2b_1^2w^2y^2-2b^2b_1^2w^2z^2-2b^2b_1^2x^2y^2+4b^2b_1^2x^2z^2-2b^2b_1^2y^2z^2-2b^2c^2c_1^2w^2-2b^2c^2c_1^2y^2+2b^2c^2w^4-2b^2c^2w^2y^2-2b^2c^2w^2z^2+2b^2c^2y^2z^2+2b^2c_1^2w^2x^2-2b^2c_1^2w^2y^2-2b^2c_1^2x^2y^2+2b^2c_1^2y^4-b_1^4x^4+2b_1^4x^2z^2-b_1^4z^4-2b_1^2c^2c_1^2x^2-2b_1^2c^2c_1^2z^2+2b_1^2c^2w^2x^2-2b_1^2c^2w^2z^2-2b_1^2c^2x^2z^2+2b_1^2c^2z^4+2b_1^2c_1^2x^4-2b_1^2c_1^2x^2y^2-2b_1^2c_1^2x^2z^2+2b_1^2c_1^2y^2z^2-c^4c_1^4+2c^4c_1^2w^2+2c^4c_1^2z^2-c^4w^4+2c^4w^2z^2-c^4z^4+2c^2c_1^4x^2+2c^2c_1^4y^2-2c^2c_1^2w^2x^2-2c^2c_1^2w^2y^2+4c^2c_1^2w^2z^2+4c^2c_1^2x^2y^2-2c^2c_1^2x^2z^2-2c^2c_1^2y^2z^2-c_1^4x^4+2c_1^4x^2y^2-c_1^4y^4=0\)

\(-2k_1yz-a^2+y^2+z^2=0\)
\(-2k_1wx-a_1^2+w^2+x^2=0\)
\(-2k_2xz-b^2+x^2+z^2=0\)
\(-2k_2wy-b_1^2+w^2+y^2=0\)
\(-2k_3xy-c_1^2+x^2+y^2=0\)
\(-2k_3wz-c_1^2+w^2+z^2=0\)

代入得到:

\(576w^4-576w^2x^2-576w^2y^2+1215x^4+288x^2y^2-2142x^2z^2+2048y^4-3808y^2z^2+2975z^4-27072w^2-20736x^2-20736y^2+746496=0\)

\(-2k_1yz+y^2+z^2-9=0\)
\(-2k_1wx+w^2+x^2-36=0\)
\(-2k_2xz+x^2+z^2-16=0\)
\(-2k_2wy+w^2+y^2-36=0\)
\(-2k_3xy+x^2+y^2-25=0\)
\(-2k_3wz+w^2+z^2-36=0\)

求解得到:

\(k_1 = -0.08125763422, k_2 = -0.3028690928, k_3 = -0.6158732730, w = 4.838906457, x = 3.176057999, y = 2.372782500, z = 1.653020392,\alpha = 94.660858238723315798^{\circ}, \beta = 107.63000925585309530^{\circ},\delta = 128.01540128977782951^{\circ}\)
页: 1 2 [3] 4 5
查看完整版本: 四面体中费马点计算