越来越笨了,dftmtx用法都老是忘记,举几个例子以免又忘了。。。
例1:
x=rand(1,Nfft);
y=dftmtx(Nfft)*x.' %和y=fft(x)是等价的
fft变换公式,n代表时域采样点,k代表频域采样点:
n=1:Nfft;
Y(k)=sum(x.*exp(-j*2*pi*n*k))
相应的,dftmtx(Nfft)产生的矩阵中,第k行,n列元素=exp(-j*2*pi*k*n/Nfft),与x.'相乘正好对应fft变换后的每个频点值。这样物理意义就清楚多了。
例2:
做信道矩阵相关时,假设信道频域相关矩阵R_ff,时域相关矩阵为R_tt,Nfft,那么
R_ff=dftmtx(Nfft)*R_tt*dftmtx(Nfft)'
时域相关矩阵R_tt=h*h',therefore
频域相关矩阵R_ff=H*H'=dftmtx(Nfft)*h*(dftmtx(Nfft)*h)'=dftmtx(Nfft)*h*h'*dftmtx(Nfft)=dftmtx(Nfft)*R_tt*dftmtx(Nfft)'
这样就不会忘了,Over