又见面了
syms x y x2 y2 w c b
k=w/c;
r=sqrt(x*x+y*y);
f=exp(i*k*(y2*y2-2*y*y2)/(2*r))
yy2=int(f,y2,-b/2,b/2)
f =
exp(1/2*i*w/c*(y2^2-2*y*y2)/(x^2+y^2)^(1/2))
yy2 =
-1/2*pi^(1/2)*exp(-1/2*i*w/c*y^2/(x^2+y^2)^(1/2))*2^(1/2)*(erf(1/4*i*2^(1/2)*w*(b+2*y)/c/(x^2+y^2)^(1/2)/(-i*w/c/(x^2+y^2)^(1/2))^(1/2))-erf(1/4*i*2^(1/2)*w*(-b+2*y)/c/(x^2+y^2)^(1/2)/(-i*w/c/(x^2+y^2)^(1/2))^(1/2)))/(-i*w/c/(x^2+y^2)^(1/2))^(1/2)
望有帮助！！！