hành vi kỳ lạ trong chức năng không thể thiếu trong MATLAB

Question

tôi tạo hai hàm fx và fy như sau:

fx = @(x) ncx2pdf(x, 10, 1); 
fy = @(y) ncx2pdf(y + 25, 10, 10); 

Sau đó, tôi xác định fs chức năng như sau:

shift = 0 
fs = @(s) fx(s) .* fy(shift - s) 

Lưu ý rằng fs luôn là dương tính (sản phẩm của hai hàm mật độ xác suất). Nếu tôi tính toán:

integral(fs, -Inf, 100) 

Tôi có được giá trị thực sự 0,0413, nhưng nếu tôi tính toán

integral(fs, -Inf, 1000) 

Tôi lấy 0. Tại sao hành vi kỳ lạ này xảy ra sử dụng chức năng không thể thiếu? Lưu ý rằng nếu tôi tính

integral(fs, -Inf, Inf) 

Tôi có giá trị thực 0.0413.

Answer 1

"Công trình bậc hai bằng cách chia tách khoảng cách xa nhau và xấp xỉ tích phân của mỗi phân đoạn. Ước tính được thực hiện bằng cách phân loại các khoảng con và sử dụng một số dạng xấp xỉ tích phân (Riemann khoản tiền ...). sẽ có một số lỗi và điều này phụ thuộc vào việc phân vùng, do đó chúng tôi là tốt nhất để giữ giới hạn của chúng ta về hội nhập gần gũi hơn với những điểm mà sẽ đóng góp nhiều nhất vào không thể thiếu "Brendan Hamm

Ref: http://www.mathworks.com/matlabcentral/answers/242910-strange-behaviour-in-integral-function-in-matlab#answer_192302