在MATLAB中生成OFDM信号通过高斯白噪声的仿真代码

以下是在MATLAB中生成OFDM信号并添加高斯白噪声的仿真代码：

clear all;close all;clc;

% OFDM参数设置
N = 64;           % 子载波数量
cp_len = 16;      % 循环前缀长度
data_len = N/2-1; % 数据载波数量
Fs = 100e3;       % 采样率
Fc = 2e3;         % 载波频率
Ts = 1/Fs;        % 采样周期

% 生成随机数据
data = randi([0 1],1,data_len);

% QPSK调制
mod_data = pskmod(data,4);

% 将数据映射到子载波上
signal = zeros(1,N);
signal(2:N/2) = mod_data(1:data_len/2);
signal(N/2+2:N) = mod_data(data_len/2+1:data_len);

% IFFT变换
ifft_signal = ifft(signal)*sqrt(N);

% 插入循环前缀
tx_signal = [ifft_signal(N-cp_len+1:N) ifft_signal];

% 发送信号
t = (0:length(tx_signal)-1)*Ts;
carrier_wave = exp(1i*2*pi*Fc*t);
tx_signal = tx_signal.*carrier_wave;

% 加入高斯白噪声
snr = 10; % 信噪比
rx_signal = awgn(tx_signal,snr,'measured');

% 接收端解调
rx_signal = rx_signal.*conj(carrier_wave);
rx_signal = rx_signal(cp_len+1:end);
rx_signal = fft(rx_signal)/sqrt(N);

% 将载波数据映射回原始数据
rx_data = zeros(1,data_len);
rx_data(1:data_len/2) = rx_signal(2:data_len/2+1);
rx_data(data_len/2+1:data_len) = rx_signal(N/2+2:N);

% QPSK解调
demod_data = pskdemod(rx_data,4);

% 统计误码率
err = sum(abs(data-demod_data));
ber = err/data_len;
fprintf('误码率: %f\n',ber);

代码中首先设置OFDM的参数，包括子载波数量、循环前缀长度、数据载波数量、采样率和载波频率等。然后生成随机数据，并进行QPSK调制。接着将数据映射到子载波上，并进行IFFT变换和插入循环前缀。然后将信号进行调制，并添加高斯白噪声。接收端对信号进行解调，并将载波数据映射回原始数据，并进行QPSK解调。最后统计误码率。

需要注意的是，代码中的误码率是理论误码率，实际误码率可能会受到信道的影响而有所偏差。