- 博客(24)
- 资源 (3)
- 收藏
- 关注
RSS订阅原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章补充 如何做正确的fft）
说明： （1）main1.m以高斯信号为例，时间为[-T/2,T/2)，做fft之前需要做fftshift （2）main2.m以高斯信号为例，时间为[0,T)，做fft之前不需要做fftshift （3）为什么main1和main2一个做了fftshift，一个没有做fftshift，结果都与理论结果一致呢？这主要是因为DFT公式是以t=0时刻为第一个数值点，所以关键点是以t=0为始进行求解，因此信号时间为[-T/2,T/2)时需要做fftshift(fftshift之后第一个信号的时间就是t=0了)，而
2020-10-22 23:13:37
483
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第四章 - 图4.4 ）
matlab代码如下： %% initial matlab workspace % 图4.4 : page 81 clc clear close all % 参数设置 Re = 6368e3; % 本地地球半径 h = 800e3; % 平台距地面高度 dr = 13.6 ; % 斜距分辨率 % 设置要扫描的所有地距 G_all = linspace( 255, 425, 101 )*1e3; Ng = length( G_all ); % 初始化结果 R0_all = zeros( 1
2020-06-15 23:21:24
691
1
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.14 和 图3.15 ）
说明： （1）本程序使用的是频域匹配滤波，书本上的可能是时域匹配滤波，所以结果有一些出入；但是趋势是一样的，能够用来理解书本内容 （2）程序添加了迭代画图的功能，能够以动画的形式得到图3.15的结果 （3）程序可以设置基带和非基带信号，改tc就行； 后面贴图给的是非基带信号的结果 matlab程序代码如下： function figure3_14_15() %% initial matlab workspace % figure 3.14 and figure 3.15 ; page 64 and pag
2020-06-14 23:04:45
1556
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.13 ）
说明： （1）本程序信号开始时间为-T/2，得到的结果与设置的时延t0是一致的，方式2和方式3都能够定位到同一个时刻 （2）与书本的结果不同；所以书本应该添加或者减少一个线性相位；不过这些不要紧，只要理解时域匹配滤波器和频域匹配滤波器就行 matlab程序代码如下： % initial matlab workspace % figure 3.13; page 63 clc clear close all %% 参数设置 % LFM脉冲宽度T 和 调频斜率K T = 7.2e-6; K = 4.1e12;
2020-06-13 23:19:53
834
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.10 ）
说明： （1）图3.10中的频域和时域不是同一组数据画出来的图；本程序改为同一组程序画的了，能说明问题就行 matlab程序代码如下： % initial matlab workspace % figure 3.10 : page 59 clc clear close all %% 参数设置 % 设置 脉冲宽度T 和 抽样点数N T = 7.2e-6; N = 128; % 设置过采样因子 alpha = 1.15; % 根据设置的参数 求其他参数 dt = T/N; %
2020-06-13 21:24:36
513
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.9 和 图3.12 ）
说明： （1）如果是非基带信号，频域匹配滤波的关键在于找到DFT一个宽度内的实际频率值 matlab代码如下： % initial matlab workspace % figure 3.9; page 55 and page 56 clc clear close all %% 参数设置 % 设置时间宽度T 和 TBP T = 7.2e-6; TBP = 42; % 设置tc 和 时延t0; 总时间为t0+T,时间范围是-T/2~t0+T/2 tc = 0e-6; t0 = 1e-6; % 是否
2020-06-13 18:09:15
802
1
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.6、图3.7 、图3.8 和 图3.11 ）
matlab代码如下： % initial matlab workspace % figure 3.6; page 55 clc clear close all % 设置时间宽度T 和 TBP T = 7.2e-6; TBP = 42; % 计算B和K B = TBP/T; K = B/T; % 设置时间序列t fs = 1.2*B; dt = 1/fs; Nt = ceil( T/dt ); Nt = 2^( ceil( log2(Nt) ) ); dt = T/Nt; fs = 1/dt; t
2020-06-10 22:43:40
563
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3A.1 和 图3.5）
matlab程序代码如下： % initial matlab workspace % figure 3A.1; page72 clc clear close all % 设置参数TBP和T TBP = 100; T = 10e-6; % 求解B和K B = TBP/T; K = B/T; % 生成时间序列t t = linspace( -T, T, 2001 ); % 计算第一部分结果p1 第二部分结果p2 完整结果a % 第一部分对T归一化了 p1 = ( T - abs(t) ) ./ T;
2020-06-07 22:57:40
476
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.4）
matlab程序代码如下： % initial matlab workspace % figure3.4 page51 clc clear close all % 图片设置 h = figure; set( h, 'position', [100,10,600,700]); sub_row = 4; sub_col = 2; sub_count = 0; % 设置时间宽度T 和 抽样频率fs T = 8e-6; fs = 17e6; % 计算抽样间隔 dt = 1/fs;
2020-06-06 21:05:05
357
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.3）
说明： （1）本程序没有相位解绕，提取的是频谱的相位 matlab程序代码如下： % initial matlab workspace % figure3.3 page50 clc clear close all % 设置调频斜率K 和 基本时间长度T0 K = 25; T0 = 1; % 图片设置 h = figure; set( h, 'position', [100,10,600,700]); sub_row = 5; sub_col = 2; sub_count = 0
2020-06-06 20:34:14
539
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.a2）
说明： （1）做图3.2的收获，以高斯信号为例,来说明如何进行正确的频域插值-频谱的幅度和相位都与理论值一致 matlab程序代码如下： % initial matlab workspace % figure3.a2 ; 附加图2; 以高斯信号为例,来说明如何进行正确的频域插值 clc clear close all % 时间序列 T = 150; Nt = 150; dt = T/Nt; t = -T/2:dt:T/2-dt; % 高斯信号 sigma = 2; ft = exp( - t.
2020-06-06 16:19:57
365
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.a1）
说明： （1）做图3.2的收获，以高斯信号为例,来验证得到频谱幅相的两种求解过程 matlab程序代码如下： % initial matlab workspace % figure3.a1 ; 附加图1; 以高斯信号为例,来验证得到频谱幅相的两种求解过程 clc clear close all % 时间序列 T = 150; Nt = 150; dt = T/Nt; t = -T/2:dt:T/2-dt; % 高斯信号 sigma = 2; ft = exp( - t.^2./2./sigma
2020-06-06 16:06:14
281
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.2）
说明： （1）以书本中的数据，使用FFT计算，得到的结果如下面第一幅图所示，除了频谱幅度与书本图片一致之外，其他均不相同 （2）将TBP改成1e5，运行结果如第二幅图所示，很明显可以发现驻定相位法得到的结果为DFT得到结果的包络，注意看红色的那根线；并且包络与书本上图3.2一致（注意相位那一张图不是连续相位，而是-pi至pi之间的值，angle值） （3）推测：书本图3.2，(a)(b)(d)三幅子图是通过驻定相位法公式得到的，©子图是通过DFT得到的； （4）如果叙述有误，请留言联系；我也希望能够找到DF
2020-06-06 09:56:34
474
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.1）
matlab程序如下： % initial matlab workspace % figure3.1 page46 clc clear close all % 设置时间宽度T 和 斜率K T = 8e-6; K = 0.66e12; % 时间序列t 实部rs 虚部is 复数信号s t = linspace( -T/2, T/2, 1001 ); rs = cos( pi*K*t.^2 ); is = sin( pi*K*t.^2 ); s = rs + 1i*is; % 频率 fre = K *
2020-06-04 22:43:17
471
2
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 - 图2.16）
matlab代码如下： % initial matlab workspace % figure2.16 page37 clc clear close all % figure设置 h = figure; set( h, 'position', [100,100,900,600]); % N为插值核的点数，flag表示是否加窗 Nall = [ 4, 8, 16, 8 ]; flag = [ 0, 0, 0, 1 ]; % 设置抽样点数sn 和 fft点数sf(sf>sn是为了频域插值 使得
2020-06-04 22:05:27
416
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 - 表2.1）
说明： （1）书本上的表格行序号好像给反了 （2）本程序也给了插值核的测试语段 matlab程序如下： % initial matlab workspace % table2.1 N核M点插值表; page37 clc clear close all % 定义N 和 M N = 2^3; M = 16; % 计算fs Ts = 1; fs = 1/Ts; % 计算delta 记为d d = 1/M; % kaiser窗函数 beta = 3; T = N; w = @(t) besse
2020-06-04 21:39:59
437
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.15）
matlab程序如下： % initial matlab workspace % figure 2.15 ; page 36 % 为了方便做子程序 把主程序封装为子程序 clc clear close all % 时间序列生成 t = -8:0.2:8; % 生成未经加权插值核 kernel1 = sinc( t ); % 经过加权插值核 kernel2 = sinc( t ) .* kaiser( length(t), 2.5 )'; % 作图 figure plot( t, kernel1, '
2020-06-03 23:35:42
508
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.11和图2.12）
说明： (1) kaiser窗的傅里叶变换为sinc型，推导需要用到盖根鲍尔多项式，等把这本书看完就专攻特殊函数; (2) 3dB展宽比为相对于矩形窗的展宽比(百分比) (3) PSLR为把频谱归一化之后 第一旁瓣电平 matlab程序如下： （为了方便贴程序，把主程序也封装为子程序了； 把以下程序段全部放在一个程序内即可） function figure2_12() % initial matlab workspace % figure 2.11 ; page 33 % 为了方便做子程序 把主程序封装为子
2020-06-03 23:14:03
914
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.9）
matlab代码如下： % initial matlab workspace % figure 2.9 ; page 30 clc clear close all % LMF信号 实信号 s(t) = cos( 2*pi*(f0.*t +0.5*k*t^2) ) % LMF信号 复信号 s(t) = exp( 1i*2*pi*(f0.*t +0.5*k*t^2) ) h = figure; set( h, 'position', [100,100,900,600] ); sub_row = 5; fo
2020-05-15 22:51:20
371
1
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.8）
说明： （1）理解这张图没什么难度，但是总想写一个程序。 matlab程序如下： % initial matlab workspace % figure 2.8 ; page 29 clc clear close all % LMF信号 实信号 s(t) = cos( 2*pi*(f0.*t +0.5*k*t^2) ) % LMF信号 复信号 s(t) = exp( 1i*2*pi*(f0.*t +0.5*k*t^2) ) h = figure; set( h, 'position', [100,100
2020-05-15 22:47:35
587
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.7）
说明： （1）这个跟书本上的图2.7有点区别，加了一些内容；主要是为了验证自己的一些猜测； matlab代码如下： % initial matlab workspace % figure2.7 page28 clc clear close all %% 验证300Hz的正弦波 以800抽样 抽样点完全符合300+k*800, k=0,正负1,正负2,正负3,正负4,... % 参数设置 f0 = 300; % 连续正弦波 频率 fs = 800; % 抽样频率 k = -2; % 计算抽样频率 T
2020-05-15 22:44:17
936
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.6）
matlab代码如下： % initial matlab workspace % figure 2.6 ; page 27 clc clear close all % 图片初始 h = figure(); set( h, 'position', [100,50,900,600] ); sub_row = 2; sub_col = 2; sub_count = 0; % LMF信号 实信号 s(t) = cos( 2*pi*(f0.*t +0.5*k*t^2) ) % LMF信号 复信号 s(t) =
2020-05-15 22:41:13
518
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 - DFT循环卷积特性 page24 ）
收获： （1）DFT的循环卷积特性（已数学证明，需要的小伙伴可以联系我，手写版就不发了） matlab程序如下： % initial matlab workspace % 验证DFT做卷积和直接卷积的异同 page25 clc clear close all % 数据设置 x = [ 1 3 -1 5 2 6 4 -2]; y = [ 1 2 3 ]; % 提取数据长度 N1 = length(x); N2 = length(y); % 直接做卷积 result1 = conv(
2020-05-15 22:37:14
421
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.2）
欢迎使用Markdown编辑器 你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。 新的改变 我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客： 全新的界面设计 ，将会带来全新的写作体验； 在创作中心设置你喜爱的代码高亮样式，Markdown 将代码片显示选择的高亮样式 进行展示； 增加了 图片
2020-05-15 22:28:06
722
2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章补充 如何做正确的fft）
2020-10-22
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人