合成孔径雷达成像——算法与实现图【2.8】仿真

T = 1;BW = 150;k = BW/T;h=figure();N_0 = length(t_0);

设置初始参数，带宽在书中给出，为150Hz，h=figure()是为了后面将所有的图在一张图中显示。

p0_0 = t_0+2;p0_max_0 = max(max(p0_0));p_0 = p0_0./p0_max_0;

信号幅度修正，这里加的2这里的数越大，产生的图像斜率越大。

dt = 1/400;t = 0:dt:T;fs = 1/dt;% 信号幅度修正p0 = t+2;p0_max = max(max(p0));p = p0./p0_max;N = length(t);% 频率f = (0:N-1)*fs/N - fs/2;

同样的，构造采样后的参数。

f0 = 25+(i-1)*50;%初始频率f0%生成信号s_0 = cos(2*pi*(f0.*t_0+0.5*k*t_0.^2)).*p_0;% 实数基带% 进行fftF_0 = fftshift(fft(s_0,N_0));F_0 = F_0./max(abs(F_0));

生成初始信号，i从1到9，与书中移动九次一致，先生成初始的信号，然后进行傅里叶变换，归一化后即可直接绘图。

%采样信号s = cos(2*pi*(f0.*t+0.5*k*t.^2)).*p;% 进行fftF = fftshift(fft(s)./ N ); F = F./max(abs(F));

生成采样信号，然后对其进行相同处理。

figure(h); hold onsubplot(9,2,2*(i-1)+1 );plot( f_0, abs(F_0) ); ylabel('幅度');grid on;title( ['FC=',num2str(f0+BW/2)])subplot(9,3,3*(i-1)+3 );plot( f, abs(F) ); ylabel('幅度');grid on;if i == 1title('采样后的频率');endend

最后绘制图形。这里可看出，如果在最开头没有加入h=figure，则会产生九张图。

附上完整代码：

%% 原始信号T = 1;%设置时间长度TBW = 150;% 带宽BWk = BW/T;%斜率h=figure();%使所有图在一张图中生成dt_0 = 1/1200; t_0 = 0:dt_0:T;% 生成原始时间fs_0 = 1/dt_0;N_0 = length(t_0);% 信号幅度修正p0_0 = t_0+2;p0_max_0 = max(max(p0_0));p_0 = p0_0./p0_max_0; %% 抽样dt = 1/400;t = 0:dt:T;fs = 1/dt;% 信号幅度修正p0 = t+2;p0_max = max(max(p0));p = p0./p0_max;N = length(t);% 频率f = (0:N-1)*fs/N - fs/2;%% 构造信号与傅里叶变换for i = 1:9f0 = 25+(i-1)*50;%初始频率f0%初始信号s_0 = cos(2*pi*(f0.*t_0+0.5*k*t_0.^2)).*p_0;% 进行fftF_0 = fftshift(fft(s_0,N_0));F_0 = F_0./max(abs(F_0));% 频率f_0 = (0:N_0-1)*fs_0/N_0 - fs_0/2;%采样信号s = cos(2*pi*(f0.*t+0.5*k*t.^2)).*p;% 进行fftF = fftshift(fft(s)./ N ); F = F./max(abs(F));% 画图figure(h); hold onsubplot(9,2,2*(i-1)+1 );plot( f_0, abs(F_0) ); ylabel('幅度');grid on;title( ['FC=',num2str(f0+BW/2)])subplot(9,3,3*(i-1)+3 );plot( f, abs(F) ); ylabel('幅度');grid on;if i == 1title('采样后的频率');endend

下面是结果图：

可在h=figure()处加上set( h, ‘position’, [100,100,900,600] );使图像更好看。

h = figure; set( h, 'position', [100,100,900,600] );

参考代码：

/zhoubanlu4/article/details/106151183