在matlab中,使用数值积分的时候,应当注意,积分边界是否会使得被积函数产生无意义的值(无穷大)。
用函数trapz举例说明。
此处,我们是求积分
∫
0
1
1
1
−
x
2
d
x
\int_0^1{1 \over \sqrt{1-x^2}}dx
∫011−x21dx
很显然,应用三角函数,我们可以算出来答案是
π
/
2
\pi/2
π/2.
但是这里得到的确实‘inf’.
因为当x=1的时候,被积函数是无穷大。因此,这里我们需要避免这种情况。
如此修改即可!