Integration of two functions, one of them is piecewise function
1 次查看(过去 30 天)
显示 更早的评论
Consider the function
p(µ, σ, t, x) = exp(-m*t)/x*sqrt(2*pi*(q.^2)*t)*exp(-1/(2*(q.^2)*t)*(log(x)-(m-(q.^2)/2)*t).^2) (0 < x < ∞)
Here let µ = .02, σ = .18 and t = .25, so p is a function on the variable x only. I created a function file pfun.m for this function, having µ, σ, t and x as input.
I need to compute int(C(x)*p(x),x,0,inf)
where C(x) = 0 if x < 1.02
x − 1.02 otherwise
0 个评论
回答(1 个)
Walter Roberson
2015-9-29
编辑:Walter Roberson
2015-9-29
general_p = @(m, q, t, x) exp(-m.*t) ./ x .*sqrt(2*pi*(q.^2).*t) .* exp(-1/(2*(q.^2).*t) .* (log(x)-(m-(q.^2)/2) .*t ).^2);
p = @(x) general_p(0.02, 0.18, 0.25, x);
C = @(x) (x - 1.02) .* (1.02 <= x);
result = integral(@(x) C(x) .* p(x), 0, inf)
3 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Assumptions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!