main.m

%% main.m
% UAV-Ship信道模型主脚本

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]; % UAV初始位置
params.Rx_pos = [600 0 1.5]; % Ship初始位置
params.phiA_SB = 0; % SB簇运动方位角
params.phiA_CA_MB = 0; % MB首簇运动方位角
params.phiA_CZ_MB = 0; % MB尾簇运动方位角
% 天线阵列与极化
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;

%% 生成信道分量（MIMO+极化+空间点积+Doppler积分）
h_los = los_channel(params, t_vec);
h_sb = sb_channel(params, t_vec);
h_mb = mb_channel(params, t_vec);

% 总信道（逐极化分量相加）
h_tot = eta1*h_los + eta2*h_sb + eta3*h_mb;

%% 可视化（默认绘制第一个Tx-Rx天线对的H-H极化分量）
figure;
plot(t_vec, abs(squeeze(h_tot(1,1,1,:))));
hold on;
plot(t_vec, abs(squeeze(h_los(1,1,1,:))), '--');
plot(t_vec, abs(squeeze(h_sb(1,1,1,:))), ':');
plot(t_vec, abs(squeeze(h_mb(1,1,1,:))), '-.');
legend('总信道','LOS','SB','MB');
xlabel('时间 (s)'); ylabel('信道幅度');
title('UAV-Ship信道模型（第1对天线H-H极化分量）');
grid on;