擅长:python、mysql、java
<p>我不知道我是对还是错,但我确实找到了解决问题的办法。在</p>
<p>我只是打印了一些错误信息(spectrum_[I]和spectrum_j]的值),我发现:</p>
<p>有时候,他们是一样的!!!在</p>
<p>(也许他们不一样,但他们太接近了,我猜)</p>
<p>所以,这里</p>
<pre><code>pa += log((spectrum[i] - spectrum[j]) *
(1. / spectrum_[j] - 1. / spectrum_[i])) + log(n_samples)
</code></pre>
<p>计算日志(0)时将报告错误。在</p>
<p>我的解决方法是把一个很小的数字1e-99加到0上,这样它就变成log(0+1e-99)</p>
<p><strong>这样您就可以将其更改为:</strong></p>
^{pr2}$