MATLAB生成BPSK水声通信发射波形从比特流到WAV文件关键词MATLAB BPSK、水声通信仿真、BPSK WAV、载波频率、采样率、符号率、通信原理课程设计水声通信的发射端最容易被低估的一步是把 0/1 比特真正变成一条可播放、可写入 WAV 的实值声学波形。这篇文章用一个可复现的 MATLAB 项目从随机比特流出发生成带根升余弦脉冲成形的 BPSK 通带波形并保存为 WAV 文件。本文只讨论发射端同步、解调、信道估计和 BER 属于接收端工作不在本文范围内。在这里插入图片描述本文配套项目ZoperIOT/underwater-acoustic-modulation1. BPSK 在水声通信中做什么BPSKBinary Phase Shift Keying二进制相移键控用两个相位表示 1 bit 信息。最常见的映射可以写成b k ∈ { 0 , 1 } ⟶ a k 2 b k − 1 ∈ { − 1 , 1 } b_k \in \{0,1\}\quad\longrightarrow\quad a_k 2b_k-1 \in \{-1,1\}bk​∈{0,1}⟶ak​2bk​−1∈{−1,1}映射后的符号先经过根升余弦RRC脉冲成形再与载波相乘得到实值通带信号。写成连续时间形式可理解为s ( t ) ℜ { ∑ k a k g ( t − k T s ) e j 2 π f c t } s(t)\Re\left\{\sum_k a_k g(t-kT_s)e^{j2\pi f_c t}\right\}s(t)ℜ{k∑​ak​g(t−kTs​)ej2πfc​t}其中g ( t ) g(t)g(t)是脉冲成形滤波器T s 1 / R s T_s1/R_sTs​1/Rs​是符号周期f c f_cfc​是载波频率。对初学者来说BPSK 是很好的发射端起点映射清楚、波形直观也便于后续扩展到 QPSK、QAM 或 OFDM。2. 环境与项目准备下载或克隆项目underwater-acoustic-modulation。在 MATLAB 中将当前目录切到项目根目录。将src加入路径。项目只使用 MATLAB 基础数值函数RRC 滤波器已经包含在仓库内不依赖 Communications Toolbox。3. 先理解这 4 个参数下面的配置是项目中适合先跑通的起点96 kHz 采样率、10 kHz 载波、200 Baud 符号率。参数含义本文取值设置时的硬约束 / 建议sampleRate/f s f_sfs​每秒采样点数96000 HzsampleRate / symbolRate必须是整数carrierFrequency/f c f_cfc​声学载波频率10000 Hz必须小于f s / 2 f_s/2fs​/2且要落在换能器可用频段symbolRate/R s R_sRs​每秒发送符号数200 Baud越高越快但带宽需求和多径敏感性通常也会上升rolloff/α \alphaαRRC 滚降系数0.5取值范围为 0 到 1影响带宽与时域拖尾本例的每符号采样点数为s p s f s R s 96000 200 480 \mathrm{sps}\frac{f_s}{R_s}\frac{96000}{200}480spsRs​fs​​20096000​480也就是说每个 BPSK 符号在离散波形中用 480 个采样点表示。项目会检查这个值是否为整数不满足时会直接报错避免生成时间轴不一致的波形。4. 一段可直接运行的 MATLAB 代码将下面代码保存为项目根目录下的examples/demo_bpsk_to_wav.m或直接在命令窗口逐段运行。%% MATLAB 生成 BPSK 水声通信发射波形并写入 WAVclear;clc;close all;% 1) 将项目源码加入 MATLAB 路径addpath(src);% 2) 固定随机种子保证每次生成的比特流可复现rng(2026);numberOfBits256;bitsrandi([01],1,numberOfBits);% 3) 设置发射端参数optionsstruct(...sampleRate,96000,...% fs96 kHzcarrierFrequency,10000,...% fc10 kHzsymbolRate,200,...% Rs200 Baudrolloff,0.5,...% RRC 滚降系数filterSpan,10,...% 滤波器跨度符号数phaseOffset,0);% 初始载波相位% 4) 比特流 - RRC 成形 - BPSK 实值通带波形[waveform,info]BPSK(bits,options);% 5) 写入标准 WAV 文件outputFilebpsk_10kHz_200baud.wav;audiowrite(outputFile,waveform,info.sampleRate);% 6) 在命令窗口查看本次生成的信息fprintf(调制方式%s\n,info.modulation);fprintf(采样率%.0f Hz\n,info.sampleRate);fprintf(载波频率%.0f Hz\n,info.carrierFrequency);fprintf(符号率%.0f Baud\n,info.symbolRate);fprintf(时长%.3f s\n,info.durationSeconds);fprintf(输出文件%s\n,outputFile);% 7) 绘制时域波形展示前 35 mst(0:numel(waveform)-1)/info.sampleRate;figure(Color,w);plot(t*1e3,waveform,b);xlim([0,min(35,t(end)*1e3)]);grid on;xlabel(时间 / ms);ylabel(幅度);title(BPSK 水声通带波形前 35 ms);% 8) 绘制归一化频谱nfft2^nextpow2(numel(waveform));% 手写 Hann 窗避免额外工具箱依赖n0:numel(waveform)-1;window0.5-0.5*cos(2*pi*n/(numel(waveform)-1));Xfft(waveform.*window,nfft);f(0:nfft/2)*info.sampleRate/nfft;magdB20*log10(abs(X(1:nfft/21))/max(abs(X))eps);figure(Color,w);plot(f/1e3,magdB,LineWidth,1.2);xlim([0,20]);ylim([-100,5]);grid on;xlabel(频率 / kHz);ylabel(归一化幅度 / dB);title(BPSK 水声通带波形频谱);运行后会在当前目录得到bpsk_10kHz_200baud.wav该 WAV 的采样率为 96 kHz样本为已归一化的双精度实值信号范围在[-1, 1]内。5. 代码内部实际发生了什么项目中的BPSK(bits, options)做了四件事检查bits是否仅包含 0 和 1用2 * bits - 1映射到-1/1进行 RRC 脉冲成形乘以复载波并取实部最后按峰值归一化为实值通带波形。归一化的作用是把波形峰值控制在[-1,1]一方面满足audiowrite的常用浮点输入范围另一方面避免后续写入 DAC 或音频链路时因幅度过大而削顶。注意归一化不等于完成了硬件幅度标定真正接功放和换能器时仍需按设备的满量程、增益和安全声压校准。6. 波形和频谱应该怎么看上图是同一组参数下的教学示意。你运行上一节的 MATLAB 脚本后会得到由当前比特流和项目 RRC 成形器生成的本机精确波形与频谱。时域图中能看到 10 kHz 附近的快速振荡而符号切换处的包络并不是生硬跳变——这是 RRC 脉冲成形带来的结果。频谱图的能量集中在 10 kHz 载波附近。对滚降系数为α \alphaα的理想升余弦系统基带单边带宽可用下面的近似关系理解B ≈ ( 1 α ) R s 2 B\approx \frac{(1\alpha)R_s}{2}B≈2(1α)Rs​​本例R s 200 R_s200Rs​200Baud、α 0.5 \alpha0.5α0.5因此主能量大致围绕 10 kHz 在数百 Hz 的窄带范围内展开。实际频谱还会受到有限长度、窗函数和滤波器截断的影响。7. 三个高频问题7.1 为什么sampleRate / symbolRate不能随便设该项目的单载波波形按“每符号固定采样点数”进行上采样和脉冲成形。因此sampleRate / symbolRate必须是整数且至少为 2。例如% 合法96000 / 200 480struct(sampleRate,96000,symbolRate,200)% 不合法96000 / 333 不是整数struct(sampleRate,96000,symbolRate,333)如果想把Rs调高建议优先从能整除fs的值选起例如在 96 kHz 下尝试 400、600、800、1200 Baud。7.2 载波频率为什么不能乱设最基本的数字信号约束是0 f c f s 2 0f_c\frac{f_s}{2}0fc​2fs​​否则会混叠。但满足奈奎斯特条件只是第一步。水声实验还要检查换能器频响、功放带宽、DAC/ADC 采样率以及水下声道在该频段的衰减。7.3 为什么不能直接用方波符号上载波可以但频谱旁瓣通常更宽。RRC 脉冲成形的价值在于控制占用带宽并为后续接收端的匹配滤波和定时恢复保留更规范的接口。对于课程设计或实验记录它也让波形从“能响”走向“可分析”。8. 下一步可以怎么扩展完成 BPSK 发射端后可以按下面顺序继续把BPSK替换为QPSK观察每符号 2 bit 时的比特分组约束比较不同rolloff和symbolRate下的频谱占用使用项目的OFDM入口生成导频辅助 QPSK-OFDM 发射波形加入预导、同步、信道、多径和接收端解调形成完整链路。项目地址再次放在这里https://github.com/ZoperIOT/underwater-acoustic-modulation。如果你是在做通信原理课程设计建议先把本节的 BPSK WAV 生成、时域图和频谱图跑通再往接收端扩展每一步都会更容易定位问题。本文代码对应项目的 BPSK 发射端实现。默认参数可在仓库src/auc/common_options.m中查看不同硬件条件下请以实际设备规格为准。