OFDM系统原理及其实现

1、通信系统综合设计报告题目：ofdm系统原理及其实现 学 部： 班 级： 姓 名： 学 号： 指导教师： 撰写日期： 目录第一章31.1要求31.2系统基本原理及基本模块31.2.1设计思路31.2.2系统基本模块4第二章52.1编程思路及框架52.1.1信道编码映射52.1.2串并/并串变换62.1.3ofdm调制解调72.1.4添加/取出循环前缀7第三章83.1 实验结果83.1.1码率计算：83.1.2试验结果8总结12附录13第一章1.1要求仿真实现ofdm调制解调，在发射端，经串/并变换和ifft变换，加上保护间隔（又称“循环前缀”），形成数字信号，通过信道到达接收端，结束端实现反变换

2、，进行误码分析。1.2系统基本原理及基本模块 1.2.1设计思路ofdm的基本原理就是把高速的数据流通过串并变换，分配到传输速率相对较低的若干个子信道中进行传输。由于每个子信道中的符号周期会相对增加，因此可以减轻由无线信道的多径时延扩展所产生的时间弥散性对系统造成的影响。并且还可以在ofdm符号之间插入保护间隔，令保护间隔大于无线信道的最大时延扩展，这样就可以最大限度地消除由于多径而带来的符号间干扰(isi)。而且，一般都采用循环前缀作为保护间隔，从而可以避免由多径带来的子载波间干扰(ici) 。 调制原理： 设ofdm系统中有n个子信道，第k个子信道采用的子载波为：ofdm信号 n路子信号之

3、和 一个码元时间 内任意两个子载波正交条件：子载频条件：子载频最小间隔：ofdm信号解调原理：根据子载波正交性质：1.2.2系统基本模块此ofmd包含4个模块，编码映射、串并/并串变换、ofdm调制/解调、添加/取出循环前缀等。第二章2.1编程思路及框架ofdm实现框架如下2.1.1信道编码映射当调制好的信号在信道里进行传输的时候，必然要受到信道的影响，导致在接收端信号解调出现错误，使系统的误码率大大增加。因而，为减少差错，提高通信系统抗干扰能力和可靠性，在一个实用的通信系统中，采用信道编码这一有效措施。方法是在发送端的信号中加入一定的冗余度，而在接收端这些冗余度可以用来检测并且纠正信号通过信

4、道后产生的错误。当然，冗余度的加入降低了系统的工作效率，但是和系统误码率的降低（即信号更加正确地传送）相比，这些代价是可以接受的。2.1.2串并/并串变换设ofdm系统的输入信号为串行二进制码元，其码元序列时间为t，先将此输入码元序列分为帧，每帧中有f个码元，即有f比特。然后将此f比特分成n组，每组中的比特数可以不同。如下图所示。2.1.3ofdm调制解调调制离散傅里叶变换离散傅里叶逆变换 ；子载波信号ofdm信号复数形式解调解调即为其逆过程，用傅里叶变换完成2.1.4添加/取出循环前缀ofdm系统中,每个并行数据支路都是窄带信号,可近似认为每个支路都经历平坦衰落,这样就减小了频率选择性衰落对

5、信号的影响。同时,每路子数据流速率的降低,减小了符号间干扰( isi) 。此外,还可以通过加保护间隔的办法完全消除符号间干扰。假设每个ofdm符号由y个样值组成,由于时延扩展,接收端将会有和信道冲激响应持续时间相对应的前l (l y) 个样值发生错误,为此,可以在发送信号前端加上m个样值,接收端收到信号时,先去掉前m个样值,然后再进行fft,只要m l就可完全消除isi。由于多径效应等因素，两个连续的分组信号之间会产生码间串扰，前一分组的信号可能会影响到后一分组信号的前面几个码元第三章3.1 实验结果3.1.1码率计算：bit_errors = find(baseband_in = baseb

6、and_out);bit_error_count = size(bit_errors,2)baseband_out_lengthk=bit_error_countb=baseband_out_lengthber=k/b（误码率）3.1.2试验结果原始数据（下图）：发送端数据（下图）：调制信号（下图）：还原出来的信号（下图）:3.2 实验分析k为误码数量：为5b为波形总量：为2000误码率ber=k/b=0.04%16qam误比特率为5.96%总结得出ofdm系统的误码率远低于16qam等，由此可以见ofdm系统具有一定优越性。总结通过本次实验，我懂得了如何仿真实现ofdm调制解调，在发射端，经

7、串/并变换和ifft变换，加上保护间隔（又称“循环前缀”），形成数字信号，通过信道到达接收端，结束端实现反变换，进行误码分析。通过本次实验，我加深了ofdm的基本原理就是把高速的数据流通过串并变换，分配到传输速率相对较低的若干个子信道中进行传输。由于每个子信道中的符号周期会相对增加，因此可以减轻由无线信道的多径时延扩展所产生的时间弥散性对系统造成的影响。并且还可以在ofdm符号之间插入保护间隔，令保护间隔大于无线信道的最大时延扩展，这样就可以最大限度地消除由于多径而带来的符号间干扰(isi)。而且，一般都采用循环前缀作为保护间隔，从而可以避免由多径带来的子载波间干扰(ici)。此次实验不但给我

8、巩固了课本的基础知识，还加强了我的动手能力。附录clear all;close all;idft_bin_length = 1024; % dft的点数carrier_count = 200; % 载波的数量 bits_per_symbol = 2; % 每个符号代表的比特数symbols_per_carrier = 50; % 每个载波使用的符号数snr = 10; % 信道中的信噪比（db） baseband_out_length=carrier_count*symbols_per_carrier*bits_per_symbol;%总比特数carriers = (1:carrier_cou

9、nt) + (floor(idft_bin_length/4) - floor(carrier_count/2);conjugate_carriers = idft_bin_length - carriers + 2;%发送端%产生随机二进制数据:baseband_out = round(rand(1,baseband_out_length);convert_matrix=reshape(baseband_out,bits_per_symbol,length(baseband_out)/bits_per_symbol);for k = 1:(length(baseband_out)/bits_

10、per_symbol) modulo_baseband(k) = 0; for i = 1:bits_per_symbolmodulo_baseband(k)=modulo_baseband(k)+convert_matrix(i,k)*2(bits_per_symbol-i);endend% 串并转换carrier_matrix=reshape(modulo_baseband,carrier_count,symbols_per_carrier);% 对每一个载波的符号进行差分编码carrier_matrix = zeros(1,carrier_count);carrier_matrix;fo

11、r i = 2:(symbols_per_carrier + 1) carrier_matrix(i,:)=rem(carrier_matrix(i,:)+carrier_matrix(i-1,:),2bits_per_symbol);end% 把差分符号代码转换成相位carrier_matrix = carrier_matrix * (2*pi)/(2bits_per_symbol);% 把相位转换成复数x,y=pol2cart(carrier_matrix,ones(size(carrier_matrix,1),size(carrier_matrix,2);complex_carrier_

12、matrix = complex(x,y);% 分配载波到指定的ifft位置idft_modulation = zeros(symbols_per_carrier + 1, idft_bin_length);idft_modulation(:,carriers) = complex_carrier_matrix;idft_modulation(:,conjugate_carriers) = conj(complex_carrier_matrix);% 画出频域中的ofdm信号代表figure (1)stem(0:idft_bin_length-1, abs(idft_modulation(2,

13、1:idft_bin_length),b*-)grid onaxis (0 idft_bin_length -0.5 1.5)ylabel(magnitude)xlabel(idft bin)title(ofdm carrier frequency magnitude)% 通过ifft将频域转化为时域，得到时域信号time_wave_matrix = ifft(idft_modulation);time_wave_matrix = time_wave_matrix;%画出一个符号周期的时域ofdm信号figure (2)plot(0:idft_bin_length-1,time_wave_ma

14、trix(2,:)grid onylabel(amplitude)xlabel(time)title(ofdm time signal, one symbol period)%画出每一个载波对应的时域信号（分离的ofdm信号）for f = 1:carrier_count temp_bins(1:idft_bin_length)=0+0j; temp_bins(carriers(f)=idft_modulation(2,carriers(f); temp_bins(conjugate_carriers(f)=idft_modulation(2,conjugate_carriers(f); te

15、mp_time = idft(temp_bins); figure(3) plot(0:idft_bin_length-1, temp_time) hold onendgrid onylabel(amplitude)xlabel(time)title(separated time waveforms carriers)for i = 1:symbols_per_carrier + 1windowed_time_wave_matrix(i,:)=real(time_wave_matrix(i,:).*hamming(idft_bin_length);windowed_time_wave_matr

16、ix(i,:) = real(time_wave_matrix(i,:);end%串并转换ofdm_modulation=reshape(windowed_time_wave_matrix,1,idft_bin_length*(symbols_per_carrier+1);% 画出整个时域ofdmtemp_time = idft_bin_length*(symbols_per_carrier+1);figure (4)plot(0:temp_time-1,ofdm_modulation)grid onylabel(amplitude (volts)xlabel(time (samples)ti

17、tle(ofdm time signal)% 上变频，这个模型中我们把经过ifft运算后ofdm直接发送tx_data = ofdm_modulation;%信道% the channel model is gaussian (awgn) +multipath(时延为1) tx_signal_power = var(tx_data);linear_snr = 10(snr/10);noise_sigma = tx_signal_power/linear_snr;noise_scale_factor = sqrt(noise_sigma);noise = randn(1, length(tx_d

18、ata)*noise_scale_factor;copy1=zeros(1,length(ofdm_modulation);for i=2:length(ofdm_modulation) copy1(i)=ofdm_modulation(i-1);endrx_data = tx_data + noise;% 根据符号长度和符号数将串行的符号转换为并行的% - 每一列是符号周期rx_data_matrix=reshape(rx_data,idft_bin_length, symbols_per_carrier + 1);%对每一列信号做fft得到频域信号rx_spectrum = dft(rx_

19、data_matrix);% 抽取接收信号中有载波的点rx_carriers = rx_spectrum(carriers,:);% 计算载波的相位% - 弧度转换为角度% - 归一化相位（0-360）rx_phase = angle(rx_carriers)*(180/pi);phase_negative = find(rx_phase 0);rx_phase(phase_negative) = rem(rx_phase(phase_negative)+360,360);% 用diff（）计算相位差rx_decoded_phase = diff(rx_phase);phase_negativ

20、e = find(rx_decoded_phase 0);rx_decoded_phase(phase_negative) = rem(rx_decoded_phase(phase_negative)+360,360);% 相位转化为符号base_phase = 360/2bits_per_symbol;delta_phase = base_phase/2;rx_decoded_symbols= zeros(size(rx_decoded_phase,1),size(rx_decoded_phase,2);for i = 1:(2bits_per_symbol - 1) center_phase = base_phase*i; plus_delta = center_phase+delta_phase; minus_delta = center_phase-delta_phase;decoded=find(rx_