clc
clear all
close all
M=7;
lambda=5.32e-4;
k=2*pi/lambda;
wo=10*lambda;
z0=k*wo^2/2;
z=0*z0;
q=1+i*z/z0;
a=1/(10*lambda);
c=10;
n=3;
fai =0:2*pi;
a1=0.02;
x=linspace(-a1,a1,M);
y=linspace(-a1,a1,M);
[x1,y1]=meshgrid(x,y);
[theta,r]=cart2pol(x1,y1);
xx=1./q.*sqrt(a*r.*(a*r-2*c*q.*exp(i*theta)));
I=1./q.*(a*r./(a*r-2*c*q.*exp(i*theta))).^n.*besselj(n,xx).*besselj(n,xx);
En=I.^(1/2).*cos(fai);
figure(471)
mesh(abs(En).^2),view(2)
ph=angle(En);
figure(288)
imshow(ph,[])
轉載請註明出處，本文鏈接：https://www.uj5u.com/shujuku/38044.html
標籤：疑難問題
上一篇：幫忙看下這個sql有什么問題