icimatrix.m

本程序求的OFDM中子载波干扰之间的相关性矩阵

%程序功能：求出子载波干扰系数矩阵C

% function c=ICIMatrix(fre)

%fre表示归一化频率偏移，fre=fd*Ts
%fd表示多普勒频移，Ts表示OFDM一个周期间隔，也就是OFDM信号一个周期的时间
clear
N=16;%表示子载波的个数
fre=0.0001
for t=-N+1:N-1
   at=sin(pi*(t+fre))
   bt=exp(-i*pi*(t+fre)*(N-1)/N)
   ct=(N*sin(pi*(t+fre)/N))
   dt(t+N+1)=at*bt/ct
end
u=1:2*N
% plot(u,abs(dt))


d=dt(N+1:2*N);



for p=1:N
    c(p,:)=d
    temp=d(N);
    for q=N:-1:2
        d(q)=d(q-1)
    end
    d(1)=temp;
end
%产生循环矩阵C

ChIR=[1.0000  0  0.7943  0  0.1259  0.1000]
ChIRf=fft(ChIR,N)

Htemp=diag(ChIRf)

H=c*Htemp

xx=1;