详细晾晒A题思路

fx0618

回复 artoriuslei 的帖子
! D1 M  }) ^! a- l  ?

    这个我们也想过了，不过因为图3明确给出内壁的数据，故认为说不通。

superDi

我也是A题用得解析式，同样的符号运算式子，MATLAB10用数分钟，07只要几秒，没有MAPLE，MATLAB的符号运算不行啊。。。
6 v$ K" I% C7 Y; o整体思路可以说和LZ很类似，但是没你分析的清晰，特别是第二问先求alpha=0时罐容表然后对比显示油降 这个没写，因为已经超过20页了。。
$ f# W. R- H. S- P- }  j7 n4 x
. T+ W1 ]: L, E$ r! W! Q残差平方和最小我们求出来是520，但参数却是悲剧的 2.07  1.91
悲剧
2 u: j: o# m# e) V0 F  ?3 y悲剧
悲剧
) H" N1 d2 S0 Q  x4 q& v5 M悲
悲

superDi

还有个问题想请教下。。我之前用的也是数值积分形式，但是求解速度特别慢，大概给30个高度算30个体积要个十几秒，更别说带600个高度做迭代了。。。楼主怎么解决的。。？

artoriuslei

回复 fx0618 的帖子
- C3 u* f: s3 R
图三是问题二的实际油罐呀，所以你发现了没有嘛，如果是问题二的话，体积通过理论计算是和实际没有偏差的。但是图四就没有给出内壁数据了，所以才有偏差呀。
( d$ v: j1 `0 ?: S9 M7 T   

artoriuslei

回复 superDi 的帖子
$ i& m2 v4 L# O" ^' @
9 @& w$ r$ F6 ?, P2 Q" ]我们在处理数值积分搜索a，b的时候，先将积分精度降低了（其实对a，b影响不大），也就是在步长取5厘米的情况下把a的范围锁定在了2到3之间，b的范围锁定在了1到5.5之间，然后再换高一点的精度（1厘米）把a的范围锁定在了2.1左右，b的范围锁定在了4.4左右。这样再用更高的精度来搜索a和b，最终得到了我们的答案。这个样子用两台电脑一起跑的话，应该一天的时间是可以跑出来的。也就是说，我们在10号的时候把一二问的模型都建好了，11号的时候算出了问题一的答案，并且把论文的框架弄好，模型写了进去，12号的时候完成了问题二的解答，再把答案填了进去，论文就写好了。
   

superDi

回复 artoriuslei 的帖子
: \- n/ ?- P3 w6 x恩，搜索方法用的跟你们一样。。
2 k! i% U& j7 a我们是花了好长时间把5个分段积分手工积出来（MATLAB用INT出来的有问题，出虚数），然后把解析式作成函数，然后用跟你们差不多的网格搜索做的，倒是很快可以算出来，唉可能积的有点问题。。。beta差别得太厉害了。。。
! @0 }2 W. N0 P" B6 f估计悲剧了
) @4 V  h8 P( u: [6 Y! H/ c悲剧了
( \' n6 \5 Y! S& s悲剧
了

artoriuslei

回复 superDi 的帖子
6 y5 @, ~) X% H3 c/ F4 sa、b可能会有一定的影响，但是我觉得这次还是整体的思路会更重要吧~
1 h7 ?" }& k2 Q: `: W9 m0 U加油~肯定没问题的~
大家都是辛苦熬出来的~肯定都没问题~
) P7 F+ W2 E4 ^  Q: C( C- f- J
: [2 @! c8 ?5 S0 h   

Tabris

回复 artoriuslei 的帖子
我们思路基本相同   A=2.40 B=3.80 不知道能不能接受了

% ?6 B9 f8 J" k/ B  e
" Y0 b: [4 T% R/ B7 U   

pku007

真是由于楼主的高水平和耐心
" J! y; T, U" V! A你的帖子讨论的最热烈
我现在理解你的a b 怎么写进表达式了
你的b 除了高度那一点点其他根本不用
所以不影响积分
% Y  Y& P9 Y( X' V尽管我还是觉得这个时候b转了a也会变
7 P, Q/ d2 r( B' H& d' j1 T6 L但我还是觉得你好强 比我们那个时候厉害多了

胡种子

大家A题 算出的 角度 和答案接近吗？