matlab未系统学过，已困惑了许久，希望得到高手、老师的帮助！积分函数如图片所示。 (;0]V+-
}z6H xB]$
;+\;^nS3d
p*G_$"KpP
我采用了dblquad函数来进行积分，将待积函数设了一个function函数 -!" 8j"pA:
7O j9~3o4
function y=E_x(r,PHI,z) oUH\SW8?
omega0=7.57e15; ?*HlAVDcFT
c=3e11; /on p<u
k=omega0/c; 0td ;Ag
A=@(phi,theta)exp(i*k*(z.*cos(theta)+r.*sin(theta).*cos(phi-PHI))).*cos(theta).*cos(phi); qW9|&GuZ$
y=dblquad(A,0,2*pi,0,1.12); .x&>H
)yJjJ:re
) PtaX|U
然后直接调用该函数： 1JY4 E2Q
lambda=2*pi/k; 1'iRx,
for z=0:1:2; &yuerNK
r=linspace(-2*lambda,2*lambda,201); mDhU wZH
HI=linspace(0,2*pi,201); x< " e
[rho,PHI ]=meshgrid(r,HI); .(krB% N
Ie=conj(E_x(rho,PHI,z)).*E_x(rho,PHI,z);
FM5$83Q
end q~Ud>{
surf(rho,PHI,Ie) =0A{z#6
}[|"db
结果运行错误，显示phi和PHI的维度不匹配，想问问该怎么办？我的解法思路对吗？如不对，可否告知编程方法，万分感谢！ l[gL(p "W
"+Kr1nW
