[求助]用matlab做光學(xué)的數(shù)值模擬 [復(fù)制鏈接]

最近在用matlab做相干合成的數(shù)值模擬，在畫圖的時(shí)候遇到些問題，不知哪位高手能指點(diǎn)一下，多謝~
]l }v  
@<P;F  
這個程序是計(jì)算遠(yuǎn)場光場中主瓣的能量，然后將主瓣能量和f的值畫成散點(diǎn)圖，最后畫圖步驟總是報(bào)錯. L7- JK3/E  
clear; B_#M)d O  
lam=500e-9;N1=4;N2=4; n<EIu  
a=2e-4;z=5; gs}&a3d7k  
k=2*3.14/lam;omega=0.00025; #U'}g *  
xm=0.4*lam*z/a;ym=0.4*lam*z/a; rG-x 3>b  
n=5; e]1)_;b*  
syms thetax; !]"M]tyv\  
d=linspace(0,1,n); DFfh!KKR$  
for i=1:n wR 2`*.O  
    B1(i,:)=(sin(0.5*N1*k*d(i)*thetax)./sin(k*0.5*d(i)*thetax)).^2 ..., k