对线性调频信号LFM的射频噪声干扰

clear ; % 清除变量 close all; % 关闭窗口 clc; % 清屏 %% 雷达基本参数 Radar.Tp = 10e-6; % 脉冲宽度 Radar.BW = 10e6; % 带宽 Radar.Kr = Radar.BW/Radar.Tp; % 调频率 Radar.Fc = 100e6; % 载频 Radar.Fs = 230e6; % 采样频率 Radar.sample = Radar.Fs*Radar.Tp; % 采样点数 t = -Radar.Tp/2:Radar.Tp/Radar.sample:Radar.Tp/2-Radar.Tp/Radar.sample; % 时间范围 f = -Radar.Fs/2:Radar.Fs/Radar.sample:Radar.Fs/2-Radar.Fs/Radar.sample; % 频率范围 %% LFM信号 Radar.return = exp(1i*(2*pi*Radar.Fc*t+pi*Radar.Kr*t.^2)); % LFM信号时域 Radar.return_f = fftshift(fft(Radar.return)); % LFM信号频谱 figure(1); % 绘制第1幅图 subplot(211); % 将窗口分割成2*1的，当前是第1个子图 plot(t*1e6,real(Radar.return)); % 绘制复信号时域波形 title('复信号时域波形'); % 标题 xlabel('时间/us'); % 横坐标 ylabel('幅度'); % 纵坐标 grid on; % 打开网格线 subplot(212); % 将窗口分割成2*1的，当前是第2个子图 plot(f/1e6,abs(Radar.return_f )); % 绘制复信号频域波形 title('复信号频域波形'); % 标题 xlabel('频率/MHz'); % 横坐标 ylabel('幅度'); % 纵坐标 grid on; % 打开网格线 %% 射频噪声干扰信号 RFNJ.wave = RFNJ2(Radar.sample,Radar.Fs,Radar.Fc,Radar.BW,t); % 射频噪声 % RFNJ.wave = RFNJ(Radar.sample,Radar.Fc,t); % 射频噪声时域 RFNJ.wave_f = fftshift(fft(RFNJ.wave)); % 射频噪声频域 figure(2); % 绘制第2幅图 subplot(211); % 将窗口分割成2*1的，当前是第1个子图 plot(t*1e6,real(RFNJ.wave)); % 绘制射频噪声时域波形 title('射频噪声时域波形'); % 标题 xlabel('时间/us'); % 横坐标 ylabel('幅度'); % 纵坐标 grid on; % 打开网格线 subplot(212); % 将窗口分割成2*1的，当前是第2个子图 plot(f/1e6,abs(RFNJ.wave_f)); % 绘制射频噪声频域波形 title('射频噪声频域波形'); % 标题 xlabel('频率/MHz'); % 横坐标 ylabel('幅度'); % 纵坐标 grid on; % 打开网格线 %% 混合信号 Radar_RFNJ.return = Radar.return + RFNJ.wave; % 混合信号时域 Radar_RFNJ.return_f = fftshift(fft(Radar_RFNJ.return)); % 混合信号频域 figure(3); % 绘制第3幅图 subplot(211); % 将窗口分割成2*1的，当前是第1个子图 plot(t*1e6,real(Radar_RFNJ.return)); % 绘制混合信号时域波形 title('混合信号时域波形'); % 标题 xlabel('时间/us'); % 横坐标 ylabel('幅度'); % 纵坐标 grid on; % 打开网格线 subplot(212); % 将窗口分割成2*1的，当前是第2个子图 plot(f/1e6,abs(Radar_RFNJ.return_f)); % 绘制混合信号频域波形 title('混合信号频域波形'); % 标题 xlabel('频率/MHz'); % 横坐标 ylabel('幅度'); % 纵坐标 grid on; % 打开网格线 %% 解调 Radar_RFNJ.return_2 = Radar_RFNJ.return.*exp(-1i*(2*pi*Radar.Fc*t)); % 解调后信号�