Matlab中的数值积分

Question

我想在Matlab中对以下等式进行数值积分。

为简单起见，我们假设只有等式的第二部分存在。

这是我到目前为止所尝试的：

for lambda=280:2500  
  numerator(lambda) = trapz(Einc(lambda:lambda+1)*St(lambda));
  denominator(lambda) = trapz(Eref(lambda:lambda+1)*St(lambda));
end
Fraction= nansum(numerator)/ nansum(denominator);

St未在整个280-2500范围内定义，这就是我最后使用nansum的原因。但总的来说，我觉得我所做的事情是错误的，并且应该有一个更加正常和正确的计算方法，但我已经尝试使用谷歌搜索并在Mathworks论坛上问过但是没有人回答我，即使这应该这是一个很常见的事情......有人可以帮我写一下这个公式吗？