关于quadgk函数积分的具体问题，求指导

mayuanzhuo

       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 by  Guass method in each dimension

以上只给出了程序的一小部分，中间使用了matlab自带的quadgk自适应积分。 由于被积函数pdfu1是一个伪概率密度函数（需归一化），函数近似对称分布，在数据点增多后会变窄峰值变尖，需要求出他的积分值然后归一化，不知道这个积分对函数特性具体要求。 我只知道quadgk可以解决奇异积分，适应性很强。不知道用在这里合适不，求高手指点

happy

相对而言quadgk是matlab自带积分函数中最强的了，它能够适应剧烈的振荡函数
至于是否适用于你的函数，那就需要自己试试了