信道估计.docVIP

重庆交通大学信息科学与工程学院 综合性设计性实验报告 业： 学 号： 姓 名： 实验所属课程： 实验室(中心)： 软件与通信实验中心 指 导 教 师 ： 2014年3月  教师评阅意见： 签名： 年 月 日 实验成绩： 一、题目 二、仿真要求要求一：系统的数据传输 ①传输的数据随机产生； ②； 要求二：性能仿真 设计仿真方案，三、仿真方案详细设计 换言之 基于DFT的信道估计算法 四、仿真结果及结论 五、总结与体会 六、主要仿真代码 %基于 LS和DFT 对信道进行估计用线性插入 clear all;clc; Nfft=32;Ng=Nfft/8;Nofdm=Nfft+Ng; pilot_loc=[1 5 9 13 17 21 25 29 32];%l插入位置 Nps=length(pilot_loc);%每个ofdm的插入序列 Nbps=4;M=2^Nbps;%调制方式 SNR=30; Xp=2*(randn(1,Nps)0)-1;%插入序列的产生 msgint=randint(1,Nfft-Nps,M);%随机产生比特流 data=qammod(msgint,M);%16QAM调制 %插入导频 X=[Xp(1),data(1:3),Xp(2),data(4:6),Xp(3),data(7:9),Xp(4),data(10:12),Xp(5),data(13:15),Xp(6),data(16:18),Xp(7),data(19:21),Xp(8),data(22:23),Xp(9)]; %OFDM 的调制 x=ifft(X,Nfft);xt=[x(Nfft-Ng+1:Nfft) x];%傅里叶逆变换和加入循环前缀 h=[(randn+1i*randn),(randn+1i*randn)/2];%产生两径的信道 %信道的频域和长度 H=fft(h,Nfft);ch_length=length(h); H_power_dB=10*log10(abs(H.*conj(H))); y_channel=conv(xt,h);%通过信道的信号 yt=awgn(y_channel,SNR,measured); y=yt(Ng+1:Nofdm);Y=fft(y);%移去 CP 和FFTLS %基于LS 信道估计 H_est1=LS_CE(Y,Xp,pilot_loc,Nfft); H_est_power_dB1=10*log10(abs(H_est1.*conj(H_est1))); %基于DFT的信道估计 h_est2=ifft(H_est1);h_DFT2=h_est2(1:ch_length); H_est2=fft(h_DFT2,Nfft); H_est_power_dB2=10*log10(abs(H_est2.*conj(H_est2))); %比较LS和DFT两种估计的性能 subplot(2,1,1); plot(1:Nfft,H_power_dB,-,1:Nfft,H_est_power_dB1,s); legend(True Channel,LS); xlabel(subcarrier index);ylabel(Power[dB]); subplot(2,1,2); plot(1:Nfft,H_power_dB,-,1:Nfft,H_est_power_dB2,s); legend(True Channel,LS with DFT); xlabel(subcarrier index);ylabel(Power[dB]); function H_LS=LS_CE(Y,Xp,pilot_loc,Nfft) %LS 信道估计函数 %输入 %Y为频域接收信号 %Xp=插入信号 %pilot_loc=插入导频的位置 %N为FFT的大小 k=1:length(pilot_loc); LS_est(k)=Y(pilot_loc(k))./Xp(k); %线性插入 H_LS=interp1(pilot_loc,LS_est,[1:Nfft],linear); clear all; close all; clc; fprintf(OFDM信道估计仿真\n\n); num_sc = 1024; %载波数目 2048简化为1024 num_symbol = 14; %OFDM符号个数 data_