最近在用matlab做相干合成的数值模拟，在画图的时候遇到些问题，不知哪位高手能指点一下，多谢~ $GV7o{"&
krxo"WgD
这个程序是计算远场光场中主瓣的能量，然后将主瓣能量和f的值画成散点图，最后画图步骤总是报错. ho{*Cjv
clear; YpHg&|Fr
lam=500e-9;N1=4;N2=4; wVXS%4|v
a=2e-4;z=5; &A/]pi-\
k=2*3.14/lam;omega=0.00025; nbp=PzZy
xm=0.4*lam*z/a;ym=0.4*lam*z/a; a~`eQ_N D
n=5; R+:yVi[F]U
syms thetax; hR
n <em
d=linspace(0,1,n); m,28u3@r
for i=1:n ZgJQ?S$ D
B1(i,:)=(sin(0.5*N1*k*d(i)*thetax)./sin(k*0.5*d(i)*thetax)).^2 ..., X56q-|
*exp(-0.5*k.^2*omega.^2*thetax.^2); T.F!+
I(i,:)=int(B1(i,:),thetax,-lam/(N1*d(i)),lam/(N1*d(i))); 5<k"K^0QS
Br1=I/16; y f)%%&
f(i,:)=2*omega/d(i); yF:1( 4
vpa(Br1,5) T~?Ff|qFC
end P{`C^W$J^
plotmatrix(f,Br1)
