>> u=17.6;s=3.4;
>> f = @(x)1/s/sqrt(2*pi)*exp((x+24-u).^2/2/s^2).*(x>0 & x<=u-12) + 1/s/sqrt(2*pi)*exp((x-u).^2/2/s^2).*(x>u-12 & x<=24);
>> x=0:0.5:24;
>> plot(x,f(x))
函数2
假设sigma d=1
>> clear;d=1;x=eps:0.05:24;
>> f =arrayfun(@(x)quad(@(p)exp(-(log(x*p/1.61/0.15)-d).^2/2/d^2),2,3)/x/d/sqrt(2*pi),x);
>> plot(x,f)
追问函数2