Integration of two functions, one of them is piecewise function
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentare
Antworten (1)
Walter Roberson
am 29 Sep. 2015
Bearbeitet: Walter Roberson
am 29 Sep. 2015
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 Kommentare
Siehe auch
Kategorien
Mehr zu Assumptions finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!