求图像梯度并加噪声

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）之间
B=awgn(A,5);%加入信噪比为5分贝的高斯白噪声
hy = fspecial('sobel');%
hx = hy';
By = imfilter(double(B), hy, 'replicate');
Bx = imfilter(double(B), hx, 'replicate');
gradmag = sqrt(Bx.^2 + By.^2);%求图像梯度
fid=fopen('C:\he.txt','wt');
[m,n]=size(gradmag);
x=1:m;
y=1:n;
k=gradmag';
for i=1:m;
    for j=1:n;
        fprintf(fid,'%d\t',x(i));
        fprintf(fid,'%d\t',y(j));
        fprintf(fid,'%d\t',k(i,j));
        fprintf(fid,'\n');
    end
end
fclose(fid);