RD-AML-Clean属性散射中心特征提取程序（可运行）

%RD-AML-Clean 属性散射中心特征提取程序
%清空环境变量
close; clc; clear All;

% Step 1: 频域数据的频率 f 和方位角 ϕ 的范围
f_min = 1e2; % 频率范围起始点，单位Hz
f_max = 10e2; % 频率范围终止点，单位Hz
num_frequencies = 150; % 频率采样点数

phi_min = 0; % 方位角范围起始点，单位度
phi_max = 1; % 方位角范围终止点，单位度
num_angles = 15; % 方位角采样点数

frequencies = linspace(f_min, f_max, num_frequencies); % 频率范围
phis = linspace(phi_min, phi_max, num_angles); % 方位角范围

%% 1.2 基于属性散射中心模型产生极坐标格式的频域复数据 D(f, ϕ)
% Step 2: 产生极坐标格式的频域复数据 D(f, ϕ)
scatter_centers = [50, 50; 65, 65; 80, 80; 95, 95; 110, 110]; % 散射中心位置 (x, y)
amplitudes = [1, 1, 1, 1, 1]; % 每个散射中心的振幅
D = zeros(num_frequencies, num_angles); % 初始化频域数据

for k = 1:length(scatter_centers)
    x0 = scatter_centers(k, 1);
    y0 = scatter_centers(k, 2);
    for m = 1:num_frequencies
        for n = 1:num_angles
            % 正确计算 r 和 phase
            r = sqrt(x0^2 + y0^2 + 1e-5); % 防止除零错误