gusucode.com > m4k20功率分配图 mimo功率优化,OFDM典型仿真程序，贪婪算法仿真程序 > code/test_16qam/calcha.m

    function h=calcha(x,y,para,fftlen2,gilen,nd)

y_remcp= girem(y,fftlen2,gilen,nd);
y1=reshape(y_remcp,para,nd);
y1=fft(y1);   	%fft : built in function
y1=reshape(y1,1,para*nd);
x1=reshape(x,para,nd);
x1=fft(x1);   	%fft : built in function
x1=reshape(x1,1,para*nd);
h=y1./x1;