当前位置：文档库 › 合成孔径雷达(SAR)的点目标仿真(附件带代码程序)

合成孔径雷达(SAR)的点目标仿真(附件带代码程序)

合成孔径雷达(SAR)的点目标仿真

一． SAR原理简介

合成孔径雷达(Synthetic Aperture Radar ，简称SAR)是一种高分辨率成像雷达技术。它利用脉冲压缩技术获得高的距离向分辨率，利用合成孔径原理获得高的方位向分辨率，从而获得大面积高分辨率雷达图像。SAR回波信号经距离向脉冲压缩后，雷达的距离分辨率由雷达发射信号带宽决定：，式中表示雷达的距离分辨率，表示雷达发射信号带宽，表示光速。同样，SAR回波信号经方位向合成孔径后，雷达的方位分辨率由雷达方位向的多谱勒带宽决定：，式中表示雷达的方位分辨率，表示雷达方位向多谱勒带宽，表示方位向SAR平台速度。

二． SAR的成像模式和空间几何关系

根据SAR波束照射的方式，SAR的典型成像模式有Stripmap(条带式)，Spotlight(聚束式)和Scan(扫描模式)，如图2.1。条带式成像是最早研究的成像模式，也是低分辨率成像最简单最有效的方式；聚束式成像是在一次飞行中，通过不同的视角对同一区域成像，因而能获得较高的分辨率；扫描模式成像较少使用，它的信号处理最复杂。

图2.1：SAR典型的成像模式

这里分析SAR点目标回波时，只讨论正侧式Stripmap SAR，正侧式表示SAR波束中心和SAR平台运动方向垂直，如图2.2，选取直角坐标系XYZ为参考坐标系，XOY平面为地平面；SAR平台距地平面高h，沿X轴正向以速度V匀速飞行；P点为SAR平台的位置矢量，设其坐标为(x,y,z)；T点为目标的位置矢量，设其坐标为；由几何关系，目标与SAR平台的斜距为：

(2.1)

由图可知：；令，其中为平台速度，s为慢时间变量（slow time），假设，其中表示SAR平台的x 坐标为的时刻；再令，表示目标与SAR的垂直斜距，重写2.1式为：

(2.2)

就表示任意时刻时，目标与雷达的斜距。一般情况下，，于是2.2式可近似写为：

(2.3)

可见，斜距是的函数，不同的目标，也不一样，但当目标距SAR较远时，在观测带内，可近似认为不变，即。

图2.2：空间几何关系(a)正视图(b)侧视图

图2.2(a)中，表示合成孔径长度，它和合成孔径时间的关系是。(b)中，为雷达天线半功率点波束角，为波束轴线与Z轴的夹角，即波束视角，为近距点距离，为远距点距离，W为测绘带宽度，它们的关系为：

(2.4)

三． SAR的回波信号模型

SAR在运动过程中，以一定的PRT(Pulse Repitition Time,脉冲重复周期)发射和接收脉冲，天线波束照射到地面上近似为一矩形区域，如图2.2（a），区域内各散射元（点）对入射波后向散射，这样，发射脉冲经目标和天线方向图的调制，携带目标和环境信息形成SAR回波。从时域来看，发射和接收的信号都是一时间序列。

图3.1：SAR发射和接收信号

图3.1表示SAR发射和接收信号的时域序列。发射序列中，为chirp信号持续时间，下标表示距离向(R ange)；PRT为脉冲重复周期；接收序列中，表示发射第个脉冲时，目标回波相对于发射序列的延时；阴影部分表示雷达接收机采样波门，采样波门的宽度要保证能罩住测绘带内所有目标的回波。

雷达发射序列的数学表达式为：

(3.1)

式中，表示矩形信号，为距离向chirp信号的调频斜率，为载频。

雷达回波信号由发射信号波形，天线方向图，斜距，目标RCS，环境等因素共同决定，若不考虑环境因素，则单点目标雷达回波信号可写成：

(3.2)

式中，为点目标的雷达散射截面, 表示点目标天线方向图双向幅度加权，表示载机发射第n个脉冲时，电磁波在雷达与目标之间传播的双程时间，，代入3.2式

(3.3)

3.3式就是单点目标回波信号模型。其中，为chirp分量，它决定距离向分辨率，为doppler分量，它决定方位向分辨率。

距离向变量远大于方位向变量t(典型相差量级)，于是一般可以假设SAR满足“停－走－停”模式，即SA R在发射和接收一个脉冲信号中间，载机未发生运动。为了理论分析方便，称为慢时间变量(slow time)，称t为快时间变量(fast time)于是，一维回波信号可以写成二维形式，正交解调去除载波后，单点目标的回波可写成：

(3.3)

图3.2：单点目标回波二维分布示意图

在方位向(慢时间域)是离散的，，其中V是SAR的速度，是0时刻目标在参考坐标系中的x坐标。为了作数字信号处理，在距离向(快时间域)也要采样，假设采样周期为Tr,则，如图3.2，方位向发射N个脉冲，距离向采样得到M个样值点，则SAR回波为一矩阵，K个理想点目标的回波经采样后的表达式为：

(3.4)

上式用Matlab语言可表示为：

%%***************************************************************************

%%Generate the raw signal data

K=Ntarget; %number of targets

N=Nslow; %number of vector in slow-time domain

M=Nfast; %number of vector in fast-time domain

T=Ptarget; %position of targets

Srnm=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

end

%%***************************************************************************

四． SAR的信号系统模型

从信号与系统的角度看，SAR回波可看作目标的散射特性通过一个二维线性系统的输出。点目标的信号与系统模型如图4.1：

图4.1：点目标信号与系统模型

模型的数学表达式为：

(4.1)

式中，表示点目标的散射特性，表示等效系统，设为发射的chirp信号，则：

(4.2)

4.2式表明只在维是线性时不变(LTI)的，在维是时变的，相同的, 不同的，响应不一样。但通常情况下可近似认为不变，即，这时，系统等效为一个二维LTI系统。

五．点目标SAR的成像处理算法仿真

SAR的回波数据不具有直观性，不经处理人无法理解它，如图5.1。从原理上讲，SAR成像处理的过程是从回波数据中提取目标区域散射系数的二维分布，本质上是一个二维相关处理过程,因此最直接的处理方法是对回波进行二维匹配滤波，但其运算量很大，再加上SAR的数据率本来就高，这使得实时处难于实现。通常，可以把二维过程分解成距离向和方位向两个一维过程，Range-Dopper Algorithm（简称RD算法）就是采用这种思想的典型算法，这里也只讨论RD算法。

图5.1：SAR回波数据(a)未经处理(b)处理后

RD算法通过距离迁移(Range Migration)矫正，消除距离和方位之间的耦合。在满足聚焦深度的前提小，将成像处理分解成两个一维的LTI系统进行相关处理，并采用频域快速相关算法提高了速度。RD算法已非

常成熟，并成为衡量其它算法优劣的标准。RD算法典型的数字处理流程如图5.2。

图5.2: Range-Dopper 算法流程

●预处理

这是对SAR回波处理的第一步，一般在SAR平台(卫星，飞机)上实时处理，包括解调和数字化。雷达信号的载频较高(～GHz)，不宜直接采样数字化处理，常常通过正交解调方式解调出基带信号，再对基带信号(～MHz)采样数字化，然后存储或传到地面做进一步处理。采样后的数据常采用矩阵形式存储，假设方位向发射(采样)N个脉冲，距离向采样得到M每个采样值(图3.2)，则待处理数据是一个的矩阵，如图5. 3。实际处理时，要在方位向上加窗截断，因此，在方位向上的开始和结束的一段数据(图中影阴区所示)是不充分的，对应的长度均为，表示SAR的合成孔径长度。仿真时，这个数据阵是程序根据3.4式产生的。

图5.3：待处理数据

●距离压缩

距离向信号是典型的Chirp信号，相关算法是在频域利用FFT进行的。Matlab语句为：

Refr=exp(j*pi*Kr*tr.^2).*(0

Sr=ifty(fty(Srnm).*(ones(N,1)*conj(fty(Refr))));

和分别是Chirp信号的调频斜率和脉冲持续时间，Refr表示参考信号，fty,ifty是对矩阵的行(对应距离向)进行FFT和IFFT运算的子程序。例如，fty的代码为：

%%FFT in row of matrix

function fs=fty(s);

fs=fftshift(fft(fftshift(s.'))).';

●距离迁移矫正

距离迁移是SAR信号处理中必然出现的现象，它的大小随系统参数不同而变化，并不总需要补偿。点目标仿真时，可以先不考虑。

●方位压缩

方位向的处理是SAR成像处理算法最核心的部分。正侧式点目标(图2.2)情况下，回波经距离压缩后在方位向也是一Chirp信号，因此其压缩处理同距离压缩处理类似，只是压缩因子不同。仿真中，调频斜率已知，因此不需要进行Doppler参数估计。

● SAR参数

SAR平台：

水平速度V=100m/s

高度H=5000m

天线等效孔径D=4m

SAR平台与测绘带的垂直斜距R0=11180m(计算结果)

发射信号：

载波频率=1GHz

Chirp信号持续时间=5us

Chirp信号调频带宽=30MHz

Chirp信号调频斜率= (计算结果)

脉冲重复频率PRF=57.6Hz(计算结果)

Doppler调频带宽=50Hz(计算结果)

Doppler调频斜率=－5.96(计算结果)

分辨率：

距离向分辨率DY=5m

方位向分辨率DX=2m

目标位置:

距离向Y=[Yc-Y0,Yc+Y0]=[9500,10500]m

方位向X=[Xmin,Xmax]=[0,50]m

目标个数Ntarget=3

目标位置矩阵: 格式[x坐标,y坐标,目标散射系数]

Ptarget=[Xmin,Yc,1

Xmin,Yc+10*DY,1

Xmin+20*DX,Yc+50*DY,1]

stripmapSAR.m程序(见附录)实现了仿真功能，图5.4到图5.7为仿真结果。运行程序，在Command Window 中列出了仿真的参数：

Parameters:

Sampling Rate in fast-time domain

3.0996

Sampling Number in fast-time domain

1024

Sampling Rate in slow-time domain

1.1525

Sampling Number in slow-time domain

512

Range Resolution

Cross-range Resolution

SAR integration length

838.5255

Position of targets

0 10000 1

0 10050 1

40 10250 1

当然，这些参数可以改变以得到不同的结果，但值得注意的是，采样点数不宜过大，否则数据量过大将导致程序运行时间过长，甚至计算机因内存耗尽而死机。本例采用的是512 1024个点。

图5.4：SAR的点目标仿真结果

图5.5:两点目标的回波仿真3D图

图5.6：两点目标距离向压缩后的3D图

图5.7：两点目标距离向和方位向压缩后的3D图

图5.8：两点目标压缩后的3dB等高线图

附录：SAR的点目标仿真Matlab程序

主程序：stripmapSAR.m

%%======================================================== 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

合成孔径雷达概述(SAR)

合成孔径雷达概述 1合成孔径雷达简介 (2) 1.1 合成孔径雷达的概念 (2) 1.2 合成孔径雷达的分类 (3) 1.3 合成孔径雷达(SAR)的特点 (4) 2合成孔径雷达的发展历史 (5) 2.1 国外合成孔径雷达的发展历程及现状 (5) 2.1.1 合成孔径雷达发展历程表 (6) 2.1.2 世界各国的SAR系统 (9) 2.2 我国的发展概况 (11) 2.2.1 我国SAR研究历程表 (11) 2.2.2 国内各单位的研究现状 (12) 2.2.2.1 电子科技大学 (12) 2.2.2.2 中科院电子所 (12) 2.2.2.3 国防科技大学 (13) 2.2.2.4 西安电子科技大学 (13) 3 合成孔径雷达的应用 (13) 4 合成孔径雷达的发展趋势 (14) 4.1 多参数SAR系统 (15) 4.2 聚束SAR (15) 4.3极化干涉SAR（POLINSAR） (16) 4.4合成孔径激光雷达（Synthetic Aperture Ladar） (16) 4.5 小型化成为星载合成孔径雷达发展的主要趋势 (17) 4.6 性能技术指标不断提高 (17) 4.7 多功能、多模式是未来星载SAR的主要特征 (18) 4.8 雷达与可见光卫星的多星组网是主要的使用模式 (18) 4.9 分布SAR成为一种很有发展潜力的星载合成孔径雷达 (18) 4.10 星载合成孔径雷达的干扰与反干扰成为电子战的重要内容 (19) 4.11 军用和民用卫星的界线越来越不明显 (19) 5 与SAR相关技术的研究动态 (20) 5.1 国内外SAR图像相干斑抑制的研究现状 (20) 5.2 合成孔径雷达干扰技术的现状和发展 (20) 5.3 SAR图像目标检测与识别 (22) 5.4 恒虚警技术的研究现状与发展动向 (25) 5.5 SAR图像变化检测方法 (27) 5.6 干涉合成孔径雷达 (31) 5.7 机载合成孔径雷达技术发展动态 (33) 5.8 SAR图像地理编码技术的发展状况 (35) 5.9 星载SAR天线方向图在轨测试的发展状况 (37) 5.10 逆合成孔径雷达的发展动态 (38) 5.11 干涉合成孔径雷达的发展简史与应用 (38)

C++图书管理系统源代码

图书管理系统系统功能： 1．借书：根据借书人提出的图书编号（id）查询该图书，如果该图书现存量（store）不为0，则提示输入借阅者的学号（num），为借书人办理借书手续，提示用户该书已被借出。 2．归书：根据借书人的学号查询该读者的信息，若有该读者，则提示输入所借书籍的编号（id）,为该读者办理还书手续，提示该书已还。 3.书籍管理：弹出书籍管理界面，输入所要执行操作的号码：（1）增加书籍：弹出注册新书的窗口，按照提示输入所增加书籍的信息，最后，提示用户该书已被注册。（2）删除书籍：弹出删除书籍的窗口，输入所要删除书籍的编号（id），输出该书的信息，确认是否删除该书，1为删除，0为放弃。（3）修改书籍：弹出修改书籍的窗口，输入所要修改书籍的编号（id），输出该书的信息，确认是否修改该书，1为修改，0为放弃。之后按照提示重新输入书籍的信息。 4.读者管理：弹出读者管理界面，输入所要执行操作的号码：（1）增加读者：弹出注册读者的窗口，按照提示输入所增加读者的信息，最后，提示用户该读者已被注册。（2）删除书籍：弹出删除读者的窗口，输入所要删除读者的学号（num），输出该读者的信息，确认是否删除该读者，1为删除，0为放弃。（3）修改书籍：弹出修改读者的窗口，输入所要修改读者的学号（num），输出该读者的信息，确认是否修改该读者，1为修改，0为放弃。之后按照提示重新输入读者的信息。 5.搜索：此搜索包括两方面的搜索，书籍搜索以及读者搜索，弹出搜索的窗口，按照提示输入所要搜索的内容，1为书籍搜索，2为读者搜索：（1）搜索书籍：弹出搜索书籍的窗口，按照提示输入所要搜索的方式，包括按<1>书名搜索， <2>书号搜索，<3>作者搜索，<4>出版社搜索，<5>出版时间搜索;根据所选方式输入相应的内容，若是该书籍存在，则输出该书籍的信息，否则，返回主界面。（2）搜索读者：弹出搜索读者的窗口，按照提示输入所要搜索的方式，包括按<1>名字搜索， <2>学号搜索；根据所选方式输入相应的内容，若是该读者存在，则输出该读者的信息，否则，返回主界面。 6.退出：退出图书管理系统。图书类设计：

雷达系统仿真matlab代码.docx

% ====================================================== =====================================% % 该程序完成16个脉冲信号的【脉压、动目标显示/动目标检测（MTI/MTD）】 % ====================================================== =====================================% % 程序中根据每个学生学号的末尾三位（依次为XYZ）来决定仿真参数，034 % 目标距离为[3000 8025 9000+(Y*10+Z)*200 8025]，4个目标 % 目标速度为[50 0 (Y*10+X+Z)*6 100] % ====================================================== =====================================% close all; %关闭所有图形 clear all; %清除所有变量 clc; % ====================================================== =============================% % 雷达参数 % % ====================================================== =============================% C=3.0e8; %光速(m/s) RF=3.140e9/2; %雷达射频 1.57GHz Lambda=C/RF;%雷达工作波长 PulseNumber=16; %回波脉冲数 BandWidth=2.0e6; %发射信号带宽带宽B=1/τ，τ是脉冲宽度TimeWidth=42.0e-6; %发射信号时宽 PRT=240e-6; % 雷达发射脉冲重复周期(s),240us对应 1/2*240*300=36000米最大无模糊距离 PRF=1/PRT; Fs=2.0e6; %采样频率

合成孔径雷达(SAR)

合成孔径雷达（SAR）合成孔径雷达（SAR）数据拥有独特的技术魅力和优势，渐成为国际上的研究热点之一，其应用领域越来越广泛。SAR数据可以全天候对研究区域进行量测、分析以及获取目标信息。高级雷达图像处理工具SARscape，能让您轻松将原始SAR数据进行处理和分析，输出SAR 图像产品、数字高程模型（DEM）和地表形变图等信息，应用永久散射体PS、短基线处理SBAS等方法快速准确地获取大范围形变信息，并可以将提取的信息与光学遥感数据、地理信息集成在一起，全面提升SAR数据应用价值。基本概念合成孔径雷达就是利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成一较大的等效天线孔径的雷达，也称综合孔径雷达。合成孔径雷达的特点是分辨率高，能全天候工作，能有效地识别伪装和穿透掩盖物。所得到的高方位分辨力相当于一个大孔径天线所能提供的方位分辨力。分类合成孔径雷达可分为聚焦型和非聚焦型两类。用在飞机上或空间飞行器上可有几种不同的工作模式,最常见的是正侧视模式,称为合成孔径侧视雷达；此外还有斜视模式、多普勒波束锐化模式和定点照射模式等。如果雷达保持相对静止，使目标运动成像,则成为逆合成孔径雷达，也称距离-多普勒成像系统。合成孔径雷达在军事侦察、测

绘、火控、制导，以及环境遥感和资源勘探等方面有广泛用途。发展概况合成孔径的概念始于50年代初期。当时，美国有些科学家想突破经典分辨力的限制，提出了一些新的设想：利用目标与雷达的相对运动所产生的多普勒频移现象来提高分辨力；用线阵天线概念证明运动着的小天线可获得高分辨力。50年代末，美国研制成第一批可供军事侦察用的机载高分辨力合成孔径雷达。60年代中期，随着遥感技术的发展，军用合成孔径雷达技术推广到民用方面，成为环境遥感的有力工具。70年代后期，卫星载合成孔径雷达和数字成像技术取得进展。美国于1978年发射的“海洋卫星”A号和80年代初发射的航天飞机都试验了合成孔径雷达的效果，证明了雷达图像的优越性。空中SAR概况 1. 1951年, Carl Wiley 首次提出利用频率分析方法改善雷达的角分辨率. 2. 1953年, 伊利诺依大学采用非聚焦方法使角度分辨率由4.13度提高到0.4度,并获得第一张SAR图像. 3. 1957年, 密西根大学采用光学处理方式, 获得了第一张全聚焦SAR图像. 4. 1978年, 美国发射了第一颗星载Seasat-1. 5. 1991年, 欧洲空间局发射了ERS-1. 6. 1995年, 加拿大发射了Radarsat-1.

附录(图书管理系统程序代码)

附录主页面（MDImain）： Private Sub addbook_Click() add_book.Show End Sub Private Sub addreader_Click() add_reader.Show End Sub Private Sub adduser_Click() add_user.Show End Sub Private Sub backbook_Click() back_book.Show End Sub Private Sub borrowbook_Click() borrow_book.Show End Sub Private Sub cuibackbook_Click() cuiback_book.Show End Sub Private Sub findbook_Click() find_book.Show End Sub Private Sub findreader_Click() find_reader.Show End Sub Private Sub MDIForm_Load() End Sub Private Sub modifybook_Click() change_book.Show End Sub Private Sub modifypwd_Click()

change_pwd.Show End Sub Private Sub modifyreader_Click() change_reader.Show End Sub Private Sub delbook_Click() del_book.Show End Sub Private Sub delreder_Click() del_reader.Show End Sub Private Sub quitsys_Click() End End Sub 登陆（login）： Option Explicit Dim cnt As Integer '记录确定次数 Private Sub Command1_Click() Dim sql As String Dim rs_login As New ADODB.Recordset If Trim(txtuser.Text) = "" Then '判断输入的用户名是否为空MsgBox "没有这个用户", vbOKOnly + vbExclamation, "" txtuser.SetFocus Else sql = "select * from 用户表where 用户名='" & txtuser.Text & "'" rs_login.Open sql, conn, adOpenKeyset, adLockPessimistic If rs_login.EOF = True Then MsgBox "没有这个用户", vbOKOnly + vbExclamation, "" txtuser.SetFocus Else '检验密码是否正确 If Trim(rs_login.Fields(1)) = Trim(txtpwd.Text) Then userID = txtuser.Text userpow = rs_login.Fields(2) rs_login.Close Unload Me MsgBox "欢迎登录到图书管理系统！", vbOKOnly + vbExclamation, "" MDImain.Show Else

图书馆管理系统程序的设计代码

1.1程序设计代码登录模块 if(username.Text.Trim()==""||password.Text.Trim()=="") MessageBox.Show("请输入用户名和密码","提示"); else { if (radioManage.Checked == true) { string strcon = "Data Source=SIMON-VAIO;Initial Catalog=lkl2;Integrated Security=True;"; //连接数据库的字符串，用于指定数据库地址，名称，账号，密码，连接方式SqlConnection sqlCon = new SqlConnection(strcon); //实例化并定义一个数据库连接 sqlCon.Open(); //打开数据库连接 string sql = "select * from login where usernum=usernum and userpassword=suerpassword"; //定义要查询sql语句 SqlCommand cmd = new SqlCommand(sql, sqlCon); //实例化并定义sql语句和数据库路径 cmd.Parameters.Add("usernum", SqlDbType.NChar, 20); //定义cmd查询命令的字段属性，loginname sqldbtype nchar（20） cmd.Parameters.Add("suerpassword", SqlDbType.NChar, 20); //同上 cmd.Parameters["usernum"].Value = username.Text; //将username中的text保存到变量loginname cmd.Parameters["suerpassword"].Value = password.Text; //同上 SqlDataReader dr = cmd.ExecuteReader(); if (dr.Read()) { this.Visible=false; Form2 Formmain = new Form2(); //应该是实例化一个主窗体的 this.Hide(); //应该是切换到主窗口的或关闭自己的 Formmain.Show(); //应该是打开一个主窗体的 dr.Close();//关闭dr的数据库连接 } else// if (dr.Read())读取失败则执行如下代码 MessageBox.Show("密码错误,请重新输入!"); //显示提示信息 } else if (radioPerson.Checked==true)

雷达信号matlab仿真

雷达系统分析大作作者：陈雪娣学号：0410420727 1. 最大不模糊距离： ,max 1252u r C R km f == 距离分辨率： 1502m c R m B ?= = 2. 天线有效面积： 22 0.07164e G A m λπ == 半功率波束宽度： 3 6.44o db G θπ == 3. 模糊函数的一般表示式为 () ()()2 2* 2 ;? ∞ ∞ -+= dt e t s t s f d f j d πττχ 对于线性调频信号 ()21 j t p p t s t ct e T T πμ??= ? ??? 则有： ()()2 21 ;Re Re p j t T j t d p p p t t f ct ct e e dt T T T πμπμτ χτ∞+-∞????+= ? ? ? ????? ? () ()()sin 1;11d p p d p d p p f T T f T f T T τπμττχττπμτ????+- ? ? ? ???????=- ? ?????+- ? ? ? ? 分别令0,0==d f τ可得()()2 2 0;,;0τχχd f ()() sin 0;d p d d p f T f f T πχπ=

()sin 1 ;01 1p p p p p T T T T T τπμττχττπμτ?? ??- ? ? ? ???????=- ? ?????- ? ?? ? 程序代码见附录1的T_3.m, 仿真结果如下:

4. 程序代码见附录1的T_4.m, 仿真结果如下:

图书馆管理系统源代码

源程序清单 1、文件名 login(login.frm) 功能说明：整个系统的登陆界面，需要输入用户名和登陆密码才能进入到系统中，进行借阅等操作。源代码： Option Explicit Dim cnt As Integer Private Sub Command1_Click() Dim sql As String Dim rs_login As New ADODB.Recordset If Trim(Combo1.Text) = "" Then MsgBox "没有这个用户", vbOKOnly + vbExclamation, "" Combo1.SetFocus Else sql = "select * from 系统管理 where 用户名='" & Combo1.Text & "'" rs_login.Open sql, conn, adOpenKeyset, adLockPessimistic If rs_login.EOF = True Then MsgBox "没有这个用户", vbOKOnly + vbExclamation, "" Combo1.SetFocus Else If Trim(rs_login.Fields(1)) = Trim(txtpwd.Text) Then userID = Combo1.Text rs_login.Close Unload Me form1.Show Else MsgBox "密码不正确", vbOKOnly + vbExclamation, "" txtpwd.SetFocus End If

End If cnt = cnt + 1 If cnt = 3 Then Unload Me End If Exit Sub End Sub Private Sub Command2_Click() Unload Me End Sub Private Sub Form_Load() Dim connectionstring As String connectionstring = "provider=Microsoft.Jet.oledb.4.0;" & _ "data source=book.mdb" conn.Open connectionstring cnt = 0 End Sub Private Sub txtuser_Change() End Sub 2、文件名 form1(form1.frm) 功能说明：整个系统的主界面，其中包括图书管理、读者管理、图书借阅管理、系统管理、关于，以及在这下面的子菜单。源代码： Private Sub add_admin_Click() frmadduser.Show End Sub Private Sub add_back_book_Click() frmbackbookinfo.Show

合成孔径雷达点目标仿真MATLAB程序

合成孔径雷达成像系统点目标仿真源程序： clc close all C=3e8; %光速 Fc=1e9; %载波频率 lambda=C/Fc; %波长 %成像区域 Xmin=0; Xmax=50; Yc=10000; Y0=500; %SAR基本参数 V=100; %雷达平台速度 H=0; %雷达平台高度 R0=sqrt(Yc^2+H^2); D=4; %天线孔径长度 Lsar=lambda*R0/D; %合成孔径长度 Tsar=Lsar/V; %合成孔径时间 Ka=-2*V^2/lambda/R0;%线性调频率 Ba=abs(Ka*Tsar); PRF=2*Ba; %脉冲重复频率 PRT=1/PRF; ds=PRT; %脉冲重复周期 Nslow=ceil((Xmax-Xmin+Lsar)/V/ds);%脉冲数 Nslow=2^nextpow2(Nslow); %量化为2的指数 sn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow); %创建时间向量PRT=(Xmax-Xmin+Lsar)/V/Nslow; %更新 PRF=1/PRT; % 更新脉冲重复频率 fa=linspace(-0.5*PRF,0.5*PRF,Nslow); Tr=5e-6; %脉冲宽度 Br=30e6; %调频信号带宽 Kr=Br/Tr; %调频率 Fsr=2*Br; %快时间域取样频率 dt=1/Fsr; %快时间域取样间隔

Rmin=sqrt((Yc-Y0)^2+H^2); Rmax=sqrt((Yc+Y0)^2+H^2+(Lsar/2)^2); Nfast=ceil(2*(Rmax-Rmin)/C/dt+Tr/dt); Nfast=2^nextpow2(Nfast); tm=linspace(2*Rmin/C,2*Rmax/C+Tr,Nfast); dt=(2*Rmax/C+Tr-2*Rmin/C)/Nfast; %更新 Fsr=1/dt; fr=linspace(-0.5*Fsr,0.5*Fsr,Nfast); DY=C/2/Br; %距离分辨率 DX=D/2; %方位分辨率 Ntarget=3; %目标数目 Ptarget=[Xmin,Yc,1 %目标位置 Xmin,Yc+10*DY,1 Xmin+20*DX,Yc+50*DY,1]; K=Ntarget; %目标数目 N=Nslow; %慢时间采样数 M=Nfast; %快时间采样数 T=Ptarget; %目标位置 %合成孔径回波仿真 Srnm=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

图书管理系统的C 代码完整版

C#代码清单共1个项目，包含5个类。项目：librarysystem 类: 主类代码： namespace librarysystem { PublicationISBN = PublicationISBN; = PublicationName; = PublicationAuthor; = PublicationType; = PublicationStatus; } } } 读者类代码： namespace librarysystem { 询图书"); ("2.借书服务"); ("3.还书服务"); ("4.新出版物入库"); ("5.录入新读者"); ("6.查询读者的借阅信息"); ("7.退出管理系统"); ("-------------------------------------------------------------------"); bool flag; do { flag = false; ("请选择："); string key = (); switch (key) { case"1": ("请输入需要查询图书的书名："); ShowPubInfo(); break; case"2": ("请输入您的姓名：");

(); ShowBorrowInfo(); break; case"3": ("请输入您的姓名："); (); ShowReturnInfo(); break; case"4": ("请按照提示逐步输入需要入库出版物的信息……"); ShowAddPublications(); break; case"5": ("请按照提示逐步输入新注册读者的信息……"); ShowAddReaders(); break; case"6": ("请输入您的姓名： "); (); ShowRedPubInfo; break; case"7": flag = false; break; default: ("无此业务，是否重新选择（y/n）？"); string answer = (); if (answer != "y") { flag = false; } else { flag = true; } break; } } while (flag); } /*显示查询到的出版物信息*/ public void ShowPubInfo() {

图书管理系统程序代码

#i n c l u d e<> #include <> #include <> #define SIZE (struct booklist *)malloc(sizeof(struct booklist)) .(Y/N) :"); ch = getchar(); if(ch == 'y' || ch == 'Y') store(head); getchar(); break; } case 2: { if(head == NULL ) { printf("请先录入图书源信息!\n"); getchar(); break; } else { borrow_head = borrow_creat(head); .(Y/N) :" ); ch = getchar(); if(ch == 'y' || ch == 'Y') borstore(borrow_head); getchar(); break; } } case 3: { if(head == NULL ) { printf("请先录入图书源信息!\n"); getchar(); break; } else { flag = pre_creat(head); if(flag == 1)

printf("没有录入预借图书信息!\n\n"); if(flag == 2) pre_head = prebor_creat(); .(Y/N) :"); ch = getchar(); if(ch == 'y' || ch == 'Y') store(head); getchar(); break; } } case 4: .(Y/N) :"); ch = getchar(); if(ch == 'y' || ch == 'Y') store(head); getchar(); break; } } case 5: .(Y/N) :"); ch = getchar(); if(ch == 'y' || ch == 'Y') store(head); getchar(); break; } } case 9: .)\n\n"); printf("请选择功能:"); scanf("%d", &n); getchar(); if(n == 0) { x = 0; } getchar(); } else { printf("请输入学号:"); gets(c); strcpy(p -> prebook_num, a); strcpy(p -> prestu_num, c); if((fp = fopen("","a+")) == NULL)

合成孔径雷达

合成孔径雷达（SAR）合成孔径雷达产生的过程为了形成一幅真实的图像增加两个关键参数：分辨率、识别能力。合成孔径打开了无限分辨能力的道路相干成像特性：以幅度和相位的形式收集信号的能力相干成像的特性可以用来进行孔径合成民用卫星接收系统SEASA T、SIR-A、SIR-B 美国军用卫星(LACROSSE) 欧洲民用卫星（ERS系列）合成孔径雷达(SAR)是利用雷达与目标的相对运动将较小的真实天线孔径用数据处理的方法合成一个较大孔径的等效天线孔径的雷达。特点：全天候、全天时、远距离、和高分辨率成像并且可以在不同频段不同极化下得到目标的高分辨率图像 SAR高分辨率成像的距离高分辨率和方位高分辨率距离分辨率取决于信号带宽方位高分辨率取决于载机与固定目标相对运动时产生的具有线性调频性质的多普勒信号带宽相干斑噪声机载合成孔径雷达是合成孔径雷达的一种极化：当一个平面将空间划分为各向同性和半无限的两个均匀介质，我们就可以定义一个电磁波的入射平面，用波矢量K来表征：该平面包含矢量K以及划分这两种介质的平面法线垂直极化（V）：无线电波的振动方向是垂直方向与水平极化（H）：无线电波的振动方向是水平方向 TE波：电场E与入射面垂直

TH波：电场E属于入射平面合成孔径雷达的应用军事上、地质和矿物资源勘探、地形测绘和制图学、海洋应用、水资源、农业和林业合成孔径雷达在军事领域的应用：战略应用、战术应用、特种应用。 SAR系统的几个发展趋势：多波段、多极化、多视角、多模式、多平台、高分辨率成像、实时成像。 SAR图像相干斑抑制的研究现状分类：成像时进行多视处理、成像后进行滤波多视处理就是对同一目标生成多幅独立的像，然后进行平均。这是最早提出的相干斑噪声去除的方法，这种技术以牺牲空间分辨率为代价来获取对斑点的抑制成像后的滤波技术成为SAR图像相干噪声抑制技术发展的主流均值滤波、中值滤波、维纳滤波用来滤去相干斑噪声，这种滤波方法能够在一定程度上减小相干斑噪声的方差合成孔径雷达理论概述合成孔径雷达是一种高分辨率成像雷达，高分辨率包含两个方面的含义：方位向的高分辨率和距离向高分辨率。它通过采用合成孔径原理提高雷达的方位分辨率，并依靠脉冲压缩技术提高距离分辨率由于SAR雷达发射信号（距离向信号）和合成孔径信号（方位信号）均具有线性调频性质，SAR成像的实质就是通过匹配滤波器对距离向和方位向具有线性调频信号的信号进行二维脉冲压缩的过程，也就是依靠脉冲压缩技术提高距离分辨率，通过合成孔径原理提高雷达的方位分辨率的过程 SAR成像处理是先利用距离向匹配滤波器，进行距离脉压，实现距离向高分辨率后，再通过方位向德匹配滤波，最终得到原始目标的高分辨图像。

图书管理系统程序代码

模块设计基础模块 SqlSession: 利用mybatis获取jdbc,其中有inter.select,update,select等方法 String resource = "mybatis-config.xml"; InputStream inputStream = Resources.getResourceAsStream(resource); SqlSessionFactory sqlSessionFactory = new SqlSessionFactoryBuilder().build(inputStream); SqlSession sesson = sqlSessionFactory.openSession(); Dao基础模块所有的dao都继承于BaseDao; 代码 /** *保存利用mybatis的xml * model 保存的对象 *id mybatis的id名 */ Public void Save(BaseModel model,String id){ sesson.inster(id model); } …… Model基础模块所有的模型都继承于BaseModel; 代码片段 /** *字段编号 */ Private int id; Public void setId(String id){ This.id = id; } Public String getId(){

Return id; } 添加书籍代码片段 Book book = new Book(); book.setId(…..); …….. Book.setTitle(….); BooksManageDao bookDao = new BooksManageDao(); bookDao .save(“addBook”,book); 删除书籍代码片段 BookDao bookDao = new BookDao(); Book book = (Book) bookDao .getModel(id); BooksManageDao bookDao = new BooksManageDao(); bookDao .delete(“delete Book”,book); 更新书籍

雷达系统建模与仿真报告

设计报告一十种随机数的产生一概述. 概论论是在已知随机变量的情况下，研究随机变量的统计特性及其参量，而随机变量的仿真正好与此相反，是在已知随机变量的统计特性及其参数的情况下研究如何在计算机上产生服从给定统计特性和参数随机变量。下面对雷达中常用的模型进行建模： ● 均匀分布 ● 高斯分布 ● 指数分布 ● 广义指数分布 ● 瑞利分布 ● 广义瑞利分布 ● Swerling 分布 ● t 分布 ● 对数一正态分布 ● 韦布尔分布二随机分布模型的产生思想及建立. 产生随机数最常用的是在(0,1)区间内均匀分布的随机数，其他分布的随机数可利用均匀分布随机数来产生。 2.1 均匀分布 1>（0，1）区间的均匀分布：用混合同余法产生 (0,1)之间均匀分布的随机数，伪随机数通常是利用递推公式产生的，所用的混和同余法的递推公式为: 1 n x =n x +C （Mod m ）

其中，C是非负整数。通过适当选取参数C可以改善随机数的统计性质。一般取作小于M的任意奇数正整数，最好使其与模M互素。其他参数的选择 (1) 的选取与计算机的字长有关。 (2) x(1)一般取为奇数。用Matlab来实现，编程语言用Matlab语言，可以用 hist 函数画出产生随机数的直方图（即统计理论概率分布的一个样本的概率密度函数），直观地看出产生随机数的有效程度。其产生程序如下： c=3;lamade=4*200+1; x(1)=11; M=2^36; for i=2:1:10000; x(i)=mod(lamade*x(i-1)+c,M); end; x=x./M; hist(x,10); mean(x) var(x) 运行结果如下：均值 = 0.4948 方差 = 0.0840 2> （a,b）区间的均匀分布：利用已产生的（0，1）均匀分布随机数的基础上采用变换法直接产生（a,b）

现役合成孔径雷达2014.11

现役全球InSAR卫星简介

1、德国TanDEM-X卫星德国TanDEM-X任务是利用两颗TerraSAR-X卫星进行编队飞行的一个高精度的雷达干涉测量系统，第一颗TerraSAR-X卫星于2007年发射升空，计划使用寿命为5年，第二颗TerraSAR-X卫星于2009年发射升空，计划使用寿命为5年，两颗卫星有三年的工作交叠期，德国预计在这三年中生成全球的高精度DEM数字高程模型，高程定位精度优于2m，DEM网格间距为12m。表2给出了以上几种星载干涉系统在不同基线的情况下高程模糊度的具体数值不同波段下系统干涉的性能比较 L、X、C波段所生成的SAR图像有其各自的特点，高程信息的精度主要取决于雷达波长和相干系数。对于同一区域的SAR图像干涉处理，L波段的图像相干性高于X、C波段的图像，但是就高程信息的敏感度，X、C波段优于L波段。

2、加拿大Radarsat-2雷达系统 Radarsat-2是加拿大第二代地球观测卫星，于2006年12月发射升空，它几乎保留了Radarsat-1的所有优点，雷达采用C波段，HH极化，数据分辨率3—100m，幅宽10—500km,设计使用寿命为7年，采用多极化工作模式，轨道定位精度15m。能够大大增加可识别地物或目标的类别，能够左视和右视，并且可以实现相互转换，主要用于测绘以及环境和自然资源的检测等方面。 3、日本ALOS观测卫星 2006年1月日本发射了先进陆地观测卫星(ALOS),它携带有L波段相控阵合成孔径雷达(PALSAR)，该卫星主要用于对全球陆地资源和环境实施全天候监测，在高分辨率模式下距离向分辨率优于2m，轨道定位精度10m。PALSAR有较高的距离向分辨率和较高的信噪比，并且在交轨方向对轨道有较好的控制。

合成孔径雷达(SAR)的点目标仿真(附件带代码程序)

最新多普勒雷达系统仿真

合成孔径雷达概述(SAR)

C++图书管理系统源代码

雷达系统仿真matlab代码.docx

合成孔径雷达(SAR)

附录(图书管理系统程序代码)

图书馆管理系统程序的设计代码

雷达信号matlab仿真

图书馆管理系统源代码

合成孔径雷达点目标仿真MATLAB程序

图书管理系统的C 代码 完整版

图书管理系统程序代码

合成孔径雷达

图书管理系统程序代码

雷达系统建模与仿真报告

现役合成孔径雷达2014.11

图书管理系统的C 代码完整版