知道特征函数(characteristic function)，请看附件的第一张图片，可以通过拉普拉斯反变换求得其密度函数，请见附件的第二，三张图片。

其中有7个参数，我的目的是对参数的不同赋值求出x的密度函数，用了下面的M-file

function [ f ] = Linverse( x0,S0,r0,b0,u0,v0,k0,t0 )
syms x s S r b u v k t
h = S^(-s)*exp(x*s)*exp((-r*s+0.5*(s*b)^2)+k*t*(exp(-(s*u)+0.5*(s*v)^2)-1));
y=1/(2*pi*i)*int(h,s,-i*inf,i*inf);
f = subs(subs(subs(subs(subs(subs(subs(subs(y,'x',x0),'S',S0),'r',r0),'b',b0),'u',u0),'v',v0),'k',k0),'t',t0)

运行，提示错误：
>> Linverse(110,100,0.05,0.1,0,0.1,3,1)
Warning: Explicit integral could not be found.
> In sym.int at 58
  In Linverse at 4
 
f =
 
-5734161139222659/36028797018963968*sqrt(-1)*int(100^(-s)*exp(110*s)*exp(-1/20*s+1/200*s^2+3*exp(1/200*s^2)-3),s = -sqrt(-1) .. sqrt(-1))
 
 
 
ans =
 
-5734161139222659/36028797018963968*sqrt(-1)*int(100^(-s)*exp(110*s)*exp(-1/20*s+1/200*s^2+3*exp(1/200*s^2)-3),s = -sqrt(-1) .. sqrt(-1))

请大家帮忙，谢谢！