在matlab中生成高斯小山图

clear all
close all
clc
%定义系统参数
B=0.2;%基线长度
alpha=pi/6;%基阵倾角
H=30;%基阵高度
r=0.01;%波长
a=50;%小山距离原点的距离
x1=0;y1=0;z1=30;%阵元1的位置
x2=-B*cos(alpha);y2=0;z2=z1+B*sin(alpha);%阵元2位置
sigma=10;%高斯函数方差
x=linspace(10,70,300);
y=linspace(-30,30,300);
[x,y]=meshgrid(x,y);
z=10000*(1/(2*pi*sigma^2)).*exp(-((x-50).^2+(y).^2)/(2*sigma^2));%高斯小山的各点的高度
p1=sqrt((x).^2+(z-z1).^2);%斜距1（阵元1到小山的距离）
p2=sqrt((x).^2+(z-z2).^2);%斜距2（阵元2到小山的距离）
PHi=((2*pi)/r).*(p1-p2);%相位差
PHi=mod(PHi,2*pi);
A=PHi;%转化到（0,2pi）之间
min=min(min(A));
max=max(max(A));
B=imnoise(A,'gaussian',0,0.01);%加入方差为0.01的高斯噪声
figure,image(B);colorbar;colormap(gray);%黑白显示
figure,imagesc(B);colorbar;colormap(jet);%彩色显示