matlab 画图 画梳状线,用matlab对求和进行画图

可这样：

clear

clc

M=16;            %行

N=16;            %列

K=32;             % number of pulse

CNR=60;

c=3*10^(8);

fr=2800;         %脉冲频率

H=6000;         % platform height

lamda=0.2;

dH=lamda/2;

dL=lamda/2;

v=140;

Range=(H+300):10:5e4;          %距离环斜距

nu_r=length(Range);            %距离环的个数

fai_r=asin(H./Range);          %求俯仰角

theta=pi/2;%主波束的方向

fai=0;%俯仰角

phi=pi/2; %锥角

S=zeros(nu_r,300);

%S=[];

for i=1:nu_r

g_fai=sum(exp(-1j*2*pi*dH/lamda*sin(fai)*(0:M-1)).*(chebwin(M,20).')*(exp(1j*2*pi*dL/lamda*sin(fai_r(i)*pi/180)*(0:M-1).')));

f_fai=sum(exp(-1j*2*pi*dH/lamda*(sin(fai)-sin(fai_r(i)))*(0:M-1)).*(chebwin(M,20).'));%列子阵发射方向图

for ii=1:300

F_dafai=sum(exp(-1j*2*pi*dH/lamda*cos(phi)*(0:N-1)').*...

exp(1j*2*pi*dH/lamda*cos(fai_r(i))*cos(ii*0.6*pi/180)*(0:N-1).').*chebwin(N,20));%行子阵的发射方向图  ii*0.6*pi/180方位角

F=g_fai*F_dafai*f_fai;            %总发射方向图

S(i,ii)=F;

%S(nu_r,ii)=F;

end

end

%x=min(fai_r):1/(nu_r-1):max(fai_r);

x1=linspace(min(fai_r),max(fai_r),300);

%y=0.6*pi/180:1/299:pi;

y1=linspace(0.6*pi/180,pi,nu_r);

[X,Y]=meshgrid(x1,y1);

[x,y]=meshgrid(x1,y1);

%z=S;

z=real(S);

Z=griddata(x,y,z,X,Y);

mesh(X,Y,Z);