【雷达通信】合成孔径雷达(SAR)的点目标仿真matlab源码
- *SAR原理简介*
用一个小天线作为单个辐射单元,将此单元沿一直线不断移动,在不同位置上接收同一地物的回波信号并进行相关解调压缩处理。一个小天线通过“运动”方式就合成一个等效“大天线”,这样可以得到较高的方位向分辨率,同时方位向分辨率与距离无关,这样SAR就可以安装在卫星平台上而可以获取较高分辨率的SAR图像。
图1 SAR成像原理示意图
1、几个参重要参数
为了更好的理解SAR和SAR图像,需要知道几个重要的参数。
- 分辨率
SAR图像分辨率包括距离向分辨率(Range Resolution)和方位向分辨率(Azimuth Resolution)。
图2 距离向和方位向示意图
- 距离向分辨率(Range Resolution)
垂直飞行方向上的分辨率,也就是侧视方向上的分辨率。距离向分辨率与雷达系统发射的脉冲信号相关,与脉冲持续时间成正比:
Res( r) = c*τ/2
其中c为光速,τ为脉冲持续时间。
- 方位向分辨率(Azimuth Resolution)
沿飞行方向上的分辨率,也称沿迹分辨率。如下为推算过程:
• 真实波束宽度: β= λ/ D
• 真实分辨率:ΔL = β*R = Ls (合成孔径长度)
• 合成波束宽度βs = λ /(2* Ls) = D / (2* R)
• 合成分辨率ΔLs = βs* R = D / 2
其中λ为波长,D为雷达孔径,R为天线与物体的距离。
从这个公式中可以看到,SAR系统使用小尺寸的天线也能得到高方位向分辨率,而且与斜距离无关(就是与遥感平台高度无关)。
图3 方位向分辨率示意图
- 极化方式
雷达发射的能量脉冲的电场矢量,可以在垂直或水平面内被偏振。无论哪个波长,雷达信号可以传送水平(H)或者垂直(V)电场矢量。接收水平(H)或者垂直(V)或者两者的返回信号。雷达遥感系统常用四种极化方式———HH、VV、HV、VH。前两者为同向极化,后两者为异向(交叉)极化。
极化是微波的一个突出特点,极化方式不同返回的图像信息也不同。返回同极化(HH或者VV)信号的基本物理过程类似准镜面反射,比如,平静的水面显示黑色。交叉极化(HV或者VH)一般返回的信号较弱,常受不同反射源影响,如粗糙表面等。
图4 HV极化示意图
- 入射角(Incidence Angle)
入射角也叫视角,是雷达波束与垂直表面直线之间的夹角(如下图中的θ)。微波与表面的相互作用是非常复杂的,不同的角度区域会产生不同的反射。低入射角通常返回较强的信号,随着入射角增加,返回信号逐渐减弱。
根据雷达距离地表高度的情况,入射角会随着近距离到远距离的改变而改变,依次影响成像几何。\
图5 SAR入射角示意图
2、SAR拍摄模式
SAR主要有三种拍摄模式:Stripmap,ScanSAR和 Spotlight。
当然最新的SAR系统拥有更多的拍摄模式,比如RADARSAT-2还用于超精细、高入射角等拍摄模式
- 条带模式-Stripmap
当运行Stripmap 模式时,雷达天线可以灵活的调整,改变入射角以获取不同的成像宽幅。
最新的SAR系统都具有这种成像模式,包括RADARSAT-1/2, ENVISAT ASAR, ALOS PALSAR,TerraSAR-X-1, COSMOSkyMed和RISAT-1。
- 扫描模式-ScanSAR
扫描模式是共享多个独立sub-swaths的操作时间,最后获取一个完整的图像覆盖区域。它能解决Stripmap模式较小的刈幅。
- 聚束模式-Spotlight
当执行聚束模式采集数据时,传感器控制天线不停向成像区域发射微波束。
它与条带模式主要区别为:
在使用相同物理天线时,聚束模式提供更好的方位分辨率;
在可能成像的以一个区域内,聚束模式在单通道上的提供更多的视角;
聚束模式可以更有效的获取多个小区域。
-
图8聚束模式
-
%%================================================================
%%Filename: stripmapSAR.m
%%Help file: stripmapSAR.doc
%%Project: Stripmap SAR Simulation using point targets and Reconstrction
%% with Range-Doppler Algorithm
%%Author: Zhihua He ,National University of Defence Tecnology ,2005/3
%%E-mail: skynismile@yahoo.com.cn
%%================================================================
clear;clc;close all;
%%================================================================
%%Parameter--constant
C=3e8; %propagation speed
%%Parameter--radar characteristics
Fc=1e9; %carrier frequency 1GHz
lambda=C/Fc; %wavelength
%%Parameter--target area
Xmin=0; %target area in azimuth is within[Xmin,Xmax]
Xmax=50;
Yc=10000; %center of imaged area
Y0=500; %target area in range is within[Yc-Y0,Yc+Y0]
%imaged width 2*Y0
%%Parameter--orbital information
V=100; %SAR velosity 100 m/s
H=5000; %height 5000 m
R0=sqrt(Yc^2+H^2);
%%Parameter--antenna
D=4; %antenna length in azimuth direction
Lsar=lambda*R0/D; %SAR integration length
Tsar=Lsar/V; %SAR integration time
%%Parameter--slow-time domain
Ka=-2*V^2/lambda/R0; %doppler frequency modulation rate
Ba=abs(Ka*Tsar); %doppler frequency modulation bandwidth
PRF=Ba; %pulse repitition frequency
PRT=1/PRF; %pulse repitition time
ds=PRT; %sample spacing in slow-time domain
Nslow=ceil((Xmax-Xmin+Lsar)/V/ds); %sample number in slow-time domain
Nslow=2^nextpow2(Nslow); %for fft
sn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);%discrete time array in slow-time domain
PRT=(Xmax-Xmin+Lsar)/V/Nslow; %refresh
PRF=1/PRT;
ds=PRT;
%%Parameter--fast-time domain
Tr=5e-6; %pulse duration 10us
Br=30e6; %chirp frequency modulation bandwidth 30MHz
Kr=Br/Tr; %chirp slope
Fsr=3*Br; %sampling frequency in fast-time domain
dt=1/Fsr; %sample spacing in fast-time domain
Rmin=sqrt((Yc-Y0)^2+H^2);
Refa=exp(j*pi*Ka*ta.^2).*(abs(ta)<Tsar/2);
Sa=iftx(ftx(Sr).*(conj(ftx(Refa)).'*ones(1,M)));
Ga=abs(Sa);
%%================================================================
%%graw the intensity image of signal
colormap(gray);
figure(1)
subplot(211);
row=tm*C/2-2008;col=sn*V-26;
imagesc(row,col,255-Gr); %intensity image of Sr
axis([Yc-Y0,Yc+Y0,Xmin-Lsar/2,Xmax+Lsar/2]);
xlabel('\rightarrow\itRange in meters'),ylabel('\itAzimuth in meters\leftarrow'),
title('Stripmap SAR after range compression'),
subplot(212);
imagesc(row,col,255-Ga); %intensity image of Sa
axis([Yc-Y0,Yc+Y0,Xmin-Lsar/2,Xmax+Lsar/2]);
xlabel('\rightarrow\itRange in meters'),ylabel('\itAzimuth in meters\leftarrow'),
title('Stripmap SAR after range and azimuth compression'),
%%================================================================
%%draw 3D picture
figure(2)
waterfall(real(Srnm((200:205),:)));axis tight
xlabel('Range'),ylabel('Azimuth'),
title('Real part of the raw signal'),
figure(3)
waterfall(Gr((200:205),(600:1000)));axis tight
xlabel('Range'),ylabel('Azimuth'),
title('Stripmap SAR after range compression'),
figure(4)
mesh(Ga((200:300),(750:860)));axis tight
xlabel('Range'),ylabel('Azimuth'),
title('Stripmap SAR after range and azimuth compression'),
%%================================================================
%%draw -3dB contour
figure(5)
a=max(max(Ga));
contour(row,col,Ga,[0.707*a,a],'b');grid on
axis([9995,10050,-20,20]),
xlabel('\rightarrow\itRange in meters'),ylabel('\itAzimuth in meters\leftarrow'),
title('Resolution Demo: -3dB contour');
%%================================================================
【雷达通信】合成孔径雷达(SAR)的点目标仿真matlab源码相关推荐
- 【卡尔曼滤波】卡尔曼滤波在雷达目标跟踪中的应用仿真matlab源码
1 模型 [摘要]目标跟踪问题的应用背景是雷达数据处理,即雷达在搜索到目标并记录目标的位置数据,对测量到的目标位置数据(称为点迹)进行处理,自动形成航迹,并对目标在下一时刻的位置进行预测.本文简要讨论 ...
- 交通流元胞自动机模拟仿真 matlab源码_元胞自动机中的时间反演
1.引言--生命游戏 1970年,英国数学家约翰·何顿·康威提出了生命游戏(Life Game).生命游戏本质是一个元胞自动机模型,每个元胞可以看作是一个细胞,细胞的产生.繁衍和死亡拥有3条演化规则. ...
- 【元胞自动机】基于元胞自动机的人口疏散仿真matlab源码
1 简介 为了消除礼堂的安全隐患,制定 行之有效的应急预案,有必要对礼堂人群疏散运动进行研究,掌握礼堂人群疏散的一般特点和规律.采用基于二维元胞自动机模型对某高校礼堂发生人群疏散运动进 行仿真,找出影 ...
- 【物理应用】大气湍流相位屏仿真matlab源码
1大气湍流相位屏模型 2 部分代码 %% example_ft_sh_phase_screen.m 利用谐波方法产生随机相位屏的例子 clc; Dx = 2; % length of one side ...
- 【图像分割】各向异性热扩散方程SAR图像分割【含Matlab源码 1879期】
⛄一.各向异性热扩散方程简介(具体理论见参考文献) 1 前言 图像分割是 图像解译的一个关键步骥. 图像分割属于图像解译的早期处理 , 目的是为更高水平的后续处理提供输人. 因此图像分割的好坏对后续 ...
- 【雷达通信】基于matlab雷达探测威力仿真【含Matlab源码 1974期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[雷达通信]基于matlab雷达探测威力仿真[含Matlab源码 1974期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付 ...
- 【雷达图像】SAR合成孔径雷达成像及处理【含Matlab源码 307期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[雷达图像]SAR合成孔径雷达成像及处理[含Matlab源码 307期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付费专 ...
- 【雷达通信】基于matlab雷达仿真模拟系统【含Matlab源码 150期】
⛄一.获取代码方式(附课程作业报告) 获取代码方式1: 完整代码已上传我的资源:[雷达通信]基于matlab雷达仿真模拟系统[含Matlab源码 150期] 点击上面蓝色字体,直接付费下载,即可. 获 ...
- 【雷达通信】基于matlab无人机FMCW毫米波高度计雷达仿真【含Matlab源码 1261期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源: [雷达通信]基于matlab无人机FMCW毫米波高度计雷达仿真[含Matlab源码 1261期] 点击上面蓝色字体,直接付费下载,即可. ...
最新文章
- 手误【删库】 == 跑路,不存在的 ——删瓦辛格
- 关于GCC Cygwin MinGW MSYS
- WTM 3.5发布,VUE来了!
- 一个大数据量表访问优化--联动下拉框查询优化
- Swagger2 最全注解说明(转载)
- 海龟交易策略要点总结
- oracle多用户导出导入用法
- .国内外主流前端开发框架对比
- BZOJ-1507 文本编辑器(Editor)
- 淘宝API开发(一) 简单的用TOPAPI来获取淘宝用户信息
- 张宇1000题高等数学 第十、十一、十二章 一元函数积分学的应用——几何应用、积分等式与积分不等式、物理应用
- Unity 内置渲染管线转URP
- 基于FPGA的电梯控制器设计(任意楼层)
- 游戏开发人员需要了解的5种ASO技术,苹果游戏aso优化
- #344 – 通过CanExecute控制按钮是可用(The CommandBinding CanExecute Determines Whether a Button is Enabled)
- QNX 7.1 交叉编译 cron
- 3D打印切片软件--cura的二次开发(2)(界面分析与汉化)(修改版增加图片)
- 回归常态啦 2020.12.08日记
- Android手机,charles安装证书时提示“键入凭据存储的密码”问题
- SQL底层执行原理详解
热门文章
- D.E.Shaw拟以8000万美元收购SunEdison光伏电站股权
- 计算机电源接通显示未充电怎么办,笔记本电脑显示“电源已接通,未充电”怎么办?...
- @Value之${}与#{}
- 软件技术功能开发思路
- 如何将图片转换base64格式?data:image/png;base64又是什么?
- 维基解密泄密狂 科技公司补漏忙
- 如何高效实现多窗口卖票
- 中国主要湖泊水库面状矢量边界图层ArcGIS使用之一篇
- 北交大计算机电子信息,北京交通大学电子信息工程学院
- 两台笔记本组建无线局域网共享上网