最近在用matlab做相干合成的数值模拟,在画图的时候遇到些问题,不知哪位高手能指点一下,多谢~ [KwwhI@3
G3t
4$3|
这个程序是计算远场光场中主瓣的能量,然后将主瓣能量和f的值画成散点图,最后画图步骤总是报错. z50f$!?
clear; K?aUIkVs
lam=500e-9;N1=4;N2=4; t8FgQ)tk
a=2e-4;z=5; +v'n[xa1v
k=2*3.14/lam;omega=0.00025; XvW
$B|
xm=0.4*lam*z/a;ym=0.4*lam*z/a; auQfWO[ u
n=5; J*O$)K%Hx
syms thetax; ,<?M/'4}G
d=linspace(0,1,n); j<*`?V^
for i=1:n >@U
lhJtW
B1(i,:)=(sin(0.5*N1*k*d(i)*thetax)./sin(k*0.5*d(i)*thetax)).^2 ..., yj{:%Km:`
*exp(-0.5*k.^2*omega.^2*thetax.^2); 5Ai$1'*p
I(i,:)=int(B1(i,:),thetax,-lam/(N1*d(i)),lam/(N1*d(i))); <0I=XsE1iX
Br1=I/16; j\8'P9~%
f(i,:)=2*omega/d(i); tc<t%]c
vpa(Br1,5) _a,XL<9 I
end YJ^TO\4WM
plotmatrix(f,Br1)