写一个雷达信号分选的MATLAB代码，从PDWdata.txt当中提取数据，PDWdata.txt文件当中的数据格式是每一行为一个数据，每一行中的顺序为序号、频率、脉冲宽度、TOA、PA，以空格隔开。要求编写基于TOA的MATLAB程序，采用SDIF算法，分选出这些信号是由几部雷达发出的，计算出每部雷达的PRI，绘制出每部雷达的频率图像

以下是一个基于TOA的MATLAB程序，用于实现SDIF算法的雷达信号分选和绘制频率图像：

% 雷达信号分选程序
clear all; clc;

% 加载数据
filename = 'PDWdata.txt';
data = load(filename);

% 设置相关参数
TOA_threshold = 5; % TOA的阈值，用于判断两个脉冲是否来自同一个雷达
PRI_threshold = 50; % PRI的阈值，用于分选不同雷达的信号
N = size(data,1); % 数据长度
num_radar = 0; % 雷达数量
selected_data = cell(1,1); % 用于存储分选后的数据
radar_PRI = zeros(1,1); % 用于存储每个雷达的PRI

% 遍历数据，进行分选
for i = 2:N
    TOA = data(i,4) - data(i-1,4); % 计算TOA差值
    if TOA > TOA_threshold % 如果TOA超过阈值，则认为两个脉冲来自不同雷达
        num_radar = num_radar + 1; % 雷达数量加1
        selected_data{num_radar} = data(i-1,:); % 将前一个脉冲存储为新雷达的第一个脉冲
    end
    selected_data{num_radar} = [selected_data{num_radar}; data(i,:)]; % 将当前脉冲存储到当前雷达的数据中
end

% 计算每个雷达的PRI
for i = 1:num_radar
    PRI = selected_data{i}(2:end,4) - selected_data{i}(1:end-1,4); % 计算每个脉冲的PRI
    mean_PRI = mean(PRI(PRI > PRI_threshold)); % 取PRI大于阈值的脉冲的平均值作为该雷达的PRI
    radar_PRI(i) = mean_PRI;
end

% 绘制每个雷达的频率图像
figure; hold on;
for i = 1:num_radar
    f = selected_data{i}(:,2); % 频率
    pw = selected_data{i}(:,3); % 脉冲宽度
    plot(f, pw./radar_PRI(i), 'o'); % 绘制频率图像
end
xlabel('频率'); ylabel('PRI/脉冲宽度');

程序首先加载PDWdata.txt文件中的数据，然后根据TOA的差值，将信号分选为不同的雷达数据，存储在selected_data数组中。接着，程序计算每个雷达的PRI，将结果存储在radar_PRI数组中。最后，程序遍历每个雷达的数据，绘制出每个雷达的频率图像，并以PRI/脉冲宽度为纵坐标，以频率为横坐标标示出来。