function [LFMPulse,targetEchoPRT,matchedFilterCoeff,pulseNumber,PRTNumber] = GenerateLFMSignal(bandWidth,pulseDuration,PRTDuration,samplingFrequency,signalPower,targetDistece,plotEnableHigh)
% 该函数用于产生线性调频信号，以及雷达的目标反射回波，仅产生单个回波
% Author:huasir 2023.9.21 @Beijing
% Input :
% * bandWidth: 信号带宽 ，参考值：2.0e6 表示2MHz
% * pulseDuration：脉冲持续时间，参考值：40.0e-6 表示40ms
% * PRTDuration：脉冲重复周期，参考值：240ms
% * samplingFrequency：采样频率，参考值：2倍的信号带宽
% * signalPower：信号能量，参考值：1
% * targetDistece：目标距离，最大无模糊距离由脉冲重复周期决定。计算公式：1/2*PRTDuration*光速
% * plotEnableHigh： 绘图控制符，1：打开绘图，0：关闭绘图
% Output :
% * LFMPulse：线性调频信号
% * targetEchoPRT： 目标反射回波
% * matchedFilterCoeff： 匹配滤波器系数
% * pulseNumber：当前采样率下线性调频信号的采样点数
% * PRTNumber：1个PRT对应的采样点数
C = 3.0e8; %光速（m/s）
BandWidth = bandWidth; %雷达发射信号带宽，带宽=B=1/tau，tau是脉冲宽度
TimeWidth = pulseDuration; %雷达发射信号的脉冲时宽
PRT = PRTDuration; %雷达发射脉冲重复周期（s），240us对应1/2*240*300=360000米最大无模糊距离
Fs = samplingFrequency; %采样频率
SampleNumber = fix(Fs*PRT);
%=========================================================================%
% 目标参数设置 %
%=========================================================================%
SigPower = signalPower; %目标功率，无量纲
TargetDistance = targetDistece; %目标距离，单位：m
DelayNumber = fix(Fs*2*TargetDistance/C); %把目标距离换算成采样点（距离门）
%=========================================================================%
% 产生线性调频信号、匹配滤波器 %
%=========================================================================%
number = fix(Fs*TimeWidth); %回波采样点数=脉压系数长度=暂态点数目+1
if rem(number,2)~=0
nember = nember + 1;
end
Chirp = zeros(1,number);
for i = -fix(number/2):fix(number/2)-1
Chirp(i+fix(number/2)+1)=exp(1j*(pi*(BandWidth/TimeWidth)*(i/Fs)^2));%产生复ChIrp信号
end
coeff = conj(fliplr(Chirp)); %把Chirp信号翻转并把复数共轭，产生脉压系数
%=========================================================================%
% 绘制线性调频信号 %
%=========================================================================%
if plotEnableHigh == 1
figure;
plot(real(Chirp)); %绘制线性调频信号
xlabel('Sampling points'); ylabel('Amplitude');title('线性调频信号实部');
end
SignalTemp = zeros(1,SampleNumber); %1个PRT
SignalTemp(DelayNumber+1:DelayNumber+number) = sqrt(SigPower)*Chirp;%将线性调频信号按照距离进行延时
if plotEnableHigh == 1
figure;
plot(real(SignalTemp)); %绘制1个完整的PRT的雷达回波信号
xlabel('Range bin'); ylabel('Amplitude');title('雷达回波的实部');
end
%=========================================================================%
% 进行脉冲压缩 %
%=========================================================================%
Echo = SignalTemp; % 目标回波
pc_time0 = conv(Echo,coeff); % 回波和滤波器卷积的结果
pc_time1 = pc_time0(number:number+SampleNumber-1); %去掉暂态点
realTargetRange = find(abs(pc_time1)==max(abs(pc_time1)))-1; %由脉压结果目标距离
fprintf('The target range bin is %d',realTargetRange);
if plotEnableHigh == 1
figure; %时域脉压结果
subplot(2,1,1);plot(abs(pc_time0),'r-');
xlabel('Range bin'); ylabel('Amplitude');title('时域脉压结果');
subplot(2,1,2);plot(abs(pc_time1),'r-');
xlabel('Range bin'); ylabel('Amplitude');title('去掉暂态点的时域脉压结果');
end
%=========================================================================%
% 返回参数 %
%=========================================================================%
LFMPulse = Chirp; %线性调频信号
targetEchoPRT = SignalTemp; %目标反射回波
matchedFilterCoeff = coeff; %匹配滤波器系数
pulseNumber = number; %线性调频信号的采样点数
PRTNumber = SampleNumber; %目标反射回波的采样点数
end
