test_functions.m

%% test_functions.m
% 各主要子函数单元测试脚本，适配空间点积+Doppler积分

clear; clc;
addpath('./');

%% 加载参数
config;

params = struct;
params.c = c;
params.lambda = lambda;
params.K_factor = K_factor;
params.K_SB = K_SB;
params.Lk = Lk;
params.N_MB = N_MB;
params.Mn = Mn;
params.v_UAV = [v_UAV 0 0];
params.v_Ship = [v_Ship 0 0];
params.v_SB = v_SB;
params.v_CA_MB = v_CA_MB;
params.v_CZ_MB = v_CZ_MB;
params.a_waveguide = a_waveguide;
params.b_waveguide = b_waveguide;
params.Tx_pos = [0 0 100];
params.Rx_pos = [600 0 1.5];
params.phiA_SB = 0;
params.phiA_CA_MB = 0;
params.phiA_CZ_MB = 0;
% 天线阵列与极化
params.P = P; params.Q = Q;
params.array_type_Tx = array_type_Tx;
params.array_type_Rx = array_type_Rx;
params.d_Tx = d_Tx; params.d_Rx = d_Rx;
params.UPA_Tx_row = UPA_Tx_row; params.UPA_Tx_col = UPA_Tx_col;
params.UPA_Rx_row = UPA_Rx_row; params.UPA_Rx_col = UPA_Rx_col;
params.polar_type = polar_type;
params.XPR = XPR;
% 空间点积与Doppler积分相关参数
params.phiA_Tx = phiA_Tx; params.phiE_Tx = phiE_Tx;
params.phiA_Rx = phiA_Rx; params.phiE_Rx = phiE_Rx;
params.phi_vec_Tx = phi_vec_Tx;
params.phi_vec_Rx = phi_vec_Rx;

%% 测试LOS信道
h_los = los_channel(params, t_vec);
h1 = squeeze(h_los(1,1,1,:));
figure; plot(t_vec, abs(h1)); title('LOS信道幅度（第1对天线H-H极化）'); xlabel('时间(s)'); ylabel('|h_{LOS}|'); grid on;

%% 测试SB信道
h_sb = sb_channel(params, t_vec);
h2 = squeeze(h_sb(1,1,1,:));
figure; plot(t_vec, abs(h2)); title('SB信道幅度（第1对天线H-H极化）'); xlabel('时间(s)'); ylabel('|h_{SB}|'); grid on;

%% 测试MB信道
h_mb = mb_channel(params, t_vec);
h3 = squeeze(h_mb(1,1,1,:));
figure; plot(t_vec, abs(h3)); title('MB信道幅度（第1对天线H-H极化）'); xlabel('时间(s)'); ylabel('|h_{MB}|'); grid on;

%% 工具函数测试（如有）
% 例如：角度向量生成、距离计算等
% 可根据后续扩展添加

disp('各主要子函数测试完毕。');