Matlab实现OFM调制一、实验目的 、进一步加深对malab的理解和使用，纯熟掌握matlb的有关库函数。 2、学习D调制的原理，通过实验加深对其理解。 、学会用matab实现BPSK调制和QSK调制。二、实验原理 FD调制原理:OFDM OFDM(Ohooal Feqeny ivion Mutiplexig)即正交频分复用技术，事实上ODM是MCM lt-CarrierMdulaion,多载波调制的一种。其重要思想是：将信道提成若干正交子信道,将高速数据信号转换成并行的低速子数据流,调制到在每个子信道上进行传播。正交信号可以通过在接受端采用有关技术来分开，这样可以减少子信道之间的互相干扰 ICI 。每个子信道上的信号带宽不不小于信道的有关带宽，因此每个子信道上的可以当作平坦性衰落，从而可以消除符号间干扰。并且由于每个子信道的带宽仅仅是原信道带宽的一小部分，信道均衡变得相对容易。在向B3G/G演进的过程中,OM是核心的技术之一,可以结合分集，时空编码，干扰和信道间干扰克制以及智能天线技术,最大限度的提高了系统性能。涉及如下类型:V-OFDM，W-ODM，FOF,IMFDM,多带-ODM。三、实验内容 实验规定:(1)不加噪声时,调制出来的信号与原信号进行对比。（2）加噪声时的误码率曲线图。用randint产生512乘以100的信号源 1、 调制原理框图 BPSK（QPSK）的数字调制加噪声添加长度为10的循环前缀512点的IFFT调制sourceMod调制信号CPOFDM调制 解调的原理框图OFDM解调去CP调制信号解调信号bitBPSK(QPSK)解调2、BPS调制：f()=,f(1)=-1代码:orilandit(1,512100,2); %源为0和的随机序列,为1行5100列；用512个子载波，做00次bsk=（-1).rigina; %PK调制，将源中的0替代为-,生成1行5200列的1和-1序列m1;iorgaeros(1,5*100); 生成零矩阵,定义变量初值od=zeros(100,1）;iodm=zers(100,52)；iifdm=ers(1,512*100);cp=zros(100,12+10);csndzros(1,(12+10)*100);i=zros(10,51);ibpsk=zeros(1,512*00）；wilm=1 ofdm(m，：)=bps(m-1)*5+1):(m512); 将b矩阵转化为100行512列的矩阵fdm fdm(m,：)ft(ofd(，:），512)； %FM调制，返回512点的逆向DF c(,:10)=ofm（m,503：51); %添加长度为10的循环前缀, p(,11:end）=ofd（m,：)； 生成矩阵p为100*（5210) cend（(（m-1)(512+1）+1）:(m*(52+1)cp（m,:);cpspnd为1*（51+10)*100） m=+1;ndcp1=csed; %没有添加噪声是调制生成的矩阵snr=2:1:5； 信噪比br=zeros（,length(nr);ight1；ort=-:1:5 m=1; %解调 wile =01 icp(m,:)=cpsend(m-1)*52+*10+1):(m*（512); %去前缀10 iofdm(,:)=t(cp(m,:),51); %将ic做返回12点的DT iiofdm(m-）12+1）：(m*512)=io(m，:); %iodm为（51*10） +1; e for n=1:(51*10) rl_ofdmrel（iif(n); imagifm=imag(idm(n）); disance=sqr(rel_iofdm1）2+magiifdm2); %根据距离来判断是还是 distcesqrt(（reaofdm1)iag_iiofdm2); if dtanc1=dsance0 ibpsk()0； ee ipsk(n)=； en e ioriginlbpsk; %计算误码率 fa=(riina=oriina); %相等时flag为1,不相等时lag为 ror=nume（fnd（flag=0）)； %记录0的个数 er(light)=eor50; %计算误码率 iht=ligt+1；endsemilogy(snr，er）BPSK调制的误码率曲线图:、K调制:f(0)= （1+i)qrt(2),f(1)(-1+i)/rt(), f(10)= （-1-i)/srt(2), f（11） =(1-)/sqrt(2）。 代码：sorc=radit(1，12*100,2)；qps=urce;qpsk=eros(1，0);m=1；=1;whle (51210+1) issk(1，m:(m1)）=0,0 qpsk(，n)=1sqrt（2)（1+i)； esef qpsk(1,m:(m+1)=0，1 ps(1,)/sqt（2)*(1); elsif qk（1,m:(m1)=， qsk(1，)/sqrt(）*(i）; es psk(1,n）=1/sr()*（1-i）; end m=+2; n=n+;ndm=1;sourcezros(1,512*00);ofdmzero(5,512）;ofdmzero(5，51)；oofdze(1,51250);c=ero（5,52+）;cpsend=zero(,(52+）*50）；cp=zro(50，12)；qqpsk=eos（,512*10);for m=1:50 ofd（m,:)=qpsk(-1)*512+）:(512)； ofm(m,：)=it(of(，:),512）; （m,1:1）=of(,503:52)； c(m，11：)=of(m,：)； cpend（(m-1)*(51+10)1):(m*(5210)）=cp(m，:）;nd=cpsnd;nr-2:1:5;br=eros（1,length（snr）;as=1;for t=-2:1：5 psend=awgn(cp1，t,msre); m=1; fr m=1：50 ccp(m,:)cpsnd(m1)*512+m*1+1):(（512+0)）； oofdm(，:)=fft（ccp（m,:），51); oofdm(（(m-1)5121):(m*51)=oofm（m，:); nd s=1; fr m=1:(5*0) real_ofdm=l(ooofm(); image_oodm=imag（oofdm（m); dsanc0sqt(real_ooo1）2+(ia_oofdm-)2); distace1qr（(real_ooodm+1)+(iage_oofdm-1)2); istaeqr(reloodm+1)2+(image_ofdm+1）2）; distance3=sqr(realooodm)2+(image_oofm+1)2); bb（1,1:4)=istnce0,istac,dsace,itace3; a(1,1）mn(bb）； f (1,)=disane0 qs(1,s:(s+1)=,0； eseif a(,1)distance1 qqpsk（1，:(s+1)=0,1; elseif (1,1)=dsnce2 qqpsk(，s:(s1)=1,0; ese aa(1，)=distace3 qpsk(1,:(1)）=1,1; end s=s+2; e ssource=qqsk； lg(sorce=source); erro=el(nd(flag=0）; be()=eror/520； ass+；dsemog(snr,br)；QPSK调制的误码率曲线图: