matlab未系统学过，已困惑了许久，希望得到高手、老师的帮助！积分函数如图片所示。 E e>j7k.G.
(p=GR#
y>|XpImZ
:g[x;Q [@
我采用了dblquad函数来进行积分，将待积函数设了一个function函数 -FV'%X$i
g/p9"eBpq
function y=E_x(r,PHI,z) lVtn$frp
omega0=7.57e15; /g'- *:a
c=3e11; ( y*X8
k=omega0/c; A'iF'< %
A=@(phi,theta)exp(i*k*(z.*cos(theta)+r.*sin(theta).*cos(phi-PHI))).*cos(theta).*cos(phi); TZNgtR{q
y=dblquad(A,0,2*pi,0,1.12); F
t/
x 5
"B3: m-'
W y*7j B
然后直接调用该函数： CL : M>(
lambda=2*pi/k; Eu.qA9,@U
for z=0:1:2; +|)1_NK
r=linspace(-2*lambda,2*lambda,201); 2 _n*u^X:_
HI=linspace(0,2*pi,201); Z[ u,1l .T
[rho,PHI ]=meshgrid(r,HI); t<,p-TM]
Ie=conj(E_x(rho,PHI,z)).*E_x(rho,PHI,z); 2Q|*xd4B^
end oh0|2 IrM
surf(rho,PHI,Ie) `yX+NRi(s
E(kpK5h{
结果运行错误，显示phi和PHI的维度不匹配，想问问该怎么办？我的解法思路对吗？如不对，可否告知编程方法，万分感谢！ ~PW}sN6ppG
% W',c u
