《OFDMmatlab源程序总结》由会员分享，可在线阅读，更多相关《OFDMmatlab源程序总结（24页珍藏版）》请在金锄头文库上搜索。
1、-比较完整的OFDM仿真，供大家学习下载。、%main_OFDM.m%一个相对完整的OFDM通信系统的仿真设计，包括编码，调制，IFFT，%上下变频，高斯信道建模，FFT，PAPR抑制，各种同步，解调和解码等模%块，并统括系统性能的仿真验证了系统设计的可靠性。clear allclose allclc%+全局变量+% seq_num 表示当前帧是第几帧% count_dds_up 上变频处的控制字的累加% count_dds_down 下变频处的控制字的累加整整% count_dds_down_tmp 下变频处的控制字的累加小数% dingshi 定时同步的定位% m_syn 记录定时同步中的
2、自相关平台global seq_numglobal count_dds_upglobal count_dds_downglobal count_dds_down_tmpglobal dingshi global m_syn%+% SNR_Pre 设定用于仿真的信噪比的初值% interval_SNR 设定用于仿真的信噪比间隔% frame_num 每一个信噪比下仿真的数据帧数% err_int_final 用于计算每帧出现的误比特数% fwc_down 设定的接收机初始载波频率控制字% fre_offset 设定接收机初始载波频率偏移调整量单位为Hz% k0 每次进入卷积编码器的信息比特数%
3、G 卷积编码的生成矩阵SNR_Pre=-5;interval_SNR=1;for SNR_System=SNR_Pre:interval_SNR:5frame_num=152;dingshi=250;err_int_final=0;fwc_down=16.050;fre_offset=0;k0=1;G=1 0 1 1 0 1 1;1 1 1 1 0 0 1 ;disp(-start-);for seq_num=1:frame_num, %frame_num 帧数%+以下为输入数据局部+datain=randint(1,90); %+%+以下为信道卷积编码局部+encodeDATA=v_encd
4、(G,k0,datain);%+%+信道交织编码+interlacedata=interlacecode(encodeDATA,8,24);%+%+以下为QPSK调制局部+QPSKdata=qpsk(interlacedata);%+%+生成训练序列+if seq_num3trainsp_temp=seq_train(); end%+%+插入导频+PILOT=(1+j);m_QPSKdata=QPSKdata;data2fft_temp=m_QPSKdata(1:8),PILOT,m_QPSKdata(9:16),PILOT,m_QPSKdata(17:24),PILOT,m_QPSKdata
5、(25:32),PILOT,m_QPSKdata(33:40),PILOT,m_QPSKdata(41:48),m_QPSKdata(49:56),PILOT,m_QPSKdata(57:64),PILOT,m_QPSKdata(65:72),PILOT,m_QPSKdata(73:80),PILOT,m_QPSKdata(81:88),PILOT,m_QPSKdata(89:end);%+trainsp_temp2=trainsp_temp,zeros(1,128);trainsp=trainsp_temp2(65:256),trainsp_temp2(1:64);%+降PAPR矩阵变换+m
6、ati*_data=nyquistimp_PS();matri*_mult=data2fft_temp*mati*_data;%+data2fft2=matri*_mult(65:128),zeros(1,128),matri*_mult(1:64);%+ifft运算+if seq_num=1 ifftin=trainsp;elseif seq_num=2 ifftin=trainsp;else ifftin=data2fft2;endIFFTdata=fft_my(conj(ifftin)/256);IFFTdata=conj(IFFTdata);% figure% plot(real(IF
7、FTdata)% *label(realIFFTdata)% figure% plot(imag(IFFTdata)% *label(imagIFFTdata) %+%+以下为插入循环前后缀，2倍升采样+data2fir=add_CYC_upsample(IFFTdata,2);% +% +fir低通滤波+guiyi_a=0.00172160.0101620.0255120.028801-0.0059219-0.060115-0.04960.0914310.296360.39560.296360.091431-0.0496-0.060115-0.00592190.0288010.0255120
8、.0101620.0017216 ;%抽样截止频率为128kHZ,通带截止频率为20kHZ，阻带截止频率为40kHZ，带纹波动小于1dB，带外衰减100dBt*FIRdatai=filter(guiyi_a,1,real(data2fir);t*FIRdataq=filter(guiyi_a,1,imag(data2fir);% +%+发射机cic滤波+CICidatai=cic_inter(t*FIRdatai,20);CICidataq=cic_inter(t*FIRdataq,20);%+%+上变频+fwc_up=16; %控制字可以选择DUCdata=up_convert_ofdm(fwc_up,CICidatai,CICidataq);%+%+高斯白噪声信道+DUCdata,datama*=guiyi_DUCdata(DUCdata);awgn_data=awgn(DUCdata,SNR_System); %+%*承受机*%+下变频+DUCdata_tmp=awgn_data;fwc_down=fwc_down+(fre_offset*128/2560000);r_fre_offset=2560000*(fwc_down-fwc_up)/128);DDCdatai,DDCdataq=down_convert_ofdm(fwc_down,DUCdata_tmp);%+%+
