% for i=1:length(bitdata)/4
% qam_data(i,:)=mapping(source(i),:);%数据映射 % end
% complex_qam_data=complex(qam_data(:,1),qam_data(:,2)); % %组合为负数形式，形成16QAM信号 % figure(1);
% plot(complex_qam_data,'*r');%绘制16QAM星座图 % title('16QAM调制后星座图'); % grid on ; figure(1);
plot(complex_carrier_matrix,'*r');%绘制16QAM星座图 title('16QAM调制后星座图'); grid on ;
% 程序功能
% （1）通过对OFDM信号各个子载波赋共轭对称的数据产生一个实OFDM符号。
% （2）给OFDM符号加循环前缀与循环后缀
% （3）给OFDM符号加窗。在程序中加入的是升余弦窗，可以通过改变升余弦窗的滚降系数，观察加入不同升余弦窗对OFDM信号频谱的的影响
% （4）信道采用嘉兴高斯白噪声信道。可以通过改变信噪比来改变信道环境，从而在接受端通过误码率或是星座图观察信道对OFDM信号传输的影响
% （5）去除循环前缀与循环后缀，对OFDM信号进行解调 % 程序代码： clear all; close all;
carrier_count=200;
%这个程序中OFDM子载波个数为512，其中400即carrier_count *2为数据符号，其余赋0值
symbols_per_carrier=20;
%每个子载波上的符号数，在这里即是OFDM符号的个数 bits_per_symbol=4;
%OFDM符号的每个子载波上传输的比特数，4bit通常采用的16QAM调制 IFFT_bin_length=512;
?T的长度，也即一个OFDM符号的子载波个数； PrefixRatio=1/4;
%循环前缀的比值，循环前缀与OFDM符号长度的壁纸，通常在1/6~1/4之间 GI=PrefixRatio * (IFFT_bin_length);%保护间隔的长度，这里为128 beta=1/32;%升余弦窗的滚降系数
GIP=beta*(IFFT_bin_length+GI);%循环后缀的长度，这里为20；
SNR=16;%本程序考虑加性高斯白噪声信道，这里信噪比为30db %===============================OFDM=====================
baseband_out_length=carrier_count*symbols_per_carrier * bits_per_symbol; %计算传输数据总的比特数，为200*20*4
%bits=16000bits。16000bits构成20个OFDM符号，每个OFDM符号200个子载波，每个子载波传输4bit信息
carriers=(1:carrier_count)+(floor(IFFT_bin_length/4)-floor(carrier_count/2)); %计算OFDM符号的在载波的序号，carriers 中存放的序号是29~228； conjugate_carriers=IFFT_bin_length-carriers +2; % conjugate_carriers=IFFT_bin_length-carriers -2;
%====================================================这个序号有问题
%计算OFDM符号子载波的序号，conjugate_carriers中存放的序号是282~481； rand('twister',0);
baseband_out=round(rand(1,baseband_out_length));
%产生16000bit待传输的二进制比特流。这里存放的是发送的二进制信号与后面解调后的二进制信号比较，可以计算误码率。
QAM调制病绘制星座图。
complex_carrier_matrix=qam16(baseband_out);
%调用子程序qam16进行16QAM调制。将base_band中的二进制比特流，每4bit转换为一个16QAM信号，即将二进制比特流每4bit转换为-3-3j，-3+3j,3-3j,3+3j,
%-1-3j,-1+3j,1-3j,1+3j,-3-j,-3+j,3-j,3+j,-1-j,-1+j,1-j,1+j中的一个。转换后complex_carrier_matrix为1*4000矩阵。
complex_carrier_matrix=reshape(complex_carrier_matrix',carrier_count,symbols_per_carrier);
%转换complex_carrier_matrix中的数据为carrier_coun6t*symbols_per_carrier矩阵，这里为20*200矩阵。
信号产生
figure(1);
plot(complex_carrier_matrix,'*r');%绘制16QAM星座图 title('16QAM调制后星座图'); grid on ;
%IFFT,即进行OFDM调制
IFFT_modulation=zeros(symbols_per_carrier,IFFT_bin_length);
%将symbols_per_carrier*IFFT_bin_length矩阵赋0值，这里讲20*512矩阵赋0值。这里512是IFFT的长度，也是OFDM符号的子载波个数
%这里明显有问题,以下矫正
% IFFT_modulation(:,carriers)=complex_carrier_matrix; IFFT_modulation(:,carriers)=complex_carrier_matrix';
%将20*200的complex_carrier_matrix的数据赋给IFFT_modulation的第29~228列，即给512个子载波中的第29~229个子载波赋值
% IFFT_modulation(:,conjugate_carriers)=conj(complex_carrier_matrix); IFFT_modulation(:,conjugate_carriers)=conj(complex_carrier_matrix');
%将20*200的complex_carrier_matrix的数据付给512个子载波中的第282~481个子载波，这段程序构造了512个在载波的OFDM符号，并且各个
%子载波上的数据是共轭对称的。这样做的目的是经过IFFT后形成的OFDM符号均为实数。另外，在512个子载波中，仅有400个子载波为数据，其余为
%0值，相当于补零。补零的目的是通常IFFT的长度应该为2的证书次幂。 signal_after_IFFT=ifft(IFFT_modulation,IFFT_bin_length,2); %IFFT实现OFDM调制
time_wave_matrix=signal_after_IFFT; figure(2);
plot(0:IFFT_bin_length-1,time_wave_matrix(2,:)); %画一个OFDM信号的时域表现 axis([0,512,-0.4,0.4]); grid on;
ylabel('Amplitude'); xlabel('Time');
title('OFDM Time Signal,One Symbol Period'); %添加循环前缀与循环后缀
XX=zeros(symbols_per_carrier,IFFT_bin_length +GI+GIP);
%IFFT_bin_length+GI+GIP为OFDM，循环前缀，循环后缀长度之和； for k=1:symbols_per_carrier; for i=1:IFFT_bin_length;
XX(k,i+GI)=signal_after_IFFT(k,i); end for i=1:GI
XX(k,i)=signal_after_IFFT(k,i+IFFT_bin_length-GI); %添加循环前缀 end for j=1:GIP
XX(k,IFFT_bin_length+GI+j)=signal_after_IFFT(k,j); %添加循环后缀 end end
time_wave_matrix_cp=XX;%带循环前缀和循环后缀的OFDM符号 figure(3)
plot(0:length(time_wave_matrix_cp)-1,time_wave_matrix_cp(2,:)); %画带循环前缀与循环后缀的OFDM信号的时域波形、 axis([0,600,-0.3,0.3]);
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库完整的OFDM系统的仿真实现(2)在线全文阅读。
相关推荐：