Question

MATLAB 里的quadgk在这种函数里怎么用？

A0=1e-6;r0=2e-3;wp=3e-4;wl=2e-4;ft=0.21;length=946e-9;
f=exp(-2*r^2*t/wp^2)/t;
E是对f积分，积分变量为t，积分区间从0到Inf;
OPD=A0*(log(r0^2/r^2)-E)+r^2/2/ft;
Integrand=exp(-1i*2*pi*OPD*exp(-r^2/(m(2,i))^2)*r/length);
求Integrand积分的值，积分变量为r，积分区间为0到0.002；
望高手不吝赐教！！

Answer 1

quadgk(@(r) exp(-1i*2*pi*OPD*exp(-r^2/(m(2,i))^2)*r/length),0,0.002)

Answer 2

最简单的
function y = myfun(x)
y = exp(x).*log(x);

Q = quadgk(@myfun,0,1)