Matlab实现OFDM调制学习资料

1、Mat l ab 实现 OFDM调制Matlab实现OFD函制一、实验目的1、进一步加深对matlab的了解和使用，熟练掌握 matlab的相关 库函数。2、学习OFD蜩制的原理，通过实验加深对其理解。3、学会用matlab实现BPS颓制和QPS调制。二、实验原理OFD蜩制原理：OFDMOFDM(Orthogonal Frequency DivisionMultiplexing) 即正交频分复用技术，实际上 OFD娓MCM MultiCarrierModulation ,多载波调制的一种。其主要思想是：将信道分 成若干正交子信道，将高速数据信号转换成并行的低速子数据流， 调制到在每个子信道上进

2、行传输。正交信号可以通过在接收端采用 相关技术来分开，这样可以减少子信道之间的相互干扰ICI 。每个子信道上的信号带宽小于信道的相关带宽，因此每个子信道上的可 以看成平坦性衰落，从而可以消除符号间干扰。而且由于每个子信 道的带宽仅仅是原信道带宽的一小部分，信道均衡变得相对容易。在向B3G/4G寅进的过程中，OFD塌关键的技术之一，可以结合分 集，时空编码，干扰和信道间干扰抑制以及智能天线技术，最大限 度的提高了系统性能。包括以下类型： V-OFDM,W-OFDM,F- OFDM,MIMO-OFDW-OFDM:、实验内容实验要求：(1)不加噪声时，调制出来的信号与原信号进行对1 、比。BPSK(

3、QPSK)解调解调信号bit2、BPS缄制:f(0)=1,f(1)=-1代码：original=randint(1,512*100,2); % 源为 0 和 1 的随机序列，为 1 行51200列；用512个子载波，做100次bpsk=(-1)Aoriginal; %BPSK 调制，将源中的0替换为-1 ,生成1行51200列的1和-1序列m=1;ioriginal=zeros(1,512*100); %生成零矩阵，定义变量初值ofdm=zeros(100,512);iofdm=zeros(100,512);iiofdm=zeros(1,512*100);cp=zeros(100,512+10

4、);cpsend=zeros(1,(512+10)*100);icp=zeros(100,512);ibpsk=zeros(1,512*100);while m=101ofdm(m,:)=bpsk(m-1)*512+1):(m*512); % 将 bpsk 矩阵转化为 100行512列的矩阵ofdmofdm(m,:)=ifft(ofdm(m,:),512); %OFDM 调制，返回 512 点的逆向 DFTcp(m,1:10)=ofdm(m,503:512); %添加长度为 10 的循环前缀，cp(m,11:end)=ofdm(m,:); % 生成矩阵 cp 为 100* (512+10)cp

5、send(m-1)*(512+10)+1):(m*(512+10)=cp(m,:);%cpspend 为 1*(512+10)*100)m=m+1;endcp1=cpsend; %没有添加噪声是调制生成的矩阵snr=-2:1:5; % 信噪比ber=zeros(1,length(snr);light=1;for t=-2:1:5m=1;% 解调while m=101icp(m,:)=cpsend(m-1)*512+m*10+1):(m*(512+10); %去前缀10iofdm(m,:)=fft(icp(m,:),512); %将 icp 做返回 512 点的 DFTiiofdm(m-1)*5

6、12+1):(m*512)=iofdm(m,:); %iiofdm为 1*(512*100)m=m+1;endfor n=1:(512*100)real_iiofdm=real(iiofdm(n);image_iiofdm=imag(iiofdm(n);distance1=sqrt(real_iiofdm+1)A2+image_iiofdmA2); %根据距离来判断是0还是1distance0=sqrt(real_iiofdm-1)A2+image_iiofdmA2);if distance1=distance0ibpsk(n)=0;elseibpsk(n)=1; end endiorigin

7、al=ibpsk; %计算误码率flag=(original=ioriginal); %相等时 flag 为 1,不相等时 flag 为0error=numel(find(flag=0); % 统计 0 的个数ber(light尸error/51200; %计算误码率light=light+1;endsemilogy(snr,ber)BPSKM制的误码率曲线图:3、QPSKW制：f(00)= (1+i)/sqrt(2),f(01)= (-1+i)/sqrt(2),f(10)= (-1-i)/sqrt(2),f(11) =(1-i)/sqrt(2)代码：source=randint(1,512*

8、100,2);sqpsk=source;qpsk=zeros(1,512*50);m=1;n=1;while m=(512*100+1)if sqpsk(1,m:(m+1)=0,0qpsk(1,n)=1/sqrt(2)*(1+i);elseif sqpsk(1,m:(m+1)=0,1qpsk(1,n)=1/sqrt(2)*(-1+i);elseif sqpsk(1,m:(m+1)=1,0qpsk(1,n)=1/sqrt(2)*(-1-i);else qpsk(1,n)=1/sqrt(2)*(1-i);endm=m+2;n=n+1;endm=1;ssource=zeros(1,512*100);

9、ofdm=zeros(50,512);oofdm=zeros(50,512);ooofdm=zeros(1,512*50);cp=zeros(50,512+10);cpsend=zeros(1,(512+10)*50);ccp=zeros(50,512);qqpsk=zeros(1,512*100);for m=1:50ofdm(m,:)=qpsk(m-1)*512+1):(m*512);ofdm(m,:)=ifft(ofdm(m,:),512);cp(m,1:10)=ofdm(m,503:512);cp(m,11:end)=ofdm(m,:);cpsend(m-1)*(512+10)+1):

10、(m*(512+10)=cp(m,:);endcp1=cpsend;snr=-2:1:5;ber=zeros(1,length(snr);as=1;for t=-2:1:5cpsend=awgn(cp1,t,measured);m=1;for m=1:50ccp(m,:)=cpsend(m-1)*512+m*10+1):(m*(512+10);oofdm(m,:)=fft(ccp(m,:),512);ooofdm(m-1)*512+1):(m*512)=oofdm(m,:);ends=1;for m=1:(512*50)real_ooofdm=real(ooofdm(m);image_ooof

11、dm=imag(ooofdm(m);distance0=sqrt(real_ooofdm-1)A2+(image_ooofdm-1)A2);distance1=sqrt(real_ooofdm+1)A2+(image_ooofdm-1)A2);distance2=sqrt(real_ooofdm+1)A2+(image_ooofdm+1)A2);distance3=sqrt(real_ooofdm-1)A2+(image_ooofdm+1)A2); bb(1,1:4)=distance0,distance1,distance2,distance3;aa(1,1)=min(bb);if aa(1,1)=distance0qqpsk(1,s:(s+1)=0,0;elseif aa(1,1)=distance1qqpsk(1,s:(s+1)=0,1;elseif aa(1,1)=distance2qqpsk(1,s:(s+1)=1,