OFDM系统注水功率分配仿真

资 源 简 介

本项目旨在开发一套基于MATLAB的正交频分复用（OFDM）系统功率分配仿真平台，核心目标是实现并验证基于香农容量最大化的注水算法（Water-filling Algorithm）。项目代码包含完整的主程序及模块化的子函数设计。功能涵盖：1. 构建OFDM通信链路基础模型，包括随机比特流生成、QAM/QPSK调制、IFFT变换、循环前缀（CP）添加以及通过频率选择性瑞利衰落信道。2. 核心算法实现中，开发了基于CSI（信道状态信息）的迭代注水算法函数，该函数利用拉格朗日乘数法计算水位线，根据各个子载波的信噪比增益动态分配发射功率，遵循“信道好则多分配，信道差则少分配或不分配”的原则，确保在总功率约束下最大化系统容量。3. 提供对比机制，同时实现等功率分配（Equal Power Allocation）算法作为性能基准。4. 仿真分析模块能够计算不同信噪比（SNR）环境下的信道容量（Capacity）和误码率（BER），并进行统计分析。5. 可视化输出功能，能够绘制子载波的功率分配直方图、信道频率响应幅度图、以及注水算法与等功率分配在容量和误码率上的性能对比曲线，直观展示优化效果。

详 情 说 明

OFDM Water-Filling Power Allocation Simulation

项目介绍

本项目实现了一套基于MATLAB的正交频分复用（OFDM）系统仿真平台，专注于研究和验证信道功率分配算法对系统性能的影响。项目核心在于对比注水算法（Water-filling Algorithm）与等功率分配（Equal Power Allocation）在频率选择性瑞利衰落信道下的性能差异。通过蒙特卡洛仿真，项目量化评估了两种策略在信道容量（Capacity）和误码率（BER）方面的表现，并提供了直观的可视化图表来展示功率分配的微观细节和宏观趋势。

功能特性

• OFDM 基础通信链路：构建了包含随机比特生成、16-QAM 调制、功率加载、多径信道传输、迫零（Zero Forcing）均衡及解调的完整链路。
• 频率选择性衰落信道：模拟具有6条多径抽头的瑞利衰落信道，每次迭代动态生成时域冲击响应并转换为频域CSI（信道状态信息）。
• 自适应功率分配：

* 等功率分配：作为基准，将总功率平均分配给所有子载波。 * 注水算法：基于CSI动态计算最优水位线，实现“优信道多分配，差信道少分配或舍弃”的策略，以最大化香农容量。

• 性能评估指标：同时计算香农信道容量（理论上限）和实际误码率（BER），支持宽范围信噪比（0-30dB）扫描。
• 多维度可视化：提供信道响应快照、功率分配直方图、水位线示意图以及容量/误码率的对比曲线。

系统要求

• MATLAB R2016b 或更高版本
• Communications Toolbox（用于 qammod, qamdemod, biterr 等函数）

---

仿真实现细节与逻辑

本项目的仿真逻辑主要集中在 main 函数及其子函数中，通过外层信噪比（SNR）循环和内层蒙特卡洛（Monte Carlo）迭代来实现统计分析。

1. 系统参数初始化

代码首先定义了核心通信参数：

• OFDM参数：64个子载波，16点循环前缀（CP）。
• 调制方式：16-QAM，每个符号携带4比特信息。
• 功率约束：总发射功率设定为30 dBm（线性值为1 Watt），并在处理过程中保持恒定。
• 信道模型：6径多径信道，并在每次生成后进行能量归一化。

2. 信道模型构建

在每次蒙特卡洛迭代中，程序生成独立的复高斯随机变量作为多径信道的时域冲击响应。通过FFT变换获得频域信道响应（H），并计算各子载波的信道功率增益（|H|^2），作为后续功率分配的依据。噪声方差根据当前信噪比和总功率反推计算，假设路径损耗归一化。

3. 核心功率分配策略

仿真中并行执行两种分配方案：

• 方案 A (等功率分配)：直接构建一个长度为64的向量，每个元素值为 P_total / 64。
• 方案 B (注水算法)：调用优化器函数，输入当前信道增益和噪声水平。算法根据干扰水平（噪声/信道增益）计算“水位线”，确保在满足总功率约束的前提下，使得分配功率与干扰水平之和为常数。对于信道条件极差的子载波，分配功率被置为0（不传输数据）。

4. 通信链路与性能计算

• 容量计算：利用香农公式，分别基于两种算法分配后的功率向量计算理论信道容量。
• 信号传输：

* 生成随机二进制比特流并进行16-QAM调制。 * 功率加载：根据分配结果，对调制符号乘以功率幅度的平方根。对于注水算法，若某子载波分配功率为0，则该载波不发送数据。 * 信道通过：信号经过多径信道并叠加高斯白噪声。 * 接收处理：接收端使用迫零检测（Zero Forcing）进行均衡，除以信道频域响应和功率因子恢复原始符号。

• 误码率统计：针对注水算法，仅统计被激活（功率>0）的子载波上的误码情况；等功率分配则统计所有子载波。结果经过多次迭代取平均值。

5. 子函数设计

• water_filling_optimizer：这是项目的核心算法函数。它首先计算每个子载波的干扰水平倒数，然后对其排序。通过迭代方式寻找满足总功率约束的最优“水位线”（Water Level or Lambda）。算法会根据水位线动态剔除信道质量最差的子载波，直到找到一个子集，使得在该子集上的注水分配满足功率限制且均为正值。
• qam_modulator / qam_demodulator：封装了MATLAB通信工具箱的调制解调功能，负责比特流与复数星座点之间的转换，并处理维度的Reshape操作以适配OFDM子载波结构。

---

结果可视化

代码执行结束后会自动生成四个图表，用于详细分析算法性能：

1. Water-Filling Power Allocation Visualization：展示特定SNR下的信道频率响应（|H|^2），并绘制经典的“注水”图示。图中灰色区域代表干扰底噪（杯底），蓝色区域代表分配的信号功率（水），直观演示了水位线如何决定功率分配。
2. Power Allocation Histogram：以直方图形式对比注水算法与等功率分配在各个子载波上的功率数值。可以清晰看到注水算法在深衰落点停止分配功率，而在信道增益高处分配更多功率。
3. Capacity Comparison：绘制信道容量随SNR变化的曲线。展示注水算法在理论容量上始终优于或等于等功率分配，特别是在低信噪比区域优势更明显。
4. BER Comparison：绘制误码率随SNR变化的曲线（半对数坐标）。展示两种算法在系统可靠性方面的差异。

使用方法

直接在MATLAB环境中运行脚本即可启动仿真。程序的控制台会实时输出当前仿真处理的SNR点、对应的两种算法的容量及误码率数据。仿真全部完成后，所有图表将一次性弹出。可以通过修改脚本头部的 SNR_dB_range 或 n_iters 来调整仿真密度和精度。