资源描述：
《基于matlab实现OFDM的编码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、clc;clearall;closeall;fprintf('OFDM系统仿真');carrier_count=input('输入系统仿真的子载波数:');%子载波数128,64,32,16symbols_per_carrier=30;%每子载波含符号数bits_per_symbol=4;%每符号含比特数,16QAM调制IFFT_bin_length=1024;%FFT点数PrefixRatio=1/4;%保护间隔与OFDM数据的比例1/6~1/4GI=PrefixRatio*IFFT_
2、bin_length;%每一个OFDM符号添加的循环前缀长度为1/4*IFFT_bin_length,即256beta=1/32;%窗函数滚降系数GIP=beta*(IFFT_bin_length+GI);%循环后缀的长度40SNR=10;%信噪比dB%================信号产生===================================baseband_out_length=carrier_count*symbols_per_carrier*bits_per_symbo
3、l;%所输入的比特数目carriers=(1:carrier_count)+(floor(IFFT_bin_length/4)-floor(carrier_count/2));%共轭对称子载波映射复数数据对应的IFFT点坐标conjugate_carriers=IFFT_bin_length-carriers+2;%共轭对称子载波映射共轭复数对应的IFFT点坐标rand('twister',0);%每次产生不相同得伪随机序列baseband_out=round(rand(1,baseband_o
4、ut_length));%产生待调制的二进制比特流figure(1);stem(baseband_out(1:50));title('二进制比特流')axis([0,50,0,1]);%==============16QAM调制====================================complex_carrier_matrix=qam16(baseband_out);%列向量complex_carrier_matrix=reshape(complex_carrier_matrix
5、',carrier_count,symbols_per_carrier)';%串并转换，转换为symbols_per_carrier*carrier_count矩阵figure(2);plot(complex_carrier_matrix,'*r');%16QAM调制后星座图title('16QAM调制后星座图')axis([-4,4,-4,4]);gridon%==========分配载波到指定的IFFT位置==========================IFFT_modulation=z
6、eros(symbols_per_carrier,IFFT_bin_length);%添0组成IFFT运算IFFT_modulation(:,carriers)=complex_carrier_matrix;%未添加导频信号，子载波映射在此处IFFT_modulation(:,conjugate_carriers)=conj(complex_carrier_matrix);%共轭复数映射figure(3);stem(0:IFFT_bin_length-1,abs(IFFT_modulation(
7、2,1:IFFT_bin_length)),'b*-')gridonaxis([0IFFT_bin_length-0.54.5]);ylabel('幅值');xlabel('频率');title('OFDM载波幅度谱');figure(4);plot(0:IFFT_bin_length-1,(180/pi)*angle(IFFT_modulation(2,1:IFFT_bin_length)),'go')holdonstem(0:carriers-1,(180/pi)*angle(IFFT_mo
8、dulation(2,1:carriers)),'b*-');stem(0:conjugate_carriers-1,(180/pi)*angle(IFFT_modulation(2,1:conjugate_carriers)),'b*-');axis([0IFFT_bin_length-200+200])gridonylabel('相位')xlabel('频率')title('OFDM载波相位谱')%========通过IFFT将频域转化为时域，得到时域信号===============sig
