雷达SAR成像仿真的应用（Matlab代码实现）

news2024/11/18 19:57:05

🍒🍒🍒欢迎关注🌈🌈🌈

📝个人主页：我爱Matlab

👍点赞➕评论➕收藏 == 养成习惯（一键三连）🌻🌻🌻

🍌希望大家多多支持🍓~一起加油 🤗

💬语录：将来的我一定会感谢现在奋斗的自己！

🍁🥬🕒摘要🕒🥬🍁

该存储库包含雷达SAR成像的软件仿真，包括线性频率调制（LFM，SAR成像中采用的基本信号），带有点目标的简单SAR成像（您可以了解SAR的工作原理），最后是真实的数据处理（在处理原始数据后获得图像。

✨🔎⚡运行结果⚡🔎✨

💂♨️👨‍🎓Matlab代码👨‍🎓♨️💂

%% 
%    线性调频信号的脉冲压缩
% 
%    介绍：对线性调频信号进行仿真，输出其时频域的相关信息，并模拟回波信号，
%    对其进行脉冲压缩和加窗处理。
% 
%    实验记录：
%     1.线性调频信号时域包络、相位；实部、虚部
%     2.线性调频信号频谱幅频、相频特性；实部、虚部
%     3.两个目标回波的时域和频域波形
%     4.信号通过匹配滤波器的输出结果（脉冲压缩）。
%     5.用Hamming窗抑制脉冲压缩结果副瓣
%%  基本参数 
    clc;clear all;close all;

T = 10e-6; % LFM周期/脉宽 10us
B = 60e6; % LFM带宽 60Mhz
fs = 100e6; % 采样率 100MHz
K = B/T;
%% 模拟发射信号
n = round(15*T*fs);
t = linspace(-10*T, 10*T,n);

lfmT = rectpuls(t,T).*exp(1j*pi*K*t.^2);
lfmF = fftshift(fft(fftshift(lfmT)));
f = linspace(-fs,fs,n);

%% 时域绘图
figure();
plot(diff(phase(lfmT)));
title('LFM信号的时间-频率变化趋势图');
xlabel('时间');
ylabel('频率');
xlim([7200,7800])
% 包络
figure();
subplot(2,2,1);
plot(t,abs(lfmT));
title('LFM信号时域包络');
xlabel('t/s');
ylabel('幅度');
xlim([-1e-5,1e-5])
ylim([-0.5,1.5])
% 相位
subplot(2,2,2);
plot(t,phase(lfmT));
title('LFM信号时域相位');
xlabel('t/s');
ylabel('相位');
xlim([-5e-6,5e-6])
% 实部
subplot(2,2,3);
plot(t,real(lfmT));
title('LFM信号时域实部');
xlabel('t/s');
ylabel('幅度');
xlim([-1.5e-6,1.5e-6]);
ylim([-1,1]);
% 虚部
subplot(2,2,4);
plot(t,imag(lfmT));
title('LFM信号时域虚部');
xlabel('t/s');
ylabel('幅度');
xlim([-1.5e-6,1.5e-6]);
ylim([-1,1]);
%% 频域绘图
figure();
subplot(2,2,1);
plot(f,abs(lfmF));
title('LFM信号幅频特性');
xlabel('Hz');
ylabel('幅度');

subplot(2,2,2);
plot(unwrap(angle(lfmF)));
title('LFM信号相频特性');
xlabel('Hz');
ylabel('相位');

subplot(2,2,3);
plot(f,real(lfmF));
title('LFM信号频谱实部');
xlabel('Hz');
ylabel('幅度');
xlim([-3e7,3e7]);

subplot(2,2,4);
plot(f,imag(lfmF));
title('LFM信号频谱虚部');
xlabel('Hz');
ylabel('幅度');
xlim([-3e7,3e7]);