矩形面天线的 三维的方向图问题求助

小弟我用MATLAB 写了这个  矩形面天线的 三维的方向图  ，出现了 维数不一致的问题，请问在那里出问题啊?

clear;
sita=meshgrid(0:pi/90:pi);
fai=meshgrid(0:2*pi/90:2*pi);
W=0.118%贴片的长度
L=0.178;%贴片的宽度
F=444;%中心频率
A1=(3.1415*W*1e6*F)/(3*1e8);
A2=(3.1415*L*1e6*F)/(3*1e8);
F=abs((sin(A1*cos(sita))*sin(sita)*cos(A2*sin(sita)*sin(fai)))/(A1*cos(sita)));%三维的方向图
Fh=abs((sin(A1*cos(sita))*sin(sita))/(A1*cos(sita)));% H 面的方向函数
Fe=abs(cos(A2*sin(fai)));% E 面的方向函数
[x1,y1,z1]=sph2cart(fai,pi/2-sita,r1);
[x2,y2,z2]=sph2cart(fai,pi/2-sita,r2);
[x3,y3,z3]=sph2cart(fai,pi/2-sita,r3);
subplot(2,2,1);
subplot(x1,y1,z1);axis([-11 -11 -11]); shading interp;
subplot(2,2,2);
subplot(x2,y2,z2);axis([-11 -11 -11]); shading interp;
subplot(2,2,3);
subplot(x3,y3,z3);axis([-11 -11 -11]); shading interp;

根据小编提供  F函数，我画出的三维图如下面所示，不知道是不是你想象的三维图
（我不知道你的贴片是如何放置的?）

theta = (0:1:180)/180*pi;
phi   = (0:2:360)/180*pi;
[Phi Theta] = meshgrid(phi, theta);
%%%这里算F
X = F.*sin(Theta).*cos(Phi);
Y = F.*sin(Theta).*sin(Phi);
Z = F.*cos(Theta);
figure;
surf(X, Y, Z);
这样才是三维方向图
你的矩阵操作好像没有带点哦..( .* 和* 的区别)