雷达信号的时域分析原始MATLAB代码，自己构造信号输入，分析输出结果包括PA 、TOA、频率 脉冲宽度 脉内调制类型和参数（脉内频率捷变、频率编码、线性调频、多项编码、非线性调频、常规脉冲、联合调制），代码尽量详细

以下是一个简单的基于 MATLAB 的雷达信号时域分析代码，其中包括了常见的参数计算和脉内调制类型和参数分析。

% 生成示例信号
fs = 100e6; % 采样率
t = 0:1/fs:10e-6; % 时域范围
f0 = 1e9; % 起始频率
f1 = 5e9; % 终止频率
B = f1 - f0; % 带宽
s = chirp(t, f0, t(end), f1); % 线性调频信号

% 分析信号
N = length(s); % 信号长度
S = fft(s); % 傅里叶变换
P = abs(S).^2/N; % 功率谱密度
f = linspace(-fs/2, fs/2, N); % 频率范围
PdB = 10*log10(P); % 转换为 dB

% 计算 TOA
[~, idx] = max(abs(s)); % 找到最大值的位置
toa = t(idx); % TOA

% 计算 PA
pa = max(abs(s)); % PA

% 计算脉冲宽度
pw = (t(find(abs(s) > 0.5*pa, 1, 'last')) - t(find(abs(s) > 0.5*pa, 1))) * 1e6; % 脉冲宽度，单位微秒

% 分析脉内调制类型和参数
fif = 10e6; % 中频
sif = s .* exp(-1i*2*pi*fif*t); % 下变频
SIF = fft(sif); % 下变频后的傅里叶变换
Pif = abs(SIF).^2/N; % 下变频后的功率谱密度
PifdB = 10*log10(Pif); % 下变频后的功率谱密度转换为 dB

% 计算脉内频率捷变
sif_fm = fmmod(sif, 1e6, fs); % 脉内频率捷变调制
demod = fmdemod(sif_fm, 1e6, fs); % 解调
fm = demod / (2*pi*fif); % 脉内频率捷变，单位 Hz

% 计算频率编码
fcode = linspace(f0, f1, N); % 频率编码

% 计算线性调频
sif_lfm = chirp(t, f0-fif, t(end), f1-fif); % 线性调频信号
lfm_bw = f1 - f0; % 线性调频带宽

% 计算多项式编码
poly_order = 3; % 多项式阶数
poly_coef = [1, 2, 3, 4]; % 多项式系数
sif_poly = polyval(poly_coef, t.^poly_order); % 多项式编码

% 计算非线性调频
sif_nlfm = chirp(t, f0-fif, t(end), f1-fif, 'quadratic'); % 非线性调频信号

% 计算常规脉冲
sif_rect = rectpuls(t, pw); % 矩形脉冲

% 计算联合调制
sif_joint = sif_lfm .* sif_poly; % 联合调制信号

以上代码生成了一个线性调频信号并计算了一些常见的参数，如 TOA、PA、脉冲宽度等。此外，还分析了脉内调制类型和参数，包括脉内频率捷变、频率编码、线性调频、多项式编码、非线性调频、常规脉冲和联合调制。