数据链Link16波形仿真走起，来学习吧(2)发波形仿真

注明：借鉴了多个AI（豆包和DeepSeek）写的程序之后，结合本人的仿真经验多次训练和改动而成，AI写的代码的改动量超过40%。

%% Link16 STDP 完整帧生成程序

% 按照7个步骤实现：

% 1. 数据打包 (210比特 + 15比特源航迹号)

% 2. CRC编码 (237,225) + 备用位插入

% 3. 报头添加 + RS编码 (16,7)和(31,15)

% 4. 码元块交织 (109个符号)

% 5. 同步头添加 (粗同步 + 精同步)

% 6. CCSK编码 (32位扩频序列循环移位)

% 7. MSK调制

lear; clc; close all;

%% ========= 初始化参数 ========

fprintf('=== Link 16 STDP 帧生成程序 ===\n');

fprintf('按照7个步骤实现完整STDP帧\n\n');

% 系统参数

M_SYMBOL = 32;                  % 调制阶数 (5比特/符号)

BITS_PER_SYMBOL = 5;            % 每符号比特数

%% ======== 第一步：数据打包 ========

fprintf('=== 第一步：数据打包 ===\n');

% 生成测试数据 (210比特)

test_data_bits = randi([0, 1], 1, 210);

source_track_bits = randi([0, 1], 1, 15);  % 15比特 航迹号

% 组合成225比特消息数据

msg_data_bits = [test_data_bits, source_track_bits];

fprintf('原始数据: %d 比特\n', length(test_data_bits));

fprintf('源航迹号: %d 比特\n', length(source_track_bits));

fprintf('消息数据: %d 比特\n', length(msg_data_bits));

%% ============ 第二步：CRC编码 ============

fprintf('\n=== 第二步：CRC编码 (237,225) ===\n');

% CRC-12多项式: x^12 + x^11 + x^3 + x^2 + x + 1（自定义）

crc_polynomial = [1 1 0 0 0 0 0 0 1 1 1 1];  % 0x80F

% 计算CRC-12校验位