我们做A题的主要思路
本帖最后由 tuotuo64 于 2010-9-15 07:17 编辑我们做的是A题，理解题后，感觉最难的就是根据实验数据估计出两个变位参数。由于在求罐中所剩油的体积的时候，就发现用matlab符号积分求不出来，也就是积分没有解析表达式。所以只好用近似法来处理。通过研究所给数据，发现每组数据取油后油面下降值都在1CM左右，所以就将所取油层体积近似为（油面上表面积+油面下表面积）/2*油层厚度，这样就将两个变位参数与每组实验数据建立了显示表达式关系，为进一步以取油体积与近似体积差的平方和为最小值为目标函数而建立非线性规划问题来求最优变位参数提供了可能，最后我们得到了最优解（也就是最小二乘法），a=2.0568,b=1.4843。在准确性的处理中，我们通过体积差的相对误差来检验，发现误差有一定的规律，也就是罐中含有三根管道所占的体积没有减掉，通过回归找到了其中规律，然后修正了我们的结果，发现误差会进一步减小，最终的相对误差的均值大概为0.1%。在稳定性分析中，我们采用了蒙特卡洛法，我们首先认为标尺的读数可有误差，故我们对所给数据中标高读数进行了随机扰动，也就是加一个均值为0，方差为某常数的正态分布的随机项，然后重新计算两个变位参数，发现变位参数的相对误差和加了随机振动后的标高产生的相对误差处在同一个数量级上，故从数值实验的角度我们认为我们的方法是比较稳定的。
顺便说一下，对无解析解的积分模型，用matlab一样可以用最小二乘法求出最优变位参数，当然这样的算法很有技巧性，一般的人就想不到了。
做B题的飘过{:3_42:}{:3_42:}{:3_42:}{:3_42:}{:3_42:} 感觉没有把问题重点搞清楚 论文写得有点乱 高手~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 终于看完了~~~ 不错。。。。。。。。。 好像看上去方法不错啊 牛人啊。。。。。。。。。。。。。。。。。。。。。 以用最小二乘法求出最优变位参数，当然这样的算法很有技巧性，一般的人就想不到了。! u7 l+ i# ?" K4 q; M7 L# |
以用最小二乘法求出最优变位参数，当然这样的算法很有技巧性，一般的人就想不到了