- 博客(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
403
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第四章 - 图4.4 ）
matlab代码如下：%% initial matlab workspace% 图4.4 : page 81clcclearclose 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
579
原创 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
1425
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.13 ）
说明：（1）本程序信号开始时间为-T/2，得到的结果与设置的时延t0是一致的，方式2和方式3都能够定位到同一个时刻（2）与书本的结果不同；所以书本应该添加或者减少一个线性相位；不过这些不要紧，只要理解时域匹配滤波器和频域匹配滤波器就行matlab程序代码如下：% initial matlab workspace% figure 3.13; page 63clcclearclose all%% 参数设置% LFM脉冲宽度T 和 调频斜率KT = 7.2e-6;K = 4.1e12;
2020-06-13 23:19:53
741
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.10 ）
说明：（1）图3.10中的频域和时域不是同一组数据画出来的图；本程序改为同一组程序画的了，能说明问题就行matlab程序代码如下：% initial matlab workspace % figure 3.10 : page 59clcclearclose all%% 参数设置% 设置 脉冲宽度T 和 抽样点数NT = 7.2e-6; N = 128; % 设置过采样因子alpha = 1.15;% 根据设置的参数 求其他参数dt = T/N; %
2020-06-13 21:24:36
427
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.9 和 图3.12 ）
说明：（1）如果是非基带信号，频域匹配滤波的关键在于找到DFT一个宽度内的实际频率值matlab代码如下：% initial matlab workspace% figure 3.9; page 55 and page 56clcclearclose all%% 参数设置% 设置时间宽度T 和 TBPT = 7.2e-6;TBP = 42;% 设置tc 和 时延t0; 总时间为t0+T,时间范围是-T/2~t0+T/2tc = 0e-6; t0 = 1e-6;% 是否
2020-06-13 18:09:15
680
1
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.6、图3.7 、图3.8 和 图3.11 ）
matlab代码如下：% initial matlab workspace% figure 3.6; page 55clcclearclose all% 设置时间宽度T 和 TBPT = 7.2e-6;TBP = 42;% 计算B和KB = TBP/T;K = B/T;% 设置时间序列tfs = 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
493
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3A.1 和 图3.5）
matlab程序代码如下：% initial matlab workspace% figure 3A.1; page72clcclearclose all% 设置参数TBP和TTBP = 100;T = 10e-6;% 求解B和KB = TBP/T;K = B/T;% 生成时间序列tt = linspace( -T, T, 2001 );% 计算第一部分结果p1 第二部分结果p2 完整结果a% 第一部分对T归一化了p1 = ( T - abs(t) ) ./ T;
2020-06-07 22:57:40
398
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.4）
matlab程序代码如下：% initial matlab workspace% figure3.4 page51clcclearclose all% 图片设置h = figure; set( h, 'position', [100,10,600,700]);sub_row = 4;sub_col = 2;sub_count = 0;% 设置时间宽度T 和 抽样频率fsT = 8e-6;fs = 17e6;% 计算抽样间隔dt = 1/fs;
2020-06-06 21:05:05
288
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.3）
说明：（1）本程序没有相位解绕，提取的是频谱的相位matlab程序代码如下：% initial matlab workspace% figure3.3 page50clcclearclose all% 设置调频斜率K 和 基本时间长度T0K = 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
466
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.a2）
说明：（1）做图3.2的收获，以高斯信号为例,来说明如何进行正确的频域插值-频谱的幅度和相位都与理论值一致matlab程序代码如下：% initial matlab workspace% figure3.a2 ; 附加图2; 以高斯信号为例,来说明如何进行正确的频域插值clcclearclose 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
297
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.a1）
说明：（1）做图3.2的收获，以高斯信号为例,来验证得到频谱幅相的两种求解过程matlab程序代码如下：% initial matlab workspace% figure3.a1 ; 附加图1; 以高斯信号为例,来验证得到频谱幅相的两种求解过程clcclearclose 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
224
原创 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
400
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第三章 - 图3.1）
matlab程序如下：% initial matlab workspace% figure3.1 page46clcclearclose all% 设置时间宽度T 和 斜率KT = 8e-6;K = 0.66e12;% 时间序列t 实部rs 虚部is 复数信号st = 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
396
2
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 - 图2.16）
matlab代码如下：% initial matlab workspace% figure2.16 page37clcclearclose 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
355
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 - 表2.1）
说明：（1）书本上的表格行序号好像给反了（2）本程序也给了插值核的测试语段matlab程序如下：% initial matlab workspace% table2.1 N核M点插值表; page37clcclearclose all% 定义N 和 MN = 2^3;M = 16;% 计算fsTs = 1;fs = 1/Ts;% 计算delta 记为dd = 1/M;% kaiser窗函数beta = 3;T = N;w = @(t) besse
2020-06-04 21:39:59
375
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.15）
matlab程序如下：% initial matlab workspace% figure 2.15 ; page 36% 为了方便做子程序 把主程序封装为子程序clcclearclose all% 时间序列生成t = -8:0.2:8;% 生成未经加权插值核kernel1 = sinc( t );% 经过加权插值核kernel2 = sinc( t ) .* kaiser( length(t), 2.5 )';% 作图figureplot( t, kernel1, '
2020-06-03 23:35:42
435
原创 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
814
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.9）
matlab代码如下：% initial matlab workspace% figure 2.9 ; page 30clcclearclose 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
288
1
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.8）
说明：（1）理解这张图没什么难度，但是总想写一个程序。matlab程序如下：% initial matlab workspace% figure 2.8 ; page 29clcclearclose 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
508
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.7）
说明：（1）这个跟书本上的图2.7有点区别，加了一些内容；主要是为了验证自己的一些猜测；matlab代码如下：% initial matlab workspace% figure2.7 page28clcclearclose 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
852
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.6）
matlab代码如下：% initial matlab workspace% figure 2.6 ; page 27clcclearclose 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
449
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 - DFT循环卷积特性 page24 ）
收获：（1）DFT的循环卷积特性（已数学证明，需要的小伙伴可以联系我，手写版就不发了）matlab程序如下：% initial matlab workspace% 验证DFT做卷积和直接卷积的异同 page25clcclearclose 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
337
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章 -图2.2）
欢迎使用Markdown编辑器你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。新的改变我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客：全新的界面设计 ，将会带来全新的写作体验；在创作中心设置你喜爱的代码高亮样式，Markdown 将代码片显示选择的高亮样式 进行展示；增加了 图片
2020-05-15 22:28:06
631
2007 合成孔径雷达成像算法与实现 各figure的matlab实现（第二章补充 如何做正确的fft）
2020-10-22
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人