数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
查看: 8732|回复: 21

【答】四棱锥藏地球之三,可能比较基础

[复制链接]
发表于 2023-5-22 07:23 | 显示全部楼层 |阅读模式
本帖最后由 dodonaomikiki 于 2023-5-30 03:02 编辑

ABCD-EFGH乃是一个单位正方体

CH,           CE,                   CF,                CG构成一个四棱锥,
现在,
在这个椎体里面塞进一个最大的球体,
求这个球体的半径

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
 楼主| 发表于 2023-5-22 08:30 | 显示全部楼层
吧其中的两个面涂色一哈
请看

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-5-23 04:56 | 显示全部楼层
剥离出来之后,
这个四棱锥,四条棱的长度来看,
倒是比较简单!
但是,
里面想建立一个最大的圆球坐标从而来计算,
感觉也甚为不易

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-5-23 04:57 | 显示全部楼层
草图只能示意用
其中两个面已经涂上黄色。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-5-24 01:04 | 显示全部楼层
制图水平非常有限,
4楼图形太差,
稍加以改进!便于解题
请见草图

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

发表于 2023-5-24 23:24 | 显示全部楼层

针对该问题可作如下分析推导:
四个平面方程为:
\[
Plane\,\,ABCD:z=0
\\
Plane\,\,ABE:\left( 0,1,0 \right) \times \left( -1,0,1 \right) \left( x-1,y,z \right) =0
\\
Plane\,\,BCE:y=1
\\
Plane\,\,CDE:x=0
\\
Plane\,\,ADE:\left( 1,0,0 \right) \times \left( 0,1,1 \right) \left( x,y,z \right) =0

\]
四棱锥内的点为满足:
\[
\left( x_0,y_0,z_0 \right) \left\{ \begin{array}{c}
        z_0>0\\
        \left( 1,0,1 \right) \left( x_0-1,y_0,z_0 \right) <0\\
        y_0<1\\
        x_0>0\\
        \left( 0,-1,1 \right) \left( x_0,y_0,z_0 \right) <0\\
\end{array} \right.

\]
锥内的点到五个平面的距离为:
\[
\left\{ \begin{array}{c}
        d_1=z_0\\
        d_2=-\left( 1,0,1 \right) \left( x_0-1,y_0,z_0 \right)\\
        d_3=1-y_0\\
        d_4=x_0\\
        d_5=-\left( 0,-1,1 \right) \left( x_0,y_0,z_0 \right)\\
\end{array} \right.

\]
求原题中的最大半径即为求\(d_i\)最小值的最大值:
\[
r=\max \left( \min d_i \right)

\]

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

发表于 2023-5-25 00:25 | 显示全部楼层
本帖最后由 NHY007 于 2023-5-25 00:28 编辑

根据楼上推导,使用数值算法可知,在原问题单位立方体中的四棱锥内,最大的内切球半径约为\(0.3317\),球心位置约为\([0.0050       \ 0.0101 \ 0.0050]\); 诚然,精度有限,并非最完美的解,只能说是接近的解。
附Matlab代码:
  1. clear; clc;
  2. N=200;
  3. x0=linspace(0,1,N);
  4. y0=x0; z0=x0;
  5. [xGrid,yGrid,zGrid]=ndgrid(x0,y0,z0);
  6. pos(:,1)=reshape(xGrid,[numel(xGrid),1]);
  7. pos(:,2)=reshape(yGrid,[numel(yGrid),1]);
  8. pos(:,3)=reshape(zGrid,[numel(zGrid),1]);
  9. d=[pos(:,3),-(pos(:,1)-1+pos(:,3)),1-pos(:,2),pos(:,1),pos(:,2)-pos(:,3)];
  10. rC=[];
  11. for index=1:5   
  12.     [rC1,~]=find(d(:,index)<=0);
  13.     rC=[rC;rC1];
  14. end
  15. rC=unique(rC);
  16. pos(rC,:)=[];
  17. d(rC,:)=[];
  18. [dMin,dMinIndex]=min(d,[],2);
  19. pos=pos(dMinIndex,:);
  20. [dMinMax,I]=max(dMin)
  21. pos(I,:)
  22. plot(1:numel(dMin),dMin,'.');
  23. figure;
  24. scatter3(pos(:,1),pos(:,2),pos(:,3))
复制代码
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-5-25 14:51 | 显示全部楼层
最大的内切球半径约为
0,3317
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-5-25 15:04 | 显示全部楼层
本帖最后由 dodonaomikiki 于 2023-5-26 22:04 编辑

圆球  与平面\(PTA\)有一个接触点,或者说切点,
那么,   她肯定有一垂足
\begin{align*}
\Longrightarrow    进而,  面(P\Omega  Q)  与面   (PTA)     与面( \Omega  TAQ )\\
之间构想出一直角三角形\\
\Longrightarrow      类似于Rt\blacktriangle     PAQ\\
\Longrightarrow    x&=\frac{PQ \bullet      AQ}{  PA+AQ+PQ   }=\frac{1 \bullet    1}{  \sqrt{2}+1+1   }=\frac{1}{2+ \sqrt{2}   }\\
&=\frac{2-  \sqrt{2}   }{    2}\\
\end{align*}

总结:  很见单!如果木有做错的话
如果四五十年不做立体几何,
再重新捡拾,
确有一定困难,
需要一定想象力

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-5-25 15:08 | 显示全部楼层
本帖最后由 dodonaomikiki 于 2023-5-25 15:10 编辑

\(=\frac{2-  \sqrt{2}   }{    2}  \approx   0.29289\)   
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2025-6-23 03:26 , Processed in 0.108510 second(s), 16 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表