OFDM 信道估计

大家好，
       请教个问题，我现在想估计OFDM的时域信道。
       OFDM的频域信道可以使用LS算法，H[K]=Y[K]/X[K]。但，OFDM符号中有部分0载波，即该子载波上未发送有效数据，举例，802.11a中一个符号有64路子载波，有12个0.所以，我们计算出的H[K]只有52个点。那么，我现在想得到时域信道，不能对这仅有的52个点做IFFT吧？
       请教大家，这种情况，如何求时域信道。谢谢

X=Fx, Y=diag(H)Fx, F:NxL Fourier transform matrix , left block .x:Lx1 vector.
Actually, circulant matrix representation is  simpler, and equivalent to this.

Y=diag(X)Fh, 最小二乘求出h. 所以说你并没有看懂。

记住一点，FFT只是一个算法，DFT那个矩阵才是最本质的。

Y=diag(X)Fh， X不是Nfft 点，而是K 点，K为导频个数,K<N. F为KXL的矩阵。

最小二乘，不需要矩阵是方阵。。

算了，跟你解释不清楚。

中间没有估计，返回时域必然有信息损失，只能使用类似差值的方法获得

放弃吧。。。

※ FROM: 140.207.196]
※ 来源:·水木社区 http://www.newsmth.net·[FROM: 140.207.196]

零点的插值得到，曲线拟合也可以吧

Y=diag(X)*F*h，其中
  - Y(Kx1), X(K,1), F(KxL), h(Lx1), Y和X中零频点和高频被去掉了。
  - L是时域脉冲响应的长度，L<=K(否则未知数个数多于方程个数)，也就是加约束
       当K<n<N时，h(n)=0。
    相当于原方阵F里的N-K行和后面N-L列以及h(n)后面的N-L个元素被去掉了。
对于11a,
  N=64;
  Y=[Y(-26),Y(-25),...,Y(-1),Y(1),Y(2),...,Y(26)].';
  X=[X(-26),X(-25),...,X(-1),X(1),X(2),...,X(26)].';
  F=[
    exp(-j*2*pi/N*(-26)*[0:L-1]);
    exp(-j*2*pi/N*(-25)*[0:L-1]);
    ...
    exp(-j*2*pi/N*(-1)*[0:L-1]);
    exp(-j*2*pi/N*(+1)*[0:L-1]);
    exp(-j*2*pi/N*(+2)*[0:L-1]);
    ...
    exp(-j*2*pi/N*(+26)*[0:L-1])];
  h=[h(0),h(1),h(2), ..., h(L-1)].';
从Y=diag(X)*F*h，也就是H=F*h (H=diag(X)^-1*Y) 得到
  h = (F'*F)^-1 * F' * H

我觉得直接对着52个点做逆傅立叶变换就可以，因为fft本身有插值的作用

fft不是随便缺哪些点都可以插值的，楼主问的就是如何在只缺少0频信息的情况下准确估
计信道响应。

为什么不能做fft，fft不含乘法，直接做的话从数学上没有问题，只是得到的时域响应没有那些0点的信息，不然就在频域先插值，然后再做

做dft不考虑0频相当于把时域响应减去或者加上了一个常数。

为什么是一个常数而不是0呢？对于频域到时域的转换，0值对应的频点应该没有贡献啊

idft：x(n)=[X(-K)*exp(j*2*pi/N*K*n)+...+X(0)*exp(j*2*pi/N*K*0)+...]/N
对任何n，X(0)*exp(j*2*pi/N*K*0)/N项都是X(0)/N

这个不是dc吗，另外如果X(0)=0则没有贡献啊

你这假设压根就不成立。对于11a，0频（dc）和正负高频部分未知，"未知" != "等于
0"。

所以才需要频域插值啊，例如dvb-t，可以每隔三个点插出一个导频，另外边界的地方导频更加密集

先做频域插值，常用的线性插值，mmse插值都可以。之后再做ifft