matlab未系统学过，已困惑了许久，希望得到高手、老师的帮助！积分函数如图片所示。 fQ5v?(
!i t orSl
DC h
!Z{I
\#,2#BmO"E
我采用了dblquad函数来进行积分，将待积函数设了一个function函数 Dz0D ^(;V
o ks;G([
function y=E_x(r,PHI,z) .aT@'a{F
omega0=7.57e15; _
dA y w
c=3e11; : 6Nb,Hh~
k=omega0/c; gkd4)\9
A=@(phi,theta)exp(i*k*(z.*cos(theta)+r.*sin(theta).*cos(phi-PHI))).*cos(theta).*cos(phi); ~3.*b% ,
y=dblquad(A,0,2*pi,0,1.12); cTq@"v di
,1{qZ(l1
Q` &#u#
然后直接调用该函数： MSmr7%g3D
lambda=2*pi/k; o4" [{LyT
for z=0:1:2; xS 1|t};
r=linspace(-2*lambda,2*lambda,201); )Jjp^U3Ub
HI=linspace(0,2*pi,201); 8H4NNj Oy
[rho,PHI ]=meshgrid(r,HI); _|u}^MLO
Ie=conj(E_x(rho,PHI,z)).*E_x(rho,PHI,z); ~s+vJvWz
end O z6$u
surf(rho,PHI,Ie) ,,u hEoH
'6 M6e(
结果运行错误，显示phi和PHI的维度不匹配，想问问该怎么办？我的解法思路对吗？如不对，可否告知编程方法，万分感谢！ WJD2(el
i3 eF_
