作者:闻雪儿1116_414 | 来源:互联网 | 2023-08-16 18:30
目录
1.算法概述
2.部分程序
3.算法部分仿真结果图
4.完整程序获取
CSDN用户:我爱C编程
CSDN主页:https://blog.csdn.net/hlayumi1234567?type=blog
擅长技术:智能优化,路径规划,通信信号,图像处理,深度学习,控制器等。
1.算法概述
信道估计器是接收机一个很重要的组成部分。在OFDM系统中,信道估计器的设计上要有两个问题:** 一是导频信息的选择,由于无线信道的时变特性,需要接收机不断对信道进行跟踪,因此导频信息也必须不断的传送: 二是既有较低的复杂度又有良好的导频跟踪能力的信道估计器的设计,在确定导频发送方式和信道估计准则条件下,寻找最佳的信道估计器结构。 **在实际设计中,导频信息的选择和最佳估计器的设计通常又是相互关联的,因为估计器的性能与导频信息的传输方式有关。
通过信道佔计算法,接收机可以得到信道的冲激响应,在现代无线通信系统中,信道的信息已经得到了充分的利用。自适应的信道均衡器利用信道估计来对抗ISI的影响。分集技术利用信道佔计,实现与接收信号最佳匹配的接收机。最大似然检测通过信道估计使得接收端错误概率最小化。此外,信道估计的个重要的好处在于它使得相干解调成为可能。因为相干解调需要知道信号的相位信息,与非相干解调相比,可以提高系统的整体性能,而信道估计技术使之成为可能。
2.部分程序
......................................................for ii = 1:length(EsNodB)disp('EsN0dB is :'); disp(EsNodB(ii));tic;ChMSE_LS = 0;ChMSE_LMMSE=0; TDMSE_LMMSE =0;TDDMSE_LMMSE=0;TDQabsMSE_LMMSE =0;for mc = 1:MC
% Random channel tapsg = randn(L,1)+1i*randn(L,1);g = g/norm(g);H = fft(g,nFFT);
% generation of symbolX = randi([0 M-1],nFFT,1); %BPSK symbolsXD = modulate(modObj,X)/sqrt(10); % normalizing symbol powerx = F'*XD;xout = [x(nFFT-nCP+1:nFFT);x];
% channel convolution and AWGNy = conv(xout,g);nt =randn(nFFT+nCP+L-1,1) + 1i*randn(nFFT+nCP+L-1,1);No = 10^(-EsNodB(ii)/10);y = y + sqrt(No/2)*nt;
% Receiver processingy = y(nCP+1:NT);Y = F*y;
% frequency doimain LS channel estimation HhatLS = Y./XD; ChMSE_LS = ChMSE_LS + ((H -HhatLS)'*(H-HhatLS))/nFFT;
% Frequency domain LMMSE estimatiOnRhh= H*H';W = Rhh/(Rhh+(beta/snr(ii))*eye(nFFT));HhatLMMSE = W*HhatLS;ChMSE_LMMSE = ChMSE_LMMSE + ((H -HhatLMMSE)'*(H-HhatLMMSE))/nFFT;
% Time domain LMMSE estimatiOnghatLS= ifft(HhatLS,nFFT);Rgg = g*g';WW = Rgg/(Rgg+(beta/snr(ii))*eye(L));ghat = WW*ghatLS(1:L);TD_HhatLMMSE = fft(ghat,nFFT);% TDMSE_LMMSE = TDMSE_LMMSE + ((H -TD_HhatLMMSE)'*(H-TD_HhatLMMSE))/nFFT; % Time domain LMMSE estimation - ignoring channel covarianceghatLS = ifft(HhatLS,nFFT);Rgg = diag(g.*conj(g));WW = Rgg/(Rgg+(beta/snr(ii))*eye(L));ghat = WW*ghatLS(1:L);TDD_HhatLMMSE = fft(ghat,nFFT);% TDDMSE_LMMSE = TDDMSE_LMMSE + ((H -TDD_HhatLMMSE)'*(H-TDD_HhatLMMSE))/nFFT; % Time domain LMMSE estimation - ignoring smoothing matrixghatLS = ifft(HhatLS,nFFT);TDQabs_HhatLMMSE = fft(ghat,nFFT);% TDQabsMSE_LMMSE = TDQabsMSE_LMMSE + ((H -TDQabs_HhatLMMSE)'*(H-TDQabs_HhatLMMSE))/nFFT; end
.......................................................
3.算法部分仿真结果图
A157
4.完整程序获取
使用版本matlab2017b
解压密码:C+123456
获得方式1(面包多下载):
获取方式2:
如果下载链接失效,加博主微信,或私信。