moving_target_influence

  • julista
    了解作者
  • matlab
    开发工具
  • 178.2KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 5 积分
    下载积分
  • 0
    下载次数
  • 2022-05-25 21:11
    上传日期
通过MATLAB展示仿真运动目标对SAR成像的影响
moving_target_influence.rar
  • rd2_3.jpg
    39.9KB
  • rd2_2.jpg
    39.2KB
  • rd2.m
    3.1KB
  • rd2_4.jpg
    39.8KB
  • rd2_1.jpg
    39.1KB
  • rd2_5.jpg
    39.8KB
内容介绍
clc;clear; c=3e8;%光速 fc=1.5e9;%载频 lambda=c/fc;%波长 %%测绘带区域 X0=200;%方位向[-X0,X0] Rc=3000; R0=150;%距离向[Rc-R0,Rc+R0] %%距离向(Range),r/t Tr=1.5e-6;%脉宽 1.5us (200m) Br=150e6; %LFM信号带宽 150MHz Kr=Br/Tr; %调频斜率 Nr=512; r=Rc+linspace(-R0,R0,Nr) t=2*r/c;%t域序列 dt=R0*4/c/Nr;%采样周期 f=linspace(-1/2/dt,1/2/dt,Nr);%f域序列 %方位向(Azimuth),x/u v=300;%SAR 平台速度 Lsar=300;%合成孔径长度 Na=1024; x=linspace(-X0,X0,Na); u=x/v;%u域序列 du=2*X0/v/Na; fu=linspace(-1/2/du,1/2/du,Na);%fu域序列 fdc=0;%Doppler调频中心 fdr=-2*v^2/lambda/Rc;%Doppler调频斜率 %%目标位置 Ntar=6;%目标个数 Ptar=[Rc , 0 ,1 %距离向坐标,方位向坐标,目标截面积RCS sigma Rc+50 , -50 ,1 Rc+50 , 50 ,1 Rc+50 , -150,1 Rc+50 , 150,1 Rc+100 , 0 ,1]; %%产生回波 s_ut=zeros(Nr,Na); U=ones(Nr,1)*u;%扩充为矩阵 T=t'*ones(1,Na); for i=1:Ntar rn=Ptar(i,1);xn=Ptar(i,2);sigma=Ptar(i,3); if i==1 R=sqrt((rn+3*U).^2+(xn+30*U-v*U).^2); else R=sqrt(rn^2+(xn-v*U).^2); end DT=T-2*R/c; phase=pi*Kr*DT.^2-2*pi/lambda*R*2; s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)<Tr/2).*(abs(v*U-xn)<Lsar/2); end; %距离压缩 p0_t=exp(j*pi*Kr*(t-2*Rc/c).^2).*(abs(t-2*Rc/c)<Tr/2);%距离向匹配函数(行向量) p0_f=fftshift(fft(fftshift(p0_t))); s_uf=fftshift(fft(fftshift(s_ut)));%对回波信号做距离向FFT src_uf=s_uf.*(conj(p0_f).'*ones(1,Na));%匹配函数扩充为矩阵,对每一列进行距离压缩 src_ut=fftshift(ifft(fftshift(src_uf)));%距离压缩后的信号 src_fut=fftshift(fft(fftshift(src_ut).')).';%方位向FFT 距离-多普勒域 %距离弯曲校正(二维去耦) src_fuf=fftshift(fft(fftshift(src_uf).')).';%距离压缩后的二维频谱 F=f'*ones(1,Na);%扩充为矩阵 FU=ones(Nr,1)*fu; p0_2f=exp(j*pi/fc^2/fdr*(FU.*F).^2+j*pi*fdc^2/fc/fdr*F-j*pi/fc/fdr*FU.^2.*F); s2rc_fuf=src_fuf.*p0_2f; s2rc_fut=fftshift(ifft(fftshift(s2rc_fuf)));%距离向IFFT 距离-多普勒域 %方位压缩 p0_2fu=exp(j*pi/fdr*FU.^2);%方位向压缩因子 s2rcac_fut=s2rc_fut.*p0_2fu;%方位压缩 s2rcac_fuf=fftshift(fft(fftshift(s2rcac_fut)));%距离方位压缩后的二维频谱 s2rcac_ut=fftshift(ifft(fftshift(s2rcac_fut).')).';%方位向IFFT figure(1) subplot(221) G=20*log10(abs(s_ut)+1e-6); gm=max(max(G)); gn=gm-40;%显示动态范围40dB G=255/(gm-gn)*(G-gn).*(G>gn); imagesc(x,r-Rc,-G),colormap(gray) grid on,axis tight, xlabel('Azimuth') ylabel('Range') title('(a)原始信号') subplot(222) G=20*log10(abs(src_fut)+1e-6); gm=max(max(G)); gn=gm-40;%显示动态范围40dB G=255/(gm-gn)*(G-gn).*(G>gn); imagesc(fu,r-Rc,-G),colormap(gray) grid on,axis tight, xlabel('Azimuth') ylabel('Range') title('(b)RD域频谱') subplot(223) G=20*log10(abs(s2rc_fut)+1e-6); gm=max(max(G)); gn=gm-40;%显示动态范围40dB G=255/(gm-gn)*(G-gn).*(G>gn); imagesc(fu,r-Rc,-G),colormap(gray) grid on,axis tight, xlabel('Azimuth') ylabel('Range') title('(c)RMC后的RD域频谱') subplot(224) G=20*log10(abs(s2rcac_ut)+1e-6); gm=max(max(G)); gn=gm-60;%显示动态范围40dB G=255/(gm-gn)*(G-gn).*(G>gn); imagesc(x,r-Rc,G),colormap(gray) grid on,axis tight, xlabel('Azimuth') ylabel('Range') title('(d)目标图象')
评论
    相关推荐
    • Matlab合集
      冈萨雷斯数字图像处理MATLAB版.中文版+数字图像处理第二版中文版(冈萨雷斯)+MATLAB-R2014a完全自学一本通+MATLAB R2016a完全自学一本通 素材文件+[模式识别与智能计算:MATLAB技术实现(第2版)].杨淑莹.扫描版
    • MATLAB教程
      MATLAB教程MATLAB教程MATLAB教程MATLAB教程MATLAB教程MATLAB教程
    • MATLAB
      MATLAB 该项目是在matlab上完成的,涉及创建和移动宇宙飞船和机器人。 太空飞船和机器人是使用Matlab中的简单几何形状创建的,并通过连续变换矩阵进行移动。 这个项目教我如何使用变换矩阵(旋转,平移等)的概念...
    • MATLAB基础
      一本学习matlab的一本好书
    • MATLAB编译器
      基于MATLAB 2018b版本介绍MATLAB编译器。介绍如何利用编译器将MATLAB代码编译为独立应用程序或组件,并在没有安装MATLAB的计算机上进行部署。
    • matlabruntime
      通过奇点容器运行您的matlab项目 可以在没有MATLAB的容器中运行matlab代码。 为此,我们首先需要通过Matlab编译器在本地构建相应的Matlab代码的独立应用程序,然后使用具有Matlab运行时( 的容器)运行该应用程序 )...
    • matlab实现
      matlab实现 matlab实现matlab实现matlab实现matlab实现
    • matlab 教程
      matlab 信号处理资料,里面包含信号处理pdf文档,一些杂乱的程序和命令等
    • matlab教程
      matlab教程,PPT格式,包含 matlab基本知识、matlab入门、matlab作图、线性规划、无约束优化、非线性规划、统计工具箱、差值、微分方程等多项知识点,并且每个知识点独立成为PPT,内还含有matlab信号处理详解等文档...
    • matlab简介
      1.MATLAB 开发环境 1.1 MATLAB 的视窗环境 进入MATLAB之后,会看到一个视窗MATLAB Command Window称为指令视窗,它是你键入指令的地方同时 MATLAB也将计算结果显示在此。 1.2 简易计算 我们先从MATLAB的...