基于matlab实现OFDM的编码.doc

clc; clear all; close all; fprintf(OFDM系统仿真\n); carrier_count=input(输入系统仿真的子载波数: \n);%子载波数128,64,32,16 symbols_per_carrier=30;%每子载波含符号数 bits_per_symbol=4;%每符号含比特数,16QAM调制 IFFT_bin_length=1024;%FFT点数 PrefixRatio=1/4;%保护间隔与OFDM数据的比例 1/6~1/4 GI=PrefixRatio*IFFT_bin_length ;%每一个OFDM符号添加的循环前缀长度为1/4*IFFT_bin_length ,即256 beta=1/32;%窗函数滚降系数 GIP=beta*(IFFT_bin_length+GI);%循环后缀的长度40 SNR=10; %信噪比dB %================信号产生=================================== baseband_out_length=carrier_count*symbols_per_carrier*bits_per_symbol;%所输入的比特数目 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_out_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,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]); grid on %==========分配载波到指定的IFFT位置========================== IFFT_modulation=zeros(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(2,1:IFFT_bin_length)),b*-) grid on axis ([0 IFFT_bin_length -0.5 4.5]); ylabel(幅值); xlabel(频率); title(OFDM载波幅度谱); figure(4); plot(0:IFFT_bin_length-1, (180/pi)*angle(IFFT_modulation(2,1:IFFT_bin_length)),go) hold on stem(0:carriers-1, (180/pi)*angle(IFFT_modulation(2,1:carriers)),b*-); stem(0:conjugate_carriers-1, (180/pi)*angle(IFFT_modulation(2,1:conjugate_carriers)),b*-); axis ([0 IFFT_bin_length -200 +200]) grid on ylab