【雷达】模拟合成孔径雷达(SAR)的点目标仿真

【雷达】模拟合成孔径雷达(SAR)的点目标仿真

1 简介

基于理想情况下聚束式SAR点目标回波信号数学模型。

2 完整代码

%%IFFT in row of matrixfunction s=ifty(fs);s=fftshift(ifft(fftshift(fs.'))).';

%%FFT in column of matrixfunction fs=ftx(s);fs=fftshift(fft(fftshift(s)));

%%IFFT in column of matrixfunction s=iftx(fs);s=fftshift(ifft(fftshift(fs)));

%%FFT in row of matrixfunction fs=fty(s);fs=fftshift(fft(fftshift(s.'))).';

%%================================================================%%Filename: stripmapSAR.m%%Help file: stripmapSAR.doc%%Project: Stripmap SAR Simulation using point targets and Reconstrction %%================================================================clear;clc;close all;%%================================================================%%Parameter--constantC=3e8; %propagation speed%%Parameter--radar characteristicsFc=1e9; %carrier frequency 1GHzlambda=C/Fc; %wavelength %%Parameter--target areaXmin=0; %target area in azimuth is within[Xmin,Xmax]Xmax=50; Yc=10000; %center of imaged areaY0=500; %target area in range is within[Yc-Y0,Yc+Y0] %imaged width 2*Y0%%Parameter--orbital informationV=100; %SAR velosity 100 m/sH=5000; %height 5000 mR0=sqrt(Yc^2+H^2);%%Parameter--antennaD=4; %antenna length in azimuth directionLsar=lambda*R0/D; %SAR integration lengthTsar=Lsar/V; %SAR integration time%%Parameter--slow-time domainKa=-2*V^2/lambda/R0; %doppler frequency modulation rateBa=abs(Ka*Tsar); %doppler frequency modulation bandwidthPRF=Ba; %pulse repitition frequencyPRT=1/PRF; %pulse repitition timeds=PRT; %sample spacing in slow-time domainNslow=ceil((Xmax-Xmin+Lsar)/V/ds); %sample number in slow-time domainNslow=2^nextpow2(Nslow); %for fftsn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);%discrete time array in slow-time domainPRT=(Xmax-Xmin+Lsar)/V/Nslow; %refreshPRF=1/PRT;ds=PRT;%%Parameter--fast-time domainTr=5e-6; %pulse duration 10usBr=30e6; %chirp frequency modulation bandwidth 30MHzKr=Br/Tr; %chirp slopeFsr=3*Br; %sampling frequency in fast-time domaindt=1/Fsr; %sample spacing in fast-time domainRmin=sqrt((Yc-Y0)^2+H^2);Rmax=sqrt((Yc+Y0)^2+H^2+(Lsar/2)^2); %SAR姝ｅ鐩爣锛屾墍浠ヤ笌鏈?繙澶勭殑璺濈涓哄悎鎴愬瓟寰勯暱搴︾殑涓?崐Nfast=ceil(2*(Rmax-Rmin)/C/dt+Tr/dt);%sample number in fast-time domainNfast=2^nextpow2(Nfast); %for ffttm=linspace(2*Rmin/C,2*Rmax/C+Tr,Nfast); %discrete time array in fast-time domaindt=(2*Rmax/C+Tr-2*Rmin/C)/Nfast; %refreshFsr=1/dt;%%Parameter--resolutionDY=C/2/Br; %range resolutionDX=D/2; %cross-range resolution%%Parameter--point targetsNtarget=1; %number of targets%format [x, y, reflectivity]Ptarget=[Xmin,Yc,1];disp('Parameters:')disp('Sampling Rate in fast-time domain');disp(Fsr/Br)disp('Sampling Number in fast-time domain');disp(Nfast)disp('Sampling Rate in slow-time domain');disp(PRF/Ba)disp('Sampling Number in slow-time domain');disp(Nslow)disp('Range Resolution');disp(DY)disp('Cross-range Resolution');disp(DX) disp('SAR integration length');disp(Lsar) disp('Position of targets');disp(Ptarget)%%================================================================%%Generate the raw signal dataK=Ntarget; %number of targetsN=Nslow; %number of vector in slow-time domainM=Nfast; %number of vector in fast-time domainT=Ptarget; %position of targetsSrnm=zeros(N,M);for k=1:1:K sigma=T(k,3); Dslow=sn*V-T(k,1); R=sqrt(Dslow.^2+T(k,2)^2+H^2); tau=2*R/C; Dfast=ones(N,1)*tm-tau'*ones(1,M); phase=pi*Kr*Dfast.^2-(4*pi/lambda)*(R'*ones(1,M)); Srnm=Srnm+sigma*exp(j*phase).*(0

3 仿真结果

4 参考文献

[1]吕守业, 龙腾, and 闫肃. "机载聚束式合成孔径雷达(SAR)点目标的仿真研究." 第十四届全国遥感技术学术交流会 0.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。