详细晾晒A题思路

ypy_solo

回复 artoriuslei 的帖子
你说得有道理，我在做搜索法的时候计算机计算的特别的慢，所以用最小二乘的时候，那个全局最优化就没有接出来，直接带了几组数据搜索的，所以出现较大误差。你们数值解是怎么实现的？是用的求积公式么？求出来的体积和高度角度的表达式有多长？
还有，就是第二问中的体积推导，你们积分上下限是不是有点问题？液面在没有到达球心的时候，上下限不是从罐体左端到右端吧，应该是球缺和油面直线的交点。
最后，我想说，你们这个结果真的很精确，我表示佩服。我觉得你们论文写得好的话，应该国奖没问题了~呵呵。

+ E( E- y7 h7 O# m8 W- l   

ypy_solo

回复 pku007 的帖子
& I( f- V' B4 o9 s( s( D  ], L6 }. v6 E我认为应该精确解，因为这道题是要解决实际问题的，如果是随机模拟解出来的话，是不是说服力就差一点。加油站要的是精确结果，解决他们的问题。

: r' |* T, a. {$ I6 @$ ]   

cdrcb

分析还是多对的！

phoenix007

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

pku007

答案出来了
1 m9 c; Y: V+ B: u  ~' w6 j. r/ w, {) E我们都没对
8 ~$ W3 w) ~$ ~+ u3 k% |* ]2 S最好的是微分
! e. y$ e2 @) W. O: y间接的说你的精确解是求不出的
* \( c6 f# Z9 P  Q. B2 D" E说我的蒙特卡罗 误差大不适用
: \6 T$ d$ Y& p4 w9 T& y' l
* k7 \" E/ R1 k" N7 G微分就是不直接计算 V 而计算 dV
- n& |! i2 F0 t3 O, p. K就是截面积 S
' G- E# Z  x( a2 i: Q; F# Z; f
7 N0 e0 ?: s: u' u我们的争论是非常好的
有位同学提出关于matlab精度的问题 非常有水平
$ _0 g  c- i6 W* z& T这个可能是问题所在 甚至评委也忽略了
' G% V* I, C! g3 l- y你说的精确解 有道理
但没有看见你的式子 我内心有点怀疑
7 S+ G* ^2 K/ q8 j4 D没什么 你已经很优秀了
最后 解答上的 想法 是个绝妙的想法
2 g0 t) z  a$ s他的特点是回避了复杂积分的运算
# R, @: j* X1 \) b! d5 [: L
# v2 }4 N8 R; g+ [) T! l0 X但我还是对蒙特卡罗 情有独钟
% A6 ?7 i7 ?/ a5 T" w* W" A正如你说的 可以先检验一下
. B" i; s+ g* c( @7 D( V% M如果误差不大 精度可以满足要求
4 ^0 F5 V  o( L5 j) D7 c8 r- b我觉得可以一试
2 s2 j6 N' R( j: W

artoriuslei

回复 ypy_solo 的帖子
首先，数值解还是利用微元法的思想取在积分上下限之间用特定步长求得函数值，再用每个函数值和步长的乘积加和，就可以了嘛。体积公式的表达式其实最终没有解析解，是由四个式子共同表示的。
第二问的积分上下限不需要用罐壁和油面直线的交点哈，只需要处理成没有油的地方横截面积为0就可以了嘛（实际方法就是当一个x对应的油面直线的y小于对应-r的时候，横截面积就取为0嘛），这样就免去了分类讨论的痛苦了呀。

   

artoriuslei

回复 pku007 的帖子
) w) q5 h. y, v$ s  b* M恩，蒙特卡洛法我也想试试编一下，主要是我们队这方面的能力还是有限，所以就没有办法在比赛的时候从这方面下手。
7 T" h5 R# ]8 R
   

artoriuslei

回复 cdrcb 的帖子
) _+ R: h' ?- v谢谢，话说，答辩的通知好久有哦？

   

1059974843

顶顶！不错！我们没有拟合出来，一个在WORD里面要写六七行的式子怎么拟也不可能！哎，悲剧了，先提前祝贺了！

李浩0123