实验-现代通信系统

《实验-现代通信系统》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、ChinaUniversityofGeosciences课程名称：现代通信系统报告题目：现代通信系统实验报告院系：机械与电子信息工程学院专业班级：071121姓名：李饪周学号：20121000471老师：陈振兴日期：2015年12月现代通信系统“实验一一、实验目的1、掌握MATLAB的编程方法。2、理解OFDM系统及基本的信道估计算法。二、实验原理对于发送的OFDM信号，间隔的插入己知信号来对传输信道进行估计，从而提高系统的差错性能。插入的已知信号分为训练序列和导频,或其他的插入方式,本次实验采用训练序列的方法。三、实验过程1、建立一个新的M文件。2、输入程序!1!实验结果LS_Ray2DOF

2、DM(100,100,8)0.0518K=1005101520253035404550K=50K=100K=200K=500K=1()OOLS_Ray2D0FDM(100,100,4)ans=0.0118K=10K=5()K=100K=200K=500K=10()0由图像可知：随着K值得不断变大，图像下降的越来越快。例程：functionSER=LS_Rey2DOFDM(SNR,K,ConstelldtionSize)NumCycle=l;Subcarrier=128;Pilotlnverval=3;NumOFDM=PilotInverval*100;NumPilot=NumOFDM/Pilo

3、tInverval;PilotSignal=randint(NumPilot,Subcarrier,[1,10]);Demapper=Mapper2D(1:ConsteHationSize,ConsteHationSize);NumSymbo丄Sum=0;forcycle=l:NumCycleNumErrorSymol=0;Datalnput=randint(1,NumOFDM*Subcarrier,[1,ConsteHationSize]);Delndex=reshape(Datalnput,NumOFDM,Subcarrier);DataFqDm=Mapper2D(Datalnput,Co

4、nstellationSize);DataFqDm=reshape(DataFqDm,NumOFDM,Subcarrier);TxDataFqDm=zeros(NumOFDM+NumPilot,Subcarrier);ReDataFqDm=zeros(NumOFDM+NumPilot,Subcarrier);a=l;forpilot=1:NumPilotTxDataFqDm(a,:)=PilotSignal(pilot,:);TxDataFqDm(a+1:a+PilotInverval,:)=DataFqDm(pilot*PilotInverval-PilotInverval+1:pilot*

5、PilotInverval,:);a=a+PilotInverval+1;endTxDataTiDm=ifft(TxDataFqDm」);%^)j6dp^E±6dTxDataTiDm=TxDataTiDm.';forb=l:NumOFDM+NumPilot%DAeAl1yDApASubTx=TxDataTiDm(b,:);L=length(SubTx);Powersignal=sum(abs(SubTx).A2/L);Sigma=sqrt(PowerSignal*10A(-SNR/10));H_s=sqrt(1/2)*(randn(1,L)+li*randn(lzL));H=sqrt(K/(K

6、+l))*ones(1,L)+sqrt(1/(K+l))*H_s;Noise=(1/sqrt(2))*(randn(1,L)+li*randn(1,L));SubTx=SubTx.*H+Sigma*Noise;ReDataTiDm(b,:)=SubTx;endTempFqDm=fft(ReDataTiDm.1);ReDataFqDm=TempFqDm.1;d=l;h=l;forc=l:PilotInverval+1:NumOFDM+NumPilotH_LS=ReDataFqDm(c,:)•/Pilotsignal(dz:);d=d+l;fore=l:PilotTnvervBlReSiLS=Re

7、DataFqDm(c+e,:)./H_LS;ReS=real(ReSiLS);ImS=imag(ReSiLS);DArray=zeros(ConstellationSizezSubcarrier);forf=l:ConstelldtionSizeDTemp=(ReS-real(Demapper(f))).A2+(ImS-imag(Demapper(f))).A2;DArray(fz:)=DTemp