关于quadgk函数积分的具体问题，求指导
syms uj;pdfu1=feval('funu',uj,m,p,d0,d,deta0,deta,z,gk,vgkc,dgkc);
eq1=simplify(vpa(pdfu1));
var=sym2poly(eq1)
eeudown=quadgk(@(uj)inter(uj,var,p),-inf,inf);
eeuup=quadgk(@(uj)inter(uj,var,p),-inf,0);
eeuu(ss)=eeuup/eeudown; % the acceptance integerated byGuass method in each dimension
以上只给出了程序的一小部分，中间使用了matlab自带的quadgk自适应积分。 由于被积函数pdfu1是一个伪概率密度函数（需归一化），函数近似对称分布，在数据点增多后会变窄峰值变尖，需要求出他的积分值然后归一化，不知道这个积分对函数特性具体要求。 我只知道quadgk可以解决奇异积分，适应性很强。不知道用在这里合适不，求高手指点{:{39}:}
相对而言quadgk是matlab自带积分函数中最强的了，它能够适应剧烈的振荡函数
至于是否适用于你的函数，那就需要自己试试了
页:
[1]