Matlab数值积分问题一个函数E（theta,phi）vl1=[0,cos(pi/4),sin(pi/4)]vl2=[0,-cos(pi/4),sin(pi/4)]vtheta=[cos(theta)*cos(phi),cos(theta)*sin(phi),-sin(theta)]vphi=[-sin(phi),cos(phi),0]vr=[sin(theta)*cos(phi),sin(th
网友回答
用符号计算求出E的表达式,再对E进行二重积分!
integrndE3 = @(theta,phi) sin(theta).*(abs(cos(phi).*(sqrt(2.0).*((exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(-1.0./2.0))+-6.123233995736766e-17-1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)+(exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(-1.0./2.0))+-6.123233995736766e-17+1i)./(sqrt(2.0).*cos(conj(theta)).*(-1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)).*(1.0./2.0)+sqrt(2.0).*((exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(1.0./2.0))+-6.123233995736766e-17-1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)-1.0)-(exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(1.0./2.0))+-6.123233995736766e-17+1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)).*(1.0./2.0))).^2+abs(sin(theta).*(sqrt(2.0).*((exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(-1.0./2.0))+-6.123233995736766e-17-1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)+(exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(-1.0./2.0))+-6.123233995736766e-17+1i)./(sqrt(2.0).*cos(conj(theta)).*(-1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)).*(1.0./2.0)-sqrt(2.0).*((exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(1.0./2.0))+-6.123233995736766e-17-1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)-1.0)-(exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(1.0./2.0))+-6.123233995736766e-17+1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)).*(1.0./2.0))+cos(theta).*sin(phi).*(sqrt(2.0).*((exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(-1.0./2.0))+-6.123233995736766e-17-1i)./(sqrt(2.0).*cos(conj(theta)).*(1.0./2.0)-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0)+(exp(pi.*(sqrt(2.0).*cos(conj(theta)).*5.0e-1i-sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*5.0e-1i).*(-1.0./2.0))+-6.123233995736766e-17+1i)./(sqrt(2.0).*cos(conj(theta)).*(-1.0./2.0)+sqrt(2.0).*sin(conj(phi)).*sin(conj(theta)).*(1.0./2.0)+1.0))