1.算法描述
在通信系统中，信道所能提供的带宽通常比传送一路信号所需的带宽要宽得多。如果一个信道只传送一路信号是非常浪费的，为了能够充分利用信道的带宽，就可以采用频分复用的方法。
OFDM主要思想是：将信道分成若干正交子信道，将高速数据信号转换成并行的低速子数据流，调制到在每个子信道上进行传输。正交信号可以通过在接收端采用相关技术来分开，这样可以减少子信道之间的相互干扰(ISI) 。每个子信道上的信号带宽小于信道的相关带宽，因此每个子信道上可以看成平坦性衰落，从而可以消除码间串扰，而且由于每个子信道的带宽仅仅是原信道带宽的一小部分，信道均衡变得相对容易。
OFDM技术是HPA联盟(HomePlug Powerline Alliance)工业规范的基础，它采用一种不连续的多音调技术，将被称为载波的不同频率中的大量信号合并成单一的信号，从而完成信号传送。由于这种技术具有在杂波干扰下传送信号的能力，因此常常会被利用在容易受外界干扰或者抵抗外界干扰能力较差的传输介质中。
在通信系统中，信道所能提供的带宽通常比传送一路信号所需的带宽要宽得多。如果一个信道只传送一路信号是非常浪费的，为了能够充分利用信道的带宽，就可以采用频分复用的方法。
OFDM主要思想是：将信道分成若干正交子信道，将高速数据信号转换成并行的低速子数据流，调制到在每个子信道上进行传输。正交信号可以通过在接收端采用相关技术来分开，这样可以减少子信道之间的相互干扰(ISI) 。每个子信道上的信号带宽小于信道的相关带宽，因此每个子信道上可以看成平坦性衰落，从而可以消除码间串扰，而且由于每个子信道的带宽仅仅是原信道带宽的一小部分，信道均衡变得相对容易。
OFDM技术是HPA联盟(HomePlug Powerline Alliance)工业规范的基础，它采用一种不连续的多音调技术，将被称为载波的不同频率中的大量信号合并成单一的信号，从而完成信号传送。由于这种技术具有在杂波干扰下传送信号的能力，因此常常会被利用在容易受外界干扰或者抵抗外界干扰能力较差的传输介质中。
一个OFDM符号之内包含多个经过相移键控（PSK）或者正交幅度调制（QAM）的子载波。
一旦要把传输的比特分配到各个子载波上，某一种调制模式则将他们映射为子载波的幅度和相位，通常采用等效基带信号来描述OFDM的输出信号：
其中信号的实部和虚部分别对应OFDM的同相和正交分量，在实际系统可以分别与对应的c o s coscos分量和s i n sinsin分量相乘。 
2.仿真效果预览
matlab2022a仿真结果如下：
3.MATLAB核心程序
pilot_sym = pilot_generator(128); % Generate the pilot symbols, all in BPSK modulation
pilot_syms=pilot_sym';
pilot_syms=pilot_syms(:)';
idx=0;
for snr = 0:4:40
idx=idx+1;
nframe = 0;
nerror = 0;
while nframe < frame_num
nframe
info_bits = round(rand(1, N_DBPF)); % Generate the information bites randomly
tail = zeros(1,6); % Tail bits for convolutional encoder
info_bits = [info_bits tail]; % Append the tail bits to the infomation bits
trl = poly2trellis(7,[133 171]); % Trellis structure of convolutional encoder,
% Constraint length K = 7,generator polynomial in octal g0 = 133, g1 = 171
code_bits = convenc(info_bits, trl); % Convolutional encoder,rate R = 1/2
TX = tx_16qam_mod(code_bits); % 16-QAM modulation
figure(1)
a=real(TX);
b=imag(TX);
plot(a,b,'bx')
pilot_num = 1; % The start number of the cyclic pilot symbols
tx = tx_fre_to_time(TX,N_SPF,pilot_syms,pilot_num); % Serial to parallel,Add pilot symbols, IFFT ( Oversampling),Add cyclic prefix
temp_channel = channel_response(:,nframe+1).'; % Channel response at the current frame
ry = channel_multipath(tx,temp_channel); % Multi-path channel
ry = awgn(ry,snr); % Add addictive white gaussian noise
[RY,R_pilot] = ry_time_to_fre(ry,N_SPF,temp_channel); % Recover the time signal to frequency symbols of each subcarrier and Channel Equalization
figure(2)
a=real(RY);
b=imag(RY);
plot(a,b,'bx')
axis([-2 2 -2 2]);
hold on
Rcode_bits = ry_16qam_demod(RY); % Demodulate the received symbols
dec_bits = ry_sovadec(Rcode_bits, trl, N_DBPF,N_DBPF); % Soft output viterbi decoding, out = ln( p(info_bits(i)=1)/p(info_bits(i)=0));
for i = 1:N_DBPF
if dec_bits(i) >= 0
hard_det(i) = 1; % hard_det refers to the hard decision
else
hard_det(i) = 0;
end
if hard_det(i) ~= info_bits(i)
nerror = nerror+1;
end
end
nframe = nframe + 1;
end
BER(idx) = nerror/(frame_num * N_DBPF);
end