f7_12.m

《Matlab时频分析技术及其应用》电子书。

clear;
%仿真合成孔径雷达信号
colormap(gray(256))
cj=sqrt(-1);
pi2=2*pi;

c=3e8;              % 传播速度
fc=200e6;           % 频率
lambda=c/fc;        % 波长
k=pi2/lambda;       % 波数
Xc=1.e3;            % 距离目标中心的距离

L=400;              % 合成孔径大小为2*L
Y0=100;             % 目标区域位于[Yc-Y0,Yc+Y0]
Yc=0;               % 距离目标中心的跨距离

theta_c=atan(Yc/Xc);  %到目标中心的斜视角
Rc=sqrt(Xc^2+Yc^2);   %到目标中心的斜视距离
kus=2*k*sin(theta_c); %Doppler频率偏移

Xcc=Xc/(cos(theta_c)^2);       
DY=(Xcc*lambda)/(4*L);         % 交叉距离分辨率
L_min=max(Y0,L);               % 补零的孔径

% 压缩信号的u域参数和阵列
du=(Xcc*lambda)/(4*(Y0+L));    % 孔径域的采样间隔
duc=(Xcc*lambda)/(4*Y0);       %压缩信号孔径域的采样间隔
mc=2*ceil(L_min/duc);            
uc=duc*(-mc/2:mc/2-1);        
dkuc=pi2/(mc*duc);                % ku域的采样间隔
kuc=dkuc*(-mc/2:mc/2-1);          % kuc阵列
%
dku=dkuc;                         % ku域的采样间隔

%合成孔径信号的u域参数和阵列
m=2*ceil(pi/(du*dku));            %样本数
du=pi2/(m*dku);
u=du*(-m/2:m/2-1);                % 合成孔径阵列
ku=dku*(-m/2:m/2-1);              % ku阵列
%
ntarget=2;             % 目标数
%目标的坐标和反射率
yn(1)=.7*Y0;         fn(1)=0.5;
yn(2)=yn(1)-4*DY;    fn(2)=1;

%测量的回波信号
s=zeros(1,mc);         
for i=1:ntarget;
 dis=sqrt(Xc^2+(Yc+yn(i)-uc).^2);
 s=s+fn(i)*exp(-cj*2*k*dis).*(abs(uc) <= L);
end;
figure(1);
plot(uc,real(s));
xlabel('合成孔径 u');
ylabel('幅值 A');
title('孔径信号的实部');
axis([uc(1) uc(mc) 1.1*min(real(s)) 1.1*max(real(s))]);

sig=s(1:50);
%计算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig',1:length(sig),256);
%显示Wigner-Ville分布等高线
figure(2);
mesh(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');

%Hough变换
[WH,rho,theta]=htl(tfr);
figure(3);
mesh(rho,theta,WH');
xlabel('极半径 rho');
ylabel('角度 theta');
zlabel('幅值 A');