齐鲁工业大学 2014 届本科毕业设计（论文）
参考文献
[1] 樊昌信,曹丽娜. 通信原理（第六版）[M].北京:国防工业出版社,2009 [2] 程相君,陈生潭. 信号与系统[M].西安:西安电子科技大学出版社,1990 [3] 曹志刚,钱压生. 现代通信原理[M].北京:清华大学出版社,1992.17-45 [4] 刘敏,魏玲. MATLAB通信仿真与应用[M].北京:国防工业出版社,2001 [5] 于万波. 基于MATLAB的图像处理（第2版）[M].清华大学出版社,2011 [6] 张晓赢．OFDM 系统信道估计研究[D]．湖南：国防科技大学,2005
[7] Rodger E.Ziemer,William H.Tranter. Principles of Communications Systems,Modulation and
Noise(Fifth Edition)[M]. Bei Jing：Higher Education Press,2000.9.1-8
[8] 刘学勇. 详解MATLAB/Simulink通信系统仿真与建模[M].北京:电子工业出版
社,2001.192-280
[9] 王文博,郑侃. 宽带无线通信OFDM技术(第二版)[M]. 内蒙古:人民邮电出版社,2007. 8-9
[10]尹长川,罗涛,乐光新. 多载波宽带无线通信技术[M]. 北京:北京邮电大学出版
社,2004.20-45
[11]周正兰. OFDM及其链路级平台的Simulink实现[J].中国数据通信,2003,(10):90-92 [12]樊昌信,徐炳祥. 通信原理[M]. 北京:清华大学出版社,2006.161-172
[13]邓华. MATIAB通信仿真及应用实例详解[M]. 北京:人民邮电出版社,2003.116-122 [14]楼顺天,姚若玉,沈俊霞.MATLAB 7.X 程序设计语言[M].西安电子科技大学出社.2007. 1-2
[15]夏细苟,胡亮,陈少平. 一种基于OFDM的多载波通信系统的设计与仿真[J]．电子工程
师,2005,(8):37-41
24
齐鲁工业大学 2014 届本科毕业设计（论文）
致 谢
本课题的设计是在孙老师的悉心指导和帮助下完成的，在设计与论文成功结束之时向孙老师表示深深的谢意。在本课题的思路、方案的设计以及在设计进行的过程中进行了正确、及时的指导建议，孙老师都投入了很多精力。
在论文完成的过程中老师给予了我最大的帮助和关心，我的毕业课题的顺利结束，包含了老师的辛勤的汗水。孙老师知识渊博，对工作认真负责，对我们严格要求，对知识精益求精，使我受到了很大的启发，对我以后走向工作岗位有很大的帮助，受益匪浅。同时，在我这个设计的过程中，孙老师一直都在关心着我课题进行的情况，随时对我的设计思路、参考文献等提出改正和参考，使我的设计能够得到合理正确的结果。在我写论文的时候也随时进行检查，对我论文的格式、内容、数据的处理等提出了正确及时的意见。对我能在规定时间内完成毕业课题及论文提供了很大的保证。因此，对孙老师再次表示深深的谢意。
感谢其他对我的课题设计提供知识帮助和实物参考的老师；感谢学校为我提供这个课题设计的机会和条件，让我对大学四年所学的知识梳理了一边，尤其是对理论知识的应用方面。
最后，谨向我尊敬的老师、参加论文评审预答辩的专家、教授致以诚挚的谢意、深深的祝福和崇高的敬意！
25
齐鲁工业大学 2014 届本科毕业设计（论文）
附 录
. .................................ofdm.m........................................................ clear all; close all; clc;
format long
...................................参数设置................................. %本次仿真载频为2GHz，带宽1MHz，子载波数128个，cp为16 %子载波间隔为7.8125kHz
%一个ofdm符号长度为128us，cp长度为16us %采用16QAM调制方式 %最大doppler频率为132Hz
%多径信道为5径，功率延迟谱服从负指数分布~exp(-t/trms),trms=(1/4)*cp时长，各径延迟取为delay=[0 2e-6 4e-6 8e-6 12e-6] pilot_inter=5;
%导频符号间隔为10,可调整，看不同导频间隔下的BER情况，和理论公式比较 pilot_symbol_bit=[0 0 0 1]; %导频为常数，对应星座点1+3*j cp_length=16; %cp长度为16 SNR_dB=[0 2 4 6 8 10 12 14 16 18 20 22]; ls_err_ber=zeros(1,length(SNR_dB)); lmmse_err_ber=zeros(1,length(SNR_dB)); lr_lmmse_err_ber=zeros(1,length(SNR_dB)); for i=1:length(SNR_dB)%每个SNR点上仿真若干次 ls_error_bit=0;
lmmse_error_bit=0; lr_lmmse_error_bit=0; total_bit_num=0;
loop_num=5; %共仿真10次 for l=1:loop_num
ofdm_symbol_num=100;%每次仿真产生100个ofdm符号,则每次仿真共有100×128个星座映射符号；16QAM调制下，1个星座映射符号包含4个bit bit_source=input_b(128,ofdm_symbol_num);%为每次仿真产生100个ofdm符号的比特个数，128为每个ofdm符号的子载波个数,相当于reshape [nbit,mbit]=size(bit_source);
total_bit_num=total_bit_num+nbit*mbit;
map_out=map_16qam(bit_source);%对一次仿真符号块进行16QAM映射 [insert_pilot_out,pilot_num,pilot_sequence]=insert_pilot(pilot_inter,pilot_symbol_bit,map_out);
26
齐鲁工业大学 2014 届本科毕业设计（论文）
%按块状导频结构，对映射后的结果插入导频序列 ofdm_modulation_out=ifft(insert_pilot_out,128); %作128点逆FFT运算，完成ofdm调制
ofdm_cp_out=insert_cp(ofdm_modulation_out,cp_length);%插入循环前缀 %******** 以下过程为ofdm符号通过频率选择性多径信道 ********** num=5;
%假设功率延迟谱服从负指数分布~exp(-t/trms),trms=(1/4)*cp时长； %t在0~cp时长上均匀分布
%若cp时长为16e-6s，可以取5径延迟如下 delay=[0 2e-6 4e-6 8e-6 12e-6]/4; trms=4e-6;
var_pow=10*log10(exp(-delay/trms)); fd=132;%最大doppler频率为132Hz t_interval=1e-6;%采样间隔为1us
counter=200000;%各径信道的采样点间隔，应该大于信道采样点数。由以上条件现在信道采样点数
count_begin=(l-1)*(5*counter);%每次仿真信道采样的开始位置 trms_1=trms/t_interval;
t_max=16e-6/t_interval;%信道采样点数，每个调制符号采一个点 passchan_ofdm_symbol=multipath_chann(ofdm_cp_out,num,var_pow,delay,fd,t_interval,counter,count_begin);
%******** 以上过程为ofdm符号通过频率选择性多径信道 ******* %******** 以下过程为ofdm符号加高斯白噪声 *********
snr=10^(SNR_dB(i)/10);
[nnl,mml]=size(passchan_ofdm_symbol); spow=0; for k=1:nnl for b=1:mml
spow=spow+real(passchan_ofdm_symbol(k,b))^2+imag(passchan_ofdm_symbol(k,b))^2;
end end
spow1=spow/(nnl*mml);
sgma=sqrt(spow1/(2*snr));
%sgma如何计算，与当前SNR和信号平均能量有关系 receive_ofdm_symbol=add_noise(sgma,passchan_ofdm_symbol);
%加入随机高斯白噪声，receive_ofdm_symbol为最终接收机收到的ofdm符号块
27
齐鲁工业大学 2014 届本科毕业设计（论文）
%***********以上过程为ofdm符号加高斯白噪声 ****************
cutcp_ofdm_symbol=cut_cp(receive_ofdm_symbol,cp_length);%去除循环前缀 ofdm_demodulation_out=fft(cutcp_ofdm_symbol,128);%作128点FFT运算，完成ofdm解调
%****** 以下就是对接收ofdm信号进行信道估计和信号检测的程****** ls_zf_detect_sig=ls_estimation(ofdm_demodulation_out,pilot_inter,pilot_sequence,pilot_num); %采用
LS
估计算法及迫零检测得到的接收信号
lmmse_zf_detect_sig=lmmse_estimation(ofdm_demodulation_out,pilot_inter,pilot_sequence,pilot_num,trms_1,t_max,snr);
%采用LMMSE估计算法及迫零检测得到的接收信号
low_rank_lmmse_sig=lr_lmmse_estimation(ofdm_demodulation_out,pilot_inter,pilot_sequence,pilot_num,trms_1,t_max,snr,cp_length);
%采用低秩LMMSE估计算法及迫零检测得到的接收信号
%** 以下就是对接收ofdm信号进行信道估计和信号检测的过程******* ls_receive_bit_sig=de_map(ls_zf_detect_sig);QAM解映射 lmmse_receive_bit_sig=de_map(lmmse_zf_detect_sig); lr_lmmse_receive_bit_sig=de_map(low_rank_lmmse_sig); %以下过程统计各种估计算法得到的接收信号中的错误比特数 ls_err_num=error_count(bit_source,ls_receive_bit_sig);
lmmse_err_num=error_count(bit_source,lmmse_receive_bit_sig);
lr_lmmse_err_num=error_count(bit_source,lr_lmmse_receive_bit_sig); ls_error_bit=ls_error_bit+ls_err_num;
lmmse_error_bit=lmmse_error_bit+lmmse_err_num;
lr_lmmse_error_bit=lr_lmmse_error_bit+lr_lmmse_err_num; end
%计算各种估计算法的误比特率
ls_err_ber(i)=ls_error_bit/total_bit_num;
lmmse_err_ber(i)=lmmse_error_bit/total_bit_num;
lr_lmmse_err_ber(i)=lr_lmmse_error_bit/total_bit_num; end
plot(SNR_dB,ls_err_ber,'b-*') hold on
plot(SNR_dB,lmmse_err_ber,'r-o') plot(SNR_dB,lr_lmmse_err_ber,'g-+') hold off
..................................................add_noise.m................................................................. function output=add_noise(sgma,input)
28