的系統的數值解得到含有術語被積的數值積分我想要計算上倍頻程（或Matlab）下面的積分： 
MATLAB/OCTAVE-從ODE
然而，我沒有H和K的明確表達式。H和K實際上是下面的微分方程的數值解。 h，k，dh/dr和dk/dr的初始條件分別爲0,1,1/2,0。 
我該如何去做這件事？我能否以某種方式直接求解它，還是需要首先找到H和K的數值解，找到多項式近似值，然後進行積分？我對Matlab和數值方法完全陌生，因此描述所有內容的詳細方法將非常值得讚賞。
如果您有能力在任何需要的點解決H和K，那麼最好的選擇是在特定的點上評估它們，稱爲「節點」，然後使用數值積分方案計算積分。基本思想是，根據整合方案，您可以得到不同準確度的結果。
一個基本的有限和將簡單地評估你的H和K在等距點，然後把dx作爲他們之間的空間。這顯然可能與你擁有的無限積分有關，但你也可以將內核重鑄成一個有限域（取1 /（x + 1）），然後你的被積函數從（0，+ inf）變成1,0） - 顯然你需要取得所得積分的負值並將（1,0）變爲（0,1））。
使用有限域，您可以使用所描述的有限和方法，或者像Gauss-Legendre Quadrature這樣的有限域上的更精確的積分方案。 http://en.wikipedia.org/wiki/Gaussian_quadrature
最後，如果你需要保持無限的領域，可以使用正交型的，更適合於無限的領域，如高斯 - 拉蓋爾正交，這將需要您精通於預乘你的內核（X）所以得到的權重（Exp（-x））不會改變你的內核。 http://en.wikipedia.org/wiki/Gauss–Laguerre_quadrature
解決微分方程需要初始條件或邊界條件，你在這裏什麼都不放。你的問題很糟糕。把第一個完整的正確的數學問題 – 2012-07-13 14:50:41
@RobertH：編輯公式，並添加所需的詳細信息..希望這有助於..請讓我知道任何其他需要更改/添加。 – 2012-07-13 15:03:20
當r = 0時你認爲會發生什麼？只要看看你的方程式並考慮一下。什麼是lambda值？當符號沒有數字時，你不能做數字。發佈完整和正確的問題。 – 2012-07-13 15:50:15