readdata1.m

一些sar数据

clear all
tp=4.00025602e-5;   %脉冲宽度
fc=9.600000000000000e+009;  %中心频率
c=3e8; 
prf=1395.0;   %重复频率
lmda=c/fc;    %波长
gama=-2.3724394e+11;  %调频率
B=gama*tp;  %带宽
fs=1.12492800e+7;  %采样率

% 对此数据进行距离匹配滤波－方位匹配滤波成像

fd1=fopen('raws.dat','rb');
s=fread(fd1,'uchar');
fclose(fd1)
s=s(1:2:2097152)+j*s(2:2:2097152);
s=reshape(s,1024,1024);
s=s-mean(mean(s));
[ran_num,azi_num]=size(s);
% 列为距离快时间  行为方位慢时间
nan=azi_num;
nrn=1024;

x=s;

for n=1:1:nan
    t=[-nrn/2:nrn/2-1]'/fs;
    y(:,n)=conv(x(:,n),exp(j*pi*gama*t.^2));
end

clear x
ka=9357;
for n=1:1:2047
    t=[-100/2:100/2-1]/prf;
    x(n,:)=conv(y(n,:),exp(j*pi*ka*t.^2));
end