求一道MATLAB程序（来自成都左右小区）

%半波\全波、1.5波，2波天线阵子方向性图k=2;%振子数量l=1/4;delta=pi/200;th=0:delta:pi;phi=0:delta:2*pi;[phi,th]=meshgrid(phi,th);F=(cos(k*l*pi.*cos(th))-cos(k*l*pi))./sin(th);X = F.*sin(th).*cos(phi);Y = F.*sin(th).*sin(phi);Z = F.*cos(th);subplot(2,2,1);surf(X,Y,Z)%axis('square')colormap(spring)shading interpl=1/2;delta=pi/200;th=0:delta:pi;phi=0:delta:2*pi;[phi,th]=meshgrid(phi,th);F=(cos(k*l*pi.*cos(th))-cos(k*l*pi))./sin(th);X = F.*sin(th).*cos(phi);Y = F.*sin(th).*sin(phi);Z = F.*cos(th);subplot(2,2,2);surf(X,Y,Z)colormap(spring)shading interpl=3/4;delta=pi/200;th=0:delta:pi;phi=0:delta:2*pi;[phi,th]=meshgrid(phi,th);F=(cos(k*l*pi.*cos(th))-cos(k*l*pi))./sin(th);X = F.*sin(th).*cos(phi);Y = F.*sin(th).*sin(phi);Z = F.*cos(th);subplot(2,2,3);surf(X,Y,Z)colormap(spring)shading interpl=1;delta=pi/200;th=0:delta:pi;phi=0:delta:2*pi;[phi,th]=meshgrid(phi,th);F=(cos(k*l*pi.*cos(th))-cos(k*l*pi))./sin(th);X = F.*sin(th).*cos(phi);Y = F.*sin(th).*sin(phi);Z = F.*cos(th);subplot(2,2,4);surf(X,Y,Z)%axis('square')colormap(spring)shading interp 收起