matlab未系统学过，已困惑了许久，希望得到高手、老师的帮助！积分函数如图片所示。 <Pf4[q&wM
Q"!GdKM
`.~*pT*u
Wj)v,v2&
我采用了dblquad函数来进行积分，将待积函数设了一个function函数 Te3 ?z
c[3x>f0
function y=E_x(r,PHI,z) [Ak 0kH >
omega0=7.57e15; ,1t|QvO
c=3e11; b!/-9 {
k=omega0/c; 3C ;;z
A=@(phi,theta)exp(i*k*(z.*cos(theta)+r.*sin(theta).*cos(phi-PHI))).*cos(theta).*cos(phi); IW}Wt{ 'm
y=dblquad(A,0,2*pi,0,1.12); [tC=P&<
cl{mRt0
Q 4L7{^[X
然后直接调用该函数： Q7zpu/5?
lambda=2*pi/k; ym` 4v5w
for z=0:1:2; qx 0F*EH|
r=linspace(-2*lambda,2*lambda,201); roi,?B_8
HI=linspace(0,2*pi,201); 6#:V3 ;
[rho,PHI ]=meshgrid(r,HI); t3v_o4`&
Ie=conj(E_x(rho,PHI,z)).*E_x(rho,PHI,z); gL@]p
end 4qd( a)NdY
surf(rho,PHI,Ie) LF{8hC[
!4 z vkJO
结果运行错误，显示phi和PHI的维度不匹配，想问问该怎么办？我的解法思路对吗？如不对，可否告知编程方法，万分感谢！ 8wkt9:
zlkW-rRk R
