计算强度分布的重心？

如何解决计算强度分布的重心？

我需要计算不对称峰的重心（强度与 x 位置）以及重心和质心之间的差异。对于我的计算，强度分布在以下数组中定义为：

array([ 2.19216735e-03,5.37355085e-04,-1.92786170e-03,-2.09720374e-04,-9.89245753e-04,3.66267317e-04,-1.08720120e-03,1.20200267e-03,1.34468683e-03,1.78588924e-03,1.49047884e-03,1.45289688e-03,-6.33576590e-04,-2.41855592e-04,-1.19240721e-03,2.47065937e-03,1.32185682e-03,1.23684367e-03,2.79990578e-04,1.03529893e-03,8.06834832e-04,8.89627813e-04,1.25917290e-03,1.78030882e-03,4.58962019e-03,2.67760495e-03,3.56203662e-03,2.54182602e-03,3.01015755e-03,1.88405979e-03,3.23527409e-05,-1.80395840e-04,1.31965758e-03,2.52311848e-03,2.69531078e-03,3.65124994e-03,3.21821725e-03,2.49144138e-03,1.00411924e-03,1.56405899e-03,1.37577960e-03,2.10396579e-03,3.15646988e-03,4.07088600e-03,5.46858630e-03,5.79245046e-03,6.65671882e-03,6.07872868e-03,7.95604246e-03,7.72638126e-03,8.18776570e-03,8.60977927e-03,8.84623384e-03,1.02519189e-02,1.02395041e-02,1.11971163e-02,1.12794920e-02,1.12670772e-02,1.21298989e-02,1.04925732e-02,1.15002320e-02,1.29607069e-02,1.46529857e-02,1.74378303e-02,2.13411330e-02,2.55848872e-02,2.74912091e-02,2.97945205e-02,3.04850263e-02,3.17121532e-02,3.43222379e-02,4.07888353e-02,4.39888358e-02,4.72722961e-02,4.99033147e-02,5.59951615e-02,6.52301195e-02,7.34484280e-02,8.20102681e-02,9.47014286e-02,1.13322540e-01,1.28020171e-01,1.47961210e-01,1.73977573e-01,2.20493987e-01,2.73113269e-01,3.31269375e-01,4.14123288e-01,4.95523055e-01,5.95621845e-01,6.84023064e-01,7.72765915e-01,8.66805026e-01,9.11172707e-01,9.30003575e-01,8.89357020e-01,8.04278855e-01,7.10927100e-01,5.99161080e-01,4.95505374e-01,4.03529525e-01,3.35915158e-01,2.82565678e-01,2.43239366e-01,2.12654801e-01,1.88088056e-01,1.67616190e-01,1.45475399e-01,1.26043389e-01,1.15869230e-01,1.05873084e-01,9.74679727e-02,8.87839791e-02,8.00927619e-02,7.29480454e-02,6.79095983e-02,6.44004537e-02,6.10575821e-02,5.49680803e-02,4.94086634e-02,4.54010625e-02,4.03070823e-02,3.71372946e-02,3.68341707e-02,3.45954005e-02,3.30989865e-02,2.92657441e-02,2.93273273e-02,2.77839588e-02,2.70974001e-02,2.24385772e-02,2.01186012e-02,1.92663279e-02,2.03701782e-02,1.80438028e-02,1.69958362e-02,1.77981371e-02,1.78357219e-02,1.79447455e-02,1.86661039e-02,1.90994701e-02,1.85390958e-02,1.92580469e-02,1.91666761e-02,1.78699674e-02,1.79132612e-02,1.76786787e-02,1.61427637e-02,1.64029383e-02,1.47868025e-02,1.36899439e-02,1.15434731e-02,1.09822217e-02,1.07784945e-02,1.08616875e-02,1.08926492e-02,1.22423917e-02,1.26381268e-02,1.18751024e-02,9.23257225e-03,9.31994413e-03,7.47262664e-03,9.28244502e-03,8.10599351e-03,9.08884165e-03,7.83636436e-03,8.57565483e-03,8.22371193e-03,5.87119049e-03,2.98790183e-03,1.40723515e-03,2.20492006e-03,3.64546037e-03,2.00096642e-03,1.45168130e-03,-8.77304855e-04,-4.91304065e-04,-1.41686139e-03,-5.24865922e-04,1.61466035e-03,1.39263532e-03,2.01518514e-03,-4.95367384e-04,-1.34921281e-03,-1.49753983e-03,-3.00597421e-04,8.12046158e-04,2.60630483e-04,8.37085358e-04])

具有相同强度数组维度的 x 数组定义为：

array([-0.01418039,-0.01402859,-0.01387682,-0.01372506,-0.01357333,-0.01342162,-0.01326992,-0.01311825,-0.01296659,-0.01281496,-0.01266335,-0.01251175,-0.01236018,-0.01220863,-0.0120571,-0.01190558,-0.01175409,-0.01160262,-0.01145117,-0.01129974,-0.01114832,-0.01099693,-0.01084556,-0.01069421,-0.01054288,-0.01039157,-0.01024028,-0.01008901,-0.00993776,-0.00978653,-0.00963532,-0.00948413,-0.00933296,-0.00918181,-0.00903069,-0.00887958,-0.00872849,-0.00857742,-0.00842637,-0.00827535,-0.00812434,-0.00797335,-0.00782238,-0.00767144,-0.00752051,-0.0073696,-0.00721872,-0.00706785,-0.00691701,-0.00676618,-0.00661538,-0.00646459,-0.00631383,-0.00616308,-0.00601236,-0.00586165,-0.00571097,-0.00556031,-0.00540966,-0.00525904,-0.00510844,-0.00495786,-0.00480729,-0.00465675,-0.00450623,-0.00435573,-0.00420525,-0.00405479,-0.00390435,-0.00375393,-0.00360352,-0.00345315,-0.00330279,-0.00315245,-0.00300213,-0.00285183,-0.00270155,-0.00255129,-0.00240105,-0.00225083,-0.00210064,-0.00195046,-0.0018003,-0.00165017,-0.00150005,-0.00134995,-0.00119988,-0.00104982,-0.00089979,-0.00074977,-0.00059978,-0.0004498,-0.00029985,-0.00014991,0.,0.00014989,0.00029977,0.00044962,0.00059945,0.00074926,0.00089906,0.00104883,0.00119858,0.00134831,0.00149802,0.00164771,0.00179738,0.00194703,0.00209666,0.00224627,0.00239586,0.00254543,0.00269498,0.0028445,0.00299401,0.0031435,0.00329297,0.00344241,0.00359184,0.00374125,0.00389063,0.00404,0.00418934,0.00433867,0.00448797,0.00463726,0.00478652,0.00493576,0.00508499,0.00523419,0.00538337,0.00553253,0.00568168,0.0058308,0.0059799,0.00612898,0.00627804,0.00642708,0.0065761,0.0067251,0.00687408,0.00702304,0.00717197,0.00732089,0.00746979,0.00761867,0.00776753,0.00791636,0.00806518,0.00821397,0.00836275,0.00851151,0.00866024,0.00880895,0.00895765,0.00910632,0.00925498,0.00940361,0.00955222,0.00970081,0.00984939,0.00999794,0.01014647,0.01029498,0.01044347,0.01059194,0.01074039,0.01088882,0.01103723,0.01118562,0.01133399,0.01148233,0.01163066,0.01177897,0.01192725,0.01207552,0.01222377,0.01237199,0.0125202,0.01266838,0.01281655,0.01296469,0.01311281,0.01326092,0.013409,0.01355706,0.0137051,0.01385312,0.01400113])

我使用了以下代码：np.sum(x*intensity)/np.sum(intensity)。你认为正确吗？图中是我得到的结果。

cgx = np.sum(x*intensity)/np.sum(intensity) 
plt.scatter(x,intensity)
plt.scatter(cgx,color='k',marker='|',s=1e6)
plt.gca().set_yticks([])  
plt.title('1 Dimensional Center of Gravity')

查看结果：enter image description here