matlab未系统学过，已困惑了许久，希望得到高手、老师的帮助！积分函数如图片所示。 G^.N$wcv
\V-N~_-H
WE\TUENac(
t]HY@@0g
我采用了dblquad函数来进行积分，将待积函数设了一个function函数 #*g=F4>t
' D 21A8*N
function y=E_x(r,PHI,z) 1 GYZ1iA
omega0=7.57e15; 6q^$}eOt
c=3e11; q %sZV>
k=omega0/c; HPp hTu}`
A=@(phi,theta)exp(i*k*(z.*cos(theta)+r.*sin(theta).*cos(phi-PHI))).*cos(theta).*cos(phi); JV~
Dly>
y=dblquad(A,0,2*pi,0,1.12); T8+[R2_
7( 84j5zb
||:> &
然后直接调用该函数： >eTbg"\
lambda=2*pi/k; sa'1hX^@
for z=0:1:2; gK h*q.
r=linspace(-2*lambda,2*lambda,201); S>oEk3zlw
HI=linspace(0,2*pi,201); `vBBJ@f4)
[rho,PHI ]=meshgrid(r,HI); Ro<779.Gn\
Ie=conj(E_x(rho,PHI,z)).*E_x(rho,PHI,z); mG1~rI
end M/ 64`lcb
surf(rho,PHI,Ie) 0q 5 J)l:
- lo?16w
结果运行错误，显示phi和PHI的维度不匹配，想问问该怎么办？我的解法思路对吗？如不对，可否告知编程方法，万分感谢！ ^5X?WA,Z99
y-hTTd"{
