- 博客(24)
- 资源 (3)
- 收藏
- 关注
原创 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
412
原创 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
601
原创 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
1436
原创 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
749
原创 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
439
原创 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
690
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
503
原创 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
408
原创 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
297
原创 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
474
原创 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
303
原创 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
229
原创 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
407
原创 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
409
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
364
原创 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
382
原创 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
443
原创 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
827
原创 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
297
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
512
原创 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
863
原创 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
457
原创 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
343
原创 2007 合成孔径雷达成像算法与实现 各figure的matlab实现(第二章 -图2.2)
欢迎使用Markdown编辑器你好! 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章,了解一下Markdown的基本语法知识。新的改变我们对Markdown编辑器进行了一些功能拓展与语法支持,除了标准的Markdown编辑器功能,我们增加了如下几点新功能,帮助你用它写博客:全新的界面设计 ,将会带来全新的写作体验;在创作中心设置你喜爱的代码高亮样式,Markdown 将代码片显示选择的高亮样式 进行展示;增加了 图片
2020-05-15 22:28:06
641
2007 合成孔径雷达成像算法与实现 各figure的matlab实现(第二章补充 如何做正确的fft)
2020-10-22
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人