(完整版)合成孔径雷达点目标仿真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

end

%距离压缩

tr=tm-2*Rmin/C;

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

F_Refr=fft((Refr));

Sr=zeros(N,M);

for k2=1:1:M

temp1=fft(Srnm(k2,:));

FSrnm=temp1.*conj(F_Refr);

Sr(k2,:)=ifft(FSrnm);

end

%方位压缩

ta=sn-(Xmin+Xmax)/2/V;

Refa=exp(j*pi*Ka*ta.^2).*(abs(ta)

Sa=zeros(N,M);

for k3=1:1:M

temp2=fft(Sr(:,k3));

F_Sa=temp2.*conj(F_Refa.');

Sa(:,k3)=fftshift(ifft(F_Sa));

end

row=tm*C/2;

col=sn*V;

%绘图

%回波雷达数据

figure(1)

subplot(211)

imagesc(abs(Srnm));title('SAR data')

subplot(212)

imagesc(angle(Srnm))

%灰度图

figure(2)

colormap(gray)

imagesc(row,col,255-abs(Sr));

title('距离压缩'),xlabel('距离向'),ylabel('方位向');

figure(3)

colormap(gray)

imagesc(row,col,255-abs(Sa));

title('方位压缩'),xlabel('距离向'),ylabel('方位向');

%轮廓图

figure(4)

Ga=abs(Sa);

a=max(max(Ga));

contour(row,col,Ga,[0.707*a,a],'b');grid on

x1=sqrt(H^2+min(abs(Ptarget(:,2))).^2)-5*DY;

x2=sqrt(H^2+max(abs(Ptarget(:,2))).^2)+5*DY;

y1=min(Ptarget(:,1))-5*DX;

y2=max(Ptarget(:,1))+5*DX;

axis([x1,x2,y1,y2])

%三维成像图

figure(5)

mesh(Ga((400:600),(200:500)));axis tight

xlabel('Range'),ylabel('Azimuth');

仿真结果图：

图4. 1 点目标原始回波数据

图4. 2 距离压缩后成像图

图4. 3 方位压缩后成像图（点阵目标成像灰度图）

图4. 4 点阵目标成像轮廓图

合成孔径雷达概述(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)

MATLAB仿真实验报告

MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件

MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays

实验二 Matlab程序设计基本方法1

实验二Matlab程序设计基本方法 覃照乘自092 电气工程学院 一、实验目的： 1、熟悉MATLAB 程序编辑与设计环境 2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计 4、了解和熟悉跨空间变量传递和赋值 二、实验基本知识： ◆for循环结构 语法：for i＝初值：增量：终值 语句1 …… 语句n end 说明：1．i＝初值：终值，则增量为1。 2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。 ◆while 循环结构 语法：while 逻辑表达式 循环体语句 end 说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达 式的值为假，则程序执行end之后的语句。 2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位 置加break语句、以便程序能正常执行。（执行循环体的次数不确定； 每一次执行循环体后，一定会改变while后面所跟关系式的值。） 3、while循环也可以嵌套、其结构如下：

while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end ◆if-else-end分支结构 if 表达式1 语句1 else if 表达式2(可选) 语句2 else(可选) 语句3 end end 说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。 2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。 3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。 ◆switch-case结构 语法：switch表达式 case常量表达式1 语句组1 case常量表达式2 语句组2 …… otherwise 语句组n end

合成孔径雷达(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.

MATLAB仿真实验全部

实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些 2、 如何判断系统稳定性 3、 系统的动态性能指标有哪些 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间范围0---Tn 。 3、),(T sys step ；表示时间范围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：) ()()()(1)(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.

实验3MATLAB程序设计

1，编写M 函数实现求一个数是否为素数，再编写一主程序（脚本文件），要求通过键盘输入一个整数，然后调用判断素数函数，从而确定它是否素数。 x=input('请输入一个整数x:'); if myprime(x) disp('您输入的整数x是一个素数。') else disp('您输入的数x不是一个素数。') end function y=myprime(x) y=1; for i=2:fix(sqrt(x)) if mod(x,i)==0 y=0; end end 2，编写M 函数统计一数值中零的个数，然后编写脚本文件，实现统计从1—2007 中零的总个数。 function num=number0(a) %统计十进制数值中0的个数 sa=num2str(a);%将数值装化为字符串 num=length(find(sa=='0'));% ));%求取字符串中'0’的个数 y=0;

for a=1:2006 num=number0(a); y=num+y; end disp(y) 504 3，编写程序计算x∈[-3,3],字长0.01：并画出曲线x = -3:0.01:3; y=zeros(size(x)); for i = 1:length(x) if -3<= x(i)& x(i)<=-1 y(i)=(-x(i).^2-4*x(i)-3)/ 2; elseif -1<= x(i) & x(i)<=1 y(i)=-x(i).^2+1; elseif 1<=x(:,i)<=3 y(i)=(-x(i).^2+4*x(i)-3)/2; end end plot(x,y) -3-2-10123

刘卫国版MATLAB程序设计与应用课后实验六八九

实验六 高层绘图操作 %第一题： 程序代码如下： x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.^2)).*cos(x); plot(x,y) 01234567 -1 -0.5 0.5 1 1.5 %第二题： %（1） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'y--'); text(4,16,'\leftarrow y1=x^2'); text(6*pi/4,-1,'\downarrow y2=cos(2*x)'); text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2');

-8 -6 -4 -2 2 4 6 8 -30-20 -10 10 20 30 40 %（2） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(1,3,1);%分区 plot(x,y1); title('y1=x^2');%设置标题 subplot(1,3,2); plot(x,y2); title('y2=cos(2*x)'); subplot(1,3,3); plot(x,y3); title('y3=x^2*cos(2*x)');

-10 10 0510 15202530 35 40y1=x 2 -10 10 -1-0.8 -0.6 -0.4-0.200.20.4 0.6 0.8 1y2=cos(2*x) -10 10 -30-20 -10 10 20 30 40 y3=x 2*cos(2*x) %（3） 程序代码如下： x=linspace(-2*pi,2*pi,20); y1=x.^2; subplot(2,2,1);%分区 bar(x,y1); title('y1=x^2的条形图');%设置标题 subplot(2,2,2); stairs(x,y1); title('y1=x^2的阶梯图'); subplot(2,2,3); stem(x,y1); title('y1=x^2的杆图'); subplot(2,2,4); fill(x,y1,'r');%如果少了'r'则会出错 title('y1=x^2的填充图'); %其他的函数照样做。

合成孔径雷达点目标仿真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

增量调制MATLAB仿真实验

增量调制MATLAB仿真实验

增量调制(DM)实验 一、实验目的 （1）进一步掌握MATLAB的应用。 （2）进一步掌握计算机仿真方法。 （3）学会用MATLAB软件进行增量调制(DM)仿真实验。 二、实验原理 增量调制是由PCM发展而来的模拟信号数字化的一种编码方式，它是PCM的一种特例。增量调制编码基本原理是指用一位编码，这一位码不是表示信号抽样值的大小，而是表示抽样幅度的增量特性，即采用一位二进制数码“1”或“0”来表示信号在抽样时刻的值相对于前一个抽样时刻的值是增大还是减小，增大则输出“1”码，减小则输出“0”码。输出的“1”，“0”只是表示信号相对于前一个时刻的增减，不表示信号的绝对值。 增量调制最主要的特点就是它所产生的二进制代码表示模拟信号前后两个抽样值的差别(增加、还是减少)而不是代表抽样值本身的大小，因此把它称为增量调制。在增量调制系统的发端调制后的二进制代码1和0只表示信号这一个抽样时刻相对于前一个抽样时刻是增加(用1码)还是减少(用0码)。收端译码器每收到一个1码，译码器的输出相对于前一个时刻的值上升一个量化阶，而收到一个0码，译码器的输出相对于前一个时刻的值下降一个量化阶。 增量调制(DM)是DPCM的一种简化形式。在增量调制方式下,采用1比特量化器,即用1位二进制码传输样值的增量信息,预测器是

一个单位延迟器,延迟一个采样时间间隔。预测滤波器的分子系数向量是[0,1],分母系数为1。当前样值与预测器输出的前一样值相比较,如果其差值大于零,则发1码,如果小于零则发0码。 三、实验内容 增量调制系统框图如图一所示,其中量化器是一个零值比较器,根据输入的电平极性,输出为 δ,预测器是一个单位延迟器,其输出为前一个采样时刻的解码样值,编码器也是一个零值比较器,若其输入为负值,则编码输出为0,否则输出为1。解码器将输入1,0符号转换为 δ,然后与预测值相加后得出解码样值输出,同时也作为预测器的输入 输入样值 e n e n =δsgn(e n ) 传输 n ) n n-1+δsgn(e n ) x n + - + + 预测输出 + n-1 + 预测输出 解码样值输出 x n-1 预测输入x n =x n-1+δsgn(e n ) 图一 增量调制原理框图 设输入信号为: x(t)=sin2π50t+0.5sin 2π150t 增量调制的采样间隔为1ms,量化阶距δ=0.4,单位延迟器初始值为0。建立仿真模型并求出前20个采样点使客商的编码输出序列以 解码 编码 二电平量化 单位延迟 单位 延迟

MATLAB程序设计实验报告

MATLAB 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法； 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 4. 通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 ： 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 （1）信号的相加和相乘：已知信号)(1t f 和)(2t f ，信号相加和相乘记为 )()(1t f t f =)(2t f +；)()(1 t f t f =)(2t f *。 （2）信号的微分和积分：对于连续时间信号，其微分运算是用diff 函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function 表示需要进行求导运算的信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；n 为求导的阶数，默认值为求一阶导数。连续信号的积分运算用int 函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function 表示需要进行被积信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；a,b 为积分上、下限，a 和b 省略时为求不定积分。 （3）信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化，可使信号压缩或扩展。)(at f 将原波形压缩a 倍，)/(a t f 将原波形扩大a 倍。 ￥ 对系统频率特性的分析

高频电子线路Matlab仿真实验

高频电子线路Matlab 仿真实验要求 1. 仿真题目 （1） 线性频谱搬移电路仿真 根据线性频谱搬移原理，仿真普通调幅波。 基本要求：载波频率为8kHz ，调制信号频率为400Hz ，调幅度为0.3；画出调制信号、载波信号、已调信号波形，以及对应的频谱图。 扩展要求1：根据你的学号更改相应参数和代码完成仿真上述仿真；载波频率改为学号的后5位，调制信号改为学号后3位，调幅度设为最后1位/10。（学号中为0的全部替换为1，例如学号2010101014，则载波为11114Hz ，调制信号频率为114，调幅度为0.4）。 扩展要求2：根据扩展要求1的条件，仿真设计相应滤波器，并获取DSB-SC 和SSB 的信号和频谱。 （2） 调频信号仿真 根据调频原理，仿真调频波。 基本要求：载波频率为30KHz ，调制信号为1KHz ，调频灵敏度32310f k π=??，仿真调制信号，瞬时角频率，瞬时相位偏移的波形。 扩展要求：调制信号改为1KHz 的方波，其它条件不变，完成上述仿真。 2. 说明 （1） 仿真的基本要求每位同学都要完成，并且记入实验基本成绩。 （2） 扩展要求可以选择完成。

1.0 >> ma = 0.3; >> omega_c = 2 * pi * 8000; >> omega = 2 * pi * 400; >> t = 0 : 5 / 400 / 1000 : 5 / 400; >> u_cm = 1; >> fc = cos(omega_c * t); >> fa = cos(omega * t); >> u_am = u_cm * (1 + fa).* fc; >> U_c =fft(fc,1024); >> U_o =fft(fa,1024); >> U_am =fft(u_am, 1024); >> figure(1); >> subplot(321);plot(t, fa, 'k');title('调制信号');grid;axis([0 2/400 -1.5 1.5]); >> subplot(323);plot(t, fc, 'k');title('高频载波');grid;axis([0 2/400 -1.5 1.5]); >> subplot(325);plot(t, u_am, 'k');title('已调信号');grid;axis([0 2/400 -3 3]); >> fs = 5000; >> w1 = (0:511)/512*(fs/2)/1000; >> subplot(322);plot(w1, abs([U_am(1:512)']),'k');title('调制信号频谱');grid;axis([0 0.7 0 500]); >> subplot(324);plot(w1, abs([U_c(1:512)']),'k');title('高频载波频谱');grid;axis([0 0.7 0 500]); >> subplot(326);plot(w1, abs([U_am(1:512)']),'k');title('已调信号频谱');grid;axis([0 0.7 0 500]); 1.1 >> ma = 0.8; >> omega_c = 2 * pi * 11138; >> omega = 2 * pi * 138; >> t = 0 : 5 / 400 / 1000 : 5 / 400; >> u_cm = 1; >> fc = cos(omega_c * t);

Matlab程序设计实验报告

实验七Matlab 程序设计 实验目的： 1、掌握建立和执行M 文件的方法； 2、掌握实现选择结构的方法； 3、掌握实现循环结构的方法。 实验内容： 1. 编写用 5 次多项式拟合函数y=sin(x), x [0, 2 ]的脚本M 文件，要求绘图观察拟合的效果。 function shiyan1 x=0:0.5:2*pi y=sin(x) p=polyfit(x,y,5) x1=0:0.2:2*pi y1=polyval(p,x1) plot(x,y, 'b' ,x1,y1, '*r' x =

Columns 1 through 9 0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 Columns 10 through 13 4.5000 5.0000 5.5000 6.0000 y = Columns 1 through 9 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411 -0.3508 -0.7568 Columns 10 through 13 -0.9775 -0.9589 -0.7055 -0.2794 p = -0.0056 0.0881 -0.3967 0.2671 0.8902 0.0029 x1 = Columns 1 through 10 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 Columns 11 through 20

2. 2.2000 2.4000 2.6000 2.8000 3.0000 3.2000 3.4000 3.6000 1.8001 Columns 21 through 30 4.0 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 5.6000 5.8000 Columns 31 through 32 6.0 6.2000 y1 = Columns 1 through 10 0.29 0.1886 0.3786 0.5585 0.7172 0.8461 0.9391 0.9926 1.0048 0.9761 Columns 11 through 20 0.9083 0.8048 0.6701 0.5098 0.3301 0.1381 -0.0590 -0.2538 -0.4389 -0.6073 Columns 21 through 30 -0.7524 -0.8685 -0.9505 -0.9949 -0.9991 -0.9626 -0.8863 -0.7732 -0.6288 -0.4606 Columns 31 through 32

合成孔径雷达

合成孔径雷达（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成像处理是先利用距离向匹配滤波器，进行距离脉压，实现距离向高分辨率后，再通过方位向德匹配滤波，最终得到原始目标的高分辨图像。

Matlab 编程方法及仿真实验

《现代机械工程基础实验》之机械工程控制基础综合实验报告 姓名 学号 班级 山东建筑大学机电工程学院 2012.06.04~06

第一部分 Matlab 编程方法及仿真实验 实验1. 三维曲面的绘制（略） 实验2. 系统零极点绘制例：求部分分式展开式和)(t g 一个线性定常系统的传递函数是 1 5422 3)(2 3 ++++= s s s s s G （1） 使用MATLAB 建立传递函数，并确定它的极点和零点，写出)(s G 的部分分式展开式并绘制 系统的脉冲响应。 实验结果：零点-0.6667 极点-0.8796 + 1.1414i -0.8796 - 1.1414i -0.2408 实验3. 系统的阶跃响应 例. )(s G 的阶跃响应 对例2中由（1）式给出的传递函数)(s G ，增加一个0=s 处的极点，使用impulse 命令绘制其拉普拉斯反变换式曲线，得到阶跃响应图。将该响应与对)(s G 使用step 命令所得到的响应比较，确定系统的DC 增益。利用初值定理和终值定理来校验结果。 实验结果：DC 增益= 2

实验4. 双输入反馈系统单位阶跃响应 考虑一个如图1所示的反馈系统，它既有参考输入也有干扰输入，其中对象和传感器的传递函数是 )12)(15.0(4)(++=s s s G p ，105.01 )(+=s s H 控制器是一个增益为80，有一个在3-=s 处的零点，极点/零点比15=α超前控制器。推导 两个独立的MATLAB 模型，其中一个模型的输入为)(s R ，另一个输入为)(s D 。使用这些模型确定闭环零点和极点，并在同一坐标系内绘制它们的阶跃响应。 D (s ) 图1 具有参考和干扰输入的反馈系统方框图 实验结果： 参考输入的CL 极点：-49.3658 -7.3336 + 7.9786i -7.3336 - 7.9786i -3.4670 参考输入的DC 增益：320 干扰输入的CL 零点：-45 干扰输入的CL 极点：-49.3658 -7.3336 + 7.9786i -7.3336 - 7.9786i -3.4670 干扰输入的DC 增益：4 -20

D实验五 M文件和MATLAB程序设计

实验五M文件和MATLAB程序设计 一、实验目的 matlab作为一种高级计算机语言，不仅可以命令行方式完成操作，也具有数据结构、控制流、输入输出等能力，本次实验通过熟悉和掌握m文件的建立与使用方法，以及函数与控制程序流程语句的使用，使学生具备一定的编程和程序调试能力。 1．掌握M文件的使用方法。 2．掌握if语句和switch语句的使用 3. 掌握循环语句的使用 4. 通过练习理解MATLAB编程方法。 二、实验原理 1．m文件 用matlab语言编写的程序，称为m文件。M文件根据调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file）。区别? 2．程序控制结构 1）顺序结构 2）选择结构 （1）if语句a) 单分支if语句b) 双分支if语句c) 多分支if语句 （2）switch 语句 （3）try语句 3）循环结构 （1）for 语句 （2）while语句 （3）break语句、continue语句、return使用,区别? 3．函数文件 function 输出形参表＝函数名（输入形参表） 注释说明部分 函数体语句 注意事项？ 三、实验要求 1．首先上机练习PPT中各种流程控制语句的有关实例。 2．然后上机练习下面的实验习题。 四、实验习题

1．数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1 运行下面的程序，按程序提示输入n=1,2,3,5,7，8,9等数来验证这一结论。 %classic "3n+1" problem from number theory. while 1 n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end a end Enter n,negative quits: 2． 编程求满足∑=>m i i 11000020的最小m 值。 a=0; i=1; while (a<100000) a=a+pow2(i); i=i+1; end m=i-1 3． 编写一个函数，计算下面函数的值，给出x 的值，调用该函数后，返回y 的值。 function [y]=myfun1(x)

现役合成孔径雷达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有较高的距离向分辨率和较高的信噪比，并且在交轨方向对轨道有较好的控制。

Matlab仿真实验教程

MATLAB的实验仿真 目录 实验一MATLAB在控制系统模型建立与仿真中的应用 (1) 实验二典型系统的时域响应分析 (13) 实验三线性控制系统的根轨迹与频域分析 (17) 实验四线性系统的校正 (22) 附录一 MATLAB6.5 控制系统工具箱函数和结构化的控制语句 (30) 附录二 SIMULINK 基本模块介绍 (34)

实验一 MATLAB 在控制系统模型建立与仿真中的应用 一、 MATLAB 基本操作与使用 1. 实验目的 1) 掌握MATLAB 仿真软件的安装及启动，熟悉 MATLAB 工作环境平台。 2) MATLAB 命令窗口，包括工具条以及菜单选项的使用；MATLAB 语言的基本规定，包括数值的表示、变量命名规定、基本运算符、预定义变量以及表达式等。 3) MATLAB 图形绘制功能、M 文件程序设计和线性控制系统传递函数模型的建立等。 2. 实验仪器 PC 计算机一台，MATLAB 软件1套 3. 实验内容 1) MATLAB 的启动 这里介绍MATLAB 装入硬盘后，如何创建MATLAB 的工作环境。 方法一 MATLAB 的工作环境由matlab.exe 创建，该程序驻留在文件夹matlab\bin\ 中。它的图标是 matlab 。只要从<我的电脑>或<资源管理器>中去找这个程序，然后双击此图标，就会自动创建如图1所示的MATLAB6.5 版的工作平台 。 图1 在英文Windows 平台上的MATLAB6.5 MATLAB 工作平台 方法二 假如经常使用MATLAB ，则可以在Windows 桌面上创建一个MATLAB 快捷方式图标。具体办法为： 把<我的电脑>中的 matlab 图标用鼠标点亮，然后直接把此图标拖到Windows 桌面上即可。此后，直接双击Windows 桌面上的matlab 图标，就可建立图1所示的 MATLAB 工作平台。 2) MATLAB 工作环境平台 桌面平台是各桌面组件的展示平台，默认设置情况下的桌面平台包括 6 个窗口，具体如下： ① MATLAB 窗口 Command Window

合成孔径雷达概述

合成孔径雷达概述 蔡 Beautyhappy521@https://www.doczj.com/doc/ab12279015.html, 二OO八年三月二十三

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