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

dodonaomikiki

1O楼  以及8楼结果不一样，希望有一个，
是正确的

dodonaomikiki

图形像一只看得更清楚些，
想方设法，
绞尽脑汁，
进行夸张制图 来一下

NHY007

受教了

NHY007

受教了

NHY007

受教了，从侧面计算确实是一个高效的算法
但对于数值算法，和算法精度有很大关系：

1. 
   
2. A=[0 0 1
   
3.     -1 0 -1
   
4.     0 -1 0
   
5.     1 0 0
   
6.     0 1 -1];
   
7. B=[0;-1;-1;0;0];
   
8. options = optimoptions('fmincon',"Algorithm","active-set","EnableFeasibilityMode",true,'PlotFcn','optimplotfval');
   
9. options.MaxFunctionEvaluations = 3.000000e+04;
   
10. options.MaxIterations = 1.000000e+04;
    
11. [pos0,fval]=fmincon(@disMin,[0.5 0.5 0.5],A,B,[],[],[0 0 0],[1 1 1],[],options)
    
12. % x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options)
    
13. function d=disMin(pos)
    
14. d=[pos(:,3),-(pos(:,1)-1+pos(:,3)),1-pos(:,2),pos(:,1),pos(:,2)-pos(:,3)];
    
15. d=-min(d);
    
16. end
    

复制代码

可求得最大半径约为 \(0.2929\)
球心坐标为：
\[
pos_0 =

    0.2929\ 0.7071\ 0.2929
\]

NHY007

结果还可由三维制图软件添加约束并测量得到
(图中立方体边长为100mm,约束每条线段垂直于一个面，线段长度相等)：

dodonaomikiki

有时候·，
数学软件，就是木有灵魂！有点像埃及木乃伊，名气很大，甚是华丽，但已经干瘪而木有灵魂

dodonaomikiki

但实际上，
反过来说，
如果一开始方法，路径正确，
即便说，
程序敲代码敲错，
木有得到正确结果，
那也很有魅力的！




大方向对啦，
是很有魅力的！
很吸引人~~~~~~至于细枝末节的计算，
慢慢改进也是很快的




如果大方向错误，
后续一切完滦

dodonaomikiki

评价

现在看来，本质上说，
题目还是太简单啦！
可以作为立体几何入门练习，可以玩耍一下



只不过，
几十年不做立体几何，
忽然捡拾，
还是感觉思绪有点繁杂的，思绪有点混乱

NHY007

dodonaomikiki 发表于 2023-5-26 22:01
有时候·，
数学软件，就是木有灵魂！有点像埃及木乃伊，名气很大，甚是华丽，但已经干瘪而木有灵魂

不太懂什么意思，能否具体些。